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.

New computational approaches for analysis of cis-regulatory networks.

The investigation and modeling of gene regulatory networks requires computational tools specific to the task. We present several locally developed software tools that have been used in support of our ongoing research into the embryogenesis of the sea urchin. These tools are especially well suited to iterative refinement of models through experimental and computational investigation. They include: BioArray, a macroarray spot processing program; SUGAR, a system to display and correlate large-BAC sequence analyses; SeqComp and FamilyRelations, programs for comparative sequence analysis; and NetBuilder, an environment for creating and analyzing models of gene networks. We also present an overview of the process used to build our model of the Strongylocentrotus purpuratus endomesoderm gene network. Several of the tools discussed in this paper are still in active development and some are available as open source.

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.

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.