Brn1/2/4, the predicted midgut regulator of the endo16 gene of the sea urchin embryo

Single-cell transcriptomics reveals evolutionary reconfiguration of embryonic cell fate specification in the sea urchin Heliocidaris erythrogramma

Altered regulatory interactions during development likely underlie a large fraction of phenotypic diversity within and between species, yet identifying specific evolutionary changes remains challenging. Analysis of single-cell developmental transcriptomes from multiple species provides a powerful framework for unbiased identification of evolutionary changes in developmental mechanisms. Here, we leverage a "natural experiment" in developmental evolution in sea urchins, where a major life history switch recently evolved in the lineage leading to Heliocidaris erythrogramma, precipitating extensive changes in early development. Comparative analyses of scRNA-seq developmental time courses from H. erythrogramma and Lytechinus variegatus (representing the derived and ancestral states respectively) reveals numerous evolutionary changes in embryonic patterning. The earliest cell fate specification events, and the primary signaling center are co-localized in the ancestral dGRN but remarkably, in H. erythrogramma they are spatially and temporally separate. Fate specification and differentiation are delayed in most embryonic cell lineages, although in some cases, these processes are conserved or even accelerated. Comparative analysis of regulator-target gene co-expression is consistent with many specific interactions being preserved but delayed in H. erythrogramma, while some otherwise widely conserved interactions have likely been lost. Finally, specific patterning events are directly correlated with evolutionary changes in larval morphology, suggesting that they are directly tied to the life history shift. Together, these findings demonstrate that comparative scRNA-seq developmental time courses can reveal a diverse set of evolutionary changes in embryonic patterning and provide an efficient way to identify likely candidate regulatory interactions for subsequent experimental validation.

Evolutionary Changes in the Chromatin Landscape Contribute to Reorganization of a Developmental Gene Network During Rapid Life History Evolution in Sea Urchins

Chromatin configuration is highly dynamic during embryonic development in animals, exerting an important point of control in transcriptional regulation. Yet there exists remarkably little information about the role of evolutionary changes in chromatin configuration to the evolution of gene expression and organismal traits. Genome-wide assays of chromatin configuration, coupled with whole-genome alignments, can help address this gap in knowledge in several ways. In this study we present a comparative analysis of regulatory element sequences and accessibility throughout embryogenesis in three sea urchin species with divergent life histories: a lecithotroph Heliocidaris erythrogramma, a closely related planktotroph H. tuberculata, and a distantly related planktotroph Lytechinus variegatus. We identified distinct epigenetic and mutational signatures of evolutionary modifications to the function of putative cis-regulatory elements in H. erythrogramma that have accumulated nonuniformly throughout the genome, suggesting selection, rather than drift, underlies many modifications associated with the derived life history. Specifically, regulatory elements composing the sea urchin developmental gene regulatory network are enriched for signatures of positive selection and accessibility changes which may function to alter binding affinity and access of developmental transcription factors to these sites. Furthermore, regulatory element changes often correlate with divergent expression patterns of genes involved in cell type specification, morphogenesis, and development of other derived traits, suggesting these evolutionary modifications have been consequential for phenotypic evolution in H. erythrogramma. Collectively, our results demonstrate that selective pressures imposed by changes in developmental life history rapidly reshape the cis-regulatory landscape of core developmental genes to generate novel traits and embryonic programs.

cis-Regulatory analysis for later phase of anterior neuroectoderm-specific foxQ2 expression in sea urchin embryos

The specification of anterior neuroectoderm is controlled by a highly conserved molecular mechanism in bilaterians. A forkhead family gene, foxQ2, is known to be one of the pivotal regulators, which is zygotically expressed in this region during embryogenesis of a broad range of bilaterians. However, what controls the expression of this essential factor has remained unclear to date. To reveal the regulatory mechanism of foxQ2, we performed cis-regulatory analysis of two foxQ2 genes, foxQ2a and foxQ2b, in a sea urchin Hemicentrotus pulcherrimus. In sea urchin embryos, foxQ2 is initially expressed in the entire animal hemisphere and subsequently shows narrower expression restricted to the anterior pole region. In this study, as a first step to understand the foxQ2 regulation, we focused on the later restricted expression and analyzed the upstream cis-regulatory sequences of foxQ2a and foxQ2b by using the constructs fused to short half-life green fluorescent protein. Based on deletion and mutation analyses of both foxQ2, we identified each of the five regulatory sequences, which were 4-9 bp long. Neither of the regulatory sequences contains any motifs for ectopic activation or spatial repression, suggesting that later mRNA localization is regulated in situ. We also suggest that the three amino acid loop extension-class homeobox gene Meis is involved in the maintenance of foxQ2b, the expression of which is dominant during embryogenesis.

Neurogenesis in the sea urchin embryo is initiated uniquely in three domains

Many marine larvae begin feeding within a day of fertilization, thus requiring rapid development of a nervous system to coordinate feeding activities. Here, we examine the patterning and specification of early neurogenesis in sea urchin embryos. Lineage analysis indicates that neurons arise locally in three regions of the embryo. Perturbation analyses showed that when patterning is disrupted, neurogenesis in the three regions is differentially affected, indicating distinct patterning requirements for each neural domain. Six transcription factors that function during proneural specification were identified and studied in detail. Perturbations of these proneural transcription factors showed that specification occurs differently in each neural domain prior to the Delta-Notch restriction signal. Though gene regulatory network state changes beyond the proneural restriction are largely unresolved, the data here show that the three neural regions already differ from each other significantly early in specification. Future studies that define the larval nervous system in the sea urchin must therefore separately characterize the three populations of neurons that enable the larva to feed, to navigate, and to move food particles through the gut.

Ca²⁺ influx-linked protein kinase C activity regulates the β-catenin localization, micromere induction signalling and the oral-aboral axis formation in early sea urchin embryos

Sea urchin embryos initiate cell specifications at the 16-cell stage by forming the mesomeres, macromeres and micromeres according to the relative position of the cells in the animal-vegetal axis. The most vegetal cells, micromeres, autonomously differentiate into skeletons and induce the neighbouring macromere cells to become mesoendoderm in the β-catenin-dependent Wnt8 signalling pathway. Although the underlying molecular mechanism for this progression is largely unknown, we have previously reported that the initial events might be triggered by the Ca2+ influxes through the egg-originated L-type Ca2+ channels distributed asymmetrically along the animal-vegetal axis and through the stretch-dependent Ca2+channels expressed specifically in the micromere at the 4th cleavage. In this communication, we have examined whether one of the earliest Ca2+ targets, protein kinase C (PKC), plays a role in cell specification upstream of β-catenin. To this end, we surveyed the expression pattern of β-catenin in early embryos in the presence or absence of the specific peptide inhibitor of Hemicentrotus pulcherrimus PKC (HpPKC-I). Unlike previous knowledge, we have found that the initial nuclear entrance of β-catenin does not take place in the micromeres, but in the macromeres at the 16-cell stage. Using the HpPKC-I, we have demonstrated further that PKC not only determines cell-specific nucleation of β-catenin, but also regulates a variety of cell specification events in the early sea urchin embryos by modulating the cell adhesion structures, actin dynamics, intracellular Ca2+ signalling, and the expression of key transcription factors.

Restricted expression of karyopherin alpha mRNA in the sea urchin suggests a role in neurogenesis

Karyopherin alpha (KAP-α) proteins are critical for the transport of many molecules into the nucleus. In this study, we identified three members of the KAP-α family in the sea urchin Lytechinus variegatus and described the developmental expression of these proteins. Although many importins are assumed to have ubiquitous expression, we found that all three genes were differentially expressed. Both LvKPNA1/5/6 and LvKPNA3/4 accumulated at high levels during cleavage, exhibiting cyclic expression as cells divided. By the blastula and gastrula stages expression decreased, remaining highest in the vegetal plate and archenteron, and by the prism/pluteus stages expression was restricted to the oral surface and gut. Expression of a third KAP-α gene, LvKPNA2/7, was examined in embryos from the mesenchyme blastula to pluteus stages. LvKPNA2/7 mRNA is present in vegetal cells of the mesenchyme blastula and, during gastrulation, it is localized to the archenteron and appears in additional groups of ectodermal cells. Prism/pluteus stage embryos expressed LvKPNA2/7 in the gut and scattered distribution of transcripts in the ciliary band resembled expression patterns of neural cells. We hypothesize that LvKPNA2/7 maintains pluripotency in the neural precursors prior to activation of neural differentiation and believe that this study is an important first step in an effort to better understand the roles of importins during embryogenesis.

Genetics of gene expression responses to temperature stress in a sea urchin gene network

Stress responses play an important role in shaping species distributions and robustness to climate change. We investigated how stress responses alter the contribution of additive genetic variation to gene expression during development of the purple sea urchin, Strongylocentrotus purpuratus, under increased temperatures that model realistic climate change scenarios. We first measured gene expression responses in the embryos by RNA-seq to characterize molecular signatures of mild, chronic temperature stress in an unbiased manner. We found that an increase from 12 to 18 °C caused widespread alterations in gene expression including in genes involved in protein folding, RNA processing and development. To understand the quantitative genetic architecture of this response, we then focused on a well-characterized gene network involved in endomesoderm and ectoderm specification. Using a breeding design with wild-caught individuals, we measured genetic and gene-environment interaction effects on 72 genes within this network. We found genetic or maternal effects in 33 of these genes and that the genetic effects were correlated in the network. Fourteen network genes also responded to higher temperatures, but we found no significant genotype-environment interactions in any of the genes. This absence may be owing to an effective buffering of the temperature perturbations within the network. In support of this hypothesis, perturbations to regulatory genes did not affect the expression of the genes that they regulate. Together, these results provide novel insights into the relationship between environmental change and developmental evolution and suggest that climate change may not expose large amounts of cryptic genetic variation to selection in this species.

Maternal Oct1/2 is required for Nodal and Vg1/Univin expression during dorsal-ventral axis specification in the sea urchin embryo

The TGFβ family member Nodal is expressed early in the presumptive ventral ectoderm of the early sea urchin embryo and its activity is crucial for dorsal-ventral (D/V) axis specification. Analysis of the nodal promoter identified a number of critical binding sites for transcription factors of different families including Sox, Oct, TCF and bZIP, but in most cases the specific factors that regulate nodal expression are not known. In this study, we report that the maternal factor Oct1/2 functions as a positive regulator of nodal and that its activity is essential for the initiation of nodal expression. Inhibition of Oct1/2 mRNA translation produced embryos with severe axial defects similar to those observed following inhibition of Nodal function. We show that perturbing Oct1/2 function specifically disrupted specification of the ventral and dorsal ectodermal regions and that these effects were caused by the failure of nodal to be expressed early in development. Furthermore, we identified the key gene vg1/univin, which is also necessary for nodal expression, as an additional factor that was completely dependent on Oct1/2 for its zygotic expression. These data demonstrate that the maternal Oct1/2 protein plays an early and essential role in D/V axis specification by initiating the expression of nodal and vg1/univin, two genes that act at the top of the D/V ectoderm gene regulatory network.

Monte Carlo analysis of an ODE Model of the Sea Urchin Endomesoderm Network

BACKGROUND

Gene Regulatory Networks (GRNs) control the differentiation, specification and function of cells at the genomic level. The levels of interactions within large GRNs are of enormous depth and complexity. Details about many GRNs are emerging, but in most cases it is unknown to what extent they control a given process, i.e. the grade of completeness is uncertain. This uncertainty stems from limited experimental data, which is the main bottleneck for creating detailed dynamical models of cellular processes. Parameter estimation for each node is often infeasible for very large GRNs. We propose a method, based on random parameter estimations through Monte-Carlo simulations to measure completeness grades of GRNs.

RESULTS

We developed a heuristic to assess the completeness of large GRNs, using ODE simulations under different conditions and randomly sampled parameter sets to detect parameter-invariant effects of perturbations. To test this heuristic, we constructed the first ODE model of the whole sea urchin endomesoderm GRN, one of the best studied large GRNs. We find that nearly 48% of the parameter-invariant effects correspond with experimental data, which is 65% of the expected optimal agreement obtained from a submodel for which kinetic parameters were estimated and used for simulations. Randomized versions of the model reproduce only 23.5% of the experimental data.

CONCLUSION

The method described in this paper enables an evaluation of network topologies of GRNs without requiring any parameter values. The benefit of this method is exemplified in the first mathematical analysis of the complete Endomesoderm Network Model. The predictions we provide deliver candidate nodes in the network that are likely to be erroneous or miss unknown connections, which may need additional experiments to improve the network topology. This mathematical model can serve as a scaffold for detailed and more realistic models. We propose that our method can be used to assess a completeness grade of any GRN. This could be especially useful for GRNs involved in human diseases, where often the amount of connectivity is unknown and/or many genes/interactions are missing.

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.

Two ParaHox genes, SpLox and SpCdx, interact to partition the posterior endoderm in the formation of a functional gut

We report the characterization of the ortholog of the Xenopus XlHbox8 ParaHox gene from the sea urchin Strongylocentrotus purpuratus, SpLox. It is expressed during embryogenesis, first appearing at late gastrulation in the posterior-most region of the endodermal tube, becoming progressively restricted to the constriction between the mid- and hindgut. The physiological effects of the absence of the activity of this gene have been analyzed through knockdown experiments using gene-specific morpholino antisense oligonucleotides. We show that blocking the translation of the SpLox mRNA reduces the capacity of the digestive tract to process food, as well as eliminating the morphological constriction normally present between the mid- and hindgut. Genetic interactions of the SpLox gene are revealed by the analysis of the expression of a set of genes involved in endoderm specification. Two such interactions have been analyzed in more detail: one involving the midgut marker gene Endo16, and another involving the other endodermally expressed ParaHox gene, SpCdx. We find that SpLox is able to bind Endo16 cis-regulatory DNA, suggesting direct repression of Endo16 expression in presumptive hindgut territories. More significantly, we provide the first evidence of interaction between ParaHox genes in establishing hindgut identity, and present a model of gene regulation involving a negative-feedback loop.

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.

Sequential logic model deciphers dynamic transcriptional control of gene expressions

Background

Cellular signaling involves a sequence of events from ligand binding to membrane receptors through transcription factors activation and the induction of mRNA expression. The transcriptional-regulatory system plays a pivotal role in the control of gene expression. A novel computational approach to the study of gene regulation circuits is presented here.

Methodology

Based on the concept of finite state machine, which provides a discrete view of gene regulation, a novel sequential logic model (SLM) is developed to decipher control mechanisms of dynamic transcriptional regulation of gene expressions. The SLM technique is also used to systematically analyze the dynamic function of transcriptional inputs, the dependency and cooperativity, such as synergy effect, among the binding sites with respect to when, how much and how fast the gene of interest is expressed.

Principal Findings

SLM is verified by a set of well studied expression data on endo16 of Strongylocentrotus purpuratus (sea urchin) during the embryonic midgut development. A dynamic regulatory mechanism for endo16 expression controlled by three binding sites, UI, R and Otx is identified and demonstrated to be consistent with experimental findings. Furthermore, we show that during transition from specification to differentiation in wild type endo16 expression profile, SLM reveals three binary activities are not sufficient to explain the transcriptional regulation of endo16 expression and additional activities of binding sites are required. Further analyses suggest detailed mechanism of R switch activity where indirect dependency occurs in between UI activity and R switch during specification to differentiation stage.

Conclusions/Significance

The sequential logic formalism allows for a simplification of regulation network dynamics going from a continuous to a discrete representation of gene activation in time. In effect our SLM is non-parametric and model-independent, yet providing rich biological insight. The demonstration of the efficacy of this approach in endo16 is a promising step for further application of the proposed method.

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.

Gene regulation: gene control network in development

Controlling the differential expression of many thousands of genes is the most fundamental task of a developing organism. It requires an enormous computational device that has the capacity to process in parallel a vast number of regulatory inputs in the various cells of the embryo and come out with regulatory outputs that are tissue specific. The regulatory genome constitutes this computational device, comprising many thousands of processing units in the form of cis-regulatory modules. The interconnected cis-regulatory modules that control regulatory gene expression create a network that is the underlying mechanism of specification. In this review we use the gene regulatory network that governs endomesoderm specification in the sea urchin embryo to demonstrate the salient features of developmental gene regulatory networks and illustrate the information processing that is done by the regulatory sequences.

Transcriptional regulation of AQP-8, a Caenorhabditis elegans aquaporin exclusively expressed in the excretory system, by the POU homeobox transcription factor CEH-6

Due to the ever changing environmental conditions in soil, regulation of osmotic homeostasis in the soil-dwelling nematode Caenorhabditis elegans is critical. AQP-8 is a C. elegans aquaporin that is expressed in the excretory cell, a renal equivalent tissue, where the protein participates in maintaining water balance. To better understand the regulation of AQP-8, we undertook a promoter analysis to identify the aqp-8 cis-regulatory elements. Using progressive 5' deletions of upstream sequence, we have mapped an essential regulatory region to roughly 300 bp upstream of the translational start site of aqp-8. Analysis of this region revealed a sequence corresponding to a known DNA functional element (octamer motif), which interacts with POU homeobox transcription factors. Phylogenetic footprinting showed that this site is perfectly conserved in four nematode species. The octamer site's function was further confirmed by deletion analyses, mutagenesis, functional studies, and electrophoretic mobility shift assays. Of the three POU homeobox proteins encoded in the C. elegans genome, CEH-6 is the only member that is expressed in the excretory cell. We show that expression of AQP-8 is regulated by CEH-6 by performing RNA interference experiments. CEH-6's mammalian ortholog, Brn1, is expressed both in the kidney and the central nervous system and binds to the same octamer consensus binding site to drive gene expression. These parallels in transcriptional control between Brn1 and CEH-6 suggest that C. elegans may well be an appropriate model for determining gene-regulatory networks in the developing vertebrate kidney.

Hindgut specification and cell-adhesion functions of Sphox11/13b in the endoderm of the sea urchin embryo

Sphox11/13b is one of the two hox genes of Strongylocentrotus purpuratus expressed in the embryo. Its dynamic pattern of expression begins during gastrulation, when the transcripts are transiently located in a ring of cells at the edge of the blastopore. After gastrulation, expression is restricted to the anus-hindgut region at the boundary between the ectoderm and the endoderm. The phenotype that results when translation of Sphox11/13b mRNA is knocked down by treatment with morpholino antisense oligonucleotides (MASO) suggests that this gene may be indirectly involved in cell adhesion functions as well as in the proper differentiation of the midgut-hindgut and midgut-foregut sphincters. The MASO experiments also reveal that Sphox11/13b negatively regulates several downstream endomesoderm genes. For some of these genes, Sphox11/13b function is required to restrict expression to the midgut by preventing ectopic expression in the hindgut. The evolutionary conservation of these functions indicates the general roles of posterior Hox genes in regulating cell-adhesion, as well as in spatial control of gene regulatory network subcircuits in the regionalizing gut.

A database of mRNA expression patterns for the sea urchin embryo

We present an initial characterization of a database that contains temporal expression profiles of sequences found in 35,282 gene predictions within the sea urchin genome. The relative RNA abundance for each sequence was determined at 5 key stages of development using high-density oligonucleotide microarrays that were hybridized with populations of polyA+ RNA sequence. These stages were two-cell, which represents maternal RNA, early blastula, the time at which major tissue territories are specified, early and late gastrula, during which important morphogenetic events occur, and the pluteus larva, which marks the culmination of pre-feeding embryogenesis. We provide evidence that the microarray reliably reports the temporal profiles for the large majority of predicted genes, as shown by comparison to data for many genes with known expression patterns. The sensitivity of this assay allows detection of mRNAs whose concentration is only several hundred copies/embryo. The temporal expression profiles indicate that 5% of the gene predictions encode mRNAs that are found only in the maternal population while 24% are embryo-specific. Further, we find that the concentration of >80% of different mRNAs is modulated by more than a factor of 3 during development. Along with the annotated sea urchin genome sequence and the whole-genome tiling array (the transcriptome, Samanta, M., Tongprasit, W., Istrrail, S., Cameron, R., Tu, Q., Davidson, E., Stolc, V., in press. A high-resolution transcriptome map of the sea urchin embryo. Science), this database proves a valuable resource for designing experiments to test the function of specific genes during development.

A genomic view of the sea urchin nervous system

The sequencing of the Strongylocentrotus purpuratus genome provides a unique opportunity to investigate the function and evolution of neural genes. The neurobiology of sea urchins is of particular interest because they have a close phylogenetic relationship with chordates, yet a distinctive pentaradiate body plan and unusual neural organization. Orthologues of transcription factors that regulate neurogenesis in other animals have been identified and several are expressed in neurogenic domains before gastrulation indicating that they may operate near the top of a conserved neural gene regulatory network. A family of genes encoding voltage-gated ion channels is present but, surprisingly, genes encoding gap junction proteins (connexins and pannexins) appear to be absent. Genes required for synapse formation and function have been identified and genes for synthesis and transport of neurotransmitters are present. There is a large family of G-protein-coupled receptors, including 874 rhodopsin-type receptors, 28 metabotropic glutamate-like receptors and a remarkably expanded group of 161 secretin receptor-like proteins. Absence of cannabinoid, lysophospholipid and melanocortin receptors indicates that this group may be unique to chordates. There are at least 37 putative G-protein-coupled peptide receptors and precursors for several neuropeptides and peptide hormones have been identified, including SALMFamides, NGFFFamide, a vasotocin-like peptide, glycoprotein hormones and insulin/insulin-like growth factors. Identification of a neurotrophin-like gene and Trk receptor in sea urchin indicates that this neural signaling system is not unique to chordates. Several hundred chemoreceptor genes have been predicted using several approaches, a number similar to that for other animals. Intriguingly, genes encoding homologues of rhodopsin, Pax6 and several other key mammalian retinal transcription factors are expressed in tube feet, suggesting tube feet function as photosensory organs. Analysis of the sea urchin genome presents a unique perspective on the evolutionary history of deuterostome nervous systems and reveals new approaches to investigate the development and neurobiology of sea urchins.

High regulatory gene use in sea urchin embryogenesis: Implications for bilaterian development and evolution

A global scan of transcription factor usage in the sea urchin embryo was carried out in the context of the Strongylocentrotus purpuratus genome sequencing project, and results from six individual studies are here considered. Transcript prevalence data were obtained for over 280 regulatory genes encoding sequence-specific transcription factors of every known family, but excluding genes encoding zinc finger proteins. This is a statistically inclusive proxy for the total "regulome" of the sea urchin genome. Close to 80% of the regulome is expressed at significant levels by the late gastrula stage. Most regulatory genes must be used repeatedly for different functions as development progresses. An evolutionary implication is that animal complexity at the stage when the regulome first evolved was far simpler than even the last common bilaterian ancestor, and is thus of deep antiquity.

Identification and characterization of homeobox transcription factor genes in Strongylocentrotus purpuratus, and their expression in embryonic development

A set of 96 homeobox transcription factors was identified in the Strongylocentrotus purpuratus genome using permissive blast searches with a large collection of authentic homeodomain sequences from mouse, human and fly. A phylogenetic tree was constructed to compare the sea urchin homeobox gene family to those of vertebrates, with the result that with the only a few exceptions, orthologs of all vertebrate homeodomain genes were uncovered by our search. QPCR time course measurements revealed that 65% of these genes are expressed within the first 48 h of development (late gastrula). For genes displaying sufficiently high levels of transcript during the first 24 h of development (late blastula), whole mount in situ hybridization was carried out up to 48 h to determine spatial patterns of expression. The results demonstrate that homeodomain transcription factors participate in multiple and diverse developmental functions, in that they are used at a range of time points and in every territory of the developing embryo.

Expression and function of blimp1/krox, an alternatively transcribed regulatory gene of the sea urchin endomesoderm network

The blimp1/krox gene of Strongylocentrotus purpuratus, formerly krox1, encodes zinc finger transcription factors which play a central role in both early and late endomesoderm specification. Here we show that there are two alternative splice forms transcribed under the control of different regulatory regions. The blimp/krox1b form was previously unknown, and is the form expressed during cleavage, beginning 6-9 h postfertilization. This form is required for the early events of endomesoderm specification. A different splice variant, blimp1/krox1a, is expressed only from gastrula stage onward. During cleavage stages the blimp1/krox gene is expressed in the large micromeres and veg2 descendents. Soon after, it is expressed in the ring of specified mesoderm cells at the vegetal pole of the blastula. Its expression is later restricted to the blastopore region and the posterior of the invaginating archenteron, and finally to the midgut and hindgut of the pluteus larva. The expression of blimp1/krox is dynamic, and involves several distinct spatial territories. A GFP recombinant BAC was created by substituting the GFP coding sequence for that of the second exon (1b), in order to distinguish the expression pattern of the early form from that of the late form. This construct closely mimics blimp1/krox1b expression during early stages of sea urchin development. To expand our knowledge of the downstream linkages of this gene, additional experiments were carried out using antisense morpholino oligos (MASO). We confirmed previously published data that blimp1/krox autoregulates its own expression, but discovered, surprisingly, that this gene represses rather than activates itself. This negative autoregulation is restricted to the mesodermal and probably skeletogenic territories during the blastula stage, as shown by in situ hybridization analysis of MASO injected embryos. The MASO perturbation analysis also revealed blimp1/krox inputs into other genes of the endomesoderm regulatory network.

Confocal quantification of cis-regulatory reporter gene expression in living sea urchin

Quantification of GFP reporter gene expression at single cell level in living sea urchin embryos can now be accomplished by a new method of confocal laser scanning microscopy (CLSM). Eggs injected with a tissue-specific GFP reporter DNA construct were grown to gastrula stage and their fluorescence recorded as a series of contiguous Z-section slices that spanned the entire embryo. To measure the depth-dependent signal decay seen in the successive slices of an image stack, the eggs were coinjected with a freely diffusible internal fluorescent standard, rhodamine dextran. The measured rhodamine fluorescence was used to generate a computational correction for the depth-dependent loss of GFP fluorescence per slice. The intensity of GFP fluorescence was converted to the number of GFP molecules using a conversion constant derived from CLSM imaging of eggs injected with a measured quantity of GFP protein. The outcome is a validated method for accurately counting GFP molecules in given cells in reporter gene transfer experiments, as we demonstrate by use of an expression construct expressed exclusively in skeletogenic cells.

cis-regulatory processing of Notch signaling input to the sea urchin glial cells missing gene during mesoderm specification

The glial cells missing regulatory gene of Strongylocentrotus purpuratus (spgcm) was proposed earlier to be the genomic target of Delta/Notch (D/N) signaling required for specification of the mesodermal precursors of pigment cells. Here, we show that microinjection of a spgcm antisense morpholino oligonucleotide results in larvae without pigment cells. Microinjection of an mRNA encoding a dominant negative form of Suppressor of Hairless (dn-Su(H)) results in reduced levels of spgcm mRNA, disruption of mesodermal founder cell specification and failure to produce pigment cells. These results confirm that this gene is required for pigment cell specification. Three cis-regulatory modules of the spgcm gene were identified, which when incorporated in a GFP expression construct recapitulate the early expression pattern of this gene. Spatial expression of this GFP expression construct is severely disrupted by co-expression of dn-Su(H) mRNA, confirming that spgcm is a direct target of canonical N signaling mediated through Su(H) inputs. cis-perturbation analysis by mutation of consensus Su(H) sites identified a conserved motif paired-site and a lone site in the middle module that function both to drive expression in SMC precursors which receive the Delta signal and to repress expression in ectopic locations which lack this signal. While these Su(H) target sites provide the cis-regulatory architecture with the core of an N signaling transcriptional response switch, both the on and off outputs from this module require additional inputs.

Regulation of spblimp1/krox1a, an alternatively transcribed isoform expressed in midgut and hindgut of the sea urchin gastrula

The sea urchin regulatory gene Spblimp1/krox produces alternatively transcribed and spliced isoforms, 1a and 1b, which have different temporal and spatial patterns of expression. Here we describe a cis-regulatory module that controls the expression of the 1a splice form in the midgut and hindgut at the beginning of gastrulation. Conserved sequence patches revealed by a comparison of the blimp1/krox locus in Strongylocentrotus purpuratus and Lytechinus variegatus genomes were tested by gene transfer, in association with GFP or CAT reporter genes. An expression construct containing a conserved sequence patch immediately 5' of exon 1a included the transcription initiation site for blimp1/krox1a. This construct displays specific mid and hindgut expression, indicating that these are the locations of endogenous blimp1/krox1a transcription during the gastrula stage. Its sequence contains binding sites for Brn1/2/4, Otx, and Blimp1/Krox itself, as predicted in a prior regulatory network analysis.