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.

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.

Genes that are upregulated by LiCl treatment of sea urchin embryos and/or downregulated by injection into the egg of mRNA encoding an internal fragment of cadherin (Cad) were detected in a differential macroarray screen. The method was that recently described by J. P. Rast et al. (2000, Dev. Biol. 228, 270-296). Almost 10(5) clones from a 12-h cDNA library were screened. Measurements on internal standards showed that the screening procedure was sufficiently sensitive to afford detection of differentially expressed mRNAs of the most rare class, those present in only a few copies per average cell. The injection of Cad mRNA, which specifically blocks nuclearization of beta-catenin, resulted in many-fold decreases in the levels of transcripts of a suite of marker genes expressed zygotically during endomesoderm specification. These measurements substantiated the use of Cad mRNA as the basis for a differential screen for discovery of new endomesodermal genes. By use of the newly developed BioArray software for analysis of macroarray screens, 1106 clones representing differentially expressed genes and yielding useful sequence were recovered. The 367 clones that gave significant BLASTX matches to known cellular proteins fell into 264 nonredundant sequence classes. Those of particular interest for this work were clones encoding DNA-binding transcription factors, signal transduction pathway components, proteases, kinases, and phosphatases. Quantitative PCR analysis of 66 such selected clones revealed that the large majority of these clones had been selected because they are upregulated by LiCl treatment, which affects the expression of a much greater diversity and number of genes than are involved in endomesoderm specification. Seven transcript species were identified that responded sharply to injection of Cad mRNA, and that are not represented in maternal mRNA. Six of those encode transcription factors. We focused on three transcription factor genes of this set that were previously unknown in sea urchin embryos. By whole-mount in situ hybridization, these genes are expressed in specific domains of the endomesodermal territory. They are: (1) Speve, an evenskipped orthologue expressed very early in all vegetal blastomeres and then gradually shifting to veg(1) derivatives by the mesenchyme blastula stage; (2) Spgcm, an orthologue of the fruit fly gene glial cells missing, which is first expressed specifically and exclusively in part of the prospective secondary mesenchyme (mesodermal) domain at late-cleavage blastula stage; and (3) Spfoxc, which is first expressed in the early blastula only in the four small micromeres, and later only expressed in that coelomic pouch which gives rise to the mesoderm of the ventral surface of the adult rudiment.

No abstract

Expression of the wnt8 gene is the key transcriptional motivator of an intercellular signaling loop which drives endomesoderm specification forward early in sea urchin embryogenesis. This gene was predicted by network perturbation analysis to be activated by inputs from the blimp1/krox gene, itself expressed zygotically in the endomesoderm during cleavage; and by a Tcf1/beta-catenin input. The implication is that zygotic expression of wnt8 is stimulated in neighboring cells by its own gene product, since reception of the Wnt8 ligand causes beta-catenin nuclearization. Here, the modular cis-regulatory system of the wnt8 gene of Strongylocentrotus purpuratus was characterized functionally, and shown to respond to blockade of both Blimp1/Krox and Tcf1/beta-catenin inputs just as does the endogenous gene. The genomic target sites for these factors were demonstrated by mutation in one of the cis-regulatory modules. The Tcf1/beta-catenin and Blimp1/Krox inputs are both necessary for normal endomesodermal expression mediated by this cis-regulatory module; thus, the genomic regulatory code underlying the predicted signaling loop thus resides in the wnt8 cis-regulatory sequence. In a second regulatory region, which initiates expression in micromere and macromere descendant cells early in cleavage, Tcf1 sites act to repress ectopic transcription in prospective ectoderm cells.