Patchy interspecific sequence similarities efficiently identify positive cis-regulatory elements in the sea urchin

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.

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.

Methods for the experimental and computational analysis of gene regulatory networks in sea urchins

The discovery of gene regulatory networks (GRNs) has opened a gate to access the genomic mechanisms controlling development. GRNs are systems of transcriptional regulatory circuits that control the differential specification of cell fates during development by regulating gene expression. The experimental analysis of GRNs involves a collection of methods, each revealing aspects of the overall control process. This review provides an overview of experimental and computational methods that have been successfully applied for solving developmental GRNs in the sea urchin embryo. The key in this approach is to obtain experimental evidence for functional interactions between transcription factors and regulatory DNA. In the second part of this review, a more generally applicable strategy is discussed that shows a path from experimental evidence to annotation of regulatory linkages to the generation of GRN models.

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.

Echinoderm development and evolution in the post-genomic era

The highly recognizable animals within the phylum Echinodermata encompass an enormous disparity of adult and larval body plans. The extensive knowledge of sea urchin development has culminated in the description of the exquisitely detailed gene regulatory network (GRN) that governs the specification of various embryonic territories. This information provides a unique opportunity for comparative studies in other echinoderm taxa to understand the evolution and developmental mechanisms underlying body plan change. This review focuses on recent work that has utilized new genomic resources and systems-level experiments to address questions of evolutionary developmental biology. In particular, we synthesize the results of several recent studies from various echinoderm classes that have explored the development and evolution of the larval skeleton, which is a major feature that distinguishes the two predominant larval subtypes within the Phylum. We specifically examine the ways in which GRNs can evolve, either through cis regulatory and/or protein-level changes in transcription factors. We also examine recent work comparing evolution across shorter time scales that occur within and between species of sea urchin, and highlight the kinds of questions that can be addressed by these comparisons. The advent of new genomic and transcriptomic datasets in additional species from all classes of echinoderm will continue to empower the use of these taxa for evolutionary developmental studies.

Unraveling the Tangled Skein: The Evolution of Transcriptional Regulatory Networks in Development

The molecular and genetic basis for the evolution of anatomical diversity is a major question that has inspired evolutionary and developmental biologists for decades. Because morphology takes form during development, a true comprehension of how anatomical structures evolve requires an understanding of the evolutionary events that alter developmental genetic programs. Vast gene regulatory networks (GRNs) that connect transcription factors to their target regulatory sequences control gene expression in time and space and therefore determine the tissue-specific genetic programs that shape morphological structures. In recent years, many new examples have greatly advanced our understanding of the genetic alterations that modify GRNs to generate newly evolved morphologies. Here, we review several aspects of GRN evolution, including their deep preservation, their mechanisms of alteration, and how they originate to generate novel developmental programs.

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.

Pattern and process during sea urchin gut morphogenesis: the regulatory landscape

The development of the endoderm is a multistage process. From the initial specification of the endodermal domain in the embryo to the final regionalization of the gut, there are multiple stages that require the involvement of complex gene regulatory networks. In one concrete case, the sea urchin embryo, some of these stages and their genetic control are (relatively) well understood. Several studies have underscored the relevance of individual transcription factor activities in the process, but very few have focused the attention on gene interactions within specific gene regulatory networks (GRNs). Sea urchins offer an ideal system to study the different factors involved in the morphogenesis of the gut. Here we review the knowledge gained over the last 10 years on the process and its regulation, from the early specification of endodermal lineages to the late events linked to the patterning of functional domains in the gut. A lesson of remarkable importance has been learnt from comparison of the mechanisms involved in gut formation in different bilaterian animals; some of these genetic mechanisms are particularly well conserved. Patterning the gut seems to involve common molecular players and shared interactions, whether we look at mammals or echinoderms. This astounding degree of conservation reveals some key aspects of deep homology that are most probably shared by all bilaterian guts.

Tools for sea urchin genomic analysis

The Sea Urchin Genome Project Web site, SpBase ( http://SpBase.org ), in association with a suite of publicly available sequence comparison tools provides a platform from which to analyze genes and genomic sequences of sea urchin. This information system is specifically designed to support laboratory bench studies in cell and molecular biology. In particular these tools and datasets have supported the description of the gene regulatory networks of the purple sea urchin S. purpuratus. This chapter details methods to undertake in the first steps to find genes and noncoding regulatory sequences for further analysis.

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.

Pantropic retroviruses as a transduction tool for sea urchin embryos

Sea urchins are an important model for experiments at the intersection of development and systems biology, and technical innovations that enhance the utility of this model are of great value. This study explores pantropic retroviruses as a transduction tool for sea urchin embryos, and demonstrates that pantropic retroviruses infect sea urchin embryos with high efficiency and genomically integrate at a copy number of one per cell. We successfully used a self-inactivation strategy to both insert a sea urchin-specific enhancer and disrupt the endogenous viral enhancer. The resulting self-inactivating viruses drive global and persistent gene expression, consistent with genomic integration during the first cell cycle. Together, these data provide substantial proof of principle for transduction technology in sea urchin embryos.

Nucleosome exclusion from the interspecies-conserved central AT-rich region of the Ars insulator

The Ars insulator is a boundary element identified in the upstream region of the arylsulfatase (HpArs) gene in the sea urchin, Hemicentrotus pulcherrimus, and possesses the ability to both block enhancer-promoter communications and protect transgenes from silent chromatin. To understand the molecular mechanism of the Ars insulator, we investigated the correlation between chromatin structure, DNA structure and insulator activity. Nuclease digestion of nuclei isolated from sea urchin embryos revealed the presence of a nuclease-hypersensitive site within the Ars insulator. Analysis of micrococcal nuclease-sensitive sites in the Ars insulator, reconstituted with nucleosomes, showed the exclusion of nucleosomes from the central AT-rich region. Furthermore, the central AT-rich region in naked DNA was sensitive to nucleotide base modification by diethylpyrocarbonate (DEPC). These observations suggest that non-B-DNA structures in the central AT-rich region may inhibit nucleosomal formation, which leads to nuclease hypersensitivity. Furthermore, comparison of nucleotide sequences between the HpArs gene and its ortholog in Strongylocentrotus purpuratus revealed that the central AT-rich region of the Ars insulator is conserved, and this conserved region showed significant enhancer blocking activity. These results suggest that the central AT-rich nucleosome-free region plays an important role in the function of the Ars insulator.

In silico characterization of the neural alpha tubulin gene promoter of the sea urchin embryo Paracentrotus lividus by phylogenetic footprinting

During Paracentrotus lividus sea urchin embryo development one alpha and one beta tubulin genes are expressed specifically in the neural cells and they are early end output of the gene regulatory network that specifies the neural commitment. In this paper we have used a comparative genomics approach to identify conserved regulatory elements in the P. lividus neural alpha tubulin gene. To this purpose, we have first isolated a genomic clone containing the entire gene plus 4.5 Kb of 5' upstream sequences. Then, we have shown by gene transfer experiments that its non-coding region drives the spatio-temporal gene expression corresponding substantially to that of the endogenous gene. In addition, we have identified by genome and EST sequence analysis the S. purpuratus alpha tubulin orthologous gene and we propose a revised annotation of some tubulin family members. Moreover, by computational techniques we delineate at least three putative regulatory regions located both in the upstream region and in the first intron containing putative binding sites for Forkhead and Nkx transcription factor families.

Uncoupling of complex regulatory patterning during evolution of larval development in echinoderms

Background

Conservation of orthologous regulatory gene expression domains, especially along the neuroectodermal anterior-posterior axis, in animals as disparate as flies and vertebrates suggests that common patterning mechanisms have been conserved since the base of Bilateria. The homology of axial patterning is far less clear for the many marine animals that undergo a radical transformation in body plan during metamorphosis. The embryos of these animals are microscopic, feeding within the plankton until they metamorphose into their adult forms.

Results

We describe here the localization of 14 transcription factors within the ectoderm during early embryogenesis in Patiria miniata, a sea star with an indirectly developing planktonic bipinnaria larva. We find that the animal-vegetal axis of this very simple embryo is surprisingly well patterned. Furthermore, the patterning that we observe throughout the ectoderm generally corresponds to that of "head/anterior brain" patterning known for hemichordates and vertebrates, which share a common ancestor with the sea star. While we suggest here that aspects of head/anterior brain patterning are generally conserved, we show that another suite of genes involved in retinal determination is absent from the ectoderm of these echinoderms and instead operates within the mesoderm.

Conclusions

Our findings therefore extend, for the first time, evidence of a conserved axial pattering to echinoderm embryos exhibiting maximal indirect development. The dissociation of head/anterior brain patterning from "retinal specification" in echinoderm blastulae might reflect modular changes to a developmental gene regulatory network within the ectoderm that facilitates the evolution of these microscopic larvae.

Gene expression analysis of Six3, Pax6, and Otx in the early development of the stalked crinoid Metacrinus rotundus

The stalked crinoid, Metacrinus rotundus, is one of the most basal extant echinoderms. Here, we show the expression patterns of Six3, Pax6, and Otx in the early development of M. rotundus. All three genes are highly expressed in stages from the gastrula to the auricularia larval stage. Ectodermal expression of MrOtx appears to be correlated with development of the ciliary band. These three genes are expressed sequentially along the embryonic body axis in the anterior and middle walls of the archenteron in the order of MrPax6, MrSix3, and MrOtx. The anterior, middle, and posterior parts of the archenteron in the late gastrula differentiate into the axo-hydrocoel, the enteric sac, and somatocoels at later stages, respectively. The three genes are expressed sequentially from the tip of the axo-hydrocoel to the bottom of enteric sac in the order of MrSix3, MrPax6, and MrOtx at the later stages. This suggests that these genes are involved in patterning of the larval endo-mesoderm in stalked crinoids. The present results suggest that radical alterations have occurred in the expression and function of homeobox genes in basal echinoderms.

Characterization and expression of a sea star otx ortholog (Protxβ1/2) in the larva of Patiriella regularis

A transcript of otx from the sea star Patiriella regularis (Protxβ1/2) was characterized and its expression in early bipinnaria larvae was documented by whole mount in situ hybridization (WMISH). The nucleotide sequence exhibited 94% identity with Amotxβ1/2 from the closely related species Patiria miniata. Protxβ1/2 was expressed strongly in the developing archenteron in the future fore and mid-gut regions. This was followed by expression of Protxβ1/2 in the developing enterocoels, mesodermal derivatives. This suggests a role for Protx in endomesoderm development. In coelom development, Protxβ1/2 was first expressed in the left coelom. Subsequently expression was evident in the right coelom, but localization was never as strong as in the left coelom. This asymmetry in Protxβ1/2 expression in the coeloms was evident up to the stage when they started to extend posteriorly. These data indicate that Protxβ1/2 may have a role in coelom development, particularly in the left coelom, a definitive adult structure.

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.

Parameters for accurate genome alignment

Background

Genome sequence alignments form the basis of much research. Genome alignment depends on various mundane but critical choices, such as how to mask repeats and which score parameters to use. Surprisingly, there has been no large-scale assessment of these choices using real genomic data. Moreover, rigorous procedures to control the rate of spurious alignment have not been employed.

Results

We have assessed 495 combinations of score parameters for alignment of animal, plant, and fungal genomes. As our gold-standard of accuracy, we used genome alignments implied by multiple alignments of proteins and of structural RNAs. We found the HOXD scoring schemes underlying alignments in the UCSC genome database to be far from optimal, and suggest better parameters. Higher values of the X-drop parameter are not always better. E-values accurately indicate the rate of spurious alignment, but only if tandem repeats are masked in a non-standard way. Finally, we show that γ-centroid (probabilistic) alignment can find highly reliable subsets of aligned bases.

Conclusions

These results enable more accurate genome alignment, with reliability measures for local alignments and for individual aligned bases. This study was made possible by our new software, LAST, which can align vertebrate genomes in a few hours http://last.cbrc.jp/.

The endoderm gene regulatory network in sea urchin embryos up to mid-blastula stage

As the result of early specification processes, sea urchin embryos eventually form various mesodermal cell lineages and a gut consisting of fore-, mid- and hindgut. The progression of specification as well as the overall spatial organization of the organism is encoded in its gene regulatory networks (GRNs). We have analyzed the GRN driving endoderm specification up to the onset of gastrulation and present in this paper the mechanisms which determine this process up to mid-blastula stage. At this stage, the embryo consists of two separate lineages of endoderm precursor cells with distinct regulatory states. One of these lineages, the veg2 cell lineage, gives rise to endoderm and mesoderm cell types. The separation of these cell fates is initiated by the spatially confined activation of the mesoderm GRN superimposed on a generally activated endoderm GRN within veg2 descendants. Here we integrate the architecture of regulatory interactions with the spatial restriction of regulatory gene expression to model the logic control of endoderm development.

Building developmental gene regulatory networks

Animal development is an elaborate process programmed by genomic regulatory instructions. Regulatory genes encode transcription factors and signal molecules, and their expression is under the control of cis-regulatory modules that define the logic of transcriptional responses to the inputs of other regulatory genes. The functional linkages among regulatory genes constitute the gene regulatory networks (GRNs) that govern cell specification and patterning in development. Constructing such networks requires identification of the regulatory genes involved and characterization of their temporal and spatial expression patterns. Interactions (activation/repression) among transcription factors or signals can be investigated by large-scale perturbation analysis, in which the function of each gene is specifically blocked. Resultant expression changes are then integrated to identify direct linkages, and to reveal the structure of the GRN. Predicted GRN linkages can be tested and verified by cis-regulatory analysis. The explanatory power of the GRN was shown in the lineage specification of sea urchin endomesoderm. Acquiring such networks is essential for a systematic and mechanistic understanding of the developmental process.

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.

A perturbation model of the gene regulatory network for oral and aboral ectoderm specification in the sea urchin embryo

The current gene regulatory network (GRN) for the sea urchin embryo pertains to pregastrular specification functions in the endomesodermal territories. Here we extend gene regulatory network analysis to the adjacent oral and aboral ectoderm territories over the same period. A large fraction of the regulatory genes predicted by the sea urchin genome project and shown in ancillary studies to be expressed in either oral or aboral ectoderm by 24 h are included, though universally expressed and pan-ectodermal regulatory genes are in general not. The loci of expression of these genes have been determined by whole mount in situ hybridization. We have carried out a global perturbation analysis in which expression of each gene was interrupted by introduction of morpholino antisense oligonucleotide, and the effects on all other genes were measured quantitatively, both by QPCR and by a new instrumental technology (NanoString Technologies nCounter Analysis System). At its current stage the network model, built in BioTapestry, includes 22 genes encoding transcription factors, 4 genes encoding known signaling ligands, and 3 genes that are yet unknown but are predicted to perform specific roles. Evidence emerged from the analysis pointing to distinctive subcircuit features observed earlier in other parts of the GRN, including a double negative transcriptional regulatory gate, and dynamic state lockdowns by feedback interactions. While much of the regulatory apparatus is downstream of Nodal signaling, as expected from previous observations, there are also cohorts of independently activated oral and aboral ectoderm regulatory genes, and we predict yet unidentified signaling interactions between oral and aboral territories.

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.

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.

Dissection of lin-11 enhancer regions in Caenorhabditis elegans and other nematodes

The Caenorhabditis elegans LIM homeobox gene lin-11 plays crucial roles in the morphogenesis of the reproductive system and differentiation of several neurons. The expression of lin-11 in different tissues is regulated by enhancer regions located upstream as well as within lin-11 introns. These regions are functionally separable suggesting that multiple regulatory inputs operate to control the spatiotemporal pattern of lin-11 expression. To further dissect apart the nature of lin-11 regulation we focused on three Caenorhabditis species C. briggsae, C. remanei, and C. brenneri that are substantially diverged from C. elegans but share almost identical vulval morphology. We show that, in these species, the 5' region of lin-11 possesses conserved sequences to activate lin-11 expression in the reproductive system. Analysis of the in vivo role of these sequences in C. elegans has led to the identification of three functionally distinct enhancers for the vulva, VC neurons, and uterine pi lineage cells. We found that the pi enhancer is regulated by FOS homolog FOS-1 and LIN-12/Notch pathway effectors, LAG-1 (Su(H)/CBF1 family) and EGL-43 (EVI1 family). These results indicate that multiple factors cooperate to regulate pi-specific expression of lin-11 and together with other findings suggest that the mechanism of lin-11 regulation by LIN-12/Notch signaling is evolutionarily conserved in Caenorhabditis species. Our work demonstrates that 4-way comparison is a powerful tool to study conserved mechanisms of gene regulation in C. elegans and other nematodes.

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

ModuleMiner - improved computational detection of cis-regulatory modules: are there different modes of gene regulation in embryonic development and adult tissues?

ModuleMiner detects cis-regulatory modules in a set of co-expressed genes in tissue-specific microarray clusters and in embryonic development datasets.

We present ModuleMiner, a novel algorithm for computationally detecting cis-regulatory modules (CRMs) in a set of co-expressed genes. ModuleMiner outperforms other methods for CRM detection on benchmark data, and successfully detects CRMs in tissue-specific microarray clusters and in embryonic development gene sets. Interestingly, CRM predictions for differentiated tissues exhibit strong enrichment close to the transcription start site, whereas CRM predictions for embryonic development gene sets are depleted in this region.

Rapid detection and curation of conserved DNA via enhanced-BLAT and EvoPrinterHD analysis

Background

Multi-genome comparative analysis has yielded important insights into the molecular details of gene regulation. We have developed EvoPrinter, a web-accessed genomics tool that provides a single uninterrupted view of conserved sequences as they appear in a species of interest. An EvoPrint reveals with near base-pair resolution those sequences that are essential for gene function.

Results

We describe here EvoPrinterHD, a 2^nd-generation comparative genomics tool that automatically generates from a single input sequence an enhanced view of sequence conservation between evolutionarily distant species. Currently available for 5 nematode, 3 mosquito, 12 Drosophila, 20 vertebrate, 17 Staphylococcus and 20 enteric bacteria genomes, EvoPrinterHD employs a modified BLAT algorithm [enhanced-BLAT (eBLAT)], which detects up to 75% more conserved bases than identified by the BLAT alignments used in the earlier EvoPrinter program. The new program also identifies conserved sequences within rearranged DNA, highlights repetitive DNA, and detects sequencing gaps. EvoPrinterHD currently holds over 112 billion bp of indexed genomes in memory and has the flexibility of selecting a subset of genomes for analysis. An EvoDifferences profile is also generated to portray conserved sequences that are uniquely lost in any one of the orthologs. Finally, EvoPrinterHD incorporates options that allow for (1) re-initiation of the analysis using a different genome's aligning region as the reference DNA to detect species-specific changes in less-conserved regions, (2) rapid extraction and curation of conserved sequences, and (3) for bacteria, identifies unique or uniquely shared sequences present in subsets of genomes.

Conclusion

EvoPrinterHD is a fast, high-resolution comparative genomics tool that automatically generates an uninterrupted species-centric view of sequence conservation and enables the discovery of conserved sequences within rearranged DNA. When combined with cis-Decoder, a program that discovers sequence elements shared among tissue specific enhancers, EvoPrinterHD facilitates the analysis of conserved sequences that are essential for coordinate gene regulation.

Collier transcription in a single Drosophila muscle lineage: the combinatorial control of muscle identity

Specification of muscle identity in Drosophila is a multistep process: early positional information defines competence groups termed promuscular clusters, from which muscle progenitors are selected, followed by asymmetric division of progenitors into muscle founder cells (FCs). Each FC seeds the formation of an individual muscle with morphological and functional properties that have been proposed to reflect the combination of transcription factors expressed by its founder. However, it is still unclear how early patterning and muscle-specific differentiation are linked. We addressed this question, using Collier (Col; also known as Knot) expression as both a determinant and read-out of DA3 muscle identity. Characterization of the col upstream region driving DA3 muscle specific expression revealed the existence of three separate phases of cis-regulation, correlating with conserved binding sites for different mesodermal transcription factors. Examination of col transcription in col and nautilus (nau) loss-of-function and gain-of-function conditions showed that both factors are required for col activation in the ;naïve' myoblasts that fuse with the DA3 FC, thereby ensuring that all DA3 myofibre nuclei express the same identity programme. Together, these results indicate that separate sets of cis-regulatory elements control the expression of identity factors in muscle progenitors and myofibre nuclei and directly support the concept of combinatorial control of muscle identity.

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.

cis-Decoder discovers constellations of conserved DNA sequences shared among tissue-specific enhancers

: The use of cis-Decoder, a new tool for discovery of conserved sequence elements that are shared between similarly regulating enhancers, suggests that enhancers use overlapping repertoires of highly conserved core elements.

A systematic approach is described for analysis of evolutionarily conserved cis-regulatory DNA using cis-Decoder, a tool for discovery of conserved sequence elements that are shared between similarly regulated enhancers. Analysis of 2,086 conserved sequence blocks (CSBs), identified from 135 characterized enhancers, reveals most CSBs consist of shorter overlapping/adjacent elements that are either enhancer type-specific or common to enhancers with divergent regulatory behaviors. Our findings suggest that enhancers employ overlapping repertoires of highly conserved core elements.

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.

The regulatory genome and the computer

The definitive feature of the many thousand cis-regulatory control modules in an animal genome is their information processing capability. These modules are "wired" together in large networks that control major processes such as development; they constitute "genomic computers." Each control module receives multiple inputs in the form of the incident transcription factors which bind to them. The functions they execute upon these inputs can be reduced to basic AND, OR and NOT logic functions, which are also the unit logic functions of electronic computers. Here we consider the operating principles of the genomic computer, the product of evolution, in comparison to those of electronic computers. For example, in the genomic computer intra-machine communication occurs by means of diffusion (of transcription factors), while in electronic computers it occurs by electron transit along pre-organized wires. There follow fundamental differences in design principle in respect to the meaning of time, speed, multiplicity of processors, memory, robustness of computation and hardware and software. The genomic computer controls spatial gene expression in the development of the body plan, and its appearance in remote evolutionary time must be considered to have been a founding requirement for animal grade life.

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.

Repression of mesodermal fate by foxa, a key endoderm regulator of the sea urchin embryo

The foxa gene is an integral component of the endoderm specification subcircuit of the endomesoderm gene regulatory network in the Strongylocentrotus purpuratus embryo. Its transcripts become confined to veg2, then veg1 endodermal territories, and, following gastrulation, throughout the gut. It is also expressed in the stomodeal ectoderm. gatae and otx genes provide input into the pregastrular regulatory system of foxa, and Foxa represses its own transcription, resulting in an oscillatory temporal expression profile. Here, we report three separate essential functions of the foxa gene: it represses mesodermal fate in the veg2 endomesoderm; it is required in postgastrular development for the expression of gut-specific genes; and it is necessary for stomodaeum formation. If its expression is reduced by a morpholino, more endomesoderm cells become pigment and other mesenchymal cell types, less gut is specified, and the larva has no mouth. Experiments in which blastomere transplantation is combined with foxa MASO treatment demonstrate that, in the normal endoderm, a crucial role of Foxa is to repress gcm expression in response to a Notch signal, and hence to repress mesodermal fate. Chimeric recombination experiments in which veg2, veg1 or ectoderm cells contained foxa MASO show which region of foxa expression controls each of the three functions. These experiments show that the foxa gene is a component of three distinct embryonic gene regulatory networks.

High regulatory gene use in sea urchin embryogenesis: Implications for bilaterian development and evolution

A global scan of transcription factor usage in the sea urchin embryo was carried out in the context of the Strongylocentrotus purpuratus genome sequencing project, and results from six individual studies are here considered. Transcript prevalence data were obtained for over 280 regulatory genes encoding sequence-specific transcription factors of every known family, but excluding genes encoding zinc finger proteins. This is a statistically inclusive proxy for the total "regulome" of the sea urchin genome. Close to 80% of the regulome is expressed at significant levels by the late gastrula stage. Most regulatory genes must be used repeatedly for different functions as development progresses. An evolutionary implication is that animal complexity at the stage when the regulome first evolved was far simpler than even the last common bilaterian ancestor, and is thus of deep antiquity.

The Transcriptome of the Sea Urchin Embryo

The sea urchin Strongylocentrotus purpuratus is a model organism for study of the genomic control circuitry underlying embryonic development. We examined the complete repertoire of genes expressed in the S. purpuratus embryo, up to late gastrula stage, by means of high-resolution custom tiling arrays covering the whole genome. We detected complete spliced structures even for genes known to be expressed at low levels in only a few cells. At least 11,000 to 12,000 genes are used in embryogenesis. These include most of the genes encoding transcription factors and signaling proteins, as well as some classes of general cytoskeletal and metabolic proteins, but only a minor fraction of genes encoding immune functions and sensory receptors. Thousands of small asymmetric transcripts of unknown function were also detected in intergenic regions throughout the genome. The tiling array data were used to correct and authenticate several thousand gene models during the genome annotation process.

Endomesoderm specification in Caenorhabditis elegans and other nematodes

The endomesoderm gene regulatory network (GRN) of C. elegans is a rich resource for studying the properties of cell-fate-specification pathways. This GRN contains both cell-autonomous and cell non-autonomous mechanisms, includes network motifs found in other GRNs, and ties maternal factors to terminal differentiation genes through a regulatory cascade. In most cases, upstream regulators and their direct downstream targets are known. With the availability of resources to study close and distant relatives of C. elegans, the molecular evolution of this network can now be examined. Within Caenorhabditis, components of the endomesoderm GRN are well conserved. A cursory examination of the preliminary genome sequences of two parasitic nematodes, Haemonchus contortus and Brugia malayi, suggests that evolution in this GRN is occurring most rapidly for the zygotic genes that specify blastomere identity.

Expression and function of blimp1/krox, an alternatively transcribed regulatory gene of the sea urchin endomesoderm network

The blimp1/krox gene of Strongylocentrotus purpuratus, formerly krox1, encodes zinc finger transcription factors which play a central role in both early and late endomesoderm specification. Here we show that there are two alternative splice forms transcribed under the control of different regulatory regions. The blimp/krox1b form was previously unknown, and is the form expressed during cleavage, beginning 6-9 h postfertilization. This form is required for the early events of endomesoderm specification. A different splice variant, blimp1/krox1a, is expressed only from gastrula stage onward. During cleavage stages the blimp1/krox gene is expressed in the large micromeres and veg2 descendents. Soon after, it is expressed in the ring of specified mesoderm cells at the vegetal pole of the blastula. Its expression is later restricted to the blastopore region and the posterior of the invaginating archenteron, and finally to the midgut and hindgut of the pluteus larva. The expression of blimp1/krox is dynamic, and involves several distinct spatial territories. A GFP recombinant BAC was created by substituting the GFP coding sequence for that of the second exon (1b), in order to distinguish the expression pattern of the early form from that of the late form. This construct closely mimics blimp1/krox1b expression during early stages of sea urchin development. To expand our knowledge of the downstream linkages of this gene, additional experiments were carried out using antisense morpholino oligos (MASO). We confirmed previously published data that blimp1/krox autoregulates its own expression, but discovered, surprisingly, that this gene represses rather than activates itself. This negative autoregulation is restricted to the mesodermal and probably skeletogenic territories during the blastula stage, as shown by in situ hybridization analysis of MASO injected embryos. The MASO perturbation analysis also revealed blimp1/krox inputs into other genes of the endomesoderm regulatory network.

Regulation of spblimp1/krox1a, an alternatively transcribed isoform expressed in midgut and hindgut of the sea urchin gastrula

The sea urchin regulatory gene Spblimp1/krox produces alternatively transcribed and spliced isoforms, 1a and 1b, which have different temporal and spatial patterns of expression. Here we describe a cis-regulatory module that controls the expression of the 1a splice form in the midgut and hindgut at the beginning of gastrulation. Conserved sequence patches revealed by a comparison of the blimp1/krox locus in Strongylocentrotus purpuratus and Lytechinus variegatus genomes were tested by gene transfer, in association with GFP or CAT reporter genes. An expression construct containing a conserved sequence patch immediately 5' of exon 1a included the transcription initiation site for blimp1/krox1a. This construct displays specific mid and hindgut expression, indicating that these are the locations of endogenous blimp1/krox1a transcription during the gastrula stage. Its sequence contains binding sites for Brn1/2/4, Otx, and Blimp1/Krox itself, as predicted in a prior regulatory network analysis.

cis-Regulatory control of cyclophilin, a member of the ETS-DRI skeletogenic gene battery in the sea urchin embryo

The Strongylocentrotus purpuratus cyclophilin1 gene (Sp-cyp1) is expressed exclusively in skeletogenic mesenchyme cells of the embryo, beginning in the micromere lineage of the early blastula stage and continuing after gastrulation during the syncytial deposition of the skeleton. This gene encodes a protein which is a member of the peptidyl-prolyl cis-trans isomerase (PPIase) family. Sp-cyp1 is among the differentiation genes activated in the skeletogenic territory as a terminal function of the endomesodermal gene regulatory network. Network perturbation analysis had predicted the skeletogenic regulators Ets1 and Deadringer (Dri) to be its driver inputs. Here, we show that a 218-bp cis-regulatory DNA fragment recapitulates skeletogenic Sp-cyp1 expression; that elimination of either Ets1 or Dri inputs severely depresses the activity of expression constructs containing this DNA fragment; and that Ets1 and Dri target sites within the 218 bp fragment are required for normal expression. This indicates that the predicted inputs are direct. Other studies indicate that the same inputs are evidently necessary for expression of several other skeletogenic differentiation genes, and these genes probably constitute a skeletogenic gene battery, defined by its Ets plus Dri regulatory inputs.

Selective constraint on noncoding regions of hominid genomes

An important challenge for human evolutionary biology is to understand the genetic basis of human–chimpanzee differences. One influential idea holds that such differences depend, to a large extent, on adaptive changes in gene expression. An important step in assessing this hypothesis involves gaining a better understanding of selective constraint on noncoding regions of hominid genomes. In noncoding sequence, functional elements are frequently small and can be separated by large nonfunctional regions. For this reason, constraint in hominid genomes is likely to be patchy. Here we use conservation in more distantly related mammals and amniotes as a way of identifying small sequence windows that are likely to be functional. We find that putatively functional noncoding elements defined in this manner are subject to significant selective constraint in hominids.

A major goal of human evolutionary biology is to understand what genetic changes make humans unique. One influential idea is that changes in gene expression are most responsible for unique human characteristics. Regulatory elements in noncoding DNA play a key role in controlling gene expression, so one approach is to study human–chimpanzee differences in these elements. Here we use conservation in more distantly related mammals and amniotes as a way of identifying small sequence windows that are likely to be functional. We find that putatively functional noncoding elements defined in this manner are subject to significant selective constraint in hominids. Contrary to some previous reports, these results argue that hominid noncoding regions are not evolving free of constraint.

Med-type GATA factors and the evolution of mesendoderm specification in nematodes

In the nematode, C. elegans, the divergent GATA-type transcription factors MED-1 and MED-2 are encoded by an unlinked, redundant pair of intronless genes. The med-1,2 genes are among the first to be activated in the embryo and are critical for the specification of the 7-cell stage MS (mesoderm) and E (endoderm) precursor cells. We have previously shown that the binding site recognized by MED-1 is a noncanonical RAGTATAC site that is not expected from the resemblance of its single C4-type zinc finger to those of other known GATA factors, which recognize the consensus HGATAR. To date, no MED-like zinc fingers have been described outside of C. elegans. In order to understand the evolution of these transcription factors, and the evolution of gene networks that specify early cell fates in Caenorhabditis, we have identified med sequence homologs in the related nematodes C. briggsae and C. remanei. While C. briggsae encodes two med-like genes similar to C. elegans, we find evidence for seven distinct med-like genes in C. remanei. Somewhat unexpectedly, the coding regions of all med genes appear to lack introns. We report that the med homologs have similar expression in their respective species. We further show that the C. briggsae homologs, and at least five of the seven C. remanei homologs, can fully complement the embryonic lethal phenotype of a C. elegans med-1,2(-) strain. We conclude that Med function and expression have been conserved over tens of millions of years of evolution, and that there may be a mechanism that selects against the acquisition of introns in these genes.

cis-Regulatory inputs of the wnt8 gene in the sea urchin endomesoderm network

Expression of the wnt8 gene is the key transcriptional motivator of an intercellular signaling loop which drives endomesoderm specification forward early in sea urchin embryogenesis. This gene was predicted by network perturbation analysis to be activated by inputs from the blimp1/krox gene, itself expressed zygotically in the endomesoderm during cleavage; and by a Tcf1/beta-catenin input. The implication is that zygotic expression of wnt8 is stimulated in neighboring cells by its own gene product, since reception of the Wnt8 ligand causes beta-catenin nuclearization. Here, the modular cis-regulatory system of the wnt8 gene of Strongylocentrotus purpuratus was characterized functionally, and shown to respond to blockade of both Blimp1/Krox and Tcf1/beta-catenin inputs just as does the endogenous gene. The genomic target sites for these factors were demonstrated by mutation in one of the cis-regulatory modules. The Tcf1/beta-catenin and Blimp1/Krox inputs are both necessary for normal endomesodermal expression mediated by this cis-regulatory module; thus, the genomic regulatory code underlying the predicted signaling loop thus resides in the wnt8 cis-regulatory sequence. In a second regulatory region, which initiates expression in micromere and macromere descendant cells early in cleavage, Tcf1 sites act to repress ectopic transcription in prospective ectoderm cells.