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.

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.

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.

Sequential Response to Multiple Developmental Network Circuits Encoded in an Intronic cis-Regulatory Module of Sea Urchin hox11/13b.

Gene expression in different spatial domains is often controlled by separate cis-regulatory modules (CRMs), but regulatory states determining CRM activity are not only distinct in space, they also change continuously during developmental time. Here, we systematically analyzed the regulatory sequences controlling hox11/13b expression and identified a single CRM required throughout embryonic gut development. We show that within this CRM, distinct sets of binding sites recognizing Ets, Tcf, and homeodomain transcription factors control the dynamic spatial expression of hox11/13b in each developmental phase. Several binding sites execute multiple, sometimes contradictory, regulatory functions, depending on the temporal and spatial regulatory context. In addition, we identified a nearby second CRM operating in inter-modular AND logic with the first CRM to control hox11/13b expression in hindgut endoderm. Our results suggest a mechanism for continuous gene expression in response to changing developmental network functions that depends on sequential combinatorial regulation of individual CRMs.

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.

Implications of Developmental Gene Regulatory Networks Inside and Outside Developmental Biology.

The insight that the genomic control of developmental process is encoded in the form of gene regulatory networks has profound impacts on many areas of modern bioscience. Most importantly, it affects developmental biology itself, as it means that a causal understanding of development requires knowledge of the architecture of regulatory network interactions. Furthermore, it follows that functional changes in developmental gene regulatory networks have to be considered as a primary mechanism for evolutionary process. We here discuss some of the recent advances in gene regulatory network biology and how they have affected our current understanding of development, evolution, and regulatory genomics.

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.

Reorganization of sea urchin gene regulatory networks at least 268 million years ago as revealed by oldest fossil cidaroid echinoid.

Echinoids, or sea urchins, are rare in the Palaeozoic fossil record, and thus the details regarding the early diversification of crown group echinoids are unclear. Here we report on the earliest probable crown group echinoid from the fossil record, recovered from Permian (Roadian-Capitanian) rocks of west Texas, which has important implications for the timing of the divergence of crown group echinoids. The presence of apophyses and rigidly sutured interambulacral areas with two columns of plates indicates this species is a cidaroid echinoid. The species, Eotiaris guadalupensis, n. sp. is therefore the earliest stem group cidaroid. The occurrence of this species in Roadian strata pushes back the divergence of cidaroids and euechinoids, the clades that comprise all living echinoids, to at least 268.8 Ma, ten million years older than the previously oldest known cidaroid. Furthermore, the genomic regulation of development in echinoids is amongst the best known, and this new species informs the timing of large-scale reorganization in echinoid gene regulatory networks that occurred at the cidaroid-euechinoid divergence, indicating that these changes took place by the Roadian stage of the Permian.

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.

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.

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.

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.

Genomics, "Discovery Science," Systems Biology, and Causal Explanation: What Really Works?

Diverse and non-coherent sets of epistemological principles currently inform research in the general area of functional genomics. Here, from the personal point of view of a scientist with over half a century of immersion in hypothesis driven scientific discovery, I compare and deconstruct the ideological bases of prominent recent alternatives, such as "discovery science," some productions of the ENCODE project, and aspects of large data set systems biology. The outputs of these types of scientific enterprise qualitatively reflect their radical definitions of scientific knowledge, and of its logical requirements. Their properties emerge in high relief when contrasted (as an example) to a recent, system-wide, predictive analysis of a developmental regulatory apparatus that was instead based directly on hypothesis-driven experimental tests of mechanism.

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.

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.

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.

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.

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.