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

A gene regulatory network (GRN) controls the process by which the endomesoderm of the sea urchin embryo is specified. In this GRN, the program of gene expression unique to the skeletogenic micromere lineage is set in train by activation of the pmar1 gene. Through a double repression system, this gene is responsible for localization of expression of downstream regulatory and signaling genes to cells of this lineage. One of these genes, delta, encodes a Notch ligand, and its expression in the right place and time is crucial to the specification of the endomesoderm. Here we report a cis-regulatory element R11 that is responsible for localizing the expression of delta by means of its response to the pmar1 repression system. R11 was identified as an evolutionarily conserved genomic sequence located about 13 kb downstream of the last exon of the delta gene. We demonstrate here that this cis-regulatory element is able to drive the expression of a reporter gene in the same cells and at the same time that the endogenous delta gene is expressed, and that temporally, spatially, and quantitatively it responds to the pmar1 repression system just as predicted for the delta gene in the endomesoderm GRN. This work illustrates the application of cis-regulatory analysis to the validation of predictions of the GRN model. In addition, we introduce new methodological tools for quantitative measurement of the output of expression constructs that promise to be of general value for cis-regulatory analysis in sea urchin embryos.     

Wnt6 activates endoderm in the sea urchin gene regulatory network

In the sea urchin, entry of β-catenin into the nuclei of the vegetal cells at 4th and 5th cleavages is necessary for activation of the endomesoderm gene regulatory network. Beyond that, little is known about how the embryo uses maternal information to initiate specification. Here, experiments establish that of the three maternal Wnts in the egg, Wnt6 is necessary for activation of endodermal genes in the endomesoderm GRN. A small region of the vegetal cortex is shown to be necessary for activation of the endomesoderm GRN. If that cortical region of the egg is removed, addition of Wnt6 rescues endoderm. At a molecular level, the vegetal cortex region contains a localized concentration of Dishevelled (Dsh) protein, a transducer of the canonical Wnt pathway; however, Wnt6 mRNA is not similarly localized. Ectopic activation of the Wnt pathway, through the expression of an activated form of β-catenin, of a dominant-negative variant of GSK-3β or of Dsh itself, rescues endomesoderm specification in eggs depleted of the vegetal cortex. Knockdown experiments in whole embryos show that absence of Wnt6 produces embryos that lack endoderm, but those embryos continue to express a number of mesoderm markers. Thus, maternal Wnt6 plus a localized vegetal cortical molecule, possibly Dsh, is necessary for endoderm specification; this has been verified in two species of sea urchin. The data also show that Wnt6 is only one of what are likely to be multiple components that are necessary for activation of the entire endomesoderm gene regulatory network.     

Gene Regulatory Network Interactions in Sea Urchin Endomesoderm Induction

A major goal of contemporary studies of embryonic development is to understand large sets of regulatory changes that accompany the phenomenon of embryonic induction. The highly resolved sea urchin pregastrular endomesoderm–gene regulatory network (EM-GRN) provides a unique framework to study the global regulatory interactions underlying endomesoderm induction. Vegetal micromeres of the sea urchin embryo constitute a classic endomesoderm signaling center, whose potential to induce archenteron formation from presumptive ectoderm was demonstrated almost a century ago. In this work, we ectopically activate the primary mesenchyme cell–GRN (PMC-GRN) that operates in micromere progeny by misexpressing the micromere determinant Pmar1 and identify the responding EM-GRN that is induced in animal blastomeres. Using localized loss-of -function analyses in conjunction with expression of endo16, the molecular definition of micromere-dependent endomesoderm specification, we show that the TGFβ cytokine, ActivinB, is an essential component of this induction in blastomeres that emit this signal, as well as in cells that respond to it. We report that normal pregastrular endomesoderm specification requires activation of the Pmar1-inducible subset of the EM-GRN by the same cytokine, strongly suggesting that early micromere-mediated endomesoderm specification, which regulates timely gastrulation in the sea urchin embryo, is also ActivinB dependent. This study unexpectedly uncovers the existence of an additional uncharacterized micromere signal to endomesoderm progenitors, significantly revising existing models. In one of the first network-level characterizations of an intercellular inductive phenomenon, we describe an important in vivo model of the requirement of ActivinB signaling in the earliest steps of embryonic endomesoderm progenitor specification.     

Evolution of the bone gene regulatory network

Current fossil, embryological and genetic data shed light on the evolution of the gene regulatory network (GRN) governing bone formation. The key proteins and genes involved in skeletogenesis are well accepted. We discuss when these essential components of the GRN evolved and propose that the Runx genes, master regulators of skeletogenesis, functioned in early cartilages well before they were co-opted to function in the making of bone. Two rounds of whole genome duplication, together with additional tandem gene duplications, created a genetic substrate for segregation of one GRN into several networks regulating the related tissues of cartilage, bone, enamel, and dentin. During this segregation, Runx2 assumed its position at the top of the bone GRN, and Sox9 was excluded from bone, retaining its ancient role in cartilage.     

Dynamic modeling of gene expression in prokaryotes: application to glucose-lactose diauxie in Escherichia coli

Coexpression of genes or, more generally, similarity in the expression profiles poses an unsurmountable obstacle to inferring the gene regulatory network (GRN) based solely on data from DNA microarray time series. Clustering of genes with similar expression profiles allows for a course-grained view of the GRN and a probabilistic determination of the connectivity among the clusters. We present a model for the temporal evolution of a gene cluster network which takes into account interactions of gene products with genes and, through a non-constant degradation rate, with other gene products. The number of model parameters is reduced by using polynomial functions to interpolate temporal data points. In this manner, the task of parameter estimation is reduced to a system of linear algebraic equations, thus making the computation time shorter by orders of magnitude. To eliminate irrelevant networks, we test each GRN for stability with respect to parameter variations, and impose restrictions on its behavior near the steady state. We apply our model and methods to DNA microarray time series' data collected on Escherichia coli during glucose-lactose diauxie and infer the most probable cluster network for different phases of the experiment. ELECTRONIC SUPPLEMENTARY MATERIAL: The online version of this article (doi:10.1007/s11693-011-9079-2) contains supplementary material, which is available to authorized users.     

From Genes to Flower Patterns and Evolution: Dynamic Models of Gene Regulatory Networks

Genes and proteins form complex dynamical systems or gene regulatory networks (GRN) that can reach several steady states (attractors). These may be associated with distinct cell types. In plants, the ABC combinatorial model establishes the necessary gene combinations for floral organ cell specification. We have developed dynamic gene regulatory network (GRN) models to understand how the combinatorial selection of gene activity is established during floral organ primordia specification as a result of the concerted action of ABC and non-ABC genes. Our analyses have shown that the floral organ specification GRN reaches six attractors with gene configurations observed in primordial cell types during early stages of flower development and four that correspond to regions of the inflorescence meristem. This suggests that it is the overall GRN dynamics rather than precise signals that underlie the ABC model. Furthermore, our analyses suggest that the steady states of the GRN are robust to random alterations of the logical functions that define the gene interactions. Here we have updated the GRN model and have systematically altered the outputs of all the logical functions and addressed in which cases the original attractors are recovered. We then reduced the original three-state GRN to a two-state (Boolean) GRN and performed the same systematic perturbation analysis. Interestingly, the Boolean GRN reaches the same number and type of attractors as reached by the three-state GRN, and it responds to perturbations in a qualitatively identical manner as the original GRN. These results suggest that a Boolean model is sufficient to capture the dynamical features of the floral network and provide additional support for the robustness of the floral GRN. These findings further support that the GRN model provides a dynamical explanation for the ABC model and that the floral GRN robustness could be behind the widespread conservation of the floral plan among eudicotyledoneous plants. Other aspects of evolution of flower organ arrangement and ABC gene expression patterns are discussed in the context of the approach proposed here. álvaro Chaos, Max Aldana and Elena Alvarez-Buylla contributed equally to this work.     

Gene regulatory network controlling embryonic specification in the sea urchin

The current state of the gene regulatory network for endomesoderm specification in sea urchin embryos is reviewed. The network was experimentally defined, and is presented as a predictive map of cis-regulatory inputs and functional regulatory gene interconnections (updated versions of the network and most of the underlying data are at ). The network illuminates the 'whys' of many aspects of zygotic control in early sea urchin development, both spatial and temporal. The network includes almost 50 genes, and these are organized in subcircuits, each of which executes a particular regulatory function.