Exclusive developmental functions of gatae cis-regulatory modules in the Strongylocentrorus purpuratus embryo

Integrative multi-omics increase resolution of the sea urchin posterior gut gene regulatory network at single-cell level

Drafting gene regulatory networks (GRNs) requires embryological knowledge pertaining to the cell type families, information on the regulatory genes, causal data from gene knockdown experiments and validations of the identified interactions by cis-regulatory analysis. We use multi-omics involving next-generation sequencing to obtain the necessary information for drafting the Strongylocentrotus purpuratus (Sp) posterior gut GRN. Here, we present an update to the GRN using: (1) a single-cell RNA-sequencing-derived cell atlas highlighting the 2 day-post-fertilization (dpf) sea urchin gastrula cell type families, as well as the genes expressed at the single-cell level; (2) a set of putative cis-regulatory modules and transcription factor-binding sites obtained from chromatin accessibility ATAC-seq data; and (3) interactions directionality obtained from differential bulk RNA sequencing following knockdown of the transcription factor Sp-Pdx1, a key regulator of gut patterning in sea urchins. Combining these datasets, we draft the GRN for the hindgut Sp-Pdx1-positive cells in the 2 dpf gastrula embryo. Overall, our data suggest the complex connectivity of the posterior gut GRN and increase the resolution of gene regulatory cascades operating within it.

Evolutionarily conserved Wnt/Sp5 signaling is critical for anterior-posterior axis patterning in sea urchin embryos

Studies across a diverse group of metazoan embryos indicate that Wnt signaling often activates the transcription factor Sp5, forming a signaling 'cassette' that plays critical roles in many developmental processes. This study explores the role of Wnt/Sp5 signaling during the specification and patterning of the primary germ layers during early anterior-posterior axis formation in the deuterostome sea urchin embryo. Our functional analyses show that Sp5 is critical for endomesoderm specification downstream of Wnt/β-catenin in posterior cells as well as anterior neuroectoderm patterning downstream of non-canonical Wnt/JNK signaling in anterior cells. Interestingly, expression and functional data comparisons show that Wnt/Sp5 signaling often plays similar roles in posterior endomesoderm as well as neuroectoderm patterning along the AP axis of several deuterostome embryos, including vertebrates. Thus, our findings provide strong support for the idea that Wnt-Sp5 signaling cassettes were critical for the establishment of early germ layers in the common deuterostome ancestor.

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.

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.

Bacterial artificial chromosomes as recombinant reporter constructs to investigate gene expression and regulation in echinoderms

Genome sequences contain all the necessary information-both coding and regulatory sequences-to construct an organism. The developmental process translates this genomic information into a three-dimensional form. One interpretation of this translation process can be described using gene regulatory network (GRN) models, which are maps of interactions among regulatory gene products in time and space. As high throughput investigations reveal increasing complexity within these GRNs, it becomes apparent that efficient methods are required to test the necessity and sufficiency of regulatory interactions. One of the most complete GRNs for early development has been described in the purple sea urchin, Strongylocentrotus purpuratus. This work has been facilitated by two resources: a well-annotated genome sequence and transgenes generated in bacterial artificial chromosome (BAC) constructs. BAC libraries played a central role in assembling the S. purpuratus genome sequence and continue to serve as platforms for generating reporter constructs for use in expression and regulatory analyses. Optically transparent echinoderm larvae are highly amenable to transgenic approaches and are therefore particularly well suited for experiments that rely on BAC-based reporter transgenes. Here, we discuss the experimental utility of BAC constructs in the context of understanding developmental processes in echinoderm embryos and larvae.

Developmental effector gene regulation: Multiplexed strategies for functional analysis

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

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.

Capsaicin Enhances Glutamatergic Synaptic Transmission to Neonatal Rat Hypoglossal Motor Neurons via a TRPV1-Independent Mechanism

We investigated whether capsaicin modulated synaptic transmission to hypoglossal motor neurons (HMNs) by acting on transient receptor potential vanilloid type 1 (TRPV1) receptors. Using whole-cell patch clamp recording from neonatal rat HMNs, we found that capsaicin increased spontaneous excitatory post-synaptic current (sEPSC) frequency and amplitude. Interestingly, the only effect of capsaicin on spontaneous inhibitory post-synaptic currents (sIPSCs) was a significant decrease in sIPSC amplitude without altering frequency, indicating a post-synaptic mechanism of action. The frequency of miniature excitatory post-synaptic currents (mEPSCs), recorded in the presence of tetrodotoxin (TTX), was also increased by capsaicin, but capsaicin did not alter mEPSC amplitude, consistent with a pre-synaptic mechanism of action. A negative shift in membrane current (I_holding) was elicited by capsaicin under both recording conditions. The effect of capsaicin on excitatory synaptic transmission remained unchanged in the presence of the TRPV1 antagonists, capsazepine or SB366791, suggesting that capsaicin acts to modulate EPSCs via a mechanism which does not require TRPV1 activation. Capsaicin, however, did not alter evoked excitatory post-synaptic currents (eEPSCs) or the paired-pulse ratio (PPR) of eEPSCs. Repetitive action potential (AP) firing in HMNs was also unaltered by capsaicin, indicating that capsaicin does not change HMN intrinsic excitability. We have demonstrated that capsaicin modulates glutamatergic excitatory, as well as glycinergic inhibitory, synaptic transmission in HMNs by differing pre- and post-synaptic mechanisms. These results expand our understanding regarding the extent to which capsaicin can modulate synaptic transmission to central neurons.

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.

Logics and properties of a genetic regulatory program that drives embryonic muscle development in an echinoderm

Evolutionary origin of muscle is a central question when discussing mesoderm evolution. Developmental mechanisms underlying somatic muscle development have mostly been studied in vertebrates and fly where multiple signals and hierarchic genetic regulatory cascades selectively specify myoblasts from a pool of naive mesodermal progenitors. However, due to the increased organismic complexity and distant phylogenetic position of the two systems, a general mechanistic understanding of myogenesis is still lacking. In this study, we propose a gene regulatory network (GRN) model that promotes myogenesis in the sea urchin embryo, an early branching deuterostome. A fibroblast growth factor signaling and four Forkhead transcription factors consist the central part of our model and appear to orchestrate the myogenic process. The topological properties of the network reveal dense gene interwiring and a multilevel transcriptional regulation of conserved and novel myogenic genes. Finally, the comparison of the myogenic network architecture among different animal groups highlights the evolutionary plasticity of developmental GRNs.

DOI:http://dx.doi.org/10.7554/eLife.07343.001

Muscles, bones, and blood vessels all develop from a tissue called the mesoderm, which forms early on in the development of an embryo. Networks of genes control which parts of the mesoderm transform into different cell types. The gene networks that control the development of muscle cells from the mesoderm have so far been investigated in flies and several species of animals with backbones. However, these species are complex, which makes it difficult to work out the general principles that control muscle cell development.

Sea urchins are often studied in developmental biology as they have many of the same genes as more complex animals, but are much simpler and easier to study in the laboratory. Andrikou et al. therefore investigated the ‘gene regulatory network’ that controls muscle development in sea urchins. This revealed that proteins called Forkhead transcription factors and a process called FGF signaling are crucial for controlling muscle development in sea urchins. These are also important factors for developing muscles in other animals.

Andrikou et al. then produced models that show the interactions between the genes that control muscle formation at three different stages of embryonic development. These models reveal several important features of the muscle development gene regulatory network. For example, the network is robust: if one gene fails, the network is connected in a way that allows it to still make muscle. This also allows the network to adapt and evolve without losing the ability to perform any of its existing roles.

Comparing the gene regulatory network that controls muscle development in sea urchins with the networks found in other animals showed that many of the same genes are used across different species, but are connected into different network structures. Investigating the similarities and differences of the regulatory networks in different species could help us to understand how muscles have evolved and could ultimately lead to a better understanding of the causes of developmental diseases.

DOI:http://dx.doi.org/10.7554/eLife.07343.002

Autoregulatory feedback controls sequential action of cis-regulatory modules at the brinker locus

cis-regulatory modules (CRMs) act sequentially to regulate temporal expression of genes, but how the switch from one to the next is accomplished is not well understood. To provide insight, here we investigate the cis-regulatory system controlling brinker (brk) expression in the Drosophila embryo. Two distally located CRMs support expression at different times, while a promoter-proximal element (PPE) is required to support their action. In the absence of Brk protein itself or upon mutagenesis of Brk binding sites within the PPE, the late-acting CRM, specifically, is delayed. This block to late-acting CRM function appears to be removed when the early-acting CRM is also deleted. These results demonstrate that autoregulatory feedback is necessary for the early-acting CRM to disengage from the promoter so that the late-acting CRM may act. Autoregulation may be a commonly used mechanism to control sequential CRM action necessary for dynamic gene expression throughout the course of development.

Employing BAC-reporter constructs in the sea anemone Nematostella vectensis

Changes in the expression and function of genes drive evolutionary change. Comparing how genes are regulated in different species is therefore becoming an important part of evo-devo studies. A key tool for investigating the regulation of genes is represented by bacterial artificial chromosomes (BAC)-reporter constructs. BACs are large insert libraries, often >100 kb, which thus capture the genomic sequences surrounding a gene of interest, including all, or nearly all, of the elements underpinning regulation. Recombinant BACs, containing a reporter gene in place of the endogenous coding sequence of genes, can be utilized to drive the expression of reporter genes under the regulatory control of the gene of interest while still embedded within its genomic context. Systematic deletions within the BAC-reporter construct can be used to identify the minimal reporter in an unbiased way, avoiding the risk of overlooking regulatory elements that may be many kilobases away from the transcription start-site. Nematostella vectensis (Edwardsiidae, Anthozoa, Cnidaria) has become an important model in regenerative biology, ecology, and especially in studies of evo-devo and gene-regulatory networks due to its interesting phylogenetic position and amenability to molecular techniques. The increasing interest in this rising model system also led to a demand for methods that can be used to study the regulation of genes in Nematostella. Here, we present our progress in employing BAC-reporter constructs to visualize gene-expression in Nematostella. Using a new Nematostella-specific recombination cassette, we made nine different BAC-reporter constructs. Although five BAC recombinants gave variable effects, three constructs, namely Nv-bra:eGFP::L10 BAC, Nv-dpp:eGFP::L10 BAC, and Nv-grm:eGFP::L10 BAC delivered promising results. We show that these three constructs express the reporter gene eGFP in 10.4-17.2% of all analyzed larvae, out of which 26.2-41.9% express GFP in a mosaic fashion within the expected domain. In addition to the expression within the known domains, we also observed cases of misexpression of eGFP and examples that could represent actual expression outside the described domain. Furthermore, we deep-sequenced and assembled five different BACs containing Nv-chordin, Nv-foxa, Nv-dpp, Nv-wnta, and Nv-wnt1, to improve assembly around these genes. The use of BAC-reporter constructs will foster cis-regulatory analyses in Nematostella and thus help to improve our understanding of the regulatory network in this cnidarian system. Ultimately, this will advance the comparison of gene-regulation across species and lead to a much better understanding of evolutionary changes and novelties.

Mesodermal gene expression in the acoel Isodiametra pulchra indicates a low number of mesodermal cell types and the endomesodermal origin of the gonads

Acoelomorphs are bilaterally symmetric small marine worms that lack a coelom and possess a digestive system with a single opening. Two alternative phylogenetic positions of this group within the animal tree are currently debated. In one view, Acoelomorpha is the sister group to all remaining Bilateria and as such, is a morphologically simple stepping stone in bilaterian evolution. In the other, the group is a lineage within the Deuterostomia, and therefore, has derived a simple morphology from a more complex ancestor. Acoels and the closely related Nemertodermatida and Xenoturbellida, which together form the Acoelomorpha, possess a very limited number of cell types. To further investigate the diversity and origin of mesodermal cell types we describe the expression pattern of 12 orthologs of bilaterian mesodermal markers including Six1/2, Twist, FoxC, GATA4/5/6, in the acoel Isodiametra pulchra. All the genes are expressed in stem cells (neoblasts), gonads, and at least subsets of the acoel musculature. Most are expressed in endomesodermal compartments of I. pulchra developing embryos similar to what has been described in cnidarians. Our molecular evidence indicates a very limited number of mesodermal cell types and suggests an endomesodermal origin of the gonads and the stem cell system. We discuss our results in light of the two prevailing phylogenetic positions of Acoelomorpha.

Diversification of oral and aboral mesodermal regulatory states in pregastrular sea urchin embryos

Specification of the non-skeletogenic mesoderm (NSM) in sea urchin embryos depends on Delta signaling. Signal reception leads to expression of regulatory genes that later contribute to the aboral NSM regulatory state. In oral NSM, this is replaced by a distinct oral regulatory state in consequence of Nodal signaling. Through regulome wide analysis we identify the homeobox gene not as an immediate Nodal target. not expression in NSM causes extinction of the aboral regulatory state in the oral NSM, and expression of a new suite of regulatory genes. All NSM specific regulatory genes are henceforth expressed exclusively, in oral or aboral domains, presaging the mesodermal cell types that will emerge. We have analyzed the regulatory linkages within the aboral NSM gene regulatory network. A linchpin of this network is gataE which as we show is a direct Gcm target and part of a feedback loop locking down the aboral regulatory state.

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.

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.

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.

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.

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.

cis-Regulatory sequences driving the expression of the Hbox12 homeobox-containing gene in the presumptive aboral ectoderm territory of the Paracentrotus lividus sea urchin embryo

Embryonic development is coordinated by networks of evolutionary conserved regulatory genes encoding transcription factors and components of cell signalling pathways. In the sea urchin embryo, a number of genes encoding transcription factors display territorial restricted expression. Among these, the zygotic Hbox12 homeobox gene is transiently transcribed in a limited number of cells of the animal-lateral half of the early Paracentrotus lividus embryo, whose descendants will constitute part of the ectoderm territory. To obtain insights on the regulation of Hbox12 expression, we have explored the cis-regulatory apparatus of the gene. In this paper, we show that the intergenic region of the tandem Hbox12 repeats drives GFP expression in the presumptive aboral ectoderm and that a 234 bp fragment, defined aboral ectoderm (AE) module, accounts for the restricted expression of the transgene. Within this module, a consensus sequence for a Sox factor and the binding of the Otx activator are both required for correct Hbox12 gene expression. Spatial restriction to the aboral ectoderm is achieved by a combination of different repressive sequence elements. Negative sequence elements necessary for repression in the endomesoderm map within the most upstream 60 bp region and nearby the Sox binding site. Strikingly, a Myb-like consensus is necessary for repression in the oral ectoderm, while down-regulation at the gastrula stage depends on a GA-rich region. These results suggest a role for Hbox12 in aboral ectoderm specification and represent our first attempt in the identification of the gene regulatory circuits involved in this process.

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.