The gene regulatory network basis of the "community effect," and analysis of a sea urchin embryo example

An RNA interference approach for functional studies in the sea urchin and its use in analysis of nodal signaling gradients

Dicer substrate interfering RNAs (DsiRNAs) destroy targeted transcripts using the RNA-Induced Silencing Complex (RISC) through a process called RNA interference (RNAi). This process is ubiquitous among eukaryotes. Here we report the utility of DsiRNA in embryos of the sea urchin Lytechinus variegatus (Lv). Specific knockdowns phenocopy known morpholino and inhibitor knockdowns, and DsiRNA offers a useful alternative to morpholinos. Methods are described for the design of specific DsiRNAs that lead to destruction of targeted mRNA. DsiRNAs directed against pks1, an enzyme necessary for pigment production, show how successful DsiRNA perturbations are monitored by RNA in situ analysis and by qPCR to determine relative destruction of targeted mRNA. DsiRNA-based knockdowns phenocopy morpholino- and drug-based inhibition of nodal and lefty. Other knockdowns demonstrate that the RISC operates early in development as well as on genes that are first transcribed hours after gastrulation is completed. Thus, DsiRNAs effectively mediate destruction of targeted mRNA in the sea urchin embryo. The approach offers significant advantages over other widely used methods in the urchin in terms of cost, and ease of procurement, and offers sizeable experimental advantages in terms of ease of handling, injection, and knockdown validation.

Combinatorial regulatory states define cell fate diversity during embryogenesis

Cell fate specification occurs along invariant species-specific trajectories that define the animal body plan. This process is controlled by gene regulatory networks that regulate the expression of the limited set of transcription factors encoded in animal genomes. Here we globally assess the spatial expression of ~90% of expressed transcription factors during sea urchin development from embryo to larva to determine the activity of gene regulatory networks and their regulatory states during cell fate specification. We show that >200 embryonically expressed transcription factors together define >70 cell fates that recapitulate the morphological and functional organization of this organism. Most cell fate-specific regulatory states consist of ~15-40 transcription factors with similarity particularly among functionally related cell types regardless of developmental origin. Temporally, regulatory states change continuously during development, indicating that progressive changes in regulatory circuit activity determine cell fate specification. We conclude that the combinatorial expression of transcription factors provides molecular definitions that suffice for the unique specification of cell states in time and space during embryogenesis.

Population-level antagonism between FGF and BMP signaling steers mesoderm differentiation in embryonic stem cells

The mesodermal precursor populations for different internal organ systems are specified during gastrulation by the combined activity of extracellular signaling systems such as BMP, Wnt, Nodal and FGF. The BMP, Wnt and Nodal signaling requirements for the differentiation of specific mesoderm subtypes in mammals have been mapped in detail, but how FGF shapes mesodermal cell type diversity is not precisely known. It is also not clear how FGF signaling integrates with the activity of other signaling systems involved in mesoderm differentiation. Here, we address these questions by analyzing the effects of targeted signaling manipulations in differentiating stem cell populations at single-cell resolution. We identify opposing functions of BMP and FGF, and map FGF-dependent and -independent mesodermal lineages. Stimulation with exogenous FGF boosts the expression of endogenous Fgf genes while repressing Bmp ligand genes. This positive autoregulation of FGF signaling, coupled with the repression of BMP signaling, may contribute to the specification of reproducible and coherent cohorts of cells with the same identity via a community effect, both in the embryo and in synthetic embryo-like systems.

Development of a larval nervous system in the sea urchin

This review reports recent findings on the specification and patterning of neurons that establish the larval nervous system of the sea urchin embryo. Neurons originate in three regions of the embryo. Perturbation analyses enabled construction of gene regulatory networks controlling the several neural cell types. Many of the mechanisms described reflect shared features of all metazoans and others are conserved among deuterostomes. This nervous system with a very small number of neurons supports the feeding and swimming behaviors of the larva until metamorphosis when an adult nervous system replaces that system.

Segmentation and patterning of the vertebrate hindbrain

During early development, the hindbrain is sub-divided into rhombomeres that underlie the organisation of neurons and adjacent craniofacial tissues. A gene regulatory network of signals and transcription factors establish and pattern segments with a distinct anteroposterior identity. Initially, the borders of segmental gene expression are imprecise, but then become sharply defined, and specialised boundary cells form. In this Review, we summarise key aspects of the conserved regulatory cascade that underlies the formation of hindbrain segments. We describe how the pattern is sharpened and stabilised through the dynamic regulation of cell identity, acting in parallel with cell segregation. Finally, we discuss evidence that boundary cells have roles in local patterning, and act as a site of neurogenesis within the hindbrain.

Deciphering and modelling the TGF-β signalling interplays specifying the dorsal-ventral axis of the sea urchin embryo

During sea urchin development, secretion of Nodal and BMP2/4 ligands and their antagonists Lefty and Chordin from a ventral organiser region specifies the ventral and dorsal territories. This process relies on a complex interplay between the Nodal and BMP pathways through numerous regulatory circuits. To decipher the interplay between these pathways, we used a combination of treatments with recombinant Nodal and BMP2/4 proteins and a computational modelling approach. We assembled a logical model focusing on cell responses to signalling inputs along the dorsal-ventral axis, which was extended to cover ligand diffusion and enable multicellular simulations. Our model simulations accurately recapitulate gene expression in wild-type embryos, accounting for the specification of ventral ectoderm, ciliary band and dorsal ectoderm. Our model simulations further recapitulate various morphant phenotypes, reveal a dominance of the BMP pathway over the Nodal pathway and stress the crucial impact of the rate of Smad activation in dorsal-ventral patterning. These results emphasise the key role of the mutual antagonism between the Nodal and BMP2/4 pathways in driving early dorsal-ventral patterning of the sea urchin embryo.

Self-organization of stem cells into embryos: A window on early mammalian development

Embryonic development is orchestrated by robust and complex regulatory mechanisms acting at different scales of organization. In vivo studies are particularly challenging for mammals after implantation, owing to the small size and inaccessibility of the embryo. The generation of stem cell models of the embryo represents a powerful system with which to dissect this complexity. Control of geometry, modulation of the physical environment, and priming with chemical signals reveal the intrinsic capacity of embryonic stem cells to make patterns. Adding the stem cells for the extraembryonic lineages generates three-dimensional models that are more autonomous from the environment and recapitulate many features of the pre- and postimplantation mouse embryo, including gastrulation. Here, we review the principles of self-organization and how they set cells in motion to create an embryo.

Causal Gene Regulatory Network Modeling and Genomics: Second-Generation Challenges

Gene regulatory network modeling has played a major role in advancing the understanding of developmental systems, by crystallizing structures of relevant extant information, by formally posing hypothetical functional relationships between network elements, and by offering clear predictive tests to improve understanding of the mechanisms driving developmental progression. Both ordinary differential equation (ODE)-based and Boolean models have also been highly successful in explaining dynamics within subcircuits of more complex processes. In a very small number of cases, gene regulatory network models of much more global scope have been proposed that successfully predict the dynamics of the processes establishing most of an embryonic body plan. Can such successes be expanded to very different developmental systems, including post-embryonic mammalian systems? This perspective discusses several problems that must be solved in more quantitative and predictive theoretical terms, to make this possible. These problems include: the effects of cellular history on chromatin state and how these affect gene accessibility; the dose dependence of activities of many transcription factors (a problem for Boolean models); stochasticity of some transcriptional outputs (a problem for simple ODE models); response timing delays due to epigenetic remodeling requirements; functionally different kinds of repression; and the regulatory syntax that governs responses of genes with multiple enhancers.

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.

Assessing regulatory information in developmental gene regulatory networks

Gene regulatory networks (GRNs) provide a transformation function between the static genomic sequence and the primary spatial specification processes operating development. The regulatory information encompassed in developmental GRNs thus goes far beyond the control of individual genes. We here address regulatory information at different levels of network organization, from single node to subcircuit to large-scale GRNs and discuss how regulatory design features such as network architecture, hierarchical organization, and cis-regulatory logic contribute to the developmental function of network circuits. Using specific subcircuits from the sea urchin endomesoderm GRN, for which both circuit design and biological function have been described, we evaluate by Boolean modeling and in silico perturbations the import of given circuit features on developmental function. The examples include subcircuits encoding positive feedback, mutual repression, and coherent feedforward, as well as signaling interaction circuitry. Within the hierarchy of the endomesoderm GRN, these subcircuits are organized in an intertwined and overlapping manner. Thus, we begin to see how regulatory information encoded at individual nodes is integrated at all levels of network organization to control developmental process.

Cell Identity Switching Regulated by Retinoic Acid Signaling Maintains Homogeneous Segments in the Hindbrain

The patterning of tissues to form subdivisions with distinct and homogeneous regional identity is potentially disrupted by cell intermingling. Transplantation studies suggest that homogeneous segmental identity in the hindbrain is maintained by identity switching of cells that intermingle into another segment. We show that switching occurs during normal development and is mediated by feedback between segment identity and the retinoic acid degrading enzymes, cyp26b1 and cyp26c1. egr2, which specifies the segmental identity of rhombomeres r3 and r5, underlies the lower expression level of cyp26b1 and cyp26c1 in r3 and r5 compared with r2, r4, and r6. Consequently, r3 or r5 cells that intermingle into adjacent segments encounter cells with higher cyp26b1/c1 expression, which we find is required for downregulation of egr2b expression. Furthermore, egr2b expression is regulated in r2, r4, and r6 by non-autonomous mechanisms that depend upon the number of neighbors that express egr2b. These findings reveal that a community regulation of retinoid signaling maintains homogeneous segmental identity.

Modeling Mammalian Gastrulation With Embryonic Stem Cells

Understanding cell fate patterning and morphogenesis in the mammalian embryo remains a formidable challenge. Recently, in vivo models based on embryonic stem cells (ESCs) have emerged as complementary methods to quantitatively dissect the physical and molecular processes that shape the embryo. Here we review recent developments in using ESCs to create both two- and three-dimensional culture models that shed light on mammalian gastrulation.

Morphogen and community effects determine cell fates in response to BMP4 signaling in human embryonic stem cells

Paracrine signals maintain developmental states and create cell fate patterns in vivo and influence differentiation outcomes in human embryonic stem cells (hESCs) in vitro Systematic investigation of morphogen signaling is hampered by the difficulty of disentangling endogenous signaling from experimentally applied ligands. Here, we grow hESCs in micropatterned colonies of 1-8 cells ('µColonies') to quantitatively investigate paracrine signaling and the response to external stimuli. We examine BMP4-mediated differentiation in µColonies and standard culture conditions and find that in µColonies, above a threshold concentration, BMP4 gives rise to only a single cell fate, contrary to its role as a morphogen in other developmental systems. Under standard culture conditions BMP4 acts as a morphogen but this requires secondary signals and particular cell densities. We find that a 'community effect' enforces a common fate within µColonies, both in the state of pluripotency and when cells are differentiated, and that this effect allows a more precise response to external signals. Using live cell imaging to correlate signaling histories with cell fates, we demonstrate that interactions between neighbors result in sustained, homogenous signaling necessary for differentiation.

A Geometric Model of Stripe Refinement

Organizing data about patterning and morphogenesis into a coherent framework remains a challenge in developmental biology. Reporting in Science, Corson et al. (2017) apply innovative analysis to an old problem of bristle patterns in Drosophila, reducing the nonlinear interactions among tens of cells to a succinct model with quantitative predictions.

Ectopic hbox12 Expression Evoked by Histone Deacetylase Inhibition Disrupts Axial Specification of the Sea Urchin Embryo

Dorsal/ventral patterning of the sea urchin embryo depends upon the establishment of a Nodal-expressing ventral organizer. Recently, we showed that spatial positioning of this organizer relies on the dorsal-specific transcription of the Hbox12 repressor. Building on these findings, we determined the influence of the epigenetic milieu on the expression of hbox12 and nodal genes. We find that Trichostatin-A, a potent and selective histone-deacetylases inhibitor, induces histone hyperacetylation in hbox12 chromatin, evoking broad ectopic expression of the gene. Transcription of nodal concomitantly drops, prejudicing dorsal/ventral polarity of the resulting larvae. Remarkably, impairing hbox12 function, either in a spatially-restricted sector or in the whole embryo, specifically rescues nodal transcription in Trichostatin-A-treated larvae. Beyond strengthen the notion that nodal expression is not allowed in the presence of functional Hbox12 in the same cells, these results highlight a critical role of histone deacetylases in regulating the spatial expression of hbox12.

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.

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

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

Specific functions of the Wnt signaling system in gene regulatory networks throughout the early sea urchin embryo

Wnt signaling affects cell-fate specification processes throughout embryonic development. Here we take advantage of the well-studied gene regulatory networks (GRNs) that control pregastrular sea urchin embryogenesis to reveal the gene regulatory functions of the entire Wnt-signaling system. Five wnt genes, three frizzled genes, two secreted frizzled-related protein 1 genes, and two Dickkopf genes are expressed in dynamic spatial patterns in the pregastrular embryo of Strongylocentrotus purpuratus. We present a comprehensive analysis of these genes in each embryonic domain. Total functions of the Wnt-signaling system in regulatory gene expression throughout the embryo were studied by use of the Porcupine inhibitor C59, which interferes with zygotic Wnt ligand secretion. Morpholino-mediated knockdown of each expressed Wnt ligand demonstrated that individual Wnt ligands are functionally distinct, despite their partially overlapping spatial expression. They target specific embryonic domains and affect particular regulatory genes. The sum of the effects of blocking expression of individual wnt genes is shown to equal C59 effects. Remarkably, zygotic Wnt-signaling inputs are required for only three general aspects of embryonic specification: the broad activation of endodermal GRNs, the regional specification of the immediately adjacent stripe of ectoderm, and the restriction of the apical neurogenic domain. All Wnt signaling in this pregastrular embryo is short range (and/or autocrine). Furthermore, we show that the transcriptional drivers of wnt genes execute important specification functions in the embryonic domains targeted by the ligands, thus connecting the expression and function of wnt genes by encoded cross-regulatory interactions within the specific regional GRNs.

Genome-Wide Mapping of Chromatin State of Mouse Forelimbs

BACKGROUND

Cell types are defined at the molecular level during embryogenesis by a process called pattern formation and created by the selective utilization of combinations of sequence specific transcription factors. Developmental programs define the sets of genes that are available to each particular cell type, and real-time biochemical signaling interactions define the extent to which these sets are used at any given time and place. Gene expression is regulated through the integrated action of many cis-regulatory elements, including core promoters, enhancers, silencers, and insulators. The chromatin state in developing body parts provides a code to cellular populations that direct their cell fates. Chromatin profiling has been a method of choice for mapping regulatory sequences in cells that go through developmental transitions.

RESULTS

We used antibodies against histone H3 lysine 4 trimethylations (H3K4me3) a modification associated with promoters and open/active chromatin, histone H3 lysine 27 trimethylations (H3K27me3) associated with Polycomb-repressed regions and RNA polymerase II (Pol2) associated with transcriptional initiation to identify the chromatin state signature of the mouse forelimb during mid-gestation, at embryonic day 12 (E12). The families of genes marked included those related to transcriptional regulation and embryogenesis. One third of the marked genes were transcriptionally active while only a small fraction were bivalent marked. Sequence specific transcription factors that were activated were involved in cell specification including bone and muscle formation.

CONCLUSION

Our results demonstrate that embryonic limb cells do not exhibit the plasticity of the ES cells but are rather programmed for a finer tuning for cell lineage specification.

Network dynamics in the tumor microenvironment

The evolutionary path from tumor initiation to metastasis can only be fully understood by considering cancer cells as part of a multi-species ecosystem within the tumor microenvironment. This paper reviews and suggests two important recent trends. Firstly, I review arguments that interactions among diverse cells in the tumor microenvironment create a distinct cellular environment that can confer growth advantages, resist interventions, and allow tumors to remain dormant for long periods. Second, I review and highlight a trend toward data-rich, molecularly detailed, computational models of the tumor microenvironment. I argue that data-driven molecularly detailed tumor microenvironment models can now be built using data from multiple emerging high-throughput technologies, and that such models can pinpoint mechanisms of dysregulation and suggest specific drug targets and follow up experiments.

Oral-aboral axis specification in the sea urchin embryo, IV: hypoxia radializes embryos by preventing the initial spatialization of nodal activity

The oral-aboral axis of the sea urchin embryo is specified conditionally via a regulated feedback circuit involving the signaling gene nodal and its antagonist lefty. In normal development nodal activity becomes localized to the prospective oral side of the blastula stage embryo, a process that requires lefty. In embryos of Strongylocentrotus purpuratus, a redox gradient established by asymmetrically distributed mitochondria provides an initial spatial input that positions the localized domain of nodal expression. This expression is perturbed by hypoxia, leading to development of radialized embryos lacking an oral-aboral axis. Here we show that this radialization is not caused by a failure to express nodal, but rather by a failure to localize nodal activity to one side of the embryo. This occurs even when embryos are removed from hypoxia at late cleavage stage when nodal is first expressed, indicating that the effect involves the initiation phase of nodal activity, rather than its positive feedback-driven amplification and maintenance. Quantitative fluorescence microscopy of MitoTracker Orange-labeled embryos expressing nodal-GFP reporter gene revealed that hypoxia abolishes the spatial correlation between mitochondrial distribution and nodal expression, suggesting that hypoxia eliminates the initial spatial bias in nodal activity normally established by the redox gradient. We propose that absent this bias, the initiation phase of nodal expression is spatially uniform, such that the ensuing Nodal-mediated community effect is not localized, and hence refractory to Lefty-mediated enforcement of localization.

New regulatory circuit controlling spatial and temporal gene expression in the sea urchin embryo oral ectoderm GRN

The sea urchin oral ectoderm gene regulatory network (GRN) model has increased in complexity as additional genes are added to it, revealing its multiple spatial regulatory state domains. The formation of the oral ectoderm begins with an oral-aboral redox gradient, which is interpreted by the cis-regulatory system of the nodal gene to cause its expression on the oral side of the embryo. Nodal signaling drives cohorts of regulatory genes within the oral ectoderm and its derived subdomains. Activation of these genes occurs sequentially, spanning the entire blastula stage. During this process the stomodeal subdomain emerges inside of the oral ectoderm, and bilateral subdomains defining the lateral portions of the future ciliary band emerge adjacent to the central oral ectoderm. Here we examine two regulatory genes encoding repressors, sip1 and ets4, which selectively prevent transcription of oral ectoderm genes until their expression is cleared from the oral ectoderm as an indirect consequence of Nodal signaling. We show that the timing of transcriptional de-repression of sip1 and ets4 targets which occurs upon their clearance explains the dynamics of oral ectoderm gene expression. In addition two other repressors, the direct Nodal target not, and the feed forward Nodal target goosecoid, repress expression of regulatory genes in the central animal oral ectoderm thereby confining their expression to the lateral domains of the animal ectoderm. These results have permitted construction of an enhanced animal ectoderm GRN model highlighting the repressive interactions providing precise temporal and spatial control of regulatory gene expression.

On the role of lateral stabilization during early patterning in the pancreas

The cell fate decision of multi-potent pancreatic progenitor cells between the exocrine and endocrine lineages is regulated by Notch signalling, mediated by cell-cell interactions. However, canonical models of Notch-mediated lateral inhibition cannot explain the scattered spatial distribution of endocrine cells and the cell-type ratio in the developing pancreas. Based on evidence from acinar-to-islet cell transdifferentiation in vitro, we propose that lateral stabilization, i.e. positive feedback between adjacent progenitor cells, acts in parallel with lateral inhibition to regulate pattern formation in the pancreas. A simple mathematical model of transcriptional regulation and cell-cell interaction reveals the existence of multi-stability of spatial patterns whose simultaneous occurrence causes scattering of endocrine cells in the presence of noise. The scattering pattern allows for control of the endocrine-to-exocrine cell-type ratio by modulation of lateral stabilization strength. These theoretical results suggest a previously unrecognized role for lateral stabilization in lineage specification, spatial patterning and cell-type ratio control in organ development.

Dpp-induced Egfr signaling triggers postembryonic wing development in Drosophila

The acquisition of flight contributed to the success of insects and winged forms are present in most orders. Key to understanding the origin of wings will be knowledge of the earliest postembryonic events promoting wing outgrowth. The Drosophila melanogaster wing is intensely studied as a model appendage, and yet little is known about the beginning of wing outgrowth. Vein (Vn) is a neuregulin-like ligand for the EGF receptor (Egfr), which is necessary for global development of the early Drosophila wing disc. vn is not expressed in the embryonic wing primordium and thus has to be induced de novo in the nascent larval wing disc. We find that Decapentaplegic (Dpp), a Bone Morphogenetic Protein (BMP) family member, provides the instructive signal for initiating vn expression. The signaling involves paracrine communication between two epithelia in the early disc. Once initiated, vn expression is amplified and maintained by autocrine signaling mediated by the E-twenty six (ETS)-factor PointedP2 (PntP2). This interplay of paracrine and autocrine signaling underlies the spatial and temporal pattern of induction of Vn/Egfr target genes and explains both body wall development and wing outgrowth. It is possible this gene regulatory network governing expression of an EGF ligand is conserved and reflects a common origin of insect wings.

Direct and indirect control of oral ectoderm regulatory gene expression by Nodal signaling in the sea urchin embryo

The Nodal signaling pathway is known from earlier work to be an essential mediator of oral ectoderm specification in the sea urchin embryo, and indirectly, of aboral ectoderm specification as well. Following expression of the Nodal ligand in the future oral ectoderm during cleavage, a sequence of regulatory gene activations occur within this territory which depend directly or indirectly on nodal gene expression. Here we describe additional regulatory genes that contribute to the oral ectoderm regulatory state during specification in Strongylocentrotus purpuratus, and show how their spatial expression changes dynamically during development. By means of system wide perturbation analyses we have significantly improved current knowledge of the epistatic relations among the regulatory genes of the oral ectoderm. From these studies there emerge diverse circuitries relating downstream regulatory genes directly and indirectly to Nodal signaling. A key intermediary regulator, the role of which had not previously been discerned, is the not gene. In addition to activating several genes earlier described as targets of Nodal signaling, the not gene product acts to repress other oral ectoderm genes, contributing crucially to the bilateral spatial organization of the embryonic oral ectoderm.

Initial deployment of the cardiogenic gene regulatory network in the basal chordate, Ciona intestinalis

The complex, partially redundant gene regulatory architecture underlying vertebrate heart formation has been difficult to characterize. Here, we dissect the primary cardiac gene regulatory network in the invertebrate chordate, Ciona intestinalis. The Ciona heart progenitor lineage is first specified by Fibroblast Growth Factor/Map Kinase (FGF/MapK) activation of the transcription factor Ets1/2 (Ets). Through microarray analysis of sorted heart progenitor cells, we identified the complete set of primary genes upregulated by FGF/Ets shortly after heart progenitor emergence. Combinatorial sequence analysis of these co-regulated genes generated a hypothetical regulatory code consisting of Ets binding sites associated with a specific co-motif, ATTA. Through extensive reporter analysis, we confirmed the functional importance of the ATTA co-motif in primary heart progenitor gene regulation. We then used the Ets/ATTA combination motif to successfully predict a number of additional heart progenitor gene regulatory elements, including an intronic element driving expression of the core conserved cardiac transcription factor, GATAa. This work significantly advances our understanding of the Ciona heart gene network. Furthermore, this work has begun to elucidate the precise regulatory architecture underlying the conserved, primary role of FGF/Ets in chordate heart lineage specification.

Axial patterning interactions in the sea urchin embryo: suppression of nodal by Wnt1 signaling

Wnt and Nodal signaling pathways are required for initial patterning of cell fates along anterior-posterior (AP) and dorsal-ventral (DV) axes, respectively, of sea urchin embryos during cleavage and early blastula stages. These mechanisms are connected because expression of nodal depends on early Wnt/β-catenin signaling. Here, we show that an important subsequent function of Wnt signaling is to control the shape of the nodal expression domain and maintain correct specification of different cell types along the axes of the embryo. In the absence of Wnt1, the posterior-ventral region of the embryo is severely altered during early gastrulation. Strikingly, at this time, nodal and its downstream target genes gsc and bra are expressed ectopically, extending posteriorly to the blastopore. They override the initial specification of posterior-ventral ectoderm and endoderm fates, eliminating the ventral contribution to the gut and displacing the ciliary band dorsally towards, and occasionally beyond, the blastopore. Consequently, in Wnt1 morphants, the blastopore is located at the border of the re-specified posterior-ventral oral ectoderm and by larval stages it is in the same plane near the stomodeum on the ventral side. In normal embryos, a Nodal-dependent process downregulates wnt1 expression in dorsal posterior cells during early gastrulation, focusing Wnt1 signaling to the posterior-ventral region where it suppresses nodal expression. These subsequent interactions between Wnt and Nodal signaling are thus mutually antagonistic, each limiting the range of the other's activity, in order to maintain and stabilize the body plan initially established by those same signaling pathways in the early embryo.

The evolution of nervous system patterning: insights from sea urchin development

Recent studies of the sea urchin embryo have elucidated the mechanisms that localize and pattern its nervous system. These studies have revealed the presence of two overlapping regions of neurogenic potential at the beginning of embryogenesis, each of which becomes progressively restricted by separate, yet linked, signals, including Wnt and subsequently Nodal and BMP. These signals act to specify and localize the embryonic neural fields - the anterior neuroectoderm and the more posterior ciliary band neuroectoderm - during development. Here, we review these conserved nervous system patterning signals and consider how the relationships between them might have changed during deuterostome evolution.

Evolution of gene regulatory networks controlling body plan development

Evolutionary change in animal morphology results from alteration of the functional organization of the gene regulatory networks (GRNs) that control development of the body plan. A major mechanism of evolutionary change in GRN structure is alteration of cis-regulatory modules that determine regulatory gene expression. Here we consider the causes and consequences of GRN evolution. Although some GRN subcircuits are of great antiquity, other aspects are highly flexible and thus in any given genome more recent. This mosaic view of the evolution of GRN structure explains major aspects of evolutionary process, such as hierarchical phylogeny and discontinuities of paleontological change.

Theoretical basis of the community effect in development

BACKGROUND

Genetically identical cells often show significant variation in gene expression profile and behaviour even in the same physiological condition. Notably, embryonic cells destined to the same tissue maintain a uniform transcriptional regulatory state and form a homogeneous cell group. One mechanism to keep the homogeneity within embryonic tissues is the so-called community effect in animal development. The community effect is an interaction among a group of many nearby precursor cells, and is necessary for them to maintain tissue-specific gene expression and differentiate in a coordinated manner. Although it has been shown that the cell-cell communication by a diffusible factor plays a crucial role, it is not immediately obvious why a community effect needs many cells.

RESULTS

In this work, we propose a model of the community effect in development, which consists in a linear gene cascade and cell-cell communication. We examined the properties of the model theoretically using a combination of stochastic and deterministic modelling methods. We have derived the analytical formula for the threshold size of a cell population that is necessary for a community effect, which is in good agreement with stochastic simulation results.

CONCLUSIONS

Our theoretical analysis indicates that a simple model with a linear gene cascade and cell-cell communication is sufficient to reproduce the community effect in development. The model explains why a community needs many cells. It suggests that the community's long-term behaviour is independent of the initial induction level, although the initiation of a community effect requires a sufficient amount of inducing signal. The mechanism of the community effect revealed by our theoretical analysis is analogous to that of quorum sensing in bacteria. The community effect may underlie the size control in animal development and also the genesis of autosomal dominant diseases including tumorigenesis.

Nodal-mediated epigenesis requires dynamin-mediated endocytosis

Nodal proteins are diffusible morphogens that drive pattern formation via short-range feedback activation coupled to long-range Lefty-mediated inhibition. In the sea urchin embryo, specification of the secondary (oral-aboral) axis occurs via zygotic expression of nodal, which is localized to the prospective oral ectoderm at early blastula stage. In mid-blastula stage embryos treated with low micromolar nickel or zinc, nodal expression expands progressively beyond the confines of this localized domain to encompass the entire equatorial circumference of the embryo, producing radialized embryos lacking an oral-aboral axis. RNAseq analysis of embryos treated with nickel, zinc, or cadmium (which does not radialize embryos) showed that several genes involved in endocytosis were similarly perturbed by nickel and zinc but not cadmium. Inhibiting dynamin, a GTPase required for receptor-mediated endocytosis, phenocopies the effects of nickel and zinc, suggesting that dynamin-mediated endocytosis is required as a sink to limit the range of Nodal signaling.

Emerging properties of animal gene regulatory networks

Gene regulatory networks (GRNs) provide system level explanations of developmental and physiological functions in the terms of the genomic regulatory code. Depending on their developmental functions, GRNs differ in their degree of hierarchy, and also in the types of modular sub-circuit of which they are composed, although there is a commonly employed sub-circuit repertoire. Mathematical modelling of some types of GRN sub-circuit has deepened biological understanding of the functions they mediate. The structural organization of various kinds of GRN reflects their roles in the life process, and causally illuminates both developmental and evolutionary process.

Oral-aboral patterning and gastrulation of sea urchin embryos depend on sulfated glycosaminoglycans

Glycosaminoglycans (GAGs) are a heavily sulfated component of the extracellular matrix (ECM) implicated in a variety of cell signaling events involved in patterning of embryos. Embryos of the sea urchin Strongylocentrotus purpuratus were exposed to several inhibitors that disrupt GAG function during development. Treatment with chlorate, a general inhibitor of sulfation that leads to undersulfated GAGs, reduced sulfation of the urchin blastocoelar ECM. It also prevented correct specification of the oral-aboral axis and mouth formation, resulting in a radialized phenotype characterized by the lack of an oral field, incomplete gastrulation and formation of multiple skeletal spicule rudiments. Oral markers were initially expressed in most of the prospective ectoderm of chlorate-treated early blastulae, but then declined as aboral markers became expressed throughout most of the ectoderm. Nodal expression in the presumptive oral field is necessary and sufficient to specify the oral-aboral axis in urchins. Several lines of evidence suggest a deregulation of Nodal signaling is involved in the radialization caused by chlorate: (1) Radial embryos resemble those in which Nodal expression was knocked down. (2) Chlorate disrupted localized nodal expression in oral ectoderm, even when applied after the oral-aboral axis is specified and expression of other oral markers is resistant to treatment. (3) Inhibition with SB-431542 of ALK-4/5/7 receptors that mediate Nodal signaling causes defects in ectodermal patterning similar to those caused by chlorate. (4) Intriguingly, treatment of embryos with a sub-threshold dose of SB-431542 rescued the radialization caused by low concentrations of chlorate. Our results indicate important roles for sulfated GAGs in Nodal signaling and oral-aboral axial patterning, and in the cellular processes necessary for archenteron extension and mouth formation during gastrulation. We propose that interaction of the Nodal ligand with sulfated GAGs limits its diffusion, and is required to specify an oral field in the urchin embryo and organize the oral-aboral axis.

TGFβ signaling positions the ciliary band and patterns neurons in the sea urchin embryo

The ciliary band is a distinct region of embryonic ectoderm that is specified between oral and aboral ectoderm. Flask-shaped ciliary cells and neurons differentiate in this region and they are patterned to form an integrated tissue that functions as the principal swimming and feeding organ of the larva. TGFβ signaling, which is known to mediate oral and aboral patterning of the ectoderm, has been implicated in ciliary band formation. We have used morpholino knockdown and ectopic expression of RNA to alter TGFβ signaling at the level of ligands, receptors, and signal transduction components and assessed the differentiation and patterning of the ciliary band cells and associated neurons. We propose that the primary effects of these signals are to position the ciliary cells, which in turn support neural differentiation. We show that Nodal signaling, which is known to be localized by Lefty, positions the oral margin of the ciliary band. Signaling from BMP through Alk3/6, affects the position of the oral and aboral margins of the ciliary band. Since both Nodal and BMP signaling produce ectoderm that does not support neurogenesis, we propose that formation of a ciliary band requires protection from these signals. Expression of BMP2/4 and Nodal suppress neural differentiation. However, the response to receptor knockdown or dominant-negative forms of signal transduction components indicate signaling is not acting directly on unspecified ectoderm cells to prevent their differentiation as neurons. Instead, it produces a restricted field of ciliary band cells that supports neurogenesis. We propose a model that incorporates spatially regulated control of Nodal and BMP signaling to determine the position and differentiation of the ciliary band, and subsequent neural patterning.

High accuracy, high-resolution prevalence measurement for the majority of locally expressed regulatory genes in early sea urchin development

Accurate measurements of transcript abundance are a prerequisite to understand gene activity in development. Using the NanoString nCounter, an RNA counting device, we measured the prevalence of 172 transcription factors and signaling molecules in early sea urchin development. These measurements show high fidelity over more than five orders of magnitude down to a few transcripts per embryo. Most of the genes included are locally restricted in their spatial expression, and contribute to the divergent regulatory states of cells in the developing embryo. In order to obtain high-resolution expression profiles from fertilization to late gastrulation samples were collected at hourly intervals. The measured time courses agree well with, and substantially extend, prior relative abundance measurements obtained by quantitative PCR. High temporal resolution permits sequences of successively activated genes to be precisely delineated providing an ancillary tool for assembling maps of gene regulatory networks. The data are available via an interactive website for quick plotting of selected time courses.

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

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

A service-oriented architecture for integrating the modeling and formal verification of genetic regulatory networks

Background

The study of biological networks has led to the development of increasingly large and detailed models. Computer tools are essential for the simulation of the dynamical behavior of the networks from the model. However, as the size of the models grows, it becomes infeasible to manually verify the predictions against experimental data or identify interesting features in a large number of simulation traces. Formal verification based on temporal logic and model checking provides promising methods to automate and scale the analysis of the models. However, a framework that tightly integrates modeling and simulation tools with model checkers is currently missing, on both the conceptual and the implementational level.

Results

We have developed a generic and modular web service, based on a service-oriented architecture, for integrating the modeling and formal verification of genetic regulatory networks. The architecture has been implemented in the context of the qualitative modeling and simulation tool GNA and the model checkers NUSMV and CADP. GNA has been extended with a verification module for the specification and checking of biological properties. The verification module also allows the display and visual inspection of the verification results.

Conclusions

The practical use of the proposed web service is illustrated by means of a scenario involving the analysis of a qualitative model of the carbon starvation response in E. coli. The service-oriented architecture allows modelers to define the model and proceed with the specification and formal verification of the biological properties by means of a unified graphical user interface. This guarantees a transparent access to formal verification technology for modelers of genetic regulatory networks.

Network design principles from the sea urchin embryo

As gene regulatory network models encompass more and more of the specification processes underlying sea urchin embryonic development, topological themes emerge that imply the existence of structural network 'building blocks'. These are subcircuits which perform given logic operations in the spatial control of gene expression. The various parts of the sea urchin gene regulatory networks offer instances of the same subcircuit topologies accomplishing the same developmental logic functions but using different genes. These subcircuits are dedicated to specific developmental functions, unlike simpler 'motifs', and may indicate a repertoire of specific devices of which developmental gene regulatory networks are composed.

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.

Monte Carlo analysis of an ODE Model of the Sea Urchin Endomesoderm Network

BACKGROUND

Gene Regulatory Networks (GRNs) control the differentiation, specification and function of cells at the genomic level. The levels of interactions within large GRNs are of enormous depth and complexity. Details about many GRNs are emerging, but in most cases it is unknown to what extent they control a given process, i.e. the grade of completeness is uncertain. This uncertainty stems from limited experimental data, which is the main bottleneck for creating detailed dynamical models of cellular processes. Parameter estimation for each node is often infeasible for very large GRNs. We propose a method, based on random parameter estimations through Monte-Carlo simulations to measure completeness grades of GRNs.

RESULTS

We developed a heuristic to assess the completeness of large GRNs, using ODE simulations under different conditions and randomly sampled parameter sets to detect parameter-invariant effects of perturbations. To test this heuristic, we constructed the first ODE model of the whole sea urchin endomesoderm GRN, one of the best studied large GRNs. We find that nearly 48% of the parameter-invariant effects correspond with experimental data, which is 65% of the expected optimal agreement obtained from a submodel for which kinetic parameters were estimated and used for simulations. Randomized versions of the model reproduce only 23.5% of the experimental data.

CONCLUSION

The method described in this paper enables an evaluation of network topologies of GRNs without requiring any parameter values. The benefit of this method is exemplified in the first mathematical analysis of the complete Endomesoderm Network Model. The predictions we provide deliver candidate nodes in the network that are likely to be erroneous or miss unknown connections, which may need additional experiments to improve the network topology. This mathematical model can serve as a scaffold for detailed and more realistic models. We propose that our method can be used to assess a completeness grade of any GRN. This could be especially useful for GRNs involved in human diseases, where often the amount of connectivity is unknown and/or many genes/interactions are missing.