Cis-regulatory logic in the endo16 gene: switching from a specification to a differentiation mode of control

Molecular patterning along the sea urchin animal-vegetal axis

The molecular regulatory mechanisms underlying primary axis formation during sea urchin development have recently been identified. Two opposing maternally inherited systems, one animalizing and one vegetalizing, set up the animal-vegetal (A-V) axis. The vegetal system relies in part on the Wnt-beta-catenin-Tcf/Lef signaling pathway and the animal system is based on a cohort of animalizing transcription factors that includes members of the Ets and Sox classes. The two systems autonomously define three zones of cell-type specification along the A-V axis. The vegetalmost zone gives rise to the skeletogenic mesenchyme lineage; the animalmost zone gives rise to ectoderm; and the zone in which the two systems overlap generates endoderm, secondary mesenchyme, and ectoderm. Patterning along the A-V also depends on cellular interactions involving Wnt, Notch, and BMP signaling. We discuss how these systems impact the formation of the second axis, the oral-aboral axis; how they connect to later developmental events; and how they lead to cell-type-specific gene expression via cis-regulatory networks associated with transcriptional control regions. We also discuss how these systems may confer on the embryo its spectacular regulatory capacity to replace missing parts.

Modeling DNA sequence-based cis-regulatory gene networks

Gene network analysis requires computationally based models which represent the functional architecture of regulatory interactions, and which provide directly testable predictions. The type of model that is useful is constrained by the particular features of developmentally active cis-regulatory systems. These systems function by processing diverse regulatory inputs, generating novel regulatory outputs. A computational model which explicitly accommodates this basic concept was developed earlier for the cis-regulatory system of the endo16 gene of the sea urchin. This model represents the genetically mandated logic functions that the system executes, but also shows how time-varying kinetic inputs are processed in different circumstances into particular kinetic outputs. The same basic design features can be utilized to construct models that connect the large number of cis-regulatory elements constituting developmental gene networks. The ultimate aim of the network models discussed here is to represent the regulatory relationships among the genomic control systems of the genes in the network, and to state their functional meaning. The target site sequences of the cis-regulatory elements of these genes constitute the physical basis of the network architecture. Useful models for developmental regulatory networks must represent the genetic logic by which the system operates, but must also be capable of explaining the real time dynamics of cis-regulatory response as kinetic input and output data become available. Most importantly, however, such models must display in a direct and transparent manner fundamental network design features such as intra- and intercellular feedback circuitry; the sources of parallel inputs into each cis-regulatory element; gene battery organization; and use of repressive spatial inputs in specification and boundary formation. Successful network models lead to direct tests of key architectural features by targeted cis-regulatory analysis.

A provisional regulatory gene network for specification of endomesoderm in the sea urchin embryo

We present the current form of a provisional DNA sequence-based regulatory gene network that explains in outline how endomesodermal specification in the sea urchin embryo is controlled. The model of the network is in a continuous process of revision and growth as new genes are added and new experimental results become available; see http://www.its.caltech.edu/~mirsky/endomeso.htm (End-mes Gene Network Update) for the latest version. The network contains over 40 genes at present, many newly uncovered in the course of this work, and most encoding DNA-binding transcriptional regulatory factors. The architecture of the network was approached initially by construction of a logic model that integrated the extensive experimental evidence now available on endomesoderm specification. The internal linkages between genes in the network have been determined functionally, by measurement of the effects of regulatory perturbations on the expression of all relevant genes in the network. Five kinds of perturbation have been applied: (1) use of morpholino antisense oligonucleotides targeted to many of the key regulatory genes in the network; (2) transformation of other regulatory factors into dominant repressors by construction of Engrailed repressor domain fusions; (3) ectopic expression of given regulatory factors, from genetic expression constructs and from injected mRNAs; (4) blockade of the beta-catenin/Tcf pathway by introduction of mRNA encoding the intracellular domain of cadherin; and (5) blockade of the Notch signaling pathway by introduction of mRNA encoding the extracellular domain of the Notch receptor. The network model predicts the cis-regulatory inputs that link each gene into the network. Therefore, its architecture is testable by cis-regulatory analysis. Strongylocentrotus purpuratus and Lytechinus variegatus genomic BAC recombinants that include a large number of the genes in the network have been sequenced and annotated. Tests of the cis-regulatory predictions of the model are greatly facilitated by interspecific computational sequence comparison, which affords a rapid identification of likely cis-regulatory elements in advance of experimental analysis. The network specifies genomically encoded regulatory processes between early cleavage and gastrula stages. These control the specification of the micromere lineage and of the initial veg(2) endomesodermal domain; the blastula-stage separation of the central veg(2) mesodermal domain (i.e., the secondary mesenchyme progenitor field) from the peripheral veg(2) endodermal domain; the stabilization of specification state within these domains; and activation of some downstream differentiation genes. Each of the temporal-spatial phases of specification is represented in a subelement of the network model, that treats regulatory events within the relevant embryonic nuclei at particular stages.

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

We demonstrate that interspecific sequence conservation can provide a systematic guide to the identification of functional cis-regulatory elements within a large expanse of genomic DNA. The test was carried out on the otx gene of Strongylocentrotus purpuratus. This gene plays a major role in the gene regulatory network that underlies endomesoderm specification in the embryo. The cis-regulatory organization of the otx gene is expected to be complex, because the gene has three different start sites (X. Li, C.-K. Chuang, C.-A. Mao, L. M. Angerer, and W. H. Klein, 1997, Dev. Biol. 187, 253-266), and it is expressed in many different spatial domains of the embryo. BAC recombinants containing the otx gene were isolated from Strongylocentrotus purpuratus and Lytechinus variegatus libraries, and the ordered sequence of these BACs was obtained and annotated. Sixty kilobases of DNA flanking the gene, and included in the BAC sequence from both species, were scanned computationally for short conserved sequence elements. For this purpose, we used a newly constructed software package assembled in our laboratory, "FamilyRelations." This tool allows detection of sequence similarities above a chosen criterion within sliding windows set at 20-50 bp. Seventeen partially conserved regions, most a few hundred base pairs long, were amplified from the S. purpuratus BAC DNA by PCR, inserted in an expression vector driving a CAT reporter, and tested for cis-regulatory activity by injection into fertilized S. purpuratus eggs. The regulatory activity of these constructs was assessed by whole-mount in situ hybridization (WMISH) using a probe against CAT mRNA. Of the 17 constructs, 11 constructs displayed spatially restricted regulatory activity, and 6 were inactive in this test. The domains within which the cis-regulatory constructs were expressed are approximately consistent with results from a WMISH study on otx expression in the embryo, in which we used probes specific for the mRNAs generated from each of the three transcription start sites. Four separate cis-regulatory elements that specifically produce endomesodermal expression were identified, as well as ubiquitously active elements, and ectoderm-specific elements. We confirm predictions from other work with respect to target sites for specific transcription factors within the elements that express in the endoderm.

Positive and negative feedback: striking a balance between necessary antagonists

Most biological regulation systems comprise feedback circuits as crucial components. Negative feedback circuits have been well understood for a very long time; indeed, their understanding has been the basis for the engineering of cybernetic machines exhibiting stable behaviour. The importance of positive feedback circuits, considered as "vicious circles", has however been underestimated. In this article, we give a demonstration based on degree theory for vector fields of the conjecture, made by René Thomas, that the presence of positive feedback circuits is a necessary condition for autonomous differential systems, covering a wide class of biologically relevant systems, to possess multiple steady states. We also show ways to derive constraints on the weights of positive and negative feedback circuits. These qualitative and quantitative results provide, respectively, structural constraints (i.e. related to the interaction graph) and numerical constraints (i.e. related to the magnitudes of the interactions) on systems exhibiting complex behaviours, and should make it easier to reverse-engineer the interaction networks animating those systems on the basis of partial, sometimes unreliable, experimental data. We illustrate these concepts on a model multistable switch, in the context of cellular differentiation, showing a requirement for sufficient cooperativity. Further developments are expected in the discovery and modelling of regulatory networks in general, and in the interpretation of bio-array hybridization and proteomics experiments in particular.

Transcriptional regulation: a genomic overview

The availability of the Arabidopsis thaliana genome sequence allows a comprehensive analysis of transcriptional regulation in plants using novel genomic approaches and methodologies. Such a genomic view of transcription first necessitates the compilation of lists of elements. Transcription factors are the most numerous of the different types of proteins involved in transcription in eukaryotes, and the Arabidopsis genome codes for more than 1,500 of them, or approximately 6% of its total number of genes. A genome-wide comparison of transcription factors across the three eukaryotic kingdoms reveals the evolutionary generation of diversity in the components of the regulatory machinery of transcription. However, as illustrated by Arabidopsis, transcription in plants follows similar basic principles and logic to those in animals and fungi. A global view and understanding of transcription at a cellular and organismal level requires the characterization of the Arabidopsis transcriptome and promoterome, as well as of the interactome, the localizome, and the phenome of the proteins involved in transcription.

A genomic regulatory network for development

Development of the body plan is controlled by large networks of regulatory genes. A gene regulatory network that controls the specification of endoderm and mesoderm in the sea urchin embryo is summarized here. The network was derived from large-scale perturbation analyses, in combination with computational methodologies, genomic data, cis-regulatory analysis, and molecular embryology. The network contains over 40 genes at present, and each node can be directly verified at the DNA sequence level by cis-regulatory analysis. Its architecture reveals specific and general aspects of development, such as how given cells generate their ordained fates in the embryo and why the process moves inexorably forward in developmental time.

A new approach to decoding life: systems biology

Systems biology studies biological systems by systematically perturbing them (biologically, genetically, or chemically); monitoring the gene, protein, and informational pathway responses; integrating these data; and ultimately, formulating mathematical models that describe the structure of the system and its response to individual perturbations. The emergence of systems biology is described, as are several examples of specific systems approaches.

The smooth muscle myosin heavy chain gene exhibits smooth muscle subtype-selective modular regulation in vivo

Previous studies in our laboratory demonstrated that the transgene consisting of the -4.2 to +11.6 kilobase (kb) region of the smooth muscle (SM) myosin heavy chain (MHC) gene was expressed in virtually all SM tissue types in vivo in transgenic mice and that the multiple CArG elements within this region were differentially required in SMC subtypes, implying that the SM-MHC gene was controlled by multiple transcriptional regulatory modules. To investigate this hypothesis, we analyzed specific regulatory regions within the SM-MHC -4.2 to +11.6 kb region by a combination of deletion analyses of various SM-MHC transgenes as well as by DNaseI hypersensitivity assays and in vivo footprinting in intact SMC tissues. The results showed that SM-MHC transgene expression depended on a large number of required regulatory modules that were widely spread over the -4.2 to +11.6 region. Moreover, the results revealed several unexpected novel features of regulation of the SM-MHC gene including: 1) unique combinations of regulatory modules were required for SM-MHC expression in different SMC-subtypes; 2) repressor modules as well as activator modules were both critical for SMC specificity of the gene; 3) certain modules were required in certain contexts but were dispensable in others within a given SMC-subtype (i.e. the net activity of the module was determined by interaction between modules not simply by the sum of module activities); and 4) we identified a highly conserved 200-base pair transcriptional regulatory module at +8 kb that was required in the large arteries but dispensable in the coronary arteries and airways in transgenic mice and contained multiple potential cis-elements that were occupied by nuclear proteins in the intact aorta based on in vivo footprinting. Taken together, the results suggest a model of complex modular control of expression of the SM-MHC gene that varies between SMC subtypes. Moreover, the studies establish the possibility of designing derivatives of the SM-MHC promoter that might be used for targeting gene expression to specific SMC subtypes in vivo.

Evidence for a mesodermal embryonic regulator of the sea urchin CyIIa gene

The CyIIa gene of the sea urchin embryo is a model for study of cis-regulation downstream of cell-type specification, as CyIIa transcription follows the specification and initial differentiation of the embryonic domains in which it is expressed. These are the skeletogenic and secondary mesenchyme and gut. We carried out a detailed structural and functional analysis of a cis-regulatory region of this gene, extending 780 bp upstream and 125 bp downstream of the transcription start site, that had been shown earlier to reproduce faithfully the complex and dynamic CyIIa pattern of expression. This analysis revealed that the overall pattern of expression of the CyIIa gene appears to be governed mainly by two independent sets of DNA elements, which are target sites for specific proteins present in blastula-stage nuclear extract. One type of element, which controls a dynamic program of expression in both skeletogenic and secondary mesenchyme cells, contains the consensus-binding site for a member of the ets transcription factor family. The other, which is responsible for the terminal or permanent phase of CyIIa expression in the gut, shares homologies with the late module of the endoderm-specific Endo16 gene (endo16 Module B). Oligonucleotides containing replicas of these two target sites fused upstream of a sea urchin basal promoter are sufficient to confer accurate mesenchyme and late gut expression of an injected GFP construct. The finding of a single protein target site that recapitulates CyIIa expression in both primary and secondary mesenchyme cells suggests the existence of a pan-mesodermal gene expression program in the sea urchin embryo.