Evolution of gene regulatory networks controlling body plan development.

Evolutionary change in animal morphology results from alteration of the functional organization of the gene regulatory networks (GRNs) that control development of the body plan. A major mechanism of evolutionary change in GRN structure is alteration of cis-regulatory modules that determine regulatory gene expression. Here we consider the causes and consequences of GRN evolution. Although some GRN subcircuits are of great antiquity, other aspects are highly flexible and thus in any given genome more recent. This mosaic view of the evolution of GRN structure explains major aspects of evolutionary process, such as hierarchical phylogeny and discontinuities of paleontological change.

Conserved regulatory state expression controlled by divergent developmental gene regulatory networks in echinoids.

Evolution of the animal body plan is driven by changes in developmental gene regulatory networks (GRNs), but how networks change to control novel developmental phenotypes remains, in most cases, unresolved. Here, we address GRN evolution by comparing the endomesoderm GRN in two echinoid sea urchins, Strongylocentrotus purpuratus and Eucidaris tribuloides, with at least 268 million years of independent evolution. We first analyzed the expression of twelve transcription factors and signaling molecules of the S. purpuratus GRN in E. tribuloides embryos, showing that orthologous regulatory genes are expressed in corresponding endomesodermal cell fates in the two species. However, perturbation of regulatory genes revealed that important regulatory circuits of the S. purpuratus GRN are significantly different in E. tribuloides For example, mesodermal Delta/Notch signaling controls exclusion of alternative cell fates in E. tribuloides but controls mesoderm induction and activation of a positive feedback circuit in S. purpuratus These results indicate that the architecture of the sea urchin endomesoderm GRN evolved by extensive gain and loss of regulatory interactions between a conserved set of regulatory factors that control endomesodermal cell fate specification.

Developmental effector gene regulation: Multiplexed strategies for functional analysis.

The staggering complexity of the genome controls for developmental processes is revealed through massively parallel cis-regulatory analysis using new methods of perturbation and readout. The choice of combinations of these new methods is tailored to the system, question and resources at hand. Our focus is on issues that include the necessity or sufficiency of given cis-regulatory modules, cis-regulatory function in the normal spatial genomic context, and easily accessible high throughput and multiplexed analysis methods. In the sea urchin embryonic model, recombineered BACs offer new opportunities for consecutive modes of cis-regulatory analyses that answer these requirements, as we here demonstrate on a diverse suite of previously unstudied sea urchin effector genes expressed in skeletogenic cells. Positively active cis-regulatory modules were located in single Nanostring experiments per BAC containing the gene of interest, by application of our previously reported "barcode" tag vectors of which> 100 can be analyzed at one time. Computational analysis of DNA sequences that drive expression, based on the known skeletogenic regulatory state, then permitted effective identification of functional target site clusters. Deletion of these sub-regions from the parent BACs revealed module necessity, as simultaneous tests of the same regions in short constructs revealed sufficiency. Predicted functional inputs were then confirmed by site mutations, all generated and tested in multiplex formats. There emerged the simple conclusion that each effector gene utilizes a small subset of inputs from the skeletogenic GRN. These inputs may function to only adjust expression levels or in some cases necessary for expression. Since we know the GRN architecture upstream of the effector genes, we could then conceptually isolate and compare the wiring of the effector gene driver sub-circuits and identify the inputs whose removal abolish expression.

Assessing regulatory information in developmental gene regulatory networks.

Gene regulatory networks (GRNs) provide a transformation function between the static genomic sequence and the primary spatial specification processes operating development. The regulatory information encompassed in developmental GRNs thus goes far beyond the control of individual genes. We here address regulatory information at different levels of network organization, from single node to subcircuit to large-scale GRNs and discuss how regulatory design features such as network architecture, hierarchical organization, and cis-regulatory logic contribute to the developmental function of network circuits. Using specific subcircuits from the sea urchin endomesoderm GRN, for which both circuit design and biological function have been described, we evaluate by Boolean modeling and in silico perturbations the import of given circuit features on developmental function. The examples include subcircuits encoding positive feedback, mutual repression, and coherent feedforward, as well as signaling interaction circuitry. Within the hierarchy of the endomesoderm GRN, these subcircuits are organized in an intertwined and overlapping manner. Thus, we begin to see how regulatory information encoded at individual nodes is integrated at all levels of network organization to control developmental process.

Mapping a multiplexed zoo of mRNA expression.

In situ hybridization methods are used across the biological sciences to map mRNA expression within intact specimens. Multiplexed experiments, in which multiple target mRNAs are mapped in a single sample, are essential for studying regulatory interactions, but remain cumbersome in most model organisms. Programmable in situ amplifiers based on the mechanism of hybridization chain reaction (HCR) overcome this longstanding challenge by operating independently within a sample, enabling multiplexed experiments to be performed with an experimental timeline independent of the number of target mRNAs. To assist biologists working across a broad spectrum of organisms, we demonstrate multiplexed in situ HCR in diverse imaging settings: bacteria, whole-mount nematode larvae, whole-mount fruit fly embryos, whole-mount sea urchin embryos, whole-mount zebrafish larvae, whole-mount chicken embryos, whole-mount mouse embryos and formalin-fixed paraffin-embedded human tissue sections. In addition to straightforward multiplexing, in situ HCR enables deep sample penetration, high contrast and subcellular resolution, providing an incisive tool for the study of interlaced and overlapping expression patterns, with implications for research communities across the biological sciences.

The uncommon roles of common gene regulatory factors in the genomes of differentiating cells.

Viewed through the lens of comparative regulatory mechanisms in developmental processes, the article of Calero-Nieto et al (2014, this issue) is of particular interest. This work uncovers the causal combinatorial subtleties of the distinct enhancer occupancy profiles displayed by ten different transcription factors, which are expressed in common in two hematopoietic cell types, a stem cell-like precursor and primary mast cells.

Genome-wide assessment of differential effector gene use in embryogenesis.

Six different populations of cells were isolated by fluorescence-activated cell sorting from disaggregated late blastula- and gastrula-stage sea urchin embryos according to the regulatory states expressed in these cells, as reported by recombineered bacterial artificial chromosomes producing fluorochromes. Transcriptomes recovered from these embryonic cell populations revealed striking, early differential expression of large cohorts of effector genes. The six cell populations were presumptive pigment cells, presumptive neurogenic cells, presumptive skeletogenic cells, cells from the stomodeal region of the oral ectoderm, ciliated band cells and cells from the endoderm/ectoderm boundary that will give rise both to hindgut and to border ectoderm. Transcriptome analysis revealed that each of these domains specifically expressed several hundred effector genes at significant levels. Annotation indicated the qualitative individuality of the functional nature of each cell population, even though they were isolated from embryos only 1-2 days old. In no case was more than a tiny fraction of the transcripts enriched in one population also enriched in any other of the six populations studied. As was particularly clear in the cases of the presumptive pigment, neurogenic and skeletogenic cells, all three of which represent precociously differentiating cell types of this embryo, most specifically expressed genes of given cell types are not significantly expressed at all in the other cell types. Thus, at the effector gene level, a dramatic, cell type-specific pattern of differential gene regulation is established well before any significant embryonic morphogenesis has occurred.

Geometric control of ciliated band regulatory states in the sea urchin embryo.

The trapezoidal ciliated band (CB) of the postgastrular sea urchin embryo surrounds the oral ectoderm, separating it from adjacent embryonic territories. Once differentiated, the CB is composed of densely arranged cells bearing long cilia that endow the larva with locomotion and feeding capability. The spatial pattern from which the CB will arise is first evidenced during pregastrular stages by expression of the pioneer gene onecut. Immediately after gastrulation, the CB consists of four separate regulatory state domains, each of which expresses a unique set of transcription factors: (1) the oral apical CB, located within the apical neurogenic field; (2) the animal lateral CB, which bilaterally separates the oral from aboral ectoderm; (3) the vegetal lateral CB, which bilaterally serves as signaling centers; and (4) the vegetal oral CB, which delineates the boundary with the underlying endoderm. Remarkably, almost all of the regulatory genes specifically expressed within these domains are downregulated by interference with SoxB1 expression, implying their common activation by this factor. Here, we show how the boundaries of the CB subdomains are established, and thus ascertain the design principle by which the geometry of this unique and complex regulatory state pattern is genomically controlled. Each of these boundaries, on either side of the CB, is defined by spatially confined transcriptional repressors, the products of regulatory genes operating across the border of each subdomain. In total this requires deployment of about ten different repressors, which we identify in this work, thus exemplifying the complexity of information required for spatial regulatory organization during embryogenesis.

Evolutionary rewiring of gene regulatory network linkages at divergence of the echinoid subclasses.

Evolution of animal body plans occurs with changes in the encoded genomic programs that direct development, by alterations in the structure of encoded developmental gene-regulatory networks (GRNs). However, study of this most fundamental of evolutionary processes requires experimentally tractable, phylogenetically divergent organisms that differ morphologically while belonging to the same monophyletic clade, plus knowledge of the relevant GRNs operating in at least one of the species. These conditions are met in the divergent embryogenesis of the two extant, morphologically distinct, echinoid (sea urchin) subclasses, Euechinoidea and Cidaroidea, which diverged from a common late Paleozoic ancestor. Here we focus on striking differences in the mode of embryonic skeletogenesis in a euechinoid, the well-known model Strongylocentrotus purpuratus (Sp), vs. the cidaroid Eucidaris tribuloides (Et). At the level of descriptive embryology, skeletogenesis in Sp and Et has long been known to occur by distinct means. The complete GRN controlling this process is known for Sp. We carried out targeted functional analyses on Et skeletogenesis to identify the presence, or demonstrate the absence, of specific regulatory linkages and subcircuits key to the operation of the Sp skeletogenic GRN. Remarkably, most of the canonical design features of the Sp skeletogenic GRN that we examined are either missing or operate differently in Et. This work directly implies a dramatic reorganization of genomic regulatory circuitry concomitant with the divergence of the euechinoids, which began before the end-Permian extinction.

Sponge grade body fossil with cellular resolution dating 60 Myr before the Cambrian.

An extraordinarily well preserved, 600-million-year (Myr)-old, three-dimensionally phosphatized fossil displaying multiple independent characters of modern adult sponges has been analyzed by SEM and synchrotron X-ray tomography. The fossilized animal (Eocyathispongia qiania gen. et sp. nov.) is slightly more than 1.2 mm wide and 1.1 mm tall, is composed of hundreds of thousands of cells, and has a gross structure consisting of three adjacent hollow tubes sharing a common base. The main tube is crowned with a large open funnel, and the others end in osculum-like openings to the exterior. The external surface is densely covered with flat tile-like cells closely resembling sponge pinacocytes, and this layer is punctuated with smaller pores. A dense patch of external structures that display the form of a lawn of sponge papillae has also survived. Within the main funnel, an area where features of the inner surface are preserved displays a regular pattern of uniform pits. Many of them are surrounded individually by distinct collars, mounted in a supporting reticulum. The possibility cannot be excluded that these pits are the remains of a field of choanocytes. The character set evinced by this specimen, ranging from general anatomy to cell type, uniquely indicates that this specimen is a fossil of probable poriferan affinity. So far, we have only this single specimen, and although its organized and complex cellular structure precludes any reasonable interpretation that its origin is abiogenic, confirmation that it is indeed a fossilized sponge will clearly require discovery of additional specimens.

cis-Regulatory control of the initial neurogenic pattern of onecut gene expression in the sea urchin embryo.

Specification of the ciliated band (CB) of echinoid embryos executes three spatial functions essential for postgastrular organization. These are establishment of a band about 5 cells wide which delimits and bounds other embryonic territories; definition of a neurogenic domain within this band; and generation within it of arrays of ciliary cells that bear the special long cilia from which the structure derives its name. In Strongylocentrotus purpuratus the spatial coordinates of the future ciliated band are initially and exactly determined by the disposition of a ring of cells that transcriptionally activate the onecut homeodomain regulatory gene, beginning in blastula stage, long before the appearance of the CB per se. Thus the cis-regulatory apparatus that governs onecut expression in the blastula directly reveals the genomic sequence code by which these aspects of the spatial organization of the embryo are initially determined. We screened the entire onecut locus and its flanking region for transcriptionally active cis-regulatory elements, and by means of BAC recombineered deletions identified three separated and required cis-regulatory modules that execute different functions. The operating logic of the crucial spatial control module accounting for the spectacularly precise and beautiful early onecut expression domain depends on spatial repression. Previously predicted oral ectoderm and aboral ectoderm repressors were identified by cis-regulatory mutation as the products of goosecoid and irxa genes respectively, while the pan-ectodermal activator SoxB1 supplies a transcriptional driver function.

General approach for in vivo recovery of cell type-specific effector gene sets.

Differentially expressed, cell type-specific effector gene sets hold the key to multiple important problems in biology, from theoretical aspects of developmental gene regulatory networks (GRNs) to various practical applications. Although individual cell types of interest have been recovered by various methods and analyzed, systematic recovery of multiple cell type-specific gene sets from whole developing organisms has remained problematic. Here we describe a general methodology using the sea urchin embryo, a material of choice because of the large-scale GRNs already solved for this model system. This method utilizes the regulatory states expressed by given cells of the embryo to define cell type and includes a fluorescence activated cell sorting (FACS) procedure that results in no perturbation of transcript representation. We have extensively validated the method by spatial and qualitative analyses of the transcriptome expressed in isolated embryonic skeletogenic cells and as a consequence, generated a prototypical cell type-specific transcriptome database.

A gene regulatory network controlling the embryonic specification of endoderm.

Specification of endoderm is the prerequisite for gut formation in the embryogenesis of bilaterian organisms. Modern lineage labelling studies have shown that in the sea urchin embryo model system, descendants of the veg1 and veg2 cell lineages produce the endoderm, and that the veg2 lineage also gives rise to mesodermal cell types. It is known that Wnt/ß-catenin signalling is required for endoderm specification and Delta/Notch signalling is required for mesoderm specification. Some direct cis-regulatory targets of these signals have been found and various phenomenological patterns of gene expression have been observed in the pre-gastrular endomesoderm. However, no comprehensive, causal explanation of endoderm specification has been conceived for sea urchins, nor for any other deuterostome. Here we propose a model, on the basis of the underlying genomic control system, that provides such an explanation, built at several levels of biological organization. The hardwired core of the control system consists of the cis-regulatory apparatus of endodermal regulatory genes, which determine the relationship between the inputs to which these genes are exposed and their outputs. The architecture of the network circuitry controlling the dynamic process of endoderm specification then explains, at the system level, a sequence of developmental logic operations, which generate the biological process. The control system initiates non-interacting endodermal and mesodermal gene regulatory networks in veg2-derived cells and extinguishes the endodermal gene regulatory network in mesodermal precursors. It also generates a cross-regulatory network that specifies future anterior endoderm in veg2 descendants and institutes a distinct network specifying posterior endoderm in veg1-derived cells. The network model provides an explanatory framework that relates endoderm specification to the genomic regulatory code.

Encoding regulatory state boundaries in the pregastrular oral ectoderm of the sea urchin embryo.

By gastrulation the ectodermal territories of the sea urchin embryo have developed an unexpectedly complex spatial pattern of sharply bounded regulatory states, organized orthogonally with respect to the animal/vegetal and oral/aboral axes of the embryo. Although much is known of the gene regulatory network (GRN) linkages that generate these regulatory states, the principles by which the boundaries between them are positioned and maintained have remained undiscovered. Here we determine the encoded genomic logic responsible for the boundaries of the oral aspect of the embryo that separate endoderm from ectoderm and ectoderm from neurogenic apical plate and that delineate the several further subdivisions into which the oral ectoderm per se is partitioned. Comprehensive regulatory state maps, including all spatially expressed oral ectoderm regulatory genes, were established. The circuitry at each boundary deploys specific repressors of regulatory states across the boundary, identified in this work, plus activation by broadly expressed positive regulators. These network linkages are integrated with previously established interactions on the oral/aboral axis to generate a GRN model encompassing the 2D organization of the regulatory state pattern in the pregastrular oral ectoderm of the embryo.

Specific functions of the Wnt signaling system in gene regulatory networks throughout the early sea urchin embryo.

Wnt signaling affects cell-fate specification processes throughout embryonic development. Here we take advantage of the well-studied gene regulatory networks (GRNs) that control pregastrular sea urchin embryogenesis to reveal the gene regulatory functions of the entire Wnt-signaling system. Five wnt genes, three frizzled genes, two secreted frizzled-related protein 1 genes, and two Dickkopf genes are expressed in dynamic spatial patterns in the pregastrular embryo of Strongylocentrotus purpuratus. We present a comprehensive analysis of these genes in each embryonic domain. Total functions of the Wnt-signaling system in regulatory gene expression throughout the embryo were studied by use of the Porcupine inhibitor C59, which interferes with zygotic Wnt ligand secretion. Morpholino-mediated knockdown of each expressed Wnt ligand demonstrated that individual Wnt ligands are functionally distinct, despite their partially overlapping spatial expression. They target specific embryonic domains and affect particular regulatory genes. The sum of the effects of blocking expression of individual wnt genes is shown to equal C59 effects. Remarkably, zygotic Wnt-signaling inputs are required for only three general aspects of embryonic specification: the broad activation of endodermal GRNs, the regional specification of the immediately adjacent stripe of ectoderm, and the restriction of the apical neurogenic domain. All Wnt signaling in this pregastrular embryo is short range (and/or autocrine). Furthermore, we show that the transcriptional drivers of wnt genes execute important specification functions in the embryonic domains targeted by the ligands, thus connecting the expression and function of wnt genes by encoded cross-regulatory interactions within the specific regional GRNs.

Juvenile skeletogenesis in anciently diverged sea urchin clades.

Mechanistic understanding of evolutionary divergence in animal body plans devolves from analysis of those developmental processes that, in forms descendant from a common ancestor, are responsible for their morphological differences. The last common ancestor of the two extant subclasses of sea urchins, i.e., euechinoids and cidaroids, existed well before the Permian/Triassic extinction (252 mya). Subsequent evolutionary divergence of these clades offers in principle a rare opportunity to solve the developmental regulatory events underlying a defined evolutionary divergence process. Thus (i) there is an excellent and fairly dense (if yet incompletely analyzed) fossil record; (ii) cladistically confined features of the skeletal structures of modern euechinoid and cidaroid sea urchins are preserved in fossils of ancestral forms; (iii) euechinoids and cidaroids are among current laboratory model systems in molecular developmental biology (here Strongylocentrotus purpuratus [Sp] and Eucidaris tribuloides [Et]); (iv) skeletogenic specification in sea urchins is uncommonly well understood at the causal level of interactions of regulatory genes with one another, and with known skeletogenic effector genes, providing a ready arsenal of available molecular tools. Here we focus on differences in test and perignathic girdle skeletal morphology that distinguish all modern euechinoid from all modern cidaroid sea urchins. We demonstrate distinct canonical test and girdle morphologies in juveniles of both species by use of SEM and X-ray microtomography. Among the sharply distinct morphological features of these clades are the internal skeletal structures of the perignathic girdle to which attach homologous muscles utilized for retraction and protraction of Aristotles׳ lantern and its teeth. We demonstrate that these structures develop de novo between one and four weeks after metamorphosis. In order to study the underlying developmental processes, a method of section whole mount in situ hybridization was adapted. This method displays current gene expression in the developing test and perignathic girdle skeletal elements of both Sp and Et juveniles. Active, specific expression of the sm37 biomineralization gene in these muscle attachment structures accompanies morphogenetic development of these clade-specific features in juveniles of both species. Skeletogenesis at these clade-specific muscle attachment structures displays molecular earmarks of the well understood embryonic skeletogenic GRN: thus the upstream regulatory gene alx1 and the gene encoding the vegfR signaling receptor are both expressed at the sites where they are formed. This work opens the way to analysis of the alternative spatial specification processes that were installed at the evolutionary divergence of the two extant subclasses of sea urchins.

Emerging properties of animal gene regulatory networks.

Gene regulatory networks (GRNs) provide system level explanations of developmental and physiological functions in the terms of the genomic regulatory code. Depending on their developmental functions, GRNs differ in their degree of hierarchy, and also in the types of modular sub-circuit of which they are composed, although there is a commonly employed sub-circuit repertoire. Mathematical modelling of some types of GRN sub-circuit has deepened biological understanding of the functions they mediate. The structural organization of various kinds of GRN reflects their roles in the life process, and causally illuminates both developmental and evolutionary process.

Quantitative developmental transcriptomes of the sea urchin Strongylocentrotus purpuratus.

Development depends on the precise control of gene expression in time and space. A critical step towards understanding the global gene regulatory networks underlying development is to obtain comprehensive information on gene expression. In this study, we measured expression profiles for the entire expressed gene set during sea urchin embryonic development. We confirmed the reliability of these profiles by comparison with NanoString measurements for a subset of genes and with literature values. The data show that ~16,500 genes have been activated by the end of embryogenesis, and for half of them the transcript abundance changes more than 10-fold during development. From this genome scale expression survey, we show that complex patterns of expression by many genes underlie embryonic development, particularly during the early stages before gastrulation. An intuitive web application for data query and visualization is presented to facilitate use of this large dataset.

Gene structure in the sea urchin Strongylocentrotus purpuratus based on transcriptome analysis.

A comprehensive transcriptome analysis has been performed on protein-coding RNAs of Strongylocentrotus purpuratus, including 10 different embryonic stages, six feeding larval and metamorphosed juvenile stages, and six adult tissues. In this study, we pooled the transcriptomes from all of these sources and focused on the insights they provide for gene structure in the genome of this recently sequenced model system. The genome had initially been annotated by use of computational gene model prediction algorithms. A large fraction of these predicted genes were recovered in the transcriptome when the reads were mapped to the genome and appropriately filtered and analyzed. However, in a manually curated subset, we discovered that more than half the computational gene model predictions were imperfect, containing errors such as missing exons, prediction of nonexistent exons, erroneous intron/exon boundaries, fusion of adjacent genes, and prediction of multiple genes from single genes. The transcriptome data have been used to provide a systematic upgrade of the gene model predictions throughout the genome, very greatly improving the research usability of the genomic sequence. We have constructed new public databases that incorporate information from the transcriptome analyses. The transcript-based gene model data were used to define average structural parameters for S. purpuratus protein-coding genes. In addition, we constructed a custom sea urchin gene ontology, and assigned about 7000 different annotated transcripts to 24 functional classes. Strong correlations became evident between given functional ontology classes and structural properties, including gene size, exon number, and exon and intron size.

The evolution of hierarchical gene regulatory networks.

Comparative developmental evidence indicates that reorganizations in developmental gene regulatory networks (GRNs) underlie evolutionary changes in animal morphology, including body plans. We argue here that the nature of the evolutionary alterations that arise from regulatory changes depends on the hierarchical position of the change within a GRN. This concept cannot be accomodated by microevolutionary nor macroevolutionary theory. It will soon be possible to investigate these ideas experimentally, by assessing the effects of GRN changes on morphological evolution.