A gene regulatory network controlling the embryonic specification of endoderm

Reprogramming of cells during embryonic transfating: overcoming a reprogramming block

Regulative development, demonstrated by many animal embryos, is the ability to replace missing cells or parts. The underlying molecular mechanism(s) of that ability is not well understood. If sea urchin micromeres (skeletogenic cell progenitors) are removed at the 16-cell stage, early endoderm initiates a sequential switch in cell fates, called transfating. Without micromeres, other mesoderm cells are absent as well, because their specification depends on signaling from micromeres. Most mesoderm cells later return by transfating, but pigment cells do not. Single-cell RNA sequencing, tracked over time, reveals the reprogramming sequence of those replacements. Beginning with an early endoderm specification state, cells progress through endomesoderm, then mesoderm, and finally distinct skeletogenic and blastocoelar cell specification states emerge, but pigment cells do not. Rescue of pigment cells was found to be a consequence of signal timing: if Delta is expressed prior to Nodal, pigment cells return. Thus, transfating operates through a series of gene regulatory state transitions, and reprogramming fails if endogenous negative signals occur prior to positive signals in the reprogramming sequence.

Chromosomal-Level Genome Assembly of the Antarctic Sea Urchin Sterechinus neumayeri: A Model for Antarctic Invertebrate Biology

The Antarctic sea urchin Sterechinus neumayeri (Echinoida; Echinidae) is routinely used as a model organism for Antarctic biology. Here, we present a high-quality genome of S. neumayeri. This chromosomal-level assembly was generated using PacBio long-read sequencing and Hi-C chromatin conformation capture sequencing. This 885.3-Mb assembly exhibits high contiguity with a scaffold length N50 of 36.7 Mb assembled into 20 chromosomal length scaffolds. These putative chromosomes exhibit a high degree of synteny compared to other sea urchin models. We used transcript evidence gene modeling combined with sequence homology to identify 21,638 gene models that capture 97.4% of BUSCO orthologs. Among these, we were able to identify and annotate conserved developmental gene regulatory network orthologs, positioning S. neumayeri as a tractable model for comparative studies on evolution and development.

Contrasting the development of larval and adult body plans during the evolution of biphasic lifecycles in sea urchins

Biphasic lifecycles are widespread among animals, but little is known about how the developmental transition between larvae and adults is regulated. Sea urchins are a unique system for studying this phenomenon because of the stark differences between their bilateral larval and pentaradial adult body plans. Here, we use single-cell RNA sequencing to analyze the development of Heliocidaris erythrogramma (He), a sea urchin species with an accelerated, non-feeding mode of larval development. The sequencing time course extends from embryogenesis to roughly a day before the onset of metamorphosis in He larvae, which is a period that has not been covered by previous datasets. We find that the non-feeding developmental strategy of He is associated with several changes in the specification of larval cell types compared to sea urchins with feeding larvae, such as the loss of a larva-specific skeletal cell population. Furthermore, the development of the larval and adult body plans in sea urchins may utilize largely different sets of regulatory genes. These findings lay the groundwork for extending existing developmental gene regulatory networks to cover additional stages of biphasic lifecycles.

Combinatorial regulatory states define cell fate diversity during embryogenesis

Cell fate specification occurs along invariant species-specific trajectories that define the animal body plan. This process is controlled by gene regulatory networks that regulate the expression of the limited set of transcription factors encoded in animal genomes. Here we globally assess the spatial expression of ~90% of expressed transcription factors during sea urchin development from embryo to larva to determine the activity of gene regulatory networks and their regulatory states during cell fate specification. We show that >200 embryonically expressed transcription factors together define >70 cell fates that recapitulate the morphological and functional organization of this organism. Most cell fate-specific regulatory states consist of ~15-40 transcription factors with similarity particularly among functionally related cell types regardless of developmental origin. Temporally, regulatory states change continuously during development, indicating that progressive changes in regulatory circuit activity determine cell fate specification. We conclude that the combinatorial expression of transcription factors provides molecular definitions that suffice for the unique specification of cell states in time and space during embryogenesis.

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

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

A neural network-based model framework for cell-fate decisions and development

Gene regulatory networks (GRNs) fulfill the essential function of maintaining the stability of cellular differentiation states by sustaining lineage-specific gene expression, while driving the progression of development. However, accounting for the relative stability of intermediate differentiation stages and their divergent trajectories remains a major challenge for models of developmental biology. Here, we develop an empirical data-based associative GRN model (AGRN) in which regulatory networks store multilineage stage-specific gene expression profiles as associative memory patterns. These networks are capable of responding to multiple instructive signals and, depending on signal timing and identity, can dynamically drive the differentiation of multipotent cells toward different cell state attractors. The AGRN dynamics can thus generate diverse lineage-committed cell populations in a robust yet flexible manner, providing an attractor-based explanation for signal-driven cell fate decisions during differentiation and offering a readily generalizable modelling tool that can be applied to a wide variety of cell specification systems.

SPREd: a simulation-supervised neural network tool for gene regulatory network reconstruction

Summary

Reconstruction of gene regulatory networks (GRNs) from expression data is a significant open problem. Common approaches train a machine learning (ML) model to predict a gene's expression using transcription factors' (TFs') expression as features and designate important features/TFs as regulators of the gene. Here, we present an entirely different paradigm, where GRN edges are directly predicted by the ML model. The new approach, named "SPREd," is a simulation-supervised neural network for GRN inference. Its inputs comprise expression relationships (e.g. correlation, mutual information) between the target gene and each TF and between pairs of TFs. The output includes binary labels indicating whether each TF regulates the target gene. We train the neural network model using synthetic expression data generated by a biophysics-inspired simulation model that incorporates linear as well as non-linear TF-gene relationships and diverse GRN configurations. We show SPREd to outperform state-of-the-art GRN reconstruction tools GENIE3, ENNET, PORTIA, and TIGRESS on synthetic datasets with high co-expression among TFs, similar to that seen in real data. A key advantage of the new approach is its robustness to relatively small numbers of conditions (columns) in the expression matrix, which is a common problem faced by existing methods. Finally, we evaluate SPREd on real data sets in yeast that represent gold-standard benchmarks of GRN reconstruction and show it to perform significantly better than or comparably to existing methods. In addition to its high accuracy and speed, SPREd marks a first step toward incorporating biophysics principles of gene regulation into ML-based approaches to GRN reconstruction.

Availability and implementation

Data and code are available from https://github.com/iiiime/SPREd.

Evolutionarily conserved Wnt/Sp5 signaling is critical for anterior-posterior axis patterning in sea urchin embryos

Studies across a diverse group of metazoan embryos indicate that Wnt signaling often activates the transcription factor Sp5, forming a signaling 'cassette' that plays critical roles in many developmental processes. This study explores the role of Wnt/Sp5 signaling during the specification and patterning of the primary germ layers during early anterior-posterior axis formation in the deuterostome sea urchin embryo. Our functional analyses show that Sp5 is critical for endomesoderm specification downstream of Wnt/β-catenin in posterior cells as well as anterior neuroectoderm patterning downstream of non-canonical Wnt/JNK signaling in anterior cells. Interestingly, expression and functional data comparisons show that Wnt/Sp5 signaling often plays similar roles in posterior endomesoderm as well as neuroectoderm patterning along the AP axis of several deuterostome embryos, including vertebrates. Thus, our findings provide strong support for the idea that Wnt-Sp5 signaling cassettes were critical for the establishment of early germ layers in the common deuterostome ancestor.

SPREd: A simulation-supervised neural network tool for gene regulatory network reconstruction

Reconstruction of gene regulatory networks (GRNs) from expression data is a significant open problem. Common approaches train a machine learning (ML) model to predict a gene's expression using transcription factors' (TFs') expression as features and designate important features/TFs as regulators of the gene. Here, we present an entirely different paradigm, where GRN edges are directly predicted by the ML model. The new approach, named "SPREd" is a simulation-supervised neural network for GRN inference. Its inputs comprise expression relationships (e.g., correlation, mutual information) between the target gene and each TF and between pairs of TFs. The output includes binary labels indicating whether each TF regulates the target gene. We train the neural network model using synthetic expression data generated by a biophysics-inspired simulation model that incorporates linear as well as non-linear TF-gene relationships and diverse GRN configurations. We show SPREd to outperform state-of-the-art GRN reconstruction tools GENIE3, ENNET, PORTIA and TIGRESS on synthetic datasets with high co-expression among TFs, similar to that seen in real data. A key advantage of the new approach is its robustness to relatively small numbers of conditions (columns) in the expression matrix, which is a common problem faced by existing methods. Finally, we evaluate SPREd on real data sets in yeast that represent gold standard benchmarks of GRN reconstruction and show it to perform significantly better than or comparably to existing methods. In addition to its high accuracy and speed, SPREd marks a first step towards incorporating biophysics principles of gene regulation into ML-based approaches to GRN reconstruction.

Hybrid Epigenomes Reveal Extensive Local Genetic Changes to Chromatin Accessibility Contribute to Divergence in Embryonic Gene Expression Between Species

Chromatin accessibility plays an important role in shaping gene expression, yet little is known about the genetic and molecular mechanisms that influence the evolution of chromatin configuration. Both local (cis) and distant (trans) genetic influences can in principle influence chromatin accessibility and are based on distinct molecular mechanisms. We, therefore, sought to characterize the role that each of these plays in altering chromatin accessibility in 2 closely related sea urchin species. Using hybrids of Heliocidaris erythrogramma and Heliocidaris tuberculata, and adapting a statistical framework previously developed for the analysis of cis and trans influences on the transcriptome, we examined how these mechanisms shape the regulatory landscape at 3 important developmental stages, and compared our results to similar analyses of the transcriptome. We found extensive cis- and trans-based influences on evolutionary changes in chromatin, with cis effects generally larger in effect. Evolutionary changes in accessibility and gene expression are correlated, especially when expression has a local genetic basis. Maternal influences appear to have more of an effect on chromatin accessibility than on gene expression, persisting well past the maternal-to-zygotic transition. Chromatin accessibility near gene regulatory network genes appears to be distinctly regulated, with trans factors appearing to play an outsized role in the configuration of chromatin near these genes. Together, our results represent the first attempt to quantify cis and trans influences on evolutionary divergence in chromatin configuration in an outbred natural study system and suggest that chromatin regulation is more genetically complex than was previously appreciated.

Widespread priming of transcriptional regulatory elements by incipient accessibility or RNA polymerase II pause in early embryos of the sea urchin Strongylocentrotus purpuratus

Transcriptional regulatory elements (TREs) are the primary nodes that control developmental gene regulatory networks. In embryo stages, larvae, and adult differentiated red spherule cells of the sea urchin Strongylocentrotus purpuratus, transcriptionally engaged TREs are detected by Precision Run-On Sequencing (PRO-seq), which maps genome-wide at base pair resolution the location of paused or elongating RNA polymerase II (Pol II). In parallel, TRE accessibility is estimated by the Assay for Transposase-Accessible Chromatin using Sequencing (ATAC-seq). Our analysis identifies surprisingly early and widespread TRE accessibility in 4-cell cleavage embryos that is not necessarily followed by concurrent or subsequent transcription. TRE transcriptional differences identified by PRO-seq provide more contrast among embryonic stages than ATAC-seq accessibility differences, in agreement with the apparent excess of accessible but inactive TREs during embryogenesis. Global TRE accessibility reaches a maximum around the 20-hour late blastula stage, which coincides with the consolidation of major embryo regionalizations and peak histone variant H2A.Z expression. A transcriptional potency model based on labile nucleosome TRE occupancy driven by DNA sequences and the prevalence of histone variants is proposed in order to explain the basal accessibility of transcriptionally inactive TREs during embryogenesis. However, our results would not reconcile well with labile nucleosome models based on simple A/T sequence enrichment. In addition, a large number of distal TREs become transcriptionally disengaged during developmental progression, in support of an early Pol II paused model for developmental gene regulation that eventually resolves in transcriptional activation or silencing. Thus, developmental potency in early embryos may be facilitated by incipient accessibility and transcriptional pause at TREs.

microRNA-124 directly suppresses Nodal and Notch to regulate mesodermal development

MicroRNAs regulate gene expression post-transcriptionally by destabilizing and/or inhibiting translation of target mRNAs in animal cells. MicroRNA-124 (miR-124) has been examined mostly in the context of neurogenesis. This study discovers a novel role of miR-124 in regulating mesodermal cell differentiation in the sea urchin embryo. The expression of miR-124 is first detectable at 12hours post fertilization at the early blastula stage, during endomesodermal specification. Mesodermally-derived immune cells come from the same progenitor cells that give rise to blastocoelar cells (BCs) and pigment cells (PCs) that must make a binary fate decision. We determined that miR-124 directly represses Nodal and Notch to regulate BC and PC differentiation. miR-124 inhibition does not impact the dorsal-ventral axis formation, but result in a significant increase in number of cells expressing BC-specific transcription factors (TFs) and a concurrent reduction of differentiated PCs. In general, removing miR-124's suppression of Nodal phenocopies miR124 inhibition. Interestingly, removing miR-124's suppression of Notch leads to an increased number of both BCs and PCs, with a subset of hybrid cells that express both BC- and PC-specific TFs in the larvae. Removal of miR-124's suppression of Notch not only affects differentiation of both BCs and PCs, but also induces cell proliferation of these cells during the first wave of Notch signaling. This study demonstrates that post-transcriptional regulation by miR-124 impacts differentiation of BCs and PCs by regulating the Nodal and Notch signaling pathways.

Feedback circuits are numerous in embryonic gene regulatory networks and offer a stabilizing influence on evolution of those networks

The developmental gene regulatory networks (dGRNs) of two sea urchin species, Lytechinus variegatus (Lv) and Strongylocentrotus purpuratus (Sp), have remained remarkably similar despite about 50 million years since a common ancestor. Hundreds of parallel experimental perturbations of transcription factors with similar outcomes support this conclusion. A recent scRNA-seq analysis suggested that the earliest expression of several genes within the dGRNs differs between Lv and Sp. Here, we present a careful reanalysis of the dGRNs in these two species, paying close attention to timing of first expression. We find that initial expression of genes critical for cell fate specification occurs during several compressed time periods in both species. Previously unrecognized feedback circuits are inferred from the temporally corrected dGRNs. Although many of these feedbacks differ in location within the respective GRNs, the overall number is similar between species. We identify several prominent differences in timing of first expression for key developmental regulatory genes; comparison with a third species indicates that these heterochronies likely originated in an unbiased manner with respect to embryonic cell lineage and evolutionary branch. Together, these results suggest that interactions can evolve even within highly conserved dGRNs and that feedback circuits may buffer the effects of heterochronies in the expression of key regulatory genes.

The evolutionary history of Brachyury genes in Hydrozoa involves duplications, divergence, and neofunctionalization

Brachyury, a member of T-box gene family, is widely known for its major role in mesoderm specification in bilaterians. It is also present in non-bilaterian metazoans, such as cnidarians, where it acts as a component of an axial patterning system. In this study, we present a phylogenetic analysis of Brachyury genes within phylum Cnidaria, investigate differential expression and address a functional framework of Brachyury paralogs in hydrozoan Dynamena pumila. Our analysis indicates two duplication events of Brachyury within the cnidarian lineage. The first duplication likely appeared in the medusozoan ancestor, resulting in two copies in medusozoans, while the second duplication arose in the hydrozoan ancestor, resulting in three copies in hydrozoans. Brachyury1 and 2 display a conservative expression pattern marking the oral pole of the body axis in D. pumila. On the contrary, Brachyury3 expression was detected in scattered presumably nerve cells of the D. pumila larva. Pharmacological modulations indicated that Brachyury3 is not under regulation of cWnt signaling in contrast to the other two Brachyury genes. Divergence in expression patterns and regulation suggest neofunctionalization of Brachyury3 in hydrozoans.

Computational approaches to understand transcription regulation in development

Gene regulatory networks (GRNs) serve as useful abstractions to understand transcriptional dynamics in developmental systems. Computational prediction of GRNs has been successfully applied to genome-wide gene expression measurements with the advent of microarrays and RNA-sequencing. However, these inferred networks are inaccurate and mostly based on correlative rather than causative interactions. In this review, we highlight three approaches that significantly impact GRN inference: (1) moving from one genome-wide functional modality, gene expression, to multi-omics, (2) single cell sequencing, to measure cell type-specific signals and predict context-specific GRNs, and (3) neural networks as flexible models. Together, these experimental and computational developments have the potential to significantly impact the quality of inferred GRNs. Ultimately, accurately modeling the regulatory interactions between transcription factors and their target genes will be essential to understand the role of transcription factors in driving developmental gene expression programs and to derive testable hypotheses for validation.

Gene regulatory divergence amongst echinoderms underlies appearance of pigment cells in sea urchin development

Larvae of the sea urchin, Strongylocentrotus purpuratus, have pigmented migratory cells implicated in immune defense and gut patterning. The transcription factor SpGcm activates the expression of many pigment cell-specific genes, including those involved in pigment biosynthesis (SpPks1 and SpFmo3) and immune related genes (e.g. SpMif5). Despite the importance of this cell type in sea urchins, pigmented cells are absent in larvae of the sea star, Patiria miniata. In this study, we tested the premises that sea stars lack genes to synthesize echinochrome pigment, that the genes are present but are not expressed in the larvae, or rather that the homologous gene expression does not contribute to echinochrome synthesis. Our results show that orthologs of sea urchin pigment cell-specific genes (PmPks1, PmFmo3-1 and PmMifL1-2) are present in the sea star genome and expressed in the larvae. Although no cell lineage homologous to migratory sea urchin pigment cells is present, dynamic gene activation accomplishes a similar spatial and temporal expression profile. The mechanisms regulating the expression of these genes, though, is highly divergent. In sea stars, PmGcm lacks the central role in pigment gene expression since it is not expressed in PmPks1 and PmFmo3-1-positive cells, and knockdown of Gcm does not abrogate pigment gene expression. Pigment genes are instead expressed in the coelomic mesoderm early in development before later being expressed in the ectoderm. These findings were supported by in situ RNA hybridization and comparative scRNA-seq analyses. We conclude that simply the coexpression of Pks1 and Fmo3 orthologs in cells of the sea star is not sufficient to underlie the emergence of the larval pigment cell in the sea urchin.

Hybrid epigenomes reveal extensive local genetic changes to chromatin accessibility contribute to divergence in embryonic gene expression between species

Chromatin accessibility plays an important role in shaping gene expression patterns across development and evolution; however, little is known about the genetic and molecular mechanisms that influence chromatin configuration itself. Because cis and trans influences can both theoretically influence the accessibility of the epigenome, we sought to better characterize the role that both of these mechanisms play in altering chromatin accessibility in two closely related sea urchin species. Using hybrids of the two species, and adapting a statistical framework previously developed for the analysis of cis and trans influences on the transcriptome, we examined how these mechanisms shape the regulatory landscape at three important developmental stages, and compared our results to similar patterns in the transcriptome. We found extensive cis- and trans-based influences on evolutionary changes in chromatin, with cis effects slightly more numerous and larger in effect. Genetic mechanisms influencing gene expression and chromatin configuration are correlated, but differ in several important ways. Maternal influences also appear to have more of an effect on chromatin accessibility than on gene expression, persisting well past the maternal-to-zygotic transition. Furthermore, chromatin accessibility near GRN genes appears to be regulated differently than the rest of the epigenome, and indicates that trans factors may play an outsized role in the configuration of chromatin near these genes. Together, our results represent the first attempt to quantify cis and trans influences on evolutionary divergence in chromatin configuration in an outbred natural study system, and suggest that the regulation of chromatin is more genetically complex than was previously appreciated.

Ethanol exposure perturbs sea urchin development and disrupts developmental timing

Ethanol is a known vertebrate teratogen that causes craniofacial defects as a component of fetal alcohol syndrome (FAS). Our results show that sea urchin embryos treated with ethanol similarly show broad skeletal patterning defects, potentially analogous to the defects associated with FAS. The sea urchin larval skeleton is a simple patterning system that involves only two cell types: the primary mesenchymal cells (PMCs) that secrete the calcium carbonate skeleton and the ectodermal cells that provide migratory, positional, and differentiation cues for the PMCs. Perturbations in RA biosynthesis and Hh signaling pathways are thought to be causal for the FAS phenotype in vertebrates. Surprisingly, our results indicate that these pathways are not functionally relevant for the teratogenic effects of ethanol in developing sea urchins. We found that developmental morphology as well as the expression of some ectodermal and PMC genes was delayed by ethanol exposure. Temporal transcriptome analysis revealed significant impacts of ethanol on signaling and metabolic gene expression, and a disruption in the timing of GRN gene expression that includes both delayed and precocious gene expression throughout the specification network. We conclude that the skeletal patterning perturbations in ethanol-treated embryos likely arise from a loss of temporal synchrony within and between the instructive and responsive tissues.

Identification and interrogation of the gene regulatory network of CEBPA-double mutant acute myeloid leukemia

Acute myeloid leukemia (AML) is a heterogeneous hematological malignancy caused by mutations in genes encoding transcriptional and epigenetic regulators together with signaling genes. It is characterized by a disturbance of differentiation and abnormal proliferation of hematopoietic progenitors. We have previously shown that each AML subtype establishes its own core gene regulatory network (GRN), consisting of transcription factors binding to their target genes and imposing a specific gene expression pattern that is required for AML maintenance. In this study, we integrate gene expression, open chromatin and ChIP data with promoter-capture Hi-C data to define a refined core GRN common to all patients with CEBPA-double mutant (CEBPAN/C) AML. These mutations disrupt the structure of a major regulator of myelopoiesis. We identify the binding sites of mutated C/EBPα proteins in primary cells, we show that C/EBPα, AP-1 factors and RUNX1 colocalize and are required for AML maintenance, and we employ single cell experiments to link important network nodes to the specific differentiation trajectory from leukemic stem to blast cells. Taken together, our study provides an important resource which predicts the specific therapeutic vulnerabilities of this AML subtype in human cells.

Homologous gene regulatory networks control development of apical organs and brains in Bilateria

Apical organs are relatively simple larval nervous systems. The extent to which apical organs are evolutionarily related to the more complex nervous systems of other animals remains unclear. To identify common developmental mechanisms, we analyzed the gene regulatory network (GRN) controlling the development of the apical organ in sea urchins. We characterized the developmental expression of 30 transcription factors and identified key regulatory functions for FoxQ2, Hbn, Delta/Notch signaling, and SoxC in the patterning of the apical organ and the specification of neurons. Almost the entire set of apical transcription factors is expressed in the nervous system of worms, flies, zebrafish, frogs, and mice. Furthermore, a regulatory module controlling the axial patterning of the vertebrate brain is expressed in the ectoderm of sea urchin embryos. We conclude that GRNs controlling the formation of bilaterian nervous systems share a common origin and that the apical GRN likely resembles an ancestral regulatory program.

The temporal dynamics of the sea urchin regulome

In this work, we used Nanostring N-counter technology, to evaluate the mRNA expression level of more than 330 regulatory genes over 34 time points covering the first three days of development of the sea urchin larvae. The hierarchical clustering of the mRNAs expression levels has identified groups corresponding to the major developmental landmarks (e.g. maternal to zygotic transition and gastrulation). Furthermore, comparison with previous experiments indicates high reproducibility of mRNA level temporal dynamics across batches. Finally, we generated an online tool to visualise gene expression during sea urchin larval development. The site can be accessed at and https://www621.lamp.le.ac.uk/nanostring_app/nanostring/.

Summary: This study uses Nanostring N-counter technology to evaluate the mRNA expression level of more than 330 regulatory genes over 34 time points covering the first three days of development of the sea urchin larvae.

A mosaic of conserved and novel modes of gene expression and morphogenesis in mesoderm and muscle formation of a larval bivalve

The mesoderm gives rise to several key morphological features of bilaterian animals including endoskeletal elements and the musculature. A number of regulatory genes involved in mesoderm and/or muscle formation (e.g., Brachyury (Bra), even-skipped (eve), Mox, myosin II heavy chain (mhc)) have been identified chiefly from chordates and the ecdysozoans Drosophila and Caenorhabditis elegans, but data for non-model protostomes, especially those belonging to the ecdysozoan sister clade, Lophotrochozoa (e.g., flatworms, annelids, mollusks), are only beginning to emerge. Within the lophotrochozoans, Mollusca constitutes the most speciose and diverse phylum. Interestingly, however, information on the morphological and molecular underpinnings of key ontogenetic processes such as mesoderm formation and myogenesis remains scarce even for prominent molluscan sublineages such as the bivalves. Here, we investigated myogenesis and developmental expression of Bra, eve, Mox, and mhc in the quagga mussel Dreissena rostriformis, an invasive freshwater bivalve and an emerging model in invertebrate evodevo. We found that all four genes are expressed during mesoderm formation, but some show additional, individual sites of expression during ontogeny. While Mox and mhc are involved in early myogenesis, eve is also expressed in the embryonic shell field and Bra is additionally present in the foregut. Comparative analysis suggests that Mox has an ancestral role in mesoderm and possibly muscle formation in bilaterians, while Bra and eve are conserved regulators of mesoderm development of nephrozoans (protostomes and deuterostomes). The fully developed Dreissena veliger larva shows a highly complex muscular architecture, supporting a muscular ground pattern of autobranch bivalve larvae that includes at least a velum muscle ring, three or four pairs of velum retractors, one or two pairs of larval retractors, two pairs of foot retractors, a pedal plexus, possibly two pairs of mantle retractors, and the muscles of the pallial line, as well as an anterior and a posterior adductor. As is typical for their molluscan kin, remodelling and loss of prominent larval features such as the velum musculature and various retractor systems appear to be also common in bivalves.

Supplementary information

The online version contains supplementary material available at 10.1007/s13127-022-00569-5.

SFRP1 Negatively Modulates Pyroptosis of Fibroblast-Like Synoviocytes in Rheumatoid Arthritis: A Review

Secreted frizzled-related protein 1 (SFRP1) is a member of secretory glycoprotein SFRP family. As a primitive gene regulating cell growth, development and transformation, SFRP1 is widely expressed in human cells, including various cancer cells and fibroblast-like synoviocytes (FLS) of rheumatoid arthritis (RA). Deletion or silencing of SFRP1 involves epigenetic and other mechanisms, and participates in biological behaviors such as cell proliferation, migration and cell pyroptosis, which leads to disease progression and poor prognosis. In this review, we discuss the role of SFRP1 in the pathogenesis of RA-FLS and summarize different experimental platforms and recent research results. These are helpful for understanding the biological characteristics of SFRP1 in RA, especially the mechanism by which SFRP1 regulates RA-FLS pyroptosis through Wnt/β-catenin and Notch signaling pathways. In addition, the epigenetic regulation of SFRP1 in RA-FLS is emphasized, which may be considered as a promising biomarker and therapeutic target of RA.

Limb Mesoderm and Head Ectomesenchyme Both Express a Core Transcriptional Program During Chondrocyte Differentiation

To explain how cartilage appeared in different parts of the vertebrate body at discrete times during evolution, we hypothesize that different embryonic populations co-opted expression of a core gene regulatory network (GRN) driving chondrocyte differentiation. To test this hypothesis, laser-capture microdissection coupled with RNA-seq was used to reveal chondrocyte transcriptomes in the developing chick humerus and ceratobranchial, which are mesoderm- and neural crest-derived, respectively. During endochondral ossification, two general types of chondrocytes differentiate. Immature chondrocytes (IMM) represent the early stages of cartilage differentiation, while mature chondrocytes (MAT) undergo additional stages of differentiation, including hypertrophy and stimulating matrix mineralization and degradation. Venn diagram analyses generally revealed a high degree of conservation between chondrocyte transcriptomes of the limb and head, including SOX9, COL2A1, and ACAN expression. Typical maturation genes, such as COL10A1, IBSP, and SPP1, were upregulated in MAT compared to IMM in both limb and head chondrocytes. Gene co-expression network (GCN) analyses of limb and head chondrocyte transcriptomes estimated the core GRN governing cartilage differentiation. Two discrete portions of the GCN contained genes that were differentially expressed in limb or head chondrocytes, but these genes were enriched for biological processes related to limb/forelimb morphogenesis or neural crest-dependent processes, respectively, perhaps simply reflecting the embryonic origin of the cells. A core GRN driving cartilage differentiation in limb and head was revealed that included typical chondrocyte differentiation and maturation markers, as well as putative novel "chondrocyte" genes. Conservation of a core transcriptional program during chondrocyte differentiation in both the limb and head suggest that the same core GRN was co-opted when cartilage appeared in different regions of the skeleton during vertebrate evolution.

Gene expression and transcriptional regulation driven by transcription factors involved in congenital heart defects

BACKGROUND

Congenital heart disease (CHD) is one of the most important birth defects caused by more than one mutated gene. Mutations in the genes could cause different types of congenital heart defects including atrial septal defect (ASD), tetralogy of Fallot (TOF), and ventricular septal defect (VSD).

OBJECTIVES

Cardiac transcription factors are key players for heart development and are actively involved in controlling stress regulation of the heart. Transcription factors are sequence-specific DNA binding proteins that control the process of transcription and work in a synergistic manner. We aim to characterize core cardiac transcription factors including NKX2-5, TBX, SRF, GATA4, and MEF2, which encode homeobox and MADS domain and play a crucial role in heart development.

METHODS

In this study, we have explored the important transcription factors involved in cardiac development and genes controlling the expression and regulation process by using the bioinformatics approach.

RESULTS

We have predicted the orthologs and homologs based on their evolutionary history, conserved protein domains, functional sites, and 3D structures for better understanding and presentation of factors responsible for causing CHD. Results showed the importance of these transcription factors for normal heart functioning and development.

CONCLUSION

Understanding the molecular pathways and genetic basis of CHD will help to open a new door for the treatment of patients with cardiac defects.

Dynamics of Drosophila endoderm specification

During early Drosophila embryogenesis, a network of gene regulatory interactions orchestrates terminal patterning, playing a critical role in the subsequent formation of the gut. We utilized CRISPR gene editing at endogenous loci to create live reporters of transcription and light-sheet microscopy to monitor the individual components of the posterior gut patterning network across 90 min prior to gastrulation. We developed a computational approach for fusing imaging datasets of the individual components into a common multivariable trajectory. Data fusion revealed low intrinsic dimensionality of posterior patterning and cell fate specification in wild-type embryos. The simple structure that we uncovered allowed us to construct a model of interactions within the posterior patterning regulatory network and make testable predictions about its dynamics at the protein level. The presented data fusion strategy is a step toward establishing a unified framework that would explore how stochastic spatiotemporal signals give rise to highly reproducible morphogenetic outcomes.

Streblospio benedicti: A genetic model for understanding the evolution of development and life-history

Investigating developmental evolution usually requires comparing differences across related species to infer how phenotypic change results from embryological modifications. However, when comparing organisms from different environments, ecologies, and evolutionary histories there can be many confounding factors to finding a genetic basis for developmental differences. In the marine annelid Streblospio benedicti, there are two distinct types of offspring with independent developmental pathways that converge on the same adult phenotype. To my knowledge, S. benedicti is the only known species that has heritable (additive) genetic variation in developmental traits that results in alternative life-history strategies. Females produce either hundreds of small, swimming and feeding larvae, or dozens of large, nonfeeding larvae. The larvae differ in their morphology, ecology, and dispersal potential. This developmental dimorphism makes S. benedicti a unique and useful model for understanding how genetic changes result in developmental modifications that ultimately lead to overall life-history differences. Because the offspring phenotypes of S. benedicti are heritable, we can use forward genetics within a single evolutionary lineage to disentangle how development evolves, and which genes and regulatory mechanisms are involved.

Extreme phenotypic divergence and the evolution of development

As analyses of developmental mechanisms extend to ever more species, it becomes important to understand not just what is conserved or altered during evolution, but why. Closely related species that exhibit extreme phenotypic divergence can be uniquely informative in this regard. A case in point is the sea urchin genus Heliocidaris, which contains species that recently evolved a life history involving nonfeeding larvae following nearly half a billion years of prior evolution with feeding larvae. The resulting shift in selective regimes produced rapid and surprisingly extensive changes in developmental mechanisms that are otherwise highly conserved among echinoderm species. The magnitude and extent of these changes challenges the notion that conservation of early development in echinoderms is largely due to internal constraints that prohibit modification and instead suggests that natural selection actively maintains stability of inherently malleable trait developmental mechanisms over immense time periods. Knowing how and why natural selection changed during the evolution of nonfeeding larvae can also reveal why developmental mechanisms do and do not change in particular ways.

Identification and prediction of developmental enhancers in sea urchin embryos

Background

The transcription of developmental regulatory genes is often controlled by multiple cis-regulatory elements. The identification and functional characterization of distal regulatory elements remains challenging, even in tractable model organisms like sea urchins.

Results

We evaluate the use of chromatin accessibility, transcription and RNA Polymerase II for their ability to predict enhancer activity of genomic regions in sea urchin embryos. ATAC-seq, PRO-seq, and Pol II ChIP-seq from early and late blastula embryos are manually contrasted with experimental cis-regulatory analyses available in sea urchin embryos, with particular attention to common developmental regulatory elements known to have enhancer and silencer functions differentially deployed among embryonic territories. Using the three functional genomic data types, machine learning models are trained and tested to classify and quantitatively predict the enhancer activity of several hundred genomic regions previously validated with reporter constructs in vivo.

Conclusions

Overall, chromatin accessibility and transcription have substantial power for predicting enhancer activity. For promoter-overlapping cis-regulatory elements in particular, the distribution of Pol II is the best predictor of enhancer activity in blastula embryos. Furthermore, ATAC- and PRO-seq predictive value is stage dependent for the promoter-overlapping subset. This suggests that the sequence of regulatory mechanisms leading to transcriptional activation have distinct relevance at different levels of the developmental gene regulatory hierarchy deployed during embryogenesis.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07936-0.

Developmental single-cell transcriptomics in the Lytechinus variegatus sea urchin embryo

Using scRNA-seq coupled with computational approaches, we studied transcriptional changes in cell states of sea urchin embryos during development to the larval stage. Eighteen closely spaced time points were taken during the first 24 h of development of Lytechinus variegatus (Lv). Developmental trajectories were constructed using Waddington-OT, a computational approach to 'stitch' together developmental time points. Skeletogenic and primordial germ cell trajectories diverged early in cleavage. Ectodermal progenitors were distinct from other lineages by the 6th cleavage, although a small percentage of ectoderm cells briefly co-expressed endoderm markers that indicated an early ecto-endoderm cell state, likely in cells originating from the equatorial region of the egg. Endomesoderm cells also originated at the 6th cleavage and this state persisted for more than two cleavages, then diverged into distinct endoderm and mesoderm fates asynchronously, with some cells retaining an intermediate specification status until gastrulation. Seventy-nine out of 80 genes (99%) examined, and included in published developmental gene regulatory networks (dGRNs), are present in the Lv-scRNA-seq dataset and are expressed in the correct lineages in which the dGRN circuits operate.

Cnidarian-bilaterian comparison reveals the ancestral regulatory logic of the β-catenin dependent axial patterning

In animals, body axis patterning is based on the concentration-dependent interpretation of graded morphogen signals, which enables correct positioning of the anatomical structures. The most ancient axis patterning system acting across animal phyla relies on β-catenin signaling, which directs gastrulation, and patterns the main body axis. However, within Bilateria, the patterning logic varies significantly between protostomes and deuterostomes. To deduce the ancestral principles of β-catenin-dependent axial patterning, we investigate the oral–aboral axis patterning in the sea anemone Nematostella—a member of the bilaterian sister group Cnidaria. Here we elucidate the regulatory logic by which more orally expressed β-catenin targets repress more aborally expressed β-catenin targets, and progressively restrict the initially global, maternally provided aboral identity. Similar regulatory logic of β-catenin-dependent patterning in Nematostella and deuterostomes suggests a common evolutionary origin of these processes and the equivalence of the cnidarian oral–aboral and the bilaterian posterior–anterior body axes.

The authors show in Nematostella that the more orally expressed β-catenin targets repress the more aborally expressed β-catenin targets, thus patterning the oral-aboral axis. This likely represents the common mechanism of β-catenin-dependent axial patterning shared by Cnidaria and Bilateria.

Composite morphogenesis during embryo development

Morphogenesis drives the formation of functional living shapes. Gene expression patterns and signaling pathways define the body plans of the animal and control the morphogenetic processes shaping the embryonic tissues. During embryogenesis, a tissue can undergo composite morphogenesis resulting from multiple concomitant shape changes. While previous studies have unraveled the mechanisms that drive simple morphogenetic processes, how a tissue can undergo multiple and simultaneous changes in shape is still not known and not much explored. In this chapter, we focus on the process of concomitant tissue folding and extension that is vital for the animal since it is key for embryo gastrulation and neurulation. Recent pioneering studies focus on this problem highlighting the roles of different spatially coordinated cell mechanisms or of the synergy between different patterns of gene expression to drive composite morphogenesis.

The tolerance to hypoxia is defined by a time-sensitive response of the gene regulatory network in sea urchin embryos

Deoxygenation, the reduction of oxygen level in the oceans induced by global warming and anthropogenic disturbances, is a major threat to marine life. This change in oxygen level could be especially harmful to marine embryos that use endogenous hypoxia and redox gradients as morphogens during normal development. Here, we show that the tolerance to hypoxic conditions changes between different developmental stages of the sea urchin embryo, possibly due to the structure of the gene regulatory networks (GRNs). We demonstrate that during normal development, the bone morphogenetic protein (BMP) pathway restricts the activity of the vascular endothelial growth factor (VEGF) pathway to two lateral domains and this restriction controls proper skeletal patterning. Hypoxia applied during early development strongly perturbs the activity of Nodal and BMP pathways that affect the VEGF pathway, dorsal-ventral (DV) and skeletogenic patterning. These pathways are largely unaffected by hypoxia applied after DV-axis formation. We propose that the use of redox and hypoxia as morphogens makes the sea urchin embryo highly sensitive to environmental hypoxia during early development, but the GRN structure provides higher tolerance to hypoxia at later stages.

Summary: The use of hypoxia and redox gradients as morphogens makes sea urchin early development sensitive to environmental hypoxia. This sensitivity decreases later, possibly due to the gene regulatory network structure.

microRNA-31 regulates skeletogenesis by direct suppression of Eve and Wnt1

microRNAs (miRNAs) play a critical role in a variety of biological processes, including embryogenesis and the physiological functions of cells. Evolutionarily conserved microRNA-31 (miR-31) has been found to be involved in cancer, bone formation, and lymphatic development. We previously discovered that, in the sea urchin, miR-31 knockdown (KD) embryos have shortened dorsoventral connecting rods, mispatterned skeletogenic primary mesenchyme cells (PMCs) and shifted and expanded Vegf3 expression domain. Vegf3 itself does not contain miR-31 binding sites; however, we identified its upstream regulators Eve and Wnt1 to be directly suppressed by miR-31. Removal of miR-31's suppression of Eve and Wnt1 resulted in skeletal and PMC patterning defects, similar to miR-31 KD phenotypes. Additionally, removal of miR-31's suppression of Eve and Wnt1 results in an expansion and anterior shift in expression of Veg1 ectodermal genes, including Vegf3 in the blastulae. This indicates that miR-31 indirectly regulates Vegf3 expression through directly suppressing Eve and Wnt1. Furthermore, removing miR-31 suppression of Eve is sufficient to cause skeletogenic defects, revealing a novel regulatory role of Eve in skeletogenesis and PMC patterning. Overall, this study provides a proposed molecular mechanism of miR-31's regulation of skeletogenesis and PMC patterning through its cross-regulation of a Wnt signaling ligand and a transcription factor of the endodermal and ectodermal gene regulatory network.

An early global role for Axin is required for correct patterning of the anterior-posterior axis in the sea urchin embryo

Activation of Wnt/β-catenin (cWnt) signaling at the future posterior end of early bilaterian embryos is a highly conserved mechanism for establishing the anterior-posterior (AP) axis. Moreover, inhibition of cWnt at the anterior end is required for development of anterior structures in many deuterostome taxa. This phenomenon, which occurs around the time of gastrulation, has been fairly well characterized, but the significance of intracellular inhibition of cWnt signaling in cleavage-stage deuterostome embryos for normal AP patterning is less well understood. To investigate this process in an invertebrate deuterostome, we defined Axin function in early sea urchin embryos. Axin is ubiquitously expressed at relatively high levels in early embryos and functional analysis revealed that Axin suppresses posterior cell fates in anterior blastomeres by blocking ectopic cWnt activation in these cells. Structure-function analysis of sea urchin Axin demonstrated that only its GSK-3β-binding domain is required for cWnt inhibition. These observations and results in other deuterostomes suggest that Axin plays a crucial conserved role in embryonic AP patterning by preventing cWnt activation in multipotent early blastomeres, thus protecting them from assuming ectopic cell fates.

Summary: Axin function is required in the early sea urchin embryo to regulate nuclear β-catenin levels and prevent ectopic cell fates in multipotent early blastomeres, and to ensure correct anterior-posterior axis patterning.

Transcriptomic analysis of Nodal - and BMP- associated genes during development to the juvenile seastar in Parvulastra exigua (Asterinidae)

The molecular mechanisms underlying development of the pentameral body of adult echinoderms are poorly understood but are important to solve with respect to evolution of a unique body plan that contrasts with the bilateral body plan of other deuterostomes. As Nodal and BMP2/4 signalling is involved in axis formation in larvae and development of the echinoderm body plan, we used the developmental transcriptome generated for the asterinid seastar Parvulastra exigua to investigate the temporal expression patterns of Nodal and BMP2/4 genes from the embryo and across metamorphosis to the juvenile. For echinoderms, the Asteroidea represents the basal-type body architecture with a distinct (separated) ray structure. Parvulastra exigua has lecithotrophic development forming the juvenile soon after gastrulation providing ready access to the developing adult stage. We identified 39 genes associated with the Nodal and BMP2/4 network in the P. exigua developmental transcriptome. Clustering analysis of these genes resulted in 6 clusters with similar temporal expression patterns across development. A co-expression analysis revealed genes that have similar expression profiles as Nodal and BMP2/4. These results indicated genes that may have a regulatory relationship in patterning morphogenesis of the juvenile seastar. Developmental RNA-seq analyses of Parvulastra exigua show changes in Nodal and BMP2/4 signalling genes across the metamorphic transition. We provide the foundation for detailed analyses of this cascade in the evolution of the unusual pentameral echinoderm body and its deuterostome affinities.

Early patterning of ABCB, ABCC, and ABCG transporters establishes unique territories of small molecule transport in embryonic mesoderm and endoderm

Directed intercellular movement of diverse small molecules, including metabolites, signal molecules and xenobiotics, is a key feature of multicellularity. Networks of small molecule transporters (SMTs), including several ATP Binding Cassette (ABC) transporters, are central to this process. While small molecule transporters are well described in differentiated organs, little is known about their patterns of expression in early embryogenesis. Here we report the pattern of ABC-type SMT expression and activity during the early development of sea urchins. Of the six major ABCs in this embryo (ABCB1, -B4, -C1, -C4, -C5 and -G2), three expression patterns were observed: 1) ABCB1 and ABCC1 are first expressed ubiquitously, and then become enriched in endoderm and ectoderm-derived structures. 2) ABCC4 and ABCC5 are restricted to a ring of mesoderm in the blastula and ABCC4 is later expressed in the coelomic pouches, the embryonic niche of the primordial germ cells. 3) ABCB4 and ABCG2 are expressed exclusively in endoderm-fated cells. Assays with fluorescent substrates and inhibitors of transporters revealed a ring of ABCC4 efflux activity emanating from ABCC4⁺ mesodermal cells. Similarly, ABCB1 and ABCB4 efflux activity was observed in the developing gut, prior to the onset of feeding. This study reveals the early establishment of unique territories of small molecule transport during embryogenesis. A pattern of ABCC4/C5 expression is consistent with signaling functions during gut invagination and germ line development, while a later pattern of ABCB1/B4 and ABCG2 is consistent with roles in the embryonic gut. This work provides a conceptual framework with which to examine the function and evolution of SMT networks and to define the specific developmental pathways that drive the expression of these genes.

How haematopoiesis research became a fertile ground for regulatory network biology as pioneered by Eric Davidson

PURPOSE OF REVIEW

This historical perspective reviews how work of Eric H. Davidson was a catalyst and exemplar for explaining haematopoietic cell fate determination through gene regulation.

RECENT FINDINGS

Researchers studying blood and immune cells pioneered many of the early mechanistic investigations of mammalian gene regulatory processes. These efforts included the characterization of complex gene regulatory sequences exemplified by the globin and T-cell/B-cell receptor gene loci, as well as the identification of many key regulatory transcription factors through the fine mapping of chromosome translocation breakpoints in leukaemia patients. As the repertoire of known regulators expanded, assembly into gene regulatory network models became increasingly important, not only to account for the truism that regulatory genes do not function in isolation but also to devise new ways of extracting biologically meaningful insights from even more complex information. Here we explore how Eric H. Davidson's pioneering studies of gene regulatory network control in nonvertebrate model organisms have had an important and lasting impact on research into blood and immune cell development.

SUMMARY

The intellectual framework developed by Davidson continues to contribute to haematopoietic research, and his insistence on demonstrating logic and causality still challenges the frontier of research today.

Gene Regulatory Network Analysis and Engineering Directs Development and Vascularization of Multilineage Human Liver Organoids

Pluripotent stem cell (PSC)-derived organoids have emerged as novel multicellular models of human tissue development but display immature phenotypes, aberrant tissue fates, and a limited subset of cells. Here, we demonstrate that integrated analysis and engineering of gene regulatory networks (GRNs) in PSC-derived multilineage human liver organoids direct maturation and vascular morphogenesis in vitro. Overexpression of PROX1 and ATF5, combined with targeted CRISPR-based transcriptional activation of endogenous CYP3A4, reprograms tissue GRNs and improves native liver functions, such as FXR signaling, CYP3A4 enzymatic activity, and stromal cell reactivity. The engineered tissues possess superior liver identity when compared with other PSC-derived liver organoids and show the presence of hepatocyte, biliary, endothelial, and stellate-like cell populations in single-cell RNA-seq analysis. Finally, they show hepatic functions when studied in vivo. Collectively, our approach provides an experimental framework to direct organogenesis in vitro by systematically probing molecular pathways and transcriptional networks that promote tissue development.

Systematic comparison of sea urchin and sea star developmental gene regulatory networks explains how novelty is incorporated in early development

The extensive array of morphological diversity among animal taxa represents the product of millions of years of evolution. Morphology is the output of development, therefore phenotypic evolution arises from changes to the topology of the gene regulatory networks (GRNs) that control the highly coordinated process of embryogenesis. A particular challenge in understanding the origins of animal diversity lies in determining how GRNs incorporate novelty while preserving the overall stability of the network, and hence, embryonic viability. Here we assemble a comprehensive GRN for endomesoderm specification in the sea star from zygote through gastrulation that corresponds to the GRN for sea urchin development of equivalent territories and stages. Comparison of the GRNs identifies how novelty is incorporated in early development. We show how the GRN is resilient to the introduction of a transcription factor, pmar1, the inclusion of which leads to a switch between two stable modes of Delta-Notch signaling. Signaling pathways can function in multiple modes and we propose that GRN changes that lead to switches between modes may be a common evolutionary mechanism for changes in embryogenesis. Our data additionally proposes a model in which evolutionarily conserved network motifs, or kernels, may function throughout development to stabilize these signaling transitions.

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.

Structure encoding in DNA

It is proposed that transposons and related long non-coding RNA define the fine structure of body parts. Although morphogens have long been known to direct the formation of many gross structures in early embryonic development, they do not have the necessary precision to define a structure down to the individual cellular level. Using the distinction between procedural and declarative knowledge in information processing as an analogy, it is hypothesized that DNA encodes fine structure in a manner that is different from the genetic code for proteins. The hypothesis states that repeated or near-repeated sequences that are in transposons and non-coding RNA define body part structures. As the cells in a body part go through the epigenetic process of differentiation, the action of methylation serves to inactivate all but the relevant structure definitions and some associated cell type genes. The transposons left active will then physically modify the DNA sequence in the heterochromatin to establish the local context in the three-dimensional body part structure. This brings the encoded definition of the cell type to the histone. The histone code for that cell type starts the regulatory cascade that turns on the genes associated with that particular type of cell, transforming it from a multipotent cell to a fully differentiated cell. This mechanism creates structures in the musculoskeletal system, the organs of the body, the major parts of the brain, and other systems.

Repressive interactions in gene regulatory networks: When you have no other choice

Tightly regulated gene expression programs, orchestrated by complex interactions between transcription factors, control cell type specification during development. Repressive interactions play a critical role in these networks, facilitating decision-making between two or more alternative cell fates. Here, we use the ventral neural tube as an example to illustrate how cross repressive interactions within a network drive pattern formation and specify cell types in response to a graded patterning signal. This and other systems serve to highlight how external signals are integrated through the cis regulatory elements controlling key genes and provide insight into the molecular underpinning of the process. Even the simplest networks can lead to counterintuitive results and we argue that a combination of experimental dissection and modeling approaches will be necessary to fully understand network behavior and the underlying design principles. Studying these gene regulatory networks as a whole ultimately allows us to extract fundamental properties applicable across systems that can expand our mechanistic understanding of how organisms develop.

Inferring Causal Gene Regulatory Networks from Coupled Single-Cell Expression Dynamics Using Scribe

Here, we present Scribe (https://github.com/aristoteleo/Scribe-py), a toolkit for detecting and visualizing causal regulatory interactions between genes and explore the potential for single-cell experiments to power network reconstruction. Scribe employs restricted directed information to determine causality by estimating the strength of information transferred from a potential regulator to its downstream target. We apply Scribe and other leading approaches for causal network reconstruction to several types of single-cell measurements and show that there is a dramatic drop in performance for "pseudotime"-ordered single-cell data compared with true time-series data. We demonstrate that performing causal inference requires temporal coupling between measurements. We show that methods such as "RNA velocity" restore some degree of coupling through an analysis of chromaffin cell fate commitment. These analyses highlight a shortcoming in experimental and computational methods for analyzing gene regulation at single-cell resolution and suggest ways of overcoming it.

The gene regulatory control of sea urchin gastrulation

The cell behaviors associated with gastrulation in sea urchins have been well described. More recently, considerable progress has been made in elucidating gene regulatory networks (GRNs) that underlie the specification of early embryonic territories in this experimental model. This review integrates information from these two avenues of work. I discuss the principal cell movements that take place during sea urchin gastrulation, with an emphasis on molecular effectors of the movements, and summarize our current understanding of the gene regulatory circuitry upstream of those effectors. A case is made that GRN biology can provide a causal explanation of gastrulation, although additional analysis is needed at several levels of biological organization in order to provide a deeper understanding of this complex morphogenetic process.

Ordered patterning of the sensory system is susceptible to stochastic features of gene expression

Sensory neuron numbers and positions are precisely organized to accurately map environmental signals in the brain. This precision emerges from biochemical processes within and between cells that are inherently stochastic. We investigated impact of stochastic gene expression on pattern formation, focusing on senseless (sens), a key determinant of sensory fate in Drosophila. Perturbing microRNA regulation or genomic location of sens produced distinct noise signatures. Noise was greatly enhanced when both sens alleles were present in homologous loci such that each allele was regulated in trans by the other allele. This led to disordered patterning. In contrast, loss of microRNA repression of sens increased protein abundance but not sensory pattern disorder. This suggests that gene expression stochasticity is a critical feature that must be constrained during development to allow rapid yet accurate cell fate resolution.

Chronic stress, physiological adaptation and developmental programming of the neuroendocrine stress system

Chronic stress undermines physical and mental health, in part via dysregulation of the neuroendocrine stress system. Key to understand this dysregulation is recognizing that the problem is not stress per se, but rather its chronicity. The optimally functioning stress system is highly dynamic, and negative feedback regulation enforces transient responses to acute stressors. Chronic stress overrides this, and adaptation to the chronicity can result in persistent dysregulation by altering sensitivity thresholds critical for control of system dynamics. Such adaptation involves plasticity within the central nervous system (CNS) as well as epigenetic regulation. When it occurs during development, it can have persistent effects on neuroendocrine regulation. Understanding how chronic stress programs development of the neuroendocrine stress system requires elucidation of stress-responsive gene regulatory networks that control CNS plasticity and development.

A biphasic role of non-canonical Wnt16 signaling during early anterior-posterior patterning and morphogenesis of the sea urchin embryo

A Wnt signaling network governs early anterior-posterior (AP) specification and patterning of the deuterostome sea urchin embryo. We have previously shown that non-canonical Fzl1/2/7 signaling antagonizes the progressive posterior-to-anterior downregulation of the anterior neuroectoderm (ANE) gene regulatory network (GRN) by canonical Wnt/β-catenin and non-canonical Wnt1/Wnt8-Fzl5/8-JNK signaling. This study focuses on the non-canonical function of the Wnt16 ligand during early AP specification and patterning. Maternally supplied wnt16 is expressed ubiquitously during cleavage and zygotic wnt16 expression is concentrated in the endoderm/mesoderm beginning at mid-blastula stage. Wnt16 antagonizes the ANE restriction mechanism and this activity depends on a functional Fzl1/2/7 receptor. Our results also show that zygotic wnt16 expression depends on both Fzl5/8 and Wnt/β-catenin signaling. Furthermore, Wnt16 is necessary for the activation and/or maintenance of key regulatory endoderm/mesoderm genes and is essential for gastrulation. Together, our data show that Wnt16 has two functions during early AP specification and patterning: (1) an initial role activating the Fzl1/2/7 pathway that antagonizes the ANE restriction mechanism; and (2) a subsequent function in activating key endoderm GRN factors and the morphogenetic movements of gastrulation.

The developmental transcriptome for Lytechinus variegatus exhibits temporally punctuated gene expression changes

Embryonic development is arguably the most complex process an organism undergoes during its lifetime, and understanding this complexity is best approached with a systems-level perspective. The sea urchin has become a highly valuable model organism for understanding developmental specification, morphogenesis, and evolution. As a non-chordate deuterostome, the sea urchin occupies an important evolutionary niche between protostomes and vertebrates. Lytechinus variegatus (Lv) is an Atlantic species that has been well studied, and which has provided important insights into signal transduction, patterning, and morphogenetic changes during embryonic and larval development. The Pacific species, Strongylocentrotus purpuratus (Sp), is another well-studied sea urchin, particularly for gene regulatory networks (GRNs) and cis-regulatory analyses. A well-annotated genome and transcriptome for Sp are available, but similar resources have not been developed for Lv. Here, we provide an analysis of the Lv transcriptome at 11 timepoints during embryonic and larval development. Temporal analysis suggests that the gene regulatory networks that underlie specification are well-conserved among sea urchin species. We show that the major transitions in variation of embryonic transcription divide the developmental time series into four distinct, temporally sequential phases. Our work shows that sea urchin development occurs via sequential intervals of relatively stable gene expression states that are punctuated by abrupt transitions.