Assessing regulatory information in developmental gene regulatory networks

The Blueprint of Logical Decisions in a NF-κB Signaling System

Nearly identical cells can exhibit substantially different responses to the same stimulus that causes phenotype diversity. Such interplay between phenotype diversity and the architecture of regulatory circuits is crucial since it determines the state of a biological cell. Here, we theoretically analyze how the circuit blueprints of NF-κB in cellular environments are formed and their role in determining the cells' metabolic state. The NF-κB is a collective name for a developmental conserved family of five different transcription factors that can form homodimers or heterodimers and often promote DNA looping to reprogram the inflammatory gene response. The NF-κB controls many biological functions, including cellular differentiation, proliferation, migration, and survival. Our model shows that nuclear localization of NF-κB differentially promotes logic operations such as AND, NAND, NOR, and OR in its regulatory network. Through the quantitative thermodynamic model of transcriptional regulation and systematic variation of promoter-enhancer interaction modes, we can account for the origin of various logic gates as formed in the NF-κB system. We further show that the interconversion or switching of logic gates yielded under systematic variations of the stimuli activity and DNA looping parameters. Such computation occurs in regulatory and signaling pathways in individual cells at a molecular scale, which one can exploit to design a biomolecular computer.

Echinobase: a resource to support the echinoderm research community

Echinobase (www.echinobase.org) is a model organism knowledgebase serving as a resource for the community that studies echinoderms, a phylum of marine invertebrates that includes sea urchins and sea stars. Echinoderms have been important experimental models for over 100 years and continue to make important contributions to environmental, evolutionary, and developmental studies, including research on developmental gene regulatory networks. As a centralized resource, Echinobase hosts genomes and collects functional genomic data, reagents, literature, and other information for the community. This third-generation site is based on the Xenbase knowledgebase design and utilizes gene-centric pages to minimize the time and effort required to access genomic information. Summary gene pages display gene symbols and names, functional data, links to the JBrowse genome browser, and orthology to other organisms and reagents, and tabs from the Summary gene page contain more detailed information concerning mRNAs, proteins, diseases, and protein-protein interactions. The gene pages also display 1:1 orthologs between the fully supported species Strongylocentrotus purpuratus (purple sea urchin), Lytechinus variegatus (green sea urchin), Patiria miniata (bat star), and Acanthaster planci (crown-of-thorns sea star). JBrowse tracks are available for visualization of functional genomic data from both fully supported species and the partially supported species Anneissia japonica (feather star), Asterias rubens (sugar star), and L. pictus (painted sea urchin). Echinobase serves a vital role by providing researchers with annotated genomes including orthology, functional genomic data aligned to the genomes, and curated reagents and data. The Echinoderm Anatomical Ontology provides a framework for standardizing developmental data across the phylum, and knowledgebase content is formatted to be findable, accessible, interoperable, and reusable by the research community.

Lhx3/4 initiates a cardiopharyngeal-specific transcriptional program in response to widespread FGF signaling

Individual signaling pathways, such as fibroblast growth factors (FGFs), can regulate a plethora of inductive events. According to current paradigms, signal-dependent transcription factors (TFs), such as FGF/MapK-activated Ets family factors, partner with lineage-determining factors to achieve regulatory specificity. However, many aspects of this model have not been rigorously investigated. One key question relates to whether lineage-determining factors dictate lineage-specific responses to inductive signals or facilitate these responses in collaboration with other inputs. We utilize the chordate model Ciona robusta to investigate mechanisms generating lineage-specific induction. Previous studies in C. robusta have shown that cardiopharyngeal progenitor cells are specified through the combined activity of FGF-activated Ets1/2.b and an inferred ATTA-binding transcriptional cofactor. Here, we show that the homeobox TF Lhx3/4 serves as the lineage-determining TF that dictates cardiopharyngeal-specific transcription in response to pleiotropic FGF signaling. Targeted knockdown of Lhx3/4 leads to loss of cardiopharyngeal gene expression. Strikingly, ectopic expression of Lhx3/4 in a neuroectodermal lineage subject to FGF-dependent specification leads to ectopic cardiopharyngeal gene expression in this lineage. Furthermore, ectopic Lhx3/4 expression disrupts neural plate morphogenesis, generating aberrant cell behaviors associated with execution of incompatible morphogenetic programs. Based on these findings, we propose that combinatorial regulation by signal-dependent and lineage-determinant factors represents a generalizable, previously uncategorized regulatory subcircuit we term "cofactor-dependent induction." Integration of this subcircuit into theoretical models will facilitate accurate predictions regarding the impact of gene regulatory network rewiring on evolutionary diversification and disease ontogeny.

Unbalanced response to growth variations reshapes the cell fate decision landscape

The global regulation of cell growth rate on gene expression perturbs the performance of gene networks, which would impose complex variations on the cell-fate decision landscape. Here we use a simple synthetic circuit of mutual repression that allows a bistable landscape to examine how such global regulation would affect the stability of phenotypic landscape and the accompanying dynamics of cell-fate determination. We show that the landscape experiences a growth-rate-induced bifurcation between monostability and bistability. Theoretical and experimental analyses reveal that this bifurcating deformation of landscape arises from the unbalanced response of gene expression to growth variations. The path of growth transition across the bifurcation would reshape cell-fate decisions. These results demonstrate the importance of growth regulation on cell-fate determination processes, regardless of specific molecular signaling or regulation.

Natural antisense transcription of presenilin in sea urchin reveals a possible role for natural antisense transcription in the general control of gene expression during development

One presenilin gene (PSEN) is expressed in the sea urchin embryo, in the vegetal pole of the gastrula and then mainly in cilia cells located around the digestive system of the pluteus, as we recently have reported. PSEN expression must be accurately regulated for correct execution of these two steps of development. While investigating PSEN expression changes in embryos after expansion of endoderm with LiCl or of ectoderm with Zn2+ by whole-mount in situ hybridization (WISH) and quantitative PCR (qPCR), we detected natural antisense transcription of PSEN. We then found that Endo16 and Wnt5, markers of endo-mesoderm, and of Hnf6 and Gsc, markers of ectoderm, are also sense and antisense transcribed. We discuss that general gene expression could depend on both sense and antisense transcription. This mechanism, together with the PSEN gene, should be included in gene regulatory networks (GRNs) that theorize diverse processes in this species. We suggest that it would also be relevant to investigate natural antisense transcription of PSEN in the field of Alzheimer's disease (AD) where the role of human PSEN1 and PSEN2 is well known.

Bioelectricity in Developmental Patterning and Size Control: Evidence and Genetically Encoded Tools in the Zebrafish Model

Developmental patterning is essential for regulating cellular events such as axial patterning, segmentation, tissue formation, and organ size determination during embryogenesis. Understanding the patterning mechanisms remains a central challenge and fundamental interest in developmental biology. Ion-channel-regulated bioelectric signals have emerged as a player of the patterning mechanism, which may interact with morphogens. Evidence from multiple model organisms reveals the roles of bioelectricity in embryonic development, regeneration, and cancers. The Zebrafish model is the second most used vertebrate model, next to the mouse model. The zebrafish model has great potential for elucidating the functions of bioelectricity due to many advantages such as external development, transparent early embryogenesis, and tractable genetics. Here, we review genetic evidence from zebrafish mutants with fin-size and pigment changes related to ion channels and bioelectricity. In addition, we review the cell membrane voltage reporting and chemogenetic tools that have already been used or have great potential to be implemented in zebrafish models. Finally, new perspectives and opportunities for bioelectricity research with zebrafish are discussed.

Decoding the IGF1 signaling gene regulatory network behind alveologenesis from a mouse model of bronchopulmonary dysplasia

Lung development is precisely controlled by underlying gene regulatory networks (GRN). Disruption of genes in the network can interrupt normal development and cause diseases such as bronchopulmonary dysplasia (BPD) - a chronic lung disease in preterm infants with morbid and sometimes lethal consequences characterized by lung immaturity and reduced alveolarization. Here, we generated a transgenic mouse exhibiting a moderate severity BPD phenotype by blocking IGF1 signaling in secondary crest myofibroblasts (SCMF) at the onset of alveologenesis. Using approaches mirroring the construction of the model GRN in sea urchin's development, we constructed the IGF1 signaling network underlying alveologenesis using this mouse model that phenocopies BPD. The constructed GRN, consisting of 43 genes, provides a bird's eye view of how the genes downstream of IGF1 are regulatorily connected. The GRN also reveals a mechanistic interpretation of how the effects of IGF1 signaling are transduced within SCMF from its specification genes to its effector genes and then from SCMF to its neighboring alveolar epithelial cells with WNT5A and FGF10 signaling as the bridge. Consistently, blocking WNT5A signaling in mice phenocopies BPD as inferred by the network. A comparative study on human samples suggests that a GRN of similar components and wiring underlies human BPD. Our network view of alveologenesis is transforming our perspective to understand and treat BPD. This new perspective calls for the construction of the full signaling GRN underlying alveologenesis, upon which targeted therapies for this neonatal chronic lung disease can be viably developed.

Dynamic Catalysis Guided by Nucleic Acid Networks and DNA Nanostructures

Nucleic acid networks conjugated to native enzymes and supramolecular DNA nanostructures modified with enzymes or DNAzymes act as functional reaction modules for guiding dynamic catalytic transformations. These systems are exemplified with the assembly of constitutional dynamic networks (CDNs) composed of nucleic acid-functionalized enzymes, as constituents, undergoing triggered structural reconfiguration, leading to dynamically switched biocatalytic cascades. By coupling two nucleic acid/enzyme networks, the intercommunicated feedback-driven dynamic biocatalytic operation of the system is demonstrated. In addition, the tailoring of a nucleic acid/enzyme reaction network driving a dissipative, transient, biocatalytic cascade is introduced as a model system for out-of-equilibrium dynamically modulated biocatalytic transformation in nature. Also, supramolecular nucleic acid machines or DNA nanostructures, modified with DNAzyme or enzyme constituents, act as functional reaction modules driving temporal dynamic catalysis. The design of dynamic supramolecular machines is exemplified with the introduction of an interlocked two-ring catenane device that is dynamically reversibly switched between two states operating two different DNAzymes, and with the tailoring of a DNA-tweezers device functionalized with enzyme/DNAzyme constituents that guides the dynamic ON/OFF operation of a biocatalytic cascade by opening and closing the molecular device. In addition, DNA origami nanostructures provide functional scaffolds for the programmed positioning of enzymes or DNAzyme for the switchable operation of catalytic transformations. This is introduced by the tailored functionalization of the edges of origami tiles with nucleic acids guiding the switchable formation of DNAzyme catalysts through the dimerization/separation of the tiles. In addition, the programmed deposition of two-enzyme/cofactor constituents on the origami raft allowed the dynamic photochemical activation of the cofactor-mediated biocatalytic cascade on the spatially biocatalytic assembly on the scaffold. Furthermore, photoinduced "mechanical" switchable and reversible unlocking and closing of nanoholes in the origami frameworks allow the "ON" and "OFF" operation of DNAzyme units in the nanoholes, confined environments. The future challenges and potential applications of dynamic nucleic acid/enzyme and DNAzyme conjugates are discussed in the conclusion paragraph.

Standardized excitable elements for scalable engineering of far-from-equilibrium chemical networks

Engineered far-from-equilibrium synthetic chemical networks that pulse or switch states in response to environmental signals could precisely regulate the kinetics of chemical synthesis or self-assembly. Currently, such networks must be extensively tuned to compensate for the different activities of and unintended reactions between a network's various chemical components. Modular elements with standardized performance could be used to rapidly construct networks with designed functions. Here we develop standardized excitable chemical regulatory elements, termed genelets, and use them to construct complex in vitro transcriptional networks. We develop a protocol for identifying >15 interchangeable genelet elements with uniform performance and minimal crosstalk. These elements can be combined to engineer feedforward and feedback modules whose dynamics match those predicted by a simple kinetic model. Modules can then be rationally integrated and organized into networks that produce tunable temporal pulses and act as multistate switchable memories. Standardized genelet elements, and the workflow to identify more, should make engineering complex far-from-equilibrium chemical dynamics routine.

Cascaded dissipative DNAzyme-driven layered networks guide transient replication of coded-strands as gene models

Dynamic, transient, out-of-equilibrium networks guide cellular genetic, metabolic or signaling processes. Designing synthetic networks emulating natural processes imposes important challenges including the ordered connectivity of transient reaction modules, engineering of the appropriate balance between production and depletion of reaction constituents, and coupling of the reaction modules with emerging chemical functions dictated by the networks. Here we introduce the assembly of three coupled reaction modules executing a cascaded dynamic process leading to the transient formation and depletion of three different Mg²⁺-ion-dependent DNAzymes. The transient operation of the DNAzyme in one layer triggers the dynamic activation of the DNAzyme in the subsequent layer, leading to a three-layer transient catalytic cascade. The kinetics of the transient cascade is computationally simulated. The cascaded network is coupled to a polymerization/nicking DNA machinery guiding transient synthesis of three coded strands acting as “gene models”, and to the rolling circle polymerization machinery leading to the transient synthesis of fluorescent Zn(II)-PPIX/G-quadruplex chains or hemin/G-quadruplex catalytic wires.

A reaction network executing a cascaded transient formation and depletion of three different catalytic strands is introduced. The system is coupled to the secondary temporal synthesis of different coded strands as gene models.

Expression Profile of Genes Related to the Th17 Pathway in Macrophages Infected by Leishmania major and Leishmania amazonensis: The Use of Gene Regulatory Networks in Modeling This Pathway

Leishmania amazonensis and Leishmania major are the causative agents of cutaneous and mucocutaneous diseases. The infections‘ outcome depends on host–parasite interactions and Th1/Th2 response, and in cutaneous form, regulation of Th17 cytokines has been reported to maintain inflammation in lesions. Despite that, the Th17 regulatory scenario remains unclear. With the aim to gain a better understanding of the transcription factors (TFs) and genes involved in Th17 induction, in this study, the role of inducing factors of the Th17 pathway in Leishmania–macrophage infection was addressed through computational modeling of gene regulatory networks (GRNs). The Th17 GRN modeling integrated experimentally validated data available in the literature and gene expression data from a time-series RNA-seq experiment (4, 24, 48, and 72 h post-infection). The generated model comprises a total of 10 TFs, 22 coding genes, and 16 cytokines related to the Th17 immune modulation. Addressing the Th17 induction in infected and uninfected macrophages, an increase of 2- to 3-fold in 4–24 h was observed in the former. However, there was a decrease in basal levels at 48–72 h for both groups. In order to evaluate the possible outcomes triggered by GRN component modulation in the Th17 pathway. The generated GRN models promoted an integrative and dynamic view of Leishmania–macrophage interaction over time that extends beyond the analysis of single-gene expression.

Differential Etv2 threshold requirement for endothelial and erythropoietic development

Endothelial and erythropoietic lineages arise from a common developmental progenitor. Etv2 is a master transcriptional regulator required for the development of both lineages. However, the mechanisms through which Etv2 initiates the gene-regulatory networks (GRNs) for endothelial and erythropoietic specification and how the two GRNs diverge downstream of Etv2 remain incompletely understood. Here, by analyzing a hypomorphic Etv2 mutant, we demonstrate different threshold requirements for initiation of the downstream GRNs for endothelial and erythropoietic development. We show that Etv2 functions directly in a coherent feedforward transcriptional network for vascular endothelial development, and a low level of Etv2 expression is sufficient to induce and sustain the endothelial GRN. In contrast, Etv2 induces the erythropoietic GRN indirectly via activation of Tal1, which requires a significantly higher threshold of Etv2 to initiate and sustain erythropoietic development. These results provide important mechanistic insight into the divergence of the endothelial and erythropoietic lineages.

Linking Molecular Mechanisms and Evolutionary Consequences of Resource Polyphenism

Resource polyphenism—the occurrence of environmentally induced, discrete, and intraspecific morphs showing differential niche use—is taxonomically widespread and fundamental to the evolution of ecological function where it has arisen. Despite longstanding appreciation for the ecological and evolutionary significance of resource polyphenism, only recently have its proximate mechanisms begun to be uncovered. Polyphenism switches, especially those influencing and influenced by trophic interactions, offer a route to integrating proximate and ultimate causation in studies of plasticity, and its potential influence on evolution more generally. Here, we use the major events in generalized polyphenic development as a scaffold for linking the molecular mechanisms of polyphenic switching with potential evolutionary outcomes of polyphenism and for discussing challenges and opportunities at each step in this process. Not only does the study of resource polyphenism uncover interesting details of discrete plasticity, it also illuminates and informs general principles at the intersection of development, ecology, and evolution.

Tenuous Transcriptional Threshold of Human Sex Determination. I. SRY and Swyer Syndrome at the Edge of Ambiguity

Male sex determination in mammals is initiated by SRY, a Y-encoded transcription factor. The protein contains a high-mobility-group (HMG) box mediating sequence-specific DNA bending. Mutations causing XY gonadal dysgenesis (Swyer syndrome) cluster in the box and ordinarily arise de novo. Rare inherited variants lead to male development in one genetic background (the father) but not another (his sterile XY daughter). De novo and inherited mutations occur at an invariant Tyr adjoining the motif's basic tail (box position 72; Y127 in SRY). In SRY-responsive cell lines CH34 and LNCaP, de novo mutations Y127H and Y127C reduced SRY activity (as assessed by transcriptional activation of principal target gene Sox9) by 5- and 8-fold, respectively. Whereas Y127H impaired testis-specific enhancer assembly, Y127C caused accelerated proteasomal proteolysis; activity was in part rescued by proteasome inhibition. Inherited variant Y127F was better tolerated: its expression was unperturbed, and activity was reduced by only twofold, a threshold similar to other inherited variants. Biochemical studies of wild-type (WT) and variant HMG boxes demonstrated similar specific DNA affinities (within a twofold range), with only subtle differences in sharp DNA bending as probed by permutation gel electrophoresis and fluorescence resonance-energy transfer (FRET); thermodynamic stabilities of the free boxes were essentially identical. Such modest perturbations are within the range of species variation. Whereas our cell-based findings rationalize the de novo genotype-phenotype relationships, a molecular understanding of inherited mutation Y127F remains elusive. Our companion study uncovers cryptic biophysical perturbations suggesting that the para-OH group of Y127 anchors a novel water-mediated DNA clamp.

Three topological features of regulatory networks control life-essential and specialized subsystems

Gene regulatory networks (GRNs) play key roles in development, phenotype plasticity, and evolution. Although graph theory has been used to explore GRNs, associations amongst topological features, transcription factors (TFs), and systems essentiality are poorly understood. Here we sought the relationship amongst the main GRN topological features that influence the control of essential and specific subsystems. We found that the K_nn, page rank, and degree are the most relevant GRN features: the ones are conserved along the evolution and are also relevant in pluripotent cells. Interestingly, life-essential subsystems are governed mainly by TFs with intermediary K_nn and high page rank or degree, whereas specialized subsystems are mainly regulated by TFs with low K_nn. Hence, we suggest that the high probability of TFs be toured by a random signal, and the high probability of the signal propagation to target genes ensures the life-essential subsystems' robustness. Gene/genome duplication is the main evolutionary process to rise K_nn as the most relevant feature. Herein, we shed light on unexplored topological GRN features to assess how they are related to subsystems and how the duplications shaped the regulatory systems along the evolution. The classification model generated can be found here: https://github.com/ivanrwolf/NoC/ .

Bioinspired Artificial Photosynthetic Systems

Mimicking photosynthesis using artificial systems, as a means for solar energy conversion and green fuel generation, is one of the holy grails of modern science. This perspective presents recent advances towards developing artificial photosynthetic systems. In one approach, native photosystems are interfaced with electrodes to yield photobioelectrochemical cells that transform light energy into electrical power. This is exemplified by interfacing photosystem I (PSI) and photosystem II (PSII) as an electrically contacted assembly mimicking the native Z-scheme, and by the assembly of an electrically wired PSI/glucose oxidase biocatalytic conjugate on an electrode support. Illumination of the functionalized electrodes led to light-induced generation of electrical power, or to the generation of photocurrents using glucose as the fuel. The second approach introduces supramolecular photosensitizer nucleic acid/electron acceptor complexes as functional modules for effective photoinduced electron transfer stimulating the subsequent biocatalyzed generation of NADPH or the Pt-nanoparticle-catalyzed evolution of molecular hydrogen. Application of the DNA machineries for scaling-up the photosystems is demonstrated. A third approach presents the integration of artificial photosynthetic modules into dynamic nucleic acid networks undergoing reversible reconfiguration or dissipative transient operation in the presence of auxiliary triggers. Control over photoinduced electron transfer reactions and photosynthetic transformations by means of the dynamic networks is demonstrated.

Towards a physical understanding of developmental patterning

The temporal coordination of events at cellular and tissue scales is essential for the proper development of organisms, and involves cell-intrinsic processes that can be coupled by local cellular signalling and instructed by global signalling, thereby creating spatial patterns of cellular states that change over time. The timing and structure of these patterns determine how an organism develops. Traditional developmental genetic methods have revealed the complex molecular circuits regulating these processes but are limited in their ability to predict and understand the emergent spatio-temporal dynamics. Increasingly, approaches from physics are now being used to help capture the dynamics of the system by providing simplified, generic descriptions. Combined with advances in imaging and computational power, such approaches aim to provide insight into timing and patterning in developing systems.

A Unifying Framework for Understanding Biological Structures and Functions Across Levels of Biological Organization

The relationship between structure and function is a major constituent of the rules of life. Structures and functions occur across all levels of biological organization. Current efforts to integrate conceptual frameworks and approaches to address new and old questions promise to allow a more holistic and robust understanding of how different biological functions are achieved across levels of biological organization. Here, we provide unifying and generalizable definitions of both structure and function that can be applied across all levels of biological organization. However, we find differences in the nature of structures at the organismal level and below as compared to above the level of the organism. We term these intrinsic and emergent structures, respectively. Intrinsic structures are directly under selection, contributing to the overall performance (fitness) of the individual organism. Emergent structures involve interactions among aggregations of organisms and are not directly under selection. Given this distinction, we argue that while the functions of many intrinsic structures remain unknown, functions of emergent structures are the result of the aggregate of processes of individual organisms. We then provide a detailed and unified framework of the structure-function relationship for intrinsic structures to explore how their unknown functions can be defined. We provide examples of how these scalable definitions applied to intrinsic structures provide a framework to address questions on structure-function relationships that can be approached simultaneously from all subdisciplines of biology. We propose that this will produce a more holistic and robust understanding of how different biological functions are achieved across levels of biological organization.

Integration of photocatalytic and dark-operating catalytic biomimetic transformations through DNA-based constitutional dynamic networks

Nucleic acid-based constitutional dynamic networks (CDNs) have recently emerged as versatile tools to control a variety of catalytic processes. A key challenge in the application of these systems is achieving intercommunication between different CDNs to mimic the complex interlinked networks found in cellular biology. In particular, the possibility to interface photochemical ‘energy-harvesting’ processes with dark-operating ‘metabolic’ processes, in a similar way to plants, represents an up to now unexplored yet enticing research direction. The present study introduces two CDNs that allow the intercommunication of photocatalytic and dark-operating catalytic functions mediated by environmental components that facilitate the dynamic coupling of the networks. The dynamic feedback-driven intercommunication of the networks is accomplished via information transfer between the two CDNs effected by hairpin fuel strands in the environment of the system, leading to the coupling of the photochemical and dark-operating modules.

Nucleic acid-based constitutional dynamic networks (CDNs) enable control of various catalytic processes, but it is challenging to achieve intercommunication between different CDNs and by that mimic complex cell biology networks. Here, the authors report two CDNs that control the integration of photochemical and dark-operating processes, and show their intercommunication afforded by environmental components.

Dynamical modules in metabolism, cell and developmental biology

Modularity is an essential feature of any adaptive complex system. Phenotypic traits are modules in the sense that they have a distinguishable structure or function, which can vary (quasi-)independently from its context. Since all phenotypic traits are the product of some underlying regulatory dynamics, the generative processes that constitute the genotype-phenotype map must also be functionally modular. Traditionally, modular processes have been identified as structural modules in regulatory networks. However, structure only constrains, but does not determine, the dynamics of a process. Here, we propose an alternative approach that decomposes the behaviour of a complex regulatory system into elementary activity-functions. Modular activities can occur in networks that show no structural modularity, making dynamical modularity more widely applicable than structural decomposition. Furthermore, the behaviour of a regulatory system closely mirrors its functional contribution to the outcome of a process, which makes dynamical modularity particularly suited for functional decomposition. We illustrate our approach with numerous examples from the study of metabolism, cellular processes, as well as development and pattern formation. We argue that dynamical modules provide a shared conceptual foundation for developmental and evolutionary biology, and serve as the foundation for a new account of process homology, which is presented in a separate contribution by DiFrisco and Jaeger to this focus issue.

Cis- and Trans-Regulatory Variations in the Domestication of the Chili Pepper Fruit

The process of domestication requires the rapid transformation of the wild morphology into the cultivated forms that humans select for. This process often takes place through changes in the regulation of genes, yet, there is no definite pattern on the role of cis- and trans-acting regulatory variations in the domestication of the fruit among crops. Using allele-specific expression and network analyses, we characterized the regulatory patterns and the inheritance of gene expression in wild and cultivated accessions of chili pepper, a crop with remarkable fruit morphological variation. We propose that gene expression differences associated to the cultivated form are best explained by cis-regulatory hubs acting through trans-regulatory cascades. We show that in cultivated chili, the expression of genes associated with fruit morphology is partially recessive with respect to those in the wild relative, consistent with the hybrid fruit phenotype. Decreased expression of fruit maturation and growth genes in cultivated chili suggest that selection for loss-of-function took place in its domestication. Trans-regulatory changes underlie the majority of the genes showing regulatory divergence and had larger effect sizes on gene expression than cis-regulatory variants. Network analysis of selected cis-regulated genes, including ARP9 and MED25, indicated their interaction with many transcription factors involved in organ growth and fruit ripening. Differentially expressed genes linked to cis-regulatory variants and their interactions with downstream trans-acting genes have the potential to drive the morphological differences observed between wild and cultivated fruits and provide an attractive mechanism of morphological transformation during the domestication of the chili pepper.

Brachyury controls Ciona notochord fate as part of a feed-forward network

The notochord is a defining feature of the chordates. The transcription factor Brachyury (Bra) is a key regulator of notochord fate but here we show that it is not a unitary master regulator in the model chordate Ciona Ectopic Bra expression only partially reprograms other cell types to a notochord-like transcriptional profile and a subset of notochord-enriched genes is unaffected by CRISPR Bra disruption. We identify Foxa.a and Mnx as potential co-regulators, and find that combinatorial cocktails are more effective at reprogramming other cell types than Bra alone. We reassess the network relationships between Bra, Foxa.a and other components of the notochord gene regulatory network, and find that Foxa.a expression in the notochord is regulated by vegetal FGF signaling. It is a direct activator of Bra expression and has a binding motif that is significantly enriched in the regulatory regions of notochord-enriched genes. These and other results indicate that Bra and Foxa.a act together in a regulatory network dominated by positive feed-forward interactions, with neither being a classically defined master regulator.

Gene Expression Changes after Parental Exposure to Metals in the Sea Urchin Affect Timing of Genetic Programme of Embryo Development

It is widely accepted that phenotypic traits can be modulated at the epigenetic level so that some conditions can affect the progeny of exposed individuals. To assess if the exposure of adult animals could result in effects on the offspring, the Mediterranean sea urchin and its well-characterized gene regulatory networks (GRNs) was chosen as a model. Adult animals were exposed to known concentrations of zinc and cadmium (both individually and in combination) for 10 days, and the resulting embryos were followed during the development. The oxidative stress occurring in parental gonads, embryo phenotypes and mortality, and the expression level of a set of selected genes, including members of the skeletogenic and endodermal GRNs, were evaluated. Increased oxidative stress at F₀, high rates of developmental aberration with impaired gastrulation, in association to deregulation of genes involved in skeletogenesis (dri, hex, sm50, p16, p19, msp130), endodermal specification (foxa, hox11/13b, wnt8) and epigenetic regulation (kat2A, hdac1, ehmt2, phf8 and UBE2a) occurred either at 24 or 48 hpf. Results strongly indicate that exposure to environmental pollutants can affect not only directly challenged animals but also their progeny (at least F₁), influencing optimal timing of genetic programme of embryo development, resulting in an overall impairment of developmental success.

Sea Urchin as a Universal Model for Studies of Gene Networks

The purple sea urchin Strongylocentrotus purpuratus has been used for over 150 years as a model organism in developmental biology. Using this model species, scientists have been able to describe, in detail, the mechanisms of cell cycle control and cell adhesion, fertilization, calcium signaling, cell differentiation, and death. Massive parallel sequencing of the sea urchin genome enabled the deciphering of the main components of gene regulatory networks during the activation of embryonic signaling pathways. This knowledge helped to extrapolate aberrations in somatic cells that may lead to diseases, including cancer in humans. Furthermore, since many, if not all, developmental signaling pathways were shown to be controlled by non-coding RNAs (ncRNAs), the sea urchin organism represents an attractive experimental model. In this review, we discuss the main discoveries in the genetics, genomics, and transcriptomics of sea urchins during embryogenesis with the main focus on the role of ncRNAs. This information may be useful for comparative studies between different organisms, and may help identify new regulatory networks controlled by ncRNAs.

A regulatory sub-circuit downstream of Wnt signaling controls developmental transitions in neural crest formation

The process of cell fate commitment involves sequential changes in the gene expression profiles of embryonic progenitors. This is exemplified in the development of the neural crest, a migratory stem cell population derived from the ectoderm of vertebrate embryos. During neural crest formation, cells transition through distinct transcriptional states in a stepwise manner. The mechanisms underpinning these shifts in cell identity are still poorly understood. Here we employ enhancer analysis to identify a genetic sub-circuit that controls developmental transitions in the nascent neural crest. This sub-circuit links Wnt target genes in an incoherent feedforward loop that controls the sequential activation of genes in the neural crest lineage. By examining the cis-regulatory apparatus of Wnt effector gene AXUD1, we found that multipotency factor SP5 directly promotes neural plate border identity, while inhibiting premature expression of specification genes. Our results highlight the importance of repressive interactions in the neural crest gene regulatory network and illustrate how genes activated by the same upstream signal become temporally segregated during progressive fate restriction.

Computational Stem Cell Biology: Open Questions and Guiding Principles

Computational biology is enabling an explosive growth in our understanding of stem cells and our ability to use them for disease modeling, regenerative medicine, and drug discovery. We discuss four topics that exemplify applications of computation to stem cell biology: cell typing, lineage tracing, trajectory inference, and regulatory networks. We use these examples to articulate principles that have guided computational biology broadly and call for renewed attention to these principles as computation becomes increasingly important in stem cell biology. We also discuss important challenges for this field with the hope that it will inspire more to join this exciting area.

Constitutional Dynamic Networks-Guided Synthesis of Programmed "Genes", Transcription of mRNAs, and Translation of Proteins

Inspired by nature, where dynamic networks control the levels of gene expression and the activities of transcribed/translated proteins, we introduce nucleic acid-based constitutional dynamic networks (CDNs) as functional modules mimicking native circuits by demonstrating CDNs-guided programmed synthesis of genes, controlled transcription of RNAs, and dictated transcription/translation synthesis of proteins. An auxiliary CDN consisting of four dynamically equilibrated constituents AA', AB', BA', and BB' is orthogonally triggered by two different inputs yielding two different compositionally reconfigured CDNs. Subjecting the parent auxiliary CDN to two hairpins, H_A and H_B, and two templates T_A and T_B and a nicking/replication machinery leads to the cleavage of the hairpins and to the activation of the nicking/replication machineries that synthesize two "genes", e.g., the histidine-dependent DNAzyme g1 and the Zn²⁺-ion-dependent DNAzyme g2. The triggered orthogonal reconfiguration of the parent CDN to the respective CDNs leads to the programmed preferred CDN-guided synthesis of g1 or g2. Similarly, the triggered reconfigured CDNs are subjected to two hairpins H_C and H_D, the templates I'/I and J'/J, and the RNA polymerase (RNAp)/NTPs machinery. While the cleavage of the hairpins by the constituents associated with the parent CDN leads to the transcription of the broccoli aptamer recognizing the DFHBI ligand and of the aptamer recognizing the malachite green (MG) ligand, the orthogonally triggered CDNs lead to the CDNs-guided enhanced transcription of either the DFHBI aptamer or the MG aptamer. In addition, subjecting the triggered reconfigured CDNs to predesigned hairpins H_E and H_F, the templates M'/M and N'/N, the RNAp/NTPs machinery, and the cell-free ribosome t-RNA machinery leads to the CDNs-guided transcription/translation of the green fluorescence protein (GFP) or red fluorescence protein (RFP).

Seven myths of how transcription factors read the cis-regulatory code

Genomics data are now being generated at large quantities, of exquisite high resolution and from single cells. They offer a unique opportunity to develop powerful machine learning algorithms, including neural networks, to uncover the rules of the cis-regulatory code. However, current modeling assumptions are often not based on state-of-the-art knowledge of the cis-regulatory code from transcription, developmental genetics, imaging and structural studies. Here I aim to fill this gap by giving a brief historical overview of the field, describing common misconceptions and providing knowledge that might help to guide computational approaches. I will describe the principles and mechanisms involved in the combinatorial requirement of transcription factor binding motifs for enhancer activity, including the role of chromatin accessibility, repressors and low-affinity motifs in the cis-regulatory code. Deciphering the cis-regulatory code would unlock an enormous amount of regulatory information in the genome and would allow us to locate cis-regulatory genetic variants involved in development and disease.

A conceptual framework of evolutionary novelty and innovation

Since 1990 the recognition of deep homologies among metazoan developmental processes and the spread of more mechanistic approaches to developmental biology have led to a resurgence of interest in evolutionary novelty and innovation. Other evolutionary biologists have proposed central roles for behaviour and phenotypic plasticity in generating the conditions for the construction of novel morphologies, or invoked the accessibility of new regions of vast sequence spaces. These approaches contrast with more traditional emphasis on the exploitation of ecological opportunities as the primary source of novelty. This definitional cornucopia reflects differing stress placed on three attributes of novelties: their radical nature, the generation of new taxa, and ecological and evolutionary impact. Such different emphasis has led to conflating four distinct issues: the origin of novel attributes (genes, developmental processes, phenotypic characters), new functions, higher clades and the ecological impact of new structures and functions. Here I distinguish novelty (the origin of new characters, deep character transformations, or new combinations) from innovation, the ecological and evolutionary success of clades. Evidence from the fossil record of macroevolutionary lags between the origin of a novelty and its ecological success demonstrates that novelty may be decoupled from innovation, and only definitions of novelty based on radicality (rather than generativity or consequentiality) can be assessed without reference to the subsequent history of the clade to which a novelty belongs. These considerations suggest a conceptual framework for novelty and innovation, involving: (i) generation of the potential for novelty; (ii) the formation of novel attributes; (iii) refinement of novelties through adaptation; (iv) exploitation of novelties by a clade, which may coincide with a new round of ecological or environmental potentiation; followed by (v) the establishment of innovations through ecological processes. This framework recognizes that there is little empirical support for either the dominance of ecological opportunity, nor abrupt discontinuities (often caricatured as 'hopeful monsters'). This general framework may be extended to aspects of cultural and social innovation.

The function of architecture and logic in developmental gene regulatory networks

An important contribution of systems biology is the insight that biological systems depend on the function of molecular interactions and not just on individual molecules. System level mechanisms are particularly important in the development of animals and plants which depends not just on transcription factors and signaling molecules, but also on regulatory circuits and gene regulatory networks (GRNs). However, since GRNs consist of transcription factors, it can be challenging to assess the function of regulatory circuits independently of the function of regulatory factors. The comparison of different GRNs offers a way to do so and leads to several observations. First, similar regulatory circuits operate in various developmental contexts and in different species, and frequently, these circuits are associated with similar developmental functions. Second, given regulatory circuits are often used at particular positions within the GRN hierarchy. Third, in some GRNs, regulatory circuits are organized in a particular order in respect to each other. And fourth, the evolution of GRNs occurs not just by co-option of regulatory genes but also by rewiring of regulatory linkages between conserved regulatory genes, indicating that the organization of interactions is important. Thus, even though in most instances the function of regulatory circuits remains to be discovered, it becomes evident that the architecture and logic of GRNs are functionally important for the control of genome activity and for the specification of the body plan.

Early Xenopus gene regulatory programs, chromatin states, and the role of maternal transcription factors

For decades, the early development of the Xenopus embryo has been an essential model system to study the gene regulatory mechanisms that govern cellular specification. At the top of the hierarchy of gene regulatory networks, maternally deposited transcription factors initiate this process and regulate the expression of zygotic genes that give rise to three distinctive germ layer cell types (ectoderm, mesoderm, and endoderm), and subsequent generation of organ precursors. The onset of germ layer specification is also closely coupled with changes associated with chromatin modifications. This review will examine the timing of maternal transcription factors initiating the zygotic genome activation, the epigenetic landscape of embryonic chromatin, and the network structure that governs the process.

Genotype Components as Predictors of Phenotype in Model Gene Regulatory Networks

Models of gene regulatory networks (GRN) have proven useful for understanding many aspects of the highly complex behavior of biological control networks. Randomly generated non-Boolean networks were used in experimental simulations to generate data on dynamic phenotypes as a function of several genotypic parameters. We found that predictive relationships between some phenotypes and quantitative genotypic parameters such as number of network genes, interaction density, and initial condition could be derived depending on the strength of the topological (positional) genotype on specific phenotypes. We quantitated the strength of the topological genotype effect (TGE) on a number of phenotypes in multi-gene networks. For phenotypes with a low influence of topological genotype, derived and empirical relationships using quantitative genotype parameters were accurate in phenotypic outcomes. We found a number of dynamic network properties, including oscillation behaviors, that were largely dependent on genotype topology, and for which no such general quantitative relationships were determinable. It remains to be determined if these results are applicable to biological gene regulatory networks.

Regulation of Cell Fate Decisions in Early Mammalian Embryos

Early embryogenesis is characterized by the segregation of cell lineages that fulfill critical roles in the establishment of pregnancy and development of the fetus. The formation of the blastocyst marks the emergence of extraembryonic precursors, needed for implantation, and of pluripotent cells, which differentiate toward the major lineages of the adult organism. The coordinated emergence of these cell types shows that these processes are broadly conserved in mammals. However, developmental heterochrony and changes in gene regulatory networks highlight unique evolutionary adaptations that may explain the diversity in placentation and in the mechanisms controlling pluripotency in mammals. The incorporation of new technologies, including single-cell omics, imaging, and gene editing, is instrumental for comparative embryology. Broadening the knowledge of mammalian embryology will provide new insights into the mechanisms driving evolution and development. This knowledge can be readily translated into biomedical and biotechnological applications in humans and livestock, respectively.

Network architecture and regulatory logic in neural crest development

The neural crest is an ectodermal cell population that gives rise to over 30 cell types during vertebrate embryogenesis. These stem cells are formed at the border of the developing central nervous system and undergo extensive migration before differentiating into components of multiple tissues and organs. Neural crest formation and differentiation is a multistep process, as these cells transition through sequential regulatory states before adopting their adult phenotype. Such changes are governed by a complex gene regulatory network (GRN) that integrates environmental and cell-intrinsic inputs to regulate cell identity. Studies of neural crest cells in a variety of vertebrate models have elucidated the function and regulation of dozens of the molecular players that are part of this network. The neural crest GRN has served as a platform to explore the molecular control of multipotency, cell differentiation, and the evolution of vertebrates. In this review, we employ this genetic program as a stepping-stone to explore the architecture and the regulatory principles of developmental GRNs. We also discuss how modern genomic approaches can further expand our understanding of genetic networks in this system and others. This article is categorized under: Physiology > Mammalian Physiology in Health and Disease Biological Mechanisms > Cell Fates Developmental Biology > Lineages Models of Systems Properties and Processes > Cellular Models.

Gastrulation in the sea urchin

Gastrulation is arguably the most important evolutionary innovation in the animal kingdom. This process provides the basic embryonic architecture, an inner layer separated from an outer layer, from which all animal forms arise. An extraordinarily simple and elegant process of gastrulation is observed in the sea urchin embryo. The cells participating in sea urchin gastrulation are specified early during cleavage. One outcome of that specification is the expression of transcription factors that control each of the many subsequent morphogenetic changes. The first of these movements is an epithelial-mesenchymal transition (EMT) of skeletogenic mesenchyme cells, then EMT of pigment cell progenitors. Shortly thereafter, invagination of the archenteron occurs. At the end of archenteron extension, a second wave of EMT occurs to release immune cells into the blastocoel and primordial germ cells that will home to the coelomic pouches. The archenteron then remodels to establish the three parts of the gut, and at the anterior end, the gut fuses with the stomodaeum to form the through-gut. As part of the anterior remodeling, mesodermal coelomic pouches bud off the lateral sides of the archenteron tip. Multiple cell biological processes conduct each of these movements and in some cases the upstream transcription factors controlling this process have been identified. Remarkably, each event seamlessly occurs at the right time to orchestrate formation of the primitive body plan. This review covers progress toward understanding many of the molecular mechanisms underlying this sequence of morphogenetic events.

Tempos and modes of collectivity in the history of life

Collective integration and processing of information have increased through the history of life, through both the formation of aggregates in which the entities may have very different properties and which jointly coarse-grained environmental variables (ranging from widely varying metabolism in microbial consortia to the ecological diversity of species on reefs) and through collectives of similar entities (such as cells within an organism or social groups). Such increases have been implicated in significant transitions in the history of life, including aspects of the origin of life, the generation of pangenomes among microbes and microbial communities such as stromatolites, multicellularity and social insects. This contribution provides a preliminary overview of the dominant modes of collective information processing in the history of life, their phylogenetic distribution and extent of convergence, and the effects of new modes for integrating and acting upon information on the tempo of evolutionary change.

Building in vitro transcriptional regulatory networks by successively integrating multiple functional circuit modules

The regulation of cellular dynamics and responses to stimuli by genetic regulatory networks suggests how in vitro chemical reaction networks might analogously direct the dynamics of synthetic materials or chemistries. A key step in developing genetic regulatory network analogues capable of this type of sophisticated regulation is the integration of multiple coordinated functions within a single network. Here, we demonstrate how such functional integration can be achieved using in vitro transcriptional genelet circuits that emulate essential features of cellular genetic regulatory networks. By successively incorporating functional genelet modules into a bistable circuit, we construct an integrated regulatory network that dynamically changes its state in response to upstream stimuli and coordinates the timing of downstream signal expression. We use quantitative models to guide module integration and develop strategies to mitigate undesired interactions between network components that arise as the size of the network increases. This approach could enable the construction of in vitro networks capable of multifaceted chemical and material regulation.

Extensive intraspecies cryptic variation in an ancient embryonic gene regulatory network

Innovations in metazoan development arise from evolutionary modification of gene regulatory networks (GRNs). We report widespread cryptic variation in the requirement for two key regulatory inputs, SKN-1/Nrf2 and MOM-2/Wnt, into the C. elegans endoderm GRN. While some natural isolates show a nearly absolute requirement for these two regulators, in others, most embryos differentiate endoderm in their absence. GWAS and analysis of recombinant inbred lines reveal multiple genetic regions underlying this broad phenotypic variation. We observe a reciprocal trend, in which genomic variants, or knockdown of endoderm regulatory genes, that result in a high SKN-1 requirement often show low MOM-2/Wnt requirement and vice-versa, suggesting that cryptic variation in the endoderm GRN may be tuned by opposing requirements for these two key regulatory inputs. These findings reveal that while the downstream components in the endoderm GRN are common across metazoan phylogeny, initiating regulatory inputs are remarkably plastic even within a single species.

How Does the Regulatory Genome Work?

The regulatory genome controls genome activity throughout the life of an organism. This requires that complex information processing functions are encoded in, and operated by, the regulatory genome. Although much remains to be learned about how the regulatory genome works, we here discuss two cases where regulatory functions have been experimentally dissected in great detail and at the systems level, and formalized by computational logic models. Both examples derive from the sea urchin embryo, but assess two distinct organizational levels of genomic information processing. The first example shows how the regulatory system of a single gene, endo16, executes logic operations through individual transcription factor binding sites and cis-regulatory modules that control the expression of this gene. The second example shows information processing at the gene regulatory network (GRN) level. The GRN controlling development of the sea urchin endomesoderm has been experimentally explored at an almost complete level. A Boolean logic model of this GRN suggests that the modular logic functions encoded at the single-gene level show compositionality and suffice to account for integrated function at the network level. We discuss these examples both from a biological-experimental point of view and from a computer science-informational point of view, as both illuminate principles of how the regulatory genome works.

Feed-forward regulation adaptively evolves via dynamics rather than topology when there is intrinsic noise

In transcriptional regulatory networks (TRNs), a canonical 3-node feed-forward loop (FFL) is hypothesized to evolve to filter out short spurious signals. We test this adaptive hypothesis against a novel null evolutionary model. Our mutational model captures the intrinsically high prevalence of weak affinity transcription factor binding sites. We also capture stochasticity and delays in gene expression that distort external signals and intrinsically generate noise. Functional FFLs evolve readily under selection for the hypothesized function but not in negative controls. Interestingly, a 4-node “diamond” motif also emerges as a short spurious signal filter. The diamond uses expression dynamics rather than path length to provide fast and slow pathways. When there is no idealized external spurious signal to filter out, but only internally generated noise, only the diamond and not the FFL evolves. While our results support the adaptive hypothesis, we also show that non-adaptive factors, including the intrinsic expression dynamics, matter.

Feed‐forward loops (FFLs) can filter out noise, but whether their overrepresentation in GRNs reflects adaptive evolution for this function is debated. Here, the authors develop a null model of regulatory evolution and find that FFLs evolve readily under selection for the noise filtering function.

Inference of Developmental Gene Regulatory Networks Beyond Classical Model Systems: New Approaches in the Post-genomic Era

The advent of high-throughput sequencing (HTS) technologies has revolutionized the way we understand the transformation of genetic information into morphological traits. Elucidating the network of interactions between genes that govern cell differentiation through development is one of the core challenges in genome research. These networks are known as developmental gene regulatory networks (dGRNs) and consist largely of the functional linkage between developmental control genes, cis-regulatory modules, and differentiation genes, which generate spatially and temporally refined patterns of gene expression. Over the last 20 years, great advances have been made in determining these gene interactions mainly in classical model systems, including human, mouse, sea urchin, fruit fly, and worm. This has brought about a radical transformation in the fields of developmental biology and evolutionary biology, allowing the generation of high-resolution gene regulatory maps to analyze cell differentiation during animal development. Such maps have enabled the identification of gene regulatory circuits and have led to the development of network inference methods that can recapitulate the differentiation of specific cell-types or developmental stages. In contrast, dGRN research in non-classical model systems has been limited to the identification of developmental control genes via the candidate gene approach and the characterization of their spatiotemporal expression patterns, as well as to the discovery of cis-regulatory modules via patterns of sequence conservation and/or predicted transcription-factor binding sites. However, thanks to the continuous advances in HTS technologies, this scenario is rapidly changing. Here, we give a historical overview on the architecture and elucidation of the dGRNs. Subsequently, we summarize the approaches available to unravel these regulatory networks, highlighting the vast range of possibilities of integrating multiple technical advances and theoretical approaches to expand our understanding on the global gene regulation during animal development in non-classical model systems. Such new knowledge will not only lead to greater insights into the evolution of molecular mechanisms underlying cell identity and animal body plans, but also into the evolution of morphological key innovations in animals.

Etv6 activates vegfa expression through positive and negative transcriptional regulatory networks in Xenopus embryos

VEGFA signaling controls physiological and pathological angiogenesis and hematopoiesis. Although many context-dependent signaling pathways downstream of VEGFA have been uncovered, vegfa transcriptional regulation in vivo remains unclear. Here, we show that the ETS transcription factor, Etv6, positively regulates vegfa expression during Xenopus blood stem cell development through multiple transcriptional inputs. In agreement with its established repressive functions, Etv6 directly inhibits expression of the repressor foxo3, to prevent Foxo3 from binding to and repressing the vegfa promoter. Etv6 also directly activates expression of the activator klf4; reflecting a genome-wide paucity in ETS-binding motifs in Etv6 genomic targets, Klf4 then recruits Etv6 to the vegfa promoter to activate its expression. These two mechanisms (double negative gate and feed-forward loop) are classic features of gene regulatory networks specifying cell fates. Thus, Etv6's dual function, as a transcriptional repressor and activator, controls a major signaling pathway involved in endothelial and blood development in vivo.

Genomic resources for the study of echinoderm development and evolution

Echinoderms are important research models for a wide range of biological questions. In particular, echinoderm embryos are exemplary models for dissecting the molecular and cellular processes that drive development and testing how these processes can be modified through evolution to produce the extensive morphological diversity observed in the phylum. Modern attempts to characterize these processes depend on some level of genomic analysis; from querying annotated gene sets to functional genomics experiments to identify candidate cis-regulatory sequences. Given how essential these data have become, it is important that researchers using available datasets or performing their own genome-scale experiments understand the nature and limitations of echinoderm genomic analyses. In this chapter we highlight the current state of echinoderm genomic data and provide methodological considerations for common approaches, including analysis of transcriptome and functional genomics datasets.

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.