Development. Built to run, not fail

How haematopoiesis research became a fertile ground for regulatory network biology as pioneered by Eric Davidson

PURPOSE OF REVIEW

This historical perspective reviews how work of Eric H. Davidson was a catalyst and exemplar for explaining haematopoietic cell fate determination through gene regulation.

RECENT FINDINGS

Researchers studying blood and immune cells pioneered many of the early mechanistic investigations of mammalian gene regulatory processes. These efforts included the characterization of complex gene regulatory sequences exemplified by the globin and T-cell/B-cell receptor gene loci, as well as the identification of many key regulatory transcription factors through the fine mapping of chromosome translocation breakpoints in leukaemia patients. As the repertoire of known regulators expanded, assembly into gene regulatory network models became increasingly important, not only to account for the truism that regulatory genes do not function in isolation but also to devise new ways of extracting biologically meaningful insights from even more complex information. Here we explore how Eric H. Davidson's pioneering studies of gene regulatory network control in nonvertebrate model organisms have had an important and lasting impact on research into blood and immune cell development.

SUMMARY

The intellectual framework developed by Davidson continues to contribute to haematopoietic research, and his insistence on demonstrating logic and causality still challenges the frontier of research today.

Modularity, criticality, and evolvability of a developmental gene regulatory network

The existence of discrete phenotypic traits suggests that the complex regulatory processes which produce them are functionally modular. These processes are usually represented by networks. Only modular networks can be partitioned into intelligible subcircuits able to evolve relatively independently. Traditionally, functional modularity is approximated by detection of modularity in network structure. However, the correlation between structure and function is loose. Many regulatory networks exhibit modular behaviour without structural modularity. Here we partition an experimentally tractable regulatory network-the gap gene system of dipteran insects-using an alternative approach. We show that this system, although not structurally modular, is composed of dynamical modules driving different aspects of whole-network behaviour. All these subcircuits share the same regulatory structure, but differ in components and sensitivity to regulatory interactions. Some subcircuits are in a state of criticality, while others are not, which explains the observed differential evolvability of the various expression features in the system.

Eric Davidson: Steps to a gene regulatory network for development

Eric Harris Davidson was a unique and creative intellectual force who grappled with the diversity of developmental processes used by animal embryos and wrestled them into an intelligible set of principles, then spent his life translating these process elements into molecularly definable terms through the architecture of gene regulatory networks. He took speculative risks in his theoretical writing but ran a highly organized, rigorous experimental program that yielded an unprecedentedly full characterization of a developing organism. His writings created logical order and a framework for mechanism from the complex phenomena at the heart of advanced multicellular organism development. This is a reminiscence of intellectual currents in his work as observed by the author through the last 30-35 years of Davidson's life.

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.

A dynamic regulatory network explains ParaHox gene control of gut patterning in the sea urchin

The anteroposterior patterning of the embryonic gut represents one of the most intriguing biological processes in development. A dynamic control of gene transcription regulation and cell movement is perfectly orchestrated to shape a functional gut in distinct specialized parts. Two ParaHox genes, Xlox and Cdx, play key roles in vertebrate and sea urchin gut patterning through molecular mechanisms that are still mostly unclear. Here, we have combined functional analysis methodologies with high-resolution imaging and RNA-seq to investigate Xlox and Cdx regulation and function. We reveal part of the regulatory machinery responsible for the onset of Xlox and Cdx transcription, uncover a Wnt10 signal that mediates Xlox repression in the intestinal cells, and provide evidence of Xlox- and Cdx-mediated control of stomach and intestine differentiation, respectively. Our findings offer a novel mechanistic explanation of how the control of transcription is linked to cell differentiation and morphogenesis for the development of a perfectly organized biological system such as the sea urchin larval gut.

Encoding regulatory state boundaries in the pregastrular oral ectoderm of the sea urchin embryo

By gastrulation the ectodermal territories of the sea urchin embryo have developed an unexpectedly complex spatial pattern of sharply bounded regulatory states, organized orthogonally with respect to the animal/vegetal and oral/aboral axes of the embryo. Although much is known of the gene regulatory network (GRN) linkages that generate these regulatory states, the principles by which the boundaries between them are positioned and maintained have remained undiscovered. Here we determine the encoded genomic logic responsible for the boundaries of the oral aspect of the embryo that separate endoderm from ectoderm and ectoderm from neurogenic apical plate and that delineate the several further subdivisions into which the oral ectoderm per se is partitioned. Comprehensive regulatory state maps, including all spatially expressed oral ectoderm regulatory genes, were established. The circuitry at each boundary deploys specific repressors of regulatory states across the boundary, identified in this work, plus activation by broadly expressed positive regulators. These network linkages are integrated with previously established interactions on the oral/aboral axis to generate a GRN model encompassing the 2D organization of the regulatory state pattern in the pregastrular oral ectoderm of the embryo.

Myogenesis in the sea urchin embryo: the molecular fingerprint of the myoblast precursors

Background

In sea urchin larvae the circumesophageal fibers form a prominent muscle system of mesodermal origin. Although the morphology and later development of this muscle system has been well-described, little is known about the molecular signature of these cells or their precise origin in the early embryo. As an invertebrate deuterostome that is more closely related to the vertebrates than other commonly used model systems in myogenesis, the sea urchin fills an important phylogenetic gap and provides a unique perspective on the evolution of muscle cell development.

Results

Here, we present a comprehensive description of the development of the sea urchin larval circumesophageal muscle lineage beginning with its mesodermal origin using high-resolution localization of the expression of several myogenic transcriptional regulators and differentiation genes. A few myoblasts are bilaterally distributed at the oral vegetal side of the tip of the archenteron and first appear at the late gastrula stage. The expression of the differentiation genes Myosin Heavy Chain, Tropomyosin I and II, as well as the regulatory genes MyoD2, FoxF, FoxC, FoxL1, Myocardin, Twist, and Tbx6 uniquely identify these cells. Interestingly, evolutionarily conserved myogenic factors such as Mef2, MyoR and Six1/2 are not expressed in sea urchin myoblasts but are found in other mesodermal domains of the tip of the archenteron. The regulatory states of these domains were characterized in detail. Moreover, using a combinatorial analysis of gene expression we followed the development of the FoxF/FoxC positive cells from the onset of expression to the end of gastrulation. Our data allowed us to build a complete map of the Non-Skeletogenic Mesoderm at the very early gastrula stage, in which specific molecular signatures identify the precursors of different cell types. Among them, a small group of cells within the FoxY domain, which also express FoxC and SoxE, have been identified as plausible myoblast precursors. Together, these data support a very early gastrula stage segregation of the myogenic lineage.

Conclusions

From this analysis, we are able to precisely define the regulatory and differentiation signatures of the circumesophageal muscle in the sea urchin embryo. Our findings have important implications in understanding the evolution of development of the muscle cell lineage at the molecular level. The data presented here suggest a high level of conservation of the myogenic specification mechanisms across wide phylogenetic distances, but also reveal clear cases of gene cooption.

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.

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.

Lineage-specific transcription factors and the evolution of gene regulatory networks

Nature is replete with examples of diverse cell types, tissues and body plans, forming very different creatures from genomes with similar gene complements. However, while the genes and the structures of proteins they encode can be highly conserved, the production of those proteins in specific cell types and at specific developmental time points might differ considerably between species. A full understanding of the factors that orchestrate gene expression will be essential to fully understand evolutionary variety. Transcription factor (TF) proteins, which form gene regulatory networks (GRNs) to act in cooperative or competitive partnerships to regulate gene expression, are key components of these unique regulatory programs. Although many TFs are conserved in structure and function, certain classes of TFs display extensive levels of species diversity. In this review, we highlight families of TFs that have expanded through gene duplication events to create species-unique repertoires in different evolutionary lineages. We discuss how the hierarchical structures of GRNs allow for flexible small to large-scale phenotypic changes. We survey evidence that explains how newly evolved TFs may be integrated into an existing GRN and how molecular changes in TFs might impact the GRNs. Finally, we review examples of traits that evolved due to lineage-specific TFs and species differences in GRNs.

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.

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.

Genomic control of patterning

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

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

The "Community Effect" denotes intra-territorial signaling amongst cells which constitute a particular tissue or embryonic progenitor field. The cells of the territory express the same transcriptional regulatory state, and the intra-territorial signaling is essential to maintenance of this specific regulatory state. The structure of the underlying gene regulatory network (GRN) subcircuitry explains the genomically wired mechanism by which community effect signaling is linked to the continuing transcriptional generation of the territorial regulatory state. A clear example is afforded by the oral ectoderm GRN of the sea urchin embryo where cis-regulatory evidence, experimental embryology, and network analysis combine to provide a complete picture. We review this example and consider less well known but similar cases in other developing systems where the same subcircuit GRN topology is present. To resolve mechanistic issues that arise in considering how community effect signaling could operate to produce its observed effects, we construct and analyze the behavior of a quantitative model of community effect signaling in the sea urchin embryo oral ectoderm. Community effect network topology could constitute part of the genomic regulatory code that defines transcriptional function in multicellular tissues composed of cells in contact, and hence may have arisen as a metazoan developmental strategy.

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.

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

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

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.

Modeling the dynamics of transcriptional gene regulatory networks for animal development

The dynamic process of cell fate specification is regulated by networks of regulatory genes. The architecture of the network defines the temporal order of specification events. To understand the dynamic control of the developmental process, the kinetics of mRNA and protein synthesis and the response of the cis-regulatory modules to transcription factor concentration must be considered. Here we review mathematical models for mRNA and protein synthesis kinetics which are based on experimental measurements of the rates of the relevant processes. The model comprises the response functions of cis-regulatory modules to their transcription factor inputs, by incorporating binding site occupancy and its dependence on biologically measurable quantities. We use this model to simulate gene expression, to distinguish between cis-regulatory execution of "AND" and "OR" logic functions, rationalize the oscillatory behavior of certain transcriptional auto-repressors and to show how linked subcircuits can be dealt with. Model simulations display the effects of mutation of binding sites, or perturbation of upstream gene expression. The model is a generally useful tool for understanding gene regulation and the dynamics of cell fate specification.

SpBase: the sea urchin genome database and web site

SpBase is a system of databases focused on the genomic information from sea urchins and related echinoderms. It is exposed to the public through a web site served with open source software (http://spbase.org/). The enterprise was undertaken to provide an easily used collection of information to directly support experimental work on these useful research models in cell and developmental biology. The information served from the databases emerges from the draft genomic sequence of the purple sea urchin, Strongylocentrotus purpuratus and includes sequence data and genomic resource descriptions for other members of the echinoderm clade which in total span 540 million years of evolutionary time. This version of the system contains two assemblies of the purple sea urchin genome, associated expressed sequences, gene annotations and accessory resources. Search mechanisms for the sequences and the gene annotations are provided. Because the system is maintained along with the Sea Urchin Genome resource, a database of sequenced clones is also provided.

A gene regulatory network armature for T lymphocyte specification

Choice of a T lymphoid fate by hematopoietic progenitor cells depends on sustained Notch-Delta signaling combined with tightly regulated activities of multiple transcription factors. To dissect the regulatory network connections that mediate this process, we have used high-resolution analysis of regulatory gene expression trajectories from the beginning to the end of specification, tests of the short-term Notch dependence of these gene expression changes, and analyses of the effects of overexpression of two essential transcription factors, namely PU.1 and GATA-3. Quantitative expression measurements of >50 transcription factor and marker genes have been used to derive the principal components of regulatory change through which T cell precursors progress from primitive multipotency to T lineage commitment. Our analyses reveal separate contributions of Notch signaling, GATA-3 activity, and down-regulation of PU.1. Using BioTapestry (www.BioTapestry.org), the results have been assembled into a draft gene regulatory network for the specification of T cell precursors and the choice of T as opposed to myeloid/dendritic or mast-cell fates. This network also accommodates effects of E proteins and mutual repression circuits of Gfi1 against Egr-2 and of TCF-1 against PU.1 as proposed elsewhere, but requires additional functions that remain unidentified. Distinctive features of this network structure include the intense dose dependence of GATA-3 effects, the gene-specific modulation of PU.1 activity based on Notch activity, the lack of direct opposition between PU.1 and GATA-3, and the need for a distinct, late-acting repressive function or functions to extinguish stem and progenitor-derived regulatory gene expression.

Properties of developmental gene regulatory networks

The modular components, or subcircuits, of developmental gene regulatory networks (GRNs) execute specific developmental functions, such as the specification of cell identity. We survey examples of such subcircuits and relate their structures to corresponding developmental functions. These relations transcend organisms and genes, as illustrated by the similar structures of the subcircuits controlling the specification of the mesectoderm in the Drosophila embryo and the endomesoderm in the sea urchin, even though the respective subcircuits are composed of nonorthologous regulatory genes.

Visualization, documentation, analysis, and communication of large-scale gene regulatory networks

Genetic regulatory networks (GRNs) are complex, large-scale, and spatially and temporally distributed. These characteristics impose challenging demands on software tools for building GRN models, and so there is a need for custom tools. In this paper, we report on our ongoing development of BioTapestry, an open source, freely available computational tool designed specifically for building GRN models. We also outline our future development plans, and give some examples of current applications of BioTapestry.

The MYB98 subcircuit of the synergid gene regulatory network includes genes directly and indirectly regulated by MYB98

The female gametophyte contains two synergid cells that play a role in many steps of the angiosperm reproductive process, including pollen tube guidance. At their micropylar poles, the synergid cells have a thickened and elaborated cell wall: the filiform apparatus that is thought to play a role in the secretion of the pollen tube attractant(s). MYB98 regulates an important subcircuit of the synergid gene regulatory network (GRN) that functions to activate the expression of genes required for pollen tube guidance and filiform apparatus formation. The MYB98 subcircuit comprises at least 83 downstream genes, including 48 genes within four gene families (CRP810, CRP3700, CRP3730 and CRP3740) that encode Cys-rich proteins. We show that the 11 CRP3700 genes, which include DD11 and DD18, are regulated by a common cis-element, GTAACNT, and that a multimer of this sequence confers MYB98-dependent synergid expression. The GTAACNT element contains the MYB98-binding site identified in vitro, suggesting that the 11 CRP3700 genes are direct targets of MYB98. We also show that five of the CRP810 genes, which include DD2, lack a functional GTAACNT element, suggesting that they are not directly regulated by MYB98. In addition, we show that the five CRP810 genes are regulated by the cis-element AACGT, and that a multimer of this sequence confers synergid expression. Together, these results suggest that the MYB98 branch of the synergid GRN is multi-tiered and, therefore, contains at least one additional downstream transcription factor.

Evolutionary design principles of modules that control cellular differentiation: consequences for hysteresis and multistationarity

MOTIVATION

Gene regulatory networks (GRNs) govern cellular differentiation processes and enable construction of multicellular organisms from single cells. Although such networks are complex, there must be evolutionary design principles that shape the network to its present form, gaining complexity from simple modules.

RESULTS

To isolate particular design principles, we have computationally evolved random regulatory networks with a preference to result either in hysteresis (switching threshold depending on current state), or in multistationarity (having multiple steady states), two commonly observed dynamical features of GRNs related to differentiation processes. We have analyzed the resulting evolved networks and compared their structures and characteristics with real GRNs reported from experiments.

CONCLUSION

We found that the artificially evolved networks have particular topologies and it was notable that these topologies share important features and similarities with the real GRNs, particularly in contrasting properties of positive and negative feedback loops. We conclude that the structures of real GRNs are consistent with selection to favor one or other of the dynamical features of multistationarity or hysteresis.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Global regulatory logic for specification of an embryonic cell lineage

Explanation of a process of development must ultimately be couched in the terms of the genomic regulatory code. Specification of an embryonic cell lineage is driven by a network of interactions among genes encoding transcription factors. Here, we present the gene regulatory network (GRN) that directs the specification of the skeletogenic micromere lineage of the sea urchin embryo. The GRN now includes all regulatory genes expressed in this lineage up to late blastula stage, as identified in a genomewide survey. The architecture of the GRN was established by a large-scale perturbation analysis in which the expression of each gene in the GRN was cut off by use of morpholinos, and the effects on all other genes were measured quantitatively. Several cis-regulatory analyses provided additional evidence. The explanatory power of the GRN suffices to provide a causal explanation for all observable developmental functions of the micromere lineage during the specification period. These functions are: (i) initial acquisition of identity through transcriptional interpretation of localized maternal cues; (ii) activation of specific regulatory genes by use of a double negative gate; (iii) dynamic stabilization of the regulatory state by activation of a feedback subcircuit; (iv) exclusion of alternative regulatory states; (v) presentation of a signal required by the micromeres themselves and of two different signals required for development of adjacent endomesodermal lineages; and (vi) lineage-specific activation of batteries of skeletogenic genes. The GRN precisely predicts gene expression responses and provides a coherent explanation of the biology of specification.

Dynamical patterning modules: physico-genetic determinants of morphological development and evolution

The shapes and forms of multicellular organisms arise by the generation of new cell states and types and changes in the numbers and rearrangements of the various kinds of cells. While morphogenesis and pattern formation in all animal species are widely recognized to be mediated by the gene products of an evolutionarily conserved 'developmental-genetic toolkit', the link between these molecular players and the physics underlying these processes has been generally ignored. This paper introduces the concept of 'dynamical patterning modules' (DPMs), units consisting of one or more products of the 'toolkit' genes that mobilize physical processes characteristic of chemically and mechanically excitable meso- to macroscopic systems such as cell aggregates: cohesion, viscoelasticity, diffusion, spatiotemporal heterogeneity based on lateral inhibition and multistable and oscillatory dynamics. We suggest that ancient toolkit gene products, most predating the emergence of multicellularity, assumed novel morphogenetic functions due to change in the scale and context inherent to multicellularity. We show that DPMs, acting individually and in concert with each other, constitute a 'pattern language' capable of generating all metazoan body plans and organ forms. The physical dimension of developmental causation implies that multicellular forms during the explosive radiation of animal body plans in the middle Cambrian, approximately 530 million years ago, could have explored an extensive morphospace without concomitant genotypic change or selection for adaptation. The morphologically plastic body plans and organ forms generated by DPMs, and their ontogenetic trajectories, would subsequently have been stabilized and consolidated by natural selection and genetic drift. This perspective also solves the apparent 'molecular homology-analogy paradox', whereby widely divergent modern animal types utilize the same molecular toolkit during development by proposing, in contrast to the Neo-Darwinian principle, that phenotypic disparity early in evolution occurred in advance of, rather than closely tracked, genotypic change.

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.

The GATA factor Serpent cross-regulates lozenge and u-shaped expression during Drosophila blood cell development

The Drosophila GATA factor Serpent interacts with the RUNX factor Lozenge to activate the crystal cell program, whereas SerpentNC binds the Friend of GATA protein U-shaped to limit crystal cell production. Here, we identified a lozenge minimal hematopoietic cis-regulatory module and showed that lozenge-lacZ reporter-gene expression was autoregulated by Serpent and Lozenge. We also showed that upregulation of u-shaped was delayed until after lozenge activation, consistent with our previous results that showed u-shaped expression in the crystal cell lineage is dependent on both Serpent and Lozenge. Together, these observations describe a feed forward regulatory motif, which controls the temporal expression of u-shaped. Finally, we showed that lozenge reporter-gene activity increased in a u-shaped mutant background and that forced expression of SerpentNC with U-shaped blocked lozenge- and u-shaped-lacZ reporter-gene activity. This is the first demonstration of GATA:FOG regulation of Runx and Fog gene expression. Moreover, these results identify components of a Serpent cross-regulatory sub-circuit that can modulate lozenge expression. Based on the sub-circuit design and the combinatorial control of crystal cell production, we present a model for the specification of a dynamic bi-potential regulatory state that contributes to the selection between a Lozenge-positive and Lozenge-negative state.

The regulatory genome and the computer

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

Logic of gene regulatory networks

Regulatory networks of transcription factors and signaling molecules lie at the heart of development. Their architecture implements logic functions whose execution propels cells from one regulatory state to the next, thus driving development forward. As an example of a subcircuit that translates transcriptional input into developmental output, we consider a particularly simple case, the regulatory processes underlying pigment cell formation in sea urchin embryos. The regulatory events in this process can be represented as elementary logic functions.

A missing link in the sea urchin embryo gene regulatory network: hesC and the double-negative specification of micromeres

Specification of sea urchin embryo micromeres occurs early in cleavage, with the establishment of a well defined regulatory state. The architecture of the gene regulatory network controlling the specification process indicates that transcription of the initial tier of control genes depends on a double-negative gate. A gene encoding a transcriptional repressor, pmar1, is activated specifically in micromeres, where it represses transcription of a second repressor that is otherwise active globally. Thus, the micromere-specific control genes, which are the target of the second repressor, are expressed exclusively in this lineage. The double-negative specification gate was logically required from the results of numerous prior experiments, but the identity of the gene encoding the second repressor remained elusive. Here we show that hesC is this gene, and we demonstrate experimentally all of its predicted functions, including global repression of micromere-specific regulatory genes. As logically required, blockade of hesC mRNA translation and global overexpression of pmar1 mRNA have the same effect, which is to cause all of the cells of the embryo to express micromere-specific genes.