R11: a cis-regulatory node of the sea urchin embryo gene network that controls early expression of SpDelta in micromeres

Identification and prediction of developmental enhancers in sea urchin embryos

Background

The transcription of developmental regulatory genes is often controlled by multiple cis-regulatory elements. The identification and functional characterization of distal regulatory elements remains challenging, even in tractable model organisms like sea urchins.

Results

We evaluate the use of chromatin accessibility, transcription and RNA Polymerase II for their ability to predict enhancer activity of genomic regions in sea urchin embryos. ATAC-seq, PRO-seq, and Pol II ChIP-seq from early and late blastula embryos are manually contrasted with experimental cis-regulatory analyses available in sea urchin embryos, with particular attention to common developmental regulatory elements known to have enhancer and silencer functions differentially deployed among embryonic territories. Using the three functional genomic data types, machine learning models are trained and tested to classify and quantitatively predict the enhancer activity of several hundred genomic regions previously validated with reporter constructs in vivo.

Conclusions

Overall, chromatin accessibility and transcription have substantial power for predicting enhancer activity. For promoter-overlapping cis-regulatory elements in particular, the distribution of Pol II is the best predictor of enhancer activity in blastula embryos. Furthermore, ATAC- and PRO-seq predictive value is stage dependent for the promoter-overlapping subset. This suggests that the sequence of regulatory mechanisms leading to transcriptional activation have distinct relevance at different levels of the developmental gene regulatory hierarchy deployed during embryogenesis.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07936-0.

CRISPR-Cas9 editing of non-coding genomic loci as a means of controlling gene expression in the sea urchin

We seek to manipulate gene function here through CRISPR-Cas9 editing of cis-regulatory sequences, rather than the more typical mutation of coding regions. This approach would minimize secondary effects of cellular responses to nonsense mediated decay pathways or to mutant protein products by premature stops. This strategy also allows for reducing gene activity in cases where a complete gene knockout would result in lethality, and it can be applied to the rapid identification of key regulatory sites essential for gene expression. We tested this strategy here with genes of known function as a proof of concept, and then applied it to examine the upstream genomic region of the germline gene Nanos2 in the sea urchin, Strongylocentrotus purpuratus. We first used CRISPR-Cas9 to target established genomic cis-regulatory regions of the skeletogenic cell transcription factor, Alx1, and the TGF-β signaling ligand, Nodal, which produce obvious developmental defects when altered in sea urchin embryos. Importantly, mutation of cis-activator sites (Alx1) and cis-repressor sites (Nodal) result in the predicted decreased and increased transcriptional output, respectively. Upon identification of efficient gRNAs by genomic mutations, we then used the same validated gRNAs to target a deadCas9-VP64 transcriptional activator to increase Nodal transcription directly. Finally, we paired these new methodologies with a more traditional, GFP reporter construct approach to further our understanding of the transcriptional regulation of Nanos2, a key gene required for germ cell identity in S. purpuratus. With a series of reporter assays, upstream Cas9-promoter targeted mutagenesis, coupled with qPCR and in situ RNA hybridization, we concluded that the promoter of Nanos2 drives strong mRNA expression in the sea urchin embryo, indicating that its primordial germ cell (PGC)-specific restriction may rely instead on post-transcriptional regulation. Overall, we present a proof-of-principle tool-kit of Cas9-mediated manipulations of promoter regions that should be applicable in most cells and embryos for which CRISPR-Cas9 is employed.

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

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

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.

Notch signaling in the division of germ layers in bilaterian embryos

Bilaterian embryos are triploblastic organisms which develop three complete germ layers (ectoderm, mesoderm, and endoderm). While the ectoderm develops mainly from the animal hemisphere, there is diversity in the location from where the endoderm and the mesoderm arise in relation to the animal-vegetal axis, ranging from endoderm being specified between the ectoderm and mesoderm in echinoderms, and the mesoderm being specified between the ectoderm and the endoderm in vertebrates. A common feature is that part of the mesoderm segregates from an ancient bipotential endomesodermal domain. The process of segregation is noisy during the initial steps but it is gradually refined. In this review, we discuss the role of the Notch pathway in the establishment and refinement of boundaries between germ layers in bilaterians, with special focus on its interaction with the Wnt/β-catenin pathway.

Genome-wide use of high- and low-affinity Tbrain transcription factor binding sites during echinoderm development

Sea stars and sea urchins are model systems for interrogating the types of deep evolutionary changes that have restructured developmental gene regulatory networks (GRNs). Although cis-regulatory DNA evolution is likely the predominant mechanism of change, it was recently shown that Tbrain, a Tbox transcription factor protein, has evolved a changed preference for a low-affinity, secondary binding motif. The primary, high-affinity motif is conserved. To date, however, no genome-wide comparisons have been performed to provide an unbiased assessment of the evolution of GRNs between these taxa, and no study has attempted to determine the interplay between transcription factor binding motif evolution and GRN topology. The study here measures genome-wide binding of Tbrain orthologs by using ChIP-sequencing and associates these orthologs with putative target genes to assess global function. Targets of both factors are enriched for other regulatory genes, although nonoverlapping sets of functional enrichments in the two datasets suggest a much diverged function. The number of low-affinity binding motifs is significantly depressed in sea urchins compared with sea star, but both motif types are associated with genes from a range of functional categories. Only a small fraction (∼10%) of genes are predicted to be orthologous targets. Collectively, these data indicate that Tbr has evolved significantly different developmental roles in these echinoderms and that the targets and the binding motifs in associated cis-regulatory sequences are dispersed throughout the hierarchy of the GRN, rather than being biased toward terminal process or discrete functional blocks, which suggests extensive evolutionary tinkering.

Comparative Study of Regulatory Circuits in Two Sea Urchin Species Reveals Tight Control of Timing and High Conservation of Expression Dynamics

Accurate temporal control of gene expression is essential for normal development and must be robust to natural genetic and environmental variation. Studying gene expression variation within and between related species can delineate the level of expression variability that development can tolerate. Here we exploit the comprehensive model of sea urchin gene regulatory networks and generate high-density expression profiles of key regulatory genes of the Mediterranean sea urchin, Paracentrotus lividus (Pl). The high resolution of our studies reveals highly reproducible gene initiation times that have lower variation than those of maximal mRNA levels between different individuals of the same species. This observation supports a threshold behavior of gene activation that is less sensitive to input concentrations. We then compare Mediterranean sea urchin gene expression profiles to those of its Pacific Ocean relative, Strongylocentrotus purpuratus (Sp). These species shared a common ancestor about 40 million years ago and show highly similar embryonic morphologies. Our comparative analyses of five regulatory circuits operating in different embryonic territories reveal a high conservation of the temporal order of gene activation but also some cases of divergence. A linear ratio of 1.3-fold between gene initiation times in Pl and Sp is partially explained by scaling of the developmental rates with temperature. Scaling the developmental rates according to the estimated Sp-Pl ratio and normalizing the expression levels reveals a striking conservation of relative dynamics of gene expression between the species. Overall, our findings demonstrate the ability of biological developmental systems to tightly control the timing of gene activation and relative dynamics and overcome expression noise induced by genetic variation and growth conditions.

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

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

Developmental cis-regulatory analysis of the cyclin D gene in the sea urchin Strongylocentrotus purpuratus

Cyclin D genes regulate the cell cycle, growth and differentiation in response to intercellular signaling. While the promoters of vertebrate cyclin D genes have been analyzed, the cis-regulatory sequences across an entire cyclin D locus have not. Doing so would increase understanding of how cyclin D genes respond to the regulatory states established by developmental gene regulatory networks, linking cell cycle and growth control to the ontogenetic program. Therefore, we conducted a cis-regulatory analysis on the cyclin D gene, SpcycD, of the sea urchin, Strongylocentrotus purpuratus, during embryogenesis, identifying upstream and intronic sequences, located within six defined regions bearing one or more cis-regulatory modules each.

Gene regulatory control in the sea urchin aboral ectoderm: spatial initiation, signaling inputs, and cell fate lockdown

The regulation of oral-aboral ectoderm specification in the sea urchin embryo has been extensively studied in recent years. The oral-aboral polarity is initially imposed downstream of a redox gradient induced by asymmetric maternal distribution of mitochondria. Two TGF-β signaling pathways, Nodal and BMP, are then respectively utilized in the generation of oral and aboral regulatory states. However, a causal understanding of the regulation of aboral ectoderm specification has been lacking. In this work control of aboral ectoderm regulatory state specification was revealed by combining detailed regulatory gene expression studies, perturbation and cis-regulatory analyses. Our analysis illuminates a dynamic system where different factors dominate at different developmental times. We found that the initial activation of aboral genes depends directly on the redox sensitive transcription factor, hypoxia inducible factor 1α (HIF-1α). Two BMP ligands, BMP2/4 and BMP5/8, then significantly enhance aboral regulatory gene transcription. Ultimately, encoded feedback wiring lockdown the aboral ectoderm regulatory state. Our study elucidates the different regulatory mechanisms that sequentially dominate the spatial localization of aboral regulatory states.

Human developmental enhancers conserved between deuterostomes and protostomes

The identification of homologies, whether morphological, molecular, or genetic, is fundamental to our understanding of common biological principles. Homologies bridging the great divide between deuterostomes and protostomes have served as the basis for current models of animal evolution and development. It is now appreciated that these two clades share a common developmental toolkit consisting of conserved transcription factors and signaling pathways. These patterning genes sometimes show common expression patterns and genetic interactions, suggesting the existence of similar or even conserved regulatory apparatus. However, previous studies have found no regulatory sequence conserved between deuterostomes and protostomes. Here we describe the first such enhancers, which we call bilaterian conserved regulatory elements (Bicores). Bicores show conservation of sequence and gene synteny. Sequence conservation of Bicores reflects conserved patterns of transcription factor binding sites. We predict that Bicores act as response elements to signaling pathways, and we show that Bicores are developmental enhancers that drive expression of transcriptional repressors in the vertebrate central nervous system. Although the small number of identified Bicores suggests extensive rewiring of cis-regulation between the protostome and deuterostome clades, additional Bicores may be revealed as our understanding of cis-regulatory logic and sample of bilaterian genomes continue to grow.

Flies and worms have long served as valuable model organisms for the study of human development and health. Despite the great morphological and evolutionary distance between them, humans, flies, and worms share many commonalities. Each develops from three major germ layers and is patterned along the two major spatial axes. At the molecular level, development in these widely diverged species is often controlled by the same signaling pathways activating members of the same transcription factor and target gene families, shared since the common ancestor of humans, flies, and worms. And yet, at the gene regulatory level, humans and flies or worms seem starkly different, with not a single regulatory region shared across the phyla. Here we discover the first two examples of developmental enhancers conserved between deuterostomes (ranging from human to sea urchins) and protostomes (a large clade that includes flies and worms). We show evidence that these ancient regulatory loci retain the capacity to respond to the same signaling pathways in these widely diverged organisms, and we show that they have been co-opted, along with the molecular pathways that control them, to pattern the vertebrate nervous systems. Our screen supports large scale regulatory rewiring, while offering the first intriguing outliers.

Barcoded DNA-tag reporters for multiplex cis-regulatory analysis

Cis-regulatory DNA sequences causally mediate patterns of gene expression, but efficient experimental analysis of these control systems has remained challenging. Here we develop a new version of "barcoded" DNA-tag reporters, "Nanotags" that permit simultaneous quantitative analysis of up to 130 distinct cis-regulatory modules (CRMs). The activities of these reporters are measured in single experiments by the NanoString RNA counting method and other quantitative procedures. We demonstrate the efficiency of the Nanotag method by simultaneously measuring hourly temporal activities of 126 CRMs from 46 genes in the developing sea urchin embryo, otherwise a virtually impossible task. Nanotags are also used in gene perturbation experiments to reveal cis-regulatory responses of many CRMs at once. Nanotag methodology can be applied to many research areas, ranging from gene regulatory networks to functional and evolutionary genomics.

Cis-regulatory logic driving glial cells missing: self-sustaining circuitry in later embryogenesis

The glial cells missing (gcm) regulatory gene of the sea urchin Strongylocentrotus purpuratus is first expressed in veg2 daughter cells as the genomic target of late cleavage stage Delta-Notch signaling from the skeletogenic mesoderm precursors. Gcm is required in veg2 progeny during late cleavages for the early phase of pigment cell precursor specification. Here we report on a later acting cis-regulatory module that assumes control of gcm expression by the early mesenchyme blastula stage and maintains it through pigment cell differentiation and dispersal. Cis-perturbation analyses reveal that the two critical elements within this late module are consensus matches to Gcm and Six1 binding sites. Significantly, six1 mRNA localizes to gcm+cells from the mesenchyme blastula stage onwards. Trans-perturbations with anti-sense morpholinos reveal a co-dependency between six1 and gcm. Six1 mRNA levels fall sharply after Gcm is depleted, while depleting Six1 leads to significant reductions in output of endogenous gcm or modular-reporters. These results support the conclusion gcm and six1 comprise a positive intergenic feedback loop in the mesodermal GRN. This often employed cross regulatory GRN feature here ensures self-sustaining gcm output in a cohort of fully specified pigment cell precursors at a relatively early developmental stage.

A comprehensive analysis of Delta signaling in pre-gastrular sea urchin embryos

In sea urchin embryos Delta signaling specifies non-skeletogenic mesoderm (NSM). Despite the identification of some direct targets, several aspects of Delta Notch (D/N) signaling remain supported only by circumstantial evidence. To obtain a detailed and more complete image of Delta function we followed a systems biology approach and evaluated the effects of D/N perturbation on expression levels of 205 genes up to gastrulation. This gene set includes virtually all transcription factors that are expressed in a localized fashion by mid-gastrulation, and which thus provide spatial regulatory information to the embryo. Also included are signaling factors and some pigment cell differentiation genes. We show that the number of pregastrular D/N signaling targets among these regulatory genes is small and is almost exclusively restricted to non-skeletogenic mesoderm genes. However, Delta signaling also activates foxY in the small micromeres. As is the early NSM, the small micromeres are in direct contact with Delta expressing skeletogenic mesoderm. In contrast, no endoderm regulatory genes are activated by Delta signaling even during the second phase of delta expression, when this gene is transcribed in NSM cells adjacent to the endoderm. During this phase Delta provides an ongoing input which continues to activate foxY expression in small micromere progeny. Disruption of the second phase of Delta expression specifically abolishes specification of late mesodermal derivatives such as the coelomic pouches to which the small micromeres contribute.

Synthetic in vivo validation of gene network circuitry

Embryonic development is controlled by networks of interacting regulatory genes. The individual linkages of gene regulatory networks (GRNs) are customarily validated by functional cis-regulatory analysis, but an additional approach to validation is to rewire GRN circuitry to test experimentally predictions derived from network structure. Here we use this synthetic method to challenge specific predictions of the sea urchin embryo endomesoderm GRN. Expression vectors generated by in vitro recombination of exogenous sequences into BACs were used to cause elements of a nonskeletogenic mesoderm GRN to be deployed in skeletogenic cells and to detect their effects. The result of reengineering the regulatory circuitry in this way was to divert the developmental program of these cells from skeletogenesis to pigment cell formation, confirming a direct prediction of the GRN. In addition, the experiment revealed previously undetected cryptic repression functions.

Precise cis-regulatory control of spatial and temporal expression of the alx-1 gene in the skeletogenic lineage of s. purpuratus

Deployment of the gene-regulatory network (GRN) responsible for skeletogenesis in the embryo of the sea urchin Strongylocentrotus purpuratus is restricted to the large micromere lineage by a double negative regulatory gate. The gate consists of a GRN subcircuit composed of the pmar1 and hesC genes, which encode repressors and are wired in tandem, plus a set of target regulatory genes under hesC control. The skeletogenic cell state is specified initially by micromere-specific expression of these regulatory genes, viz. alx1, ets1, tbrain and tel, plus the gene encoding the Notch ligand Delta. Here we use a recently developed high throughput methodology for experimental cis-regulatory analysis to elucidate the genomic regulatory system controlling alx1 expression in time and embryonic space. The results entirely confirm the double negative gate control system at the cis-regulatory level, including definition of the functional HesC target sites, and add the crucial new information that the drivers of alx1 expression are initially Ets1, and then Alx1 itself plus Ets1. Cis-regulatory analysis demonstrates that these inputs quantitatively account for the magnitude of alx1 expression. Furthermore, the Alx1 gene product not only performs an auto-regulatory role, promoting a fast rise in alx1 expression, but also, when at high levels, it behaves as an auto-repressor. A synthetic experiment indicates that this behavior is probably due to dimerization. In summary, the results we report provide the sequence level basis for control of alx1 spatial expression by the double negative gate GRN architecture, and explain the rising, then falling temporal expression profile of the alx1 gene in terms of its auto-regulatory genetic wiring.

Information processing at the foxa node of the sea urchin endomesoderm specification network

The foxa regulatory gene is of central importance for endoderm specification across Bilateria, and this gene lies at an essential node of the well-characterized sea urchin endomesoderm gene regulatory network (GRN). Here we experimentally dissect the cis-regulatory system that controls the complex pattern of foxa expression in these embryos. Four separate cis-regulatory modules (CRMs) cooperate to control foxa expression in different spatial domains of the endomesoderm, and at different times. A detailed mutational analysis revealed the inputs to each of these cis-regulatory modules. The complex and dynamic expression of foxa is regulated by a combination of repressors, a permissive switch, and multiple activators. A mathematical kinetic model was applied to study the dynamic response of foxa cis-regulatory modules to transient inputs. This study shed light on the mesoderm-endoderm fate decision and provides a functional explanation, in terms of the genomic regulatory code, for the spatial and temporal expression of a key developmental control gene.

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.

Functional cis-regulatory genomics for systems biology

Gene expression is controlled by interactions between trans-regulatory factors and cis-regulatory DNA sequences, and these interactions constitute the essential functional linkages of gene regulatory networks (GRNs). Validation of GRN models requires experimental cis-regulatory tests of predicted linkages to authenticate their identities and proposed functions. However, cis-regulatory analysis is, at present, at a severe bottleneck in genomic system biology because of the demanding experimental methodologies currently in use for discovering cis-regulatory modules (CRMs), in the genome, and for measuring their activities. Here we demonstrate a high-throughput approach to both discovery and quantitative characterization of CRMs. The unique aspect is use of DNA sequence tags to "barcode" CRM expression constructs, which can then be mixed, injected together into sea urchin eggs, and subsequently deconvolved. This method has increased the rate of cis-regulatory analysis by >100-fold compared with conventional one-by-one reporter assays. The utility of the DNA-tag reporters was demonstrated by the rapid discovery of 81 active CRMs from 37 previously unexplored sea urchin genes. We then obtained simultaneous high-resolution temporal characterization of the regulatory activities of more than 80 CRMs. On average 2-3 CRMs were discovered per gene. Comparison of endogenous gene expression profiles with those of the CRMs recovered from each gene showed that, for most cases, at least one CRM is active in each phase of endogenous expression, suggesting that CRM recovery was comprehensive. This approach will qualitatively alter the practice of GRN construction as well as validation, and will impact many additional areas of regulatory system biology.

Cis-regulatory analysis of the sea urchin pigment cell gene polyketide synthase

The Strongylocentrotus purpuratus polyketide synthase gene (SpPks) encodes an enzyme required for the biosynthesis of the larval pigment echinochrome. SpPks is expressed exclusively in pigment cells and their precursors starting at blastula stage. The 7th-9th cleavage Delta-Notch signaling, required for pigment cell development, positively regulates SpPks. In previous studies, the transcription factors glial cell missing (SpGcm), SpGatae and kruppel-like (SpKrl/z13) have been shown to positively regulate SpPks. To uncover the structure of the Gene Regulatory Network (GRN) regulating the specification and differentiation processes of pigment cells, we experimentally analyzed the putative SpPks cis-regulatory region. We established that the -1.5kb region is sufficient to recapitulate the correct spatial and temporal expression of SpPks. Predicted DNA-binding sites for SpGcm, SpGataE and SpKrl are located within this region. The mutagenesis of these DNA-binding sites indicated that SpGcm, SpGataE and SpKrl are direct positive regulators of SpPks. These results demonstrate that the sea urchin GRN for pigment cell development is quite shallow, which is typical of type I embryo development.

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.

Practical computational methods for regulatory genomics: a cisGRN-Lexicon and cisGRN-browser for gene regulatory networks

The CYRENE Project focuses on the study of cis-regulatory genomics and gene regulatory networks (GRN) and has three components: a cisGRN-Lexicon, a cisGRN-Browser, and the Virtual Sea Urchin software system. The project has been done in collaboration with Eric Davidson and is deeply inspired by his experimental work in genomic regulatory systems and gene regulatory networks. The current CYRENE cisGRN-Lexicon contains the regulatory architecture of 200 transcription factors encoding genes and 100 other regulatory genes in eight species: human, mouse, fruit fly, sea urchin, nematode, rat, chicken, and zebrafish, with higher priority on the first five species. The only regulatory genes included in the cisGRN-Lexicon (CYRENE genes) are those whose regulatory architecture is validated by what we call the Davidson Criterion: they contain functionally authenticated sites by site-specific mutagenesis, conducted in vivo, and followed by gene transfer and functional test. This is recognized as the most stringent experimental validation criterion to date for such a genomic regulatory architecture. The CYRENE cisGRN-Browser is a full genome browser tailored for cis-regulatory annotation and investigation. It began as a branch of the Celera Genome Browser (available as open source at http://sourceforge.net/projects/celeragb /) and has been transformed to a genome browser fully devoted to regulatory genomics. Its access paradigm for genomic data is zoom-to-the-DNA-base in real time. A more recent component of the CYRENE project is the Virtual Sea Urchin system (VSU), an interactive visualization tool that provides a four-dimensional (spatial and temporal) map of the gene regulatory networks of the sea urchin embryo.

Regulative recovery in the sea urchin embryo and the stabilizing role of fail-safe gene network wiring

Design features that ensure reproducible and invariant embryonic processes are major characteristics of current gene regulatory network models. New cis-regulatory studies on a gene regulatory network subcircuit activated early in the development of the sea urchin embryo reveal a sequence of encoded "fail-safe" regulatory devices. These ensure the maintenance of fate separation between skeletogenic and nonskeletogenic mesoderm lineages. An unexpected consequence of the network design revealed in the course of these experiments is that it enables the embryo to "recover" from regulatory interference that has catastrophic effects if this feature is disarmed. A reengineered regulatory system inserted into the embryo was used to prove how this system operates in vivo. Genomically encoded backup control circuitry thus provides the mechanism underlying a specific example of the regulative development for which the sea urchin embryo has long been famous.

Genomic control of patterning

The development of multicellular organisms involves the partitioning of the organism into territories of cells of specific structure and function. The information for spatial patterning processes is directly encoded in the genome. The genome determines its own usage depending on stage and position, by means of interactions that constitute gene regulatory networks (GRNs). The GRN driving endomesoderm development in sea urchin embryos illustrates different regulatory strategies by which developmental programs are initiated, orchestrated, stabilized or excluded to define the pattern of specified territories in the developing embryo.

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.

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.

Functional studies of regulatory genes in the sea urchin embryo

Sea urchin embryos are characterized by an extremely simple mode of development, rapid cleavage, high transparency, and well-defined cell lineage. Although they are not suitable for genetic studies, other approaches are successfully used to unravel mechanisms and molecules involved in cell fate specification and morphogenesis. Microinjection is the elective method to study gene function in sea urchin embryos. It is used to deliver precise amounts of DNA, RNA, oligonucleotides, peptides, or antibodies into the eggs or even into blastomeres. Here we describe microinjection as it is currently applied in our laboratory and show how it has been used in gene perturbation analyses and dissection of cis-regulatory DNA elements.

The evolution of hierarchical gene regulatory networks

Comparative developmental evidence indicates that reorganizations in developmental gene regulatory networks (GRNs) underlie evolutionary changes in animal morphology, including body plans. We argue here that the nature of the evolutionary alterations that arise from regulatory changes depends on the hierarchical position of the change within a GRN. This concept cannot be accomodated by microevolutionary nor macroevolutionary theory. It will soon be possible to investigate these ideas experimentally, by assessing the effects of GRN changes on morphological evolution.

Evolution of gene regulatory network architectures: examples of subcircuit conservation and plasticity between classes of echinoderms

Developmental gene regulatory networks (GRNs) explain how regulatory states are established in particular cells during development and how these states then determine the final form of the embryo. Evolutionary changes to the sequence of the genome will direct reorganization of GRN architectures, which in turn will lead to the alteration of developmental programs. A comparison of GRN architectures must consequently reveal the molecular basis for the evolution of developmental programs among different organisms. This review highlights some of the important findings that have emerged from the most extensive direct comparison of GRN architectures to date. Comparison of the orthologous GRNs for endomesodermal specification in the sea urchin and sea star, provides examples of several discrete, functional GRN subcircuits and shows that they are subject to diverse selective pressures. This demonstrates that different regulatory linkages may be more or less amenable to evolutionary change. One of the more surprising findings from this comparison is that GRN-level functions may be maintained while the factors performing the functions have changed, suggesting that GRNs have a high capacity for compensatory changes involving transcription factor binding to cis regulatory modules.

Gene regulatory network subcircuit controlling a dynamic spatial pattern of signaling in the sea urchin embryo

We dissect the transcriptional regulatory relationships coordinating the dynamic expression patterns of two signaling genes, wnt8 and delta, which are central to specification of the sea urchin embryo endomesoderm. cis-Regulatory analysis shows that transcription of the gene encoding the Notch ligand Delta is activated by the widely expressed Runx transcription factor, but spatially restricted by HesC-mediated repression through a site in the delta 5'UTR. Spatial transcription of the hesC gene, however, is controlled by Blimp1 repression. Blimp1 thus represses the repressor of delta, thereby permitting its transcription. The blimp1 gene is itself linked into a feedback circuit that includes the wnt8 signaling ligand gene, and we showed earlier that this circuit generates an expanding torus of blimp1 and wnt8 expression. The finding that delta expression is also controlled at the cis-regulatory level by the blimp1-wnt8 torus-generating subcircuit now explains the progression of Notch signaling from the mesoderm to the endoderm of the developing embryo. Thus the specific cis-regulatory linkages of the gene regulatory network encode the coordinated spatial expression of Wnt and Notch signaling as they sweep outward across the vegetal plate of the embryo.

A protocol describing the principles of cis-regulatory analysis in the sea urchin

cis-Regulatory analysis (CRA) is the precise identification of the cis-acting genomic sequences regulating gene transcription. As such, CRA provides essential mechanistic insight into key biological processes such as development. The first phase of this protocol involves identification of a large (approximately 100 kb) clone of genomic material surrounding the gene of interest and use of this clone to establish a reliable and unambiguous reporter assay. In the second phase, phylogenetic footprinting is used to identify candidate regulatory modules; these genomic sequences are then recursively tested for reporter activity. In the final phase, potential transcription factor binding sites are identified and disrupted in reporter constructs for individual testing. The strengths of this method reflect the use of large clones containing all relevant genomic regulatory sequences to establish a reporter assay with high fidelity. Given these foundational elements, predicted or suspected regulatory inputs can be rigorously tested and novel regulatory inputs identified. Although the expected time line varies greatly with the depth of information required, results may be obtained in as little as 4-6 months, but more detailed analyses will require several years to complete.

The surprising complexity of the transcriptional regulation of the spdri gene reveals the existence of new linkages inside sea urchin's PMC and Oral Ectoderm Gene Regulatory Networks

During sea urchin embryogenesis the spdri gene participates in two separate Gene Regulatory Networks (GRNs): the Primary Mesenchyme Cells' (PMCs) and the Oral Ectoderm's one. In both cases, activation of the gene follows initial specification events [Amore, G., Yavrouian, R., Peterson, K., Ransick, A., McClay, D., Davidson, E., 2003. Spdeadringer, a sea urchin embryo gene required separately in skeletogenic and oral ectoderm gene regulatory networks. Dev. Biol. 261, 55-81.]. We identified a portion of genomic DNA ("4.7IL" -3456;+389) which is sufficient to replicate sdpri's expression pattern in experiments of transgenesis, using a GFP reporter. In our experiments, the activation kinetic of 4.7IL-GFP was similar to that of the endogenous gene and the reporter responded to known spdri's transcriptional regulators (Ets1, Alx1, Gsc and Dri). Here we present a dissection of this regulatory region and a description of the modules involved in spdri's transcriptional regulation. Both in the PMCs' and Oral Ectoderm's expression phases, activation of spdri is obtained through the integration of three kinds of inputs: positive and globally distributed ones; negative ones (that prevent ectopic expression); positive and tissue-specific ones. Our results allow to expand the map of the regulatory connections at the spdri node, both in the PMCs and in the Oral Ectoderm Gene Regulatory Networks (GRNs).

A new method, using cis-regulatory control, for blocking embryonic gene expression

Many genes, and particularly regulatory genes, are utilized multiple times in unrelated phases of development. For studies of gene function during embryogenesis, there is often need of a method for interfering with expression only at a specific developmental time or place. Here we show that in sea urchin embryos cis-regulatory control systems which operate only at specific times and places can be used to drive expression of short designed sequences targeting given primary transcripts, thereby effectively taking out the function of the target genes. The active sequences are designed to be complementary to intronic sequences of the primary transcript of the target genes. In this work, the target genes were the transcription factors alx1 and ets1, both required for skeletogenesis, and the regulatory drivers were from the sm30 and tbr genes. The sm30 gene is expressed only after skeletogenic cell ingression. When its regulatory apparatus was used as driver, the alx1 and ets1 repression constructs had the effect of preventing postgastrular skeletogenesis, while not interfering with earlier alx1 and ets1 function in promoting skeletogenic mesenchyme ingression. In contrast, repression constructs using the tbr driver, which is active in blastula stage, block ingression. This method thus provides the opportunity to study regulatory requirements of skeletogenesis after ingression, and may be similarly useful in many other developmental contexts.

Evolutionary plasticity of developmental gene regulatory network architecture

Sea stars and sea urchins evolved from a last common ancestor that lived at the end of the Cambrian, approximately half a billion years ago. In a previous comparative study of the gene regulatory networks (GRNs) that embody the genomic program for embryogenesis in these animals, we discovered an almost perfectly conserved five-gene network subcircuit required for endoderm specification. We show here that the GRN structure upstream and downstream of the conserved network kernel has, by contrast, diverged extensively. Mesoderm specification is accomplished quite differently; the Delta-Notch signaling system is used in radically distinct ways; and various regulatory genes have been coopted to different functions. The conservation of the conserved kernel is thus the more remarkable. The results indicate types of network linkage subject to evolutionary change. An emergent theme is that subcircuit design may be preserved even while the identity of genes performing given roles changes because of alteration in their cis-regulatory control systems.

A spatially dynamic cohort of regulatory genes in the endomesodermal gene network of the sea urchin embryo

A gene regulatory network subcircuit comprising the otx, wnt8, and blimp1 genes accounts for a moving torus of gene expression that sweeps concentrically across the vegetal domain of the sea urchin embryo. Here we confirm by mutation the inputs into the blimp1cis-regulatory module predicted by network analysis. Its essential design feature is that it includes both activation and autorepression sites. The wnt8 gene is functionally linked into the subcircuit in that cells receiving this ligand generate a beta-catenin/Tcf input required for blimp1 expression, while the wnt8 gene in turn requires a Blimp1 input. Their torus-like spatial expression patterns and gene regulatory analysis indicate that the genes even-skipped and hox11/13b are also entrained by this subcircuit. We verify the cis-regulatory inputs of even-skipped predicted by network analysis. These include activation by beta-catenin/Tcf and Blimp1, repression within the torus by Hox11/13b, and repression outside the torus by Tcf in the absence of Wnt8 signal input. Thus even-skipped and hox11/13b, along with blimp1 and wnt8, are members of a cohort of torus genes with similar regulatory inputs and similar, though slightly out-of-phase, expression patterns, which reflect differences in cis-regulatory design.

Analysis of cis-regulatory elements controlling spatio-temporal expression of T-brain gene in sea urchin, Hemicentrotus pulcherrimus

In sea urchin development, micromere descendants play important roles in skeletogenesis and induction of gastrulation. We previously reported that the T-brain homolog of sea urchin Hemicentrotus pulcherrimus, HpTb expresses specifically in micromere descendants and is required for induction of gastrulation and skeletogenesis. Thus, HpTb is thought to play important roles in the function of micromere-lineage cells. To identify cis-regulatory regions responsible for spatio-temporal gene expression of HpTb, we isolated approximately 7kb genomic region of HpTb gene and showed that GFP expression driven by this region exhibits the spatio-temporal pattern corresponding substantially to that of endogenous HpTb expression. Deletion of interspecifically conserved C2 and C4 regions resulted in an increase of ectopic expression. Mutations in Hairy family and Snail family consensus sequences in C1 and C2 regions also increased ectopic expression. Furthermore, we demonstrated that C4 region functions as enhancer, and that three Ets family consensus sequences are involved in this activity but not in spatial regulation. Therefore, we concluded that expression of HpTb gene is regulated by multiple cis-regulatory elements.

Cis-regulatory analysis of nodal and maternal control of dorsal-ventral axis formation by Univin, a TGF-β related to Vg1

The TGF-β family member Nodal is essential for specification of the dorsal-ventral axis of the sea urchin embryo, but the molecular factors regulating its expression are not known. Analysis of the nodalpromoter is an excellent entry point to identify these factors and to dissect the regulatory logic driving dorsal-ventral axis specification. Using phylogenetic footprinting, we delineated two regulatory regions located in the 5′ region of the nodal promoter and in the intron that are required for correct spatial expression and for autoregulation. The 5′regulatory region contains essential binding sites for homeodomain, bZIP, Oct,Tcf/Lef, Sox and Smad transcription factors, and a binding site for an unidentified spatial repressor possibly related to Myb. Soon after its initiation, nodal expression critically requires autoregulation by Nodal and signaling by the maternal TGF-β Univin. We show that Univin is related to Vg1, that both Nodal and Univin signal through Alk4/5/7, and that zygotic expression of univin, like that of nodal, is dependent on SoxB1 function and Tcf/β-catenin signaling. This work shows that Tcf, SoxB1 and Univin play essential roles in the regulation of nodal expression in the sea urchin and suggests that some of the regulatory interactions controlling nodal expression predate the chordates. The data are consistent with a model of nodal regulation in which a maternal TGF-β acts in synergy with maternal transcription factors and with spatial repressors to establish the dorsal-ventral axis of the sea urchin embryo.

Caught in the evolutionary act: precise cis-regulatory basis of difference in the organization of gene networks of sea stars and sea urchins

The regulatory control of otxbeta1/2 in the sea urchin Strongylocentrotus purpuratus and the sea star Asterina miniata provides an exceptional opportunity to determine the genomic basis of evolutionary change in gene regulatory network (GRN) architectures. Network perturbation analyses in both taxa show that Otx regulates the transcription factors gatae and krox/blimp1 and both of these transcription factors also feed back and regulate otx. The otx gene also autoregulates. This three way interaction is an example of a GRN kernel. It has been conserved for 500 million years since these two taxa last shared a common ancestor. Amid this high level of conservation we show here one significant regulatory change. Tbrain is required for correct otxbeta1/2 expression in the sea star but not in the sea urchin. In sea urchin, tbrain is not co-expressed with otxbeta1/2 and instead has an essential role in specification of the embryonic skeleton. Tbrain in these echinoderms is thus a perfect example of an orthologous gene co-opted for entirely different developmental processes. We isolate and test the sea star otxbeta1/2 cis-regulatory module and demonstrate functional binding sites for each of the predicted inputs, including Tbrain. We compare it to the logic processing operating in the sea urchin otxbeta1/2 cis-regulatory module and present an evolutionary scenario of the change in Tbrain dependence. Finally, inter-specific gene transfer experiments confirm this scenario and demonstrate evolution occurring at the level of sequence changes to the cis-regulatory module.

Gene regulation: gene control network in development

Controlling the differential expression of many thousands of genes is the most fundamental task of a developing organism. It requires an enormous computational device that has the capacity to process in parallel a vast number of regulatory inputs in the various cells of the embryo and come out with regulatory outputs that are tissue specific. The regulatory genome constitutes this computational device, comprising many thousands of processing units in the form of cis-regulatory modules. The interconnected cis-regulatory modules that control regulatory gene expression create a network that is the underlying mechanism of specification. In this review we use the gene regulatory network that governs endomesoderm specification in the sea urchin embryo to demonstrate the salient features of developmental gene regulatory networks and illustrate the information processing that is done by the regulatory sequences.

Progressive activation of Delta-Notch signaling from around the blastopore is required to set up a functional caudal lobe in the spiderAchaearanea tepidariorum

In the development of most arthropods, the caudal region of the elongating germ band (the growth zone) sequentially produces new segments. Previous work with the spider Cupiennius salei suggested involvement of Delta-Notch signaling in segmentation. Here, we report that, in the spider Achaearanea tepidariorum, the same signaling pathway exerts a different function in the presumptive caudal region before initiation of segmentation. In the developing spider embryo, the growth zone becomes morphologically apparent as a caudal lobe around the closed blastopore. We found that, preceding caudal lobe formation, transcripts of a Delta homolog, At-Delta,are expressed in evenly spaced cells in a small area covering the closing blastopore and then in a progressively wider area of the germ disc epithelium. Cells with high At-Delta expression are likely to be prospective mesoderm cells, which later express a twist homolog, At-twist, and individually internalize. Cells remaining at the surface begin to express a caudal homolog, At-caudal, to differentiate as caudal ectoderm. Knockdown of At-Delta by parental RNA interference results in overproduction of At-twist-expressing mesoderm cells at the expense of At-caudal-expressing ectoderm cells. This condition gives rise to a disorganized caudal region that fails to pattern the opisthosoma. In addition, knockdown of Notch and Suppressor of Hairless homologs produces similar phenotypes. We suggest that, in the spider, progressive activation of Delta-Notch signaling from around the blastopore leads to stochastic cell fate decisions between mesoderm and caudal ectoderm through a process of lateral inhibition to set up a functional caudal lobe.

Cis-regulatory control of the nodal gene, initiator of the sea urchin oral ectoderm gene network

Expression of the nodal gene initiates the gene regulatory network which establishes the transcriptional specification of the oral ectoderm in the sea urchin embryo. This gene encodes a TGFbeta ligand, and in Strongylocentrotus purpuratus its transcription is activated in the presumptive oral ectoderm at about the 30-cell stage. Thereafter Nodal signaling occurs among all cells of the oral ectoderm territory, and nodal expression is required for expression of oral ectoderm regulatory genes. The cis-regulatory system of the nodal gene transduces anisotropically distributed cytoplasmic cues that distinguish the future oral and aboral domains of the early embryo. Here we establish the genomic basis for the initiation and maintenance of nodal gene expression in the oral ectoderm. Functional cis-regulatory control modules of the nodal gene were identified by interspecific sequence conservation. A 5' cis-regulatory module functions both to initiate expression of the nodal gene and to maintain its expression by means of feedback input from the Nodal signal transduction system. These functions are mediated respectively by target sites for bZIP transcription factors, and by SMAD target sites. At least one SMAD site is also needed for the initiation of expression. An intron module also contains SMAD sites which respond to Nodal feedback, and in addition acts to repress vegetal expression. These observations explain the main features of nodal expression in the oral ectoderm: since the activity of bZIP factors is redox sensitive, and the initial polarization of oral vs. aboral fate is manifested in a redox differential, the bZIP sites account for the activation of nodal on the oral side; and since the immediate early signal transduction response factors for Nodal are SMAD factors, the SMAD sites account for the feedback maintenance of nodal gene expression.

A genomic view of the sea urchin nervous system

The sequencing of the Strongylocentrotus purpuratus genome provides a unique opportunity to investigate the function and evolution of neural genes. The neurobiology of sea urchins is of particular interest because they have a close phylogenetic relationship with chordates, yet a distinctive pentaradiate body plan and unusual neural organization. Orthologues of transcription factors that regulate neurogenesis in other animals have been identified and several are expressed in neurogenic domains before gastrulation indicating that they may operate near the top of a conserved neural gene regulatory network. A family of genes encoding voltage-gated ion channels is present but, surprisingly, genes encoding gap junction proteins (connexins and pannexins) appear to be absent. Genes required for synapse formation and function have been identified and genes for synthesis and transport of neurotransmitters are present. There is a large family of G-protein-coupled receptors, including 874 rhodopsin-type receptors, 28 metabotropic glutamate-like receptors and a remarkably expanded group of 161 secretin receptor-like proteins. Absence of cannabinoid, lysophospholipid and melanocortin receptors indicates that this group may be unique to chordates. There are at least 37 putative G-protein-coupled peptide receptors and precursors for several neuropeptides and peptide hormones have been identified, including SALMFamides, NGFFFamide, a vasotocin-like peptide, glycoprotein hormones and insulin/insulin-like growth factors. Identification of a neurotrophin-like gene and Trk receptor in sea urchin indicates that this neural signaling system is not unique to chordates. Several hundred chemoreceptor genes have been predicted using several approaches, a number similar to that for other animals. Intriguingly, genes encoding homologues of rhodopsin, Pax6 and several other key mammalian retinal transcription factors are expressed in tube feet, suggesting tube feet function as photosensory organs. Analysis of the sea urchin genome presents a unique perspective on the evolutionary history of deuterostome nervous systems and reveals new approaches to investigate the development and neurobiology of sea urchins.

Deciphering the Underlying Mechanism of Specification and Differentiation: The Sea Urchin Gene Regulatory Network

The regulatory genome functions as a vast information processor through development. It processes the initial conditions that are set by asymmetric distributions of cellular components in the egg and translates them into the onset of spatially localized specification states. It regulates the timely differential activation of signaling molecules and transcription factors that divide the emerging domains into subdomains. It also governs the activation of groups of differentiation genes, the genes that encode, at the protein level, the functional and the structural properties of a cell type. The sea urchin endomesoderm gene regulatory network provides a window into the different levels of the regulatory apparatus. It demonstrates how the static physical genomic components define functional connections between the various regulatory genes that act together to conduct the dynamical developmental program.