A provisional regulatory gene network for specification of endomesoderm in the sea urchin embryo

An RNA interference approach for functional studies in the sea urchin and its use in analysis of nodal signaling gradients

Dicer substrate interfering RNAs (DsiRNAs) destroy targeted transcripts using the RNA-Induced Silencing Complex (RISC) through a process called RNA interference (RNAi). This process is ubiquitous among eukaryotes. Here we report the utility of DsiRNA in embryos of the sea urchin Lytechinus variegatus (Lv). Specific knockdowns phenocopy known morpholino and inhibitor knockdowns, and DsiRNA offers a useful alternative to morpholinos. Methods are described for the design of specific DsiRNAs that lead to destruction of targeted mRNA. DsiRNAs directed against pks1, an enzyme necessary for pigment production, show how successful DsiRNA perturbations are monitored by RNA in situ analysis and by qPCR to determine relative destruction of targeted mRNA. DsiRNA-based knockdowns phenocopy morpholino- and drug-based inhibition of nodal and lefty. Other knockdowns demonstrate that the RISC operates early in development as well as on genes that are first transcribed hours after gastrulation is completed. Thus, DsiRNAs effectively mediate destruction of targeted mRNA in the sea urchin embryo. The approach offers significant advantages over other widely used methods in the urchin in terms of cost, and ease of procurement, and offers sizeable experimental advantages in terms of ease of handling, injection, and knockdown validation.

An RNA interference approach for functional studies in the sea urchin and its use in analysis of Nodal signaling gradients

Dicer substrate interfering RNAs (DsiRNAs) destroy targeted transcripts using the RNA-Induced Silencing Complex (RISC) through a process called RNA interference (RNAi). This process is ubiquitous among eukaryotes. Here we report the utility of DsiRNA in embryos of the sea urchin Lytechinus variagatus (Lv). Specific knockdowns phenocopy known morpholino and inhibitor knockdowns, and DsiRNA offers a useful alternative to morpholinos. Methods for designing and obtaining specific DsiRNAs that lead to destruction of targeted mRNA are described. DsiRNAs directed against pks1, an enzyme necessary for pigment production, show how successful DsiRNA perturbations are monitored by RNA in situ analysis and by qPCR to determine relative destruction of targeted mRNA. DsiRNA-based knockdowns phenocopy morpholino- and drug-based inhibition of nodal and lefty. Other knockdowns demonstrate that the RISC operates early in development as well as on genes that are first transcribed hours after gastrulation is completed. Thus, DsiRNAs effectively mediate destruction of targeted mRNA in the sea urchin embryo. The approach offers significant advantages over other widely used methods in the urchin in terms of cost, and ease of procurement, and offers sizeable experimental advantages in terms of ease of handling, injection, and knockdown validation.

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.

The roles of ABCB1/P-glycoprotein drug transporters in regulating gut microbes and inflammation: insights from animal models, old and new

Commensal enteric bacteria have evolved systems that enable growth in the ecologic niche of the host gastrointestinal tract. Animals evolved parallel mechanisms to survive the constant exposure to bacteria and their metabolic by-products. We propose that drug transporters encompass a crucial system to managing the gut microbiome. Drug transporters are present in the apical surface of gut epithelia. They detoxify cells from small molecules and toxins (xenobiotics) in the lumen. Here, we review what is known about commensal structure in the absence of the transporter ABCB1/P-glycoprotein in mammalian models. Knockout or low-activity alleles of ABCB1 lead to dysbiosis, Crohn's disease and ulcerative colitis in mammals. However, the exact function of ABCB1 in these contexts remain unclear. We highlight emerging models-the zebrafish Danio rerio and sea urchin Lytechinus pictus-that are poised to help dissect the fundamental mechanisms of ATP-binding cassette (ABC) transporters in the tolerance of commensal and pathogenic communities in the gut. We and others hypothesize that ABCB1 plays a direct role in exporting inflammatory bacterial products from host epithelia. Interdisciplinary work in this research area will lend novel insight to the transporter-mediated pathways that impact microbiome community structure and accelerate the pathogenesis of inflammatory bowel disease when perturbed. This article is part of the theme issue 'Sculpting the microbiome: how host factors determine and respond to microbial colonization'.

Semi-automated, high-content imaging of drug transporter knockout sea urchin (Lytechinus pictus) embryos

A defining feature of sea urchins is their extreme fecundity. Urchins produce millions of transparent, synchronously developing embryos, ideal for spatial and temporal analysis of development. This biological feature has been effectively utilized for ensemble measurement of biochemical changes. However, it has been underutilized in imaging studies, where single embryo measurements are used. Here we present an example of how stable genetics and high content imaging, along with machine learning-based image analysis, can be used to exploit the fecundity and synchrony of sea urchins in imaging-based drug screens. Building upon our recently created sea urchin ABCB1 knockout line, we developed a high-throughput assay to probe the role of this drug transporter in embryos. We used high content imaging to compare accumulation and toxicity of canonical substrates and inhibitors of the transporter, including fluorescent molecules and antimitotic cancer drugs, in homozygous knockout and wildtype embryos. To measure responses from the resulting image data, we used a nested convolutional neural network, which rapidly classified embryos according to fluorescence or cell division. This approach identified sea urchin embryos with 99.8% accuracy and determined two-cell and aberrant embryos with 96.3% and 89.1% accuracy, respectively. The results revealed that ABCB1 knockout embryos accumulated the transporter substrate calcein 3.09 times faster than wildtypes. Similarly, knockouts were 4.71 and 3.07 times more sensitive to the mitotic poisons vinblastine and taxol. This study paves the way for large scale pharmacological screens in the sea urchin embryo.

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

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

CRISPR-Cas9-Mediated Gene Knockout in a Non-Model Sea Urchin, Heliocidaris crassispina

Sea urchins have been used as model organisms in developmental biology research and the genomes of several sea urchin species have been sequenced. Recently, genome editing technologies have become available for sea urchins, and methods for gene knockout using the CRISPRCas9 system have been established. Heliocidaris crassispina is an important marine fishery resource with edible gonads. Although H. crassispina has been used as a biological research material, its genome has not yet been published, and it is a non-model sea urchin for molecular biology research. However, as recent advances in genome editing technology have facilitated genome modification in non-model organisms, we applied genome editing using the CRISPR-Cas9 system to H. crassispina. In this study, we targeted genes encoding ETS transcription factor (HcEts) and pigmentation-related polyketide synthase (HcPks1). Gene fragments were isolated using primers designed by inter-specific sequence comparisons within Echinoidea. When Ets gene was targeted using two sgRNAs, one successfully introduced mutations and impaired skeletogenesis. In the Pks1 gene knockout, when two sgRNAs targeting the close vicinity of the site corresponding to the target site that showed 100% mutagenesis efficiency of the Pks1 gene in Hemicentrotus pulcherrimus, mutagenesis was not observed. However, two other sgRNAs targeting distant sites efficiently introduced mutations. In addition, Pks1 knockout H. crassispina exhibited an albino phenotype in the pluteus larvae and adult sea urchins after metamorphosis. This indicates that the CRISPRCas9 system can be used to modify the genome of the non-model sea urchin H. crassispina.

microRNA-1 regulates sea urchin skeletogenesis by directly targeting skeletogenic genes and modulating components of signaling pathways

microRNAs are evolutionarily conserved non-coding RNAs that direct post-transcriptional regulation of target transcripts. In vertebrates, microRNA-1 (miR-1) is expressed in muscle and has been found to play critical regulatory roles in vertebrate angiogenesis, a process that has been proposed to be analogous to sea urchin skeletogenesis. Results indicate that both miR-1 inhibitor and miR-1 mimic-injected larvae have significantly less F-actin enriched circumpharyngeal muscle fibers and fewer gut contractions. In addition, miR-1 regulates the positioning of skeletogenic primary mesenchyme cells (PMCs) and skeletogenesis of the sea urchin embryo. Interestingly, the gain-of-function of miR-1 leads to more severe PMC patterning and skeletal branching defects than its loss-of-function. The results suggest that miR-1 directly suppresses Ets1/2, Tbr, and VegfR7 of the skeletogenic gene regulatory network, and Nodal, and Wnt1 signaling components. This study identifies potential targets of miR-1 that impacts skeletogenesis and muscle formation and contributes to a deeper understanding of miR-1's function during development.

New insights into mesoderm and endoderm development, and the nature of the onychophoran blastopore

BACKGROUND

Early during onychophoran development and prior to the formation of the germ band, a posterior tissue thickening forms the posterior pit. Anterior to this thickening forms a groove, the embryonic slit, that marks the anterior-posterior orientation of the developing embryo. This slit is by some authors considered the blastopore, and thus the origin of the endoderm, while others argue that the posterior pit represents the blastopore. This controversy is of evolutionary significance because if the slit represents the blastopore, then this would support the amphistomy hypothesis that suggests that a slit-like blastopore in the bilaterian ancestor evolved into protostomy and deuterostomy.

RESULTS

In this paper, we summarize our current knowledge about endoderm and mesoderm development in onychophorans and provide additional data on early endoderm- and mesoderm-determining marker genes such as Blimp, Mox, and the T-box genes.

CONCLUSION

We come to the conclusion that the endoderm of onychophorans forms prior to the development of the embryonic slit, and thus that the slit is not the primary origin of the endoderm. It is thus unlikely that the embryonic slit represents the blastopore. We suggest instead that the posterior pit indeed represents the lips of the blastopore, and that the embryonic slit (and surrounding tissue) represents a morphologically superficial archenteron-like structure. We conclude further that both endoderm and mesoderm development are under control of conserved gene regulatory networks, and that many of the features found in arthropods including the model Drosophila melanogaster are likely derived.

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.

The crucial role of CTCF in mitotic progression during early development of sea urchin

CCCTC-binding factor (CTCF), an insulator protein with 11 zinc fingers, is enriched at the boundaries of topologically associated domains (TADs) in eukaryotic genomes. In this study, we isolated and analyzed the cDNAs encoding HpCTCF, the CTCF homolog in the sea urchin Hemicentrotus pulcherrimus, to investigate its expression patterns and functions during the early development of sea urchin. HpCTCF contains nine zinc fingers corresponding to fingers 2-10 of the vertebrate CTCF. Expression pattern analysis revealed that HpCTCF mRNA was detected at all developmental stages and in the entire embryo. Upon expressing the HpCTCF-GFP fusion protein in early embryos, we observed its uniform distribution within interphase nuclei. However, during mitosis, it disappeared from the chromosomes and subsequently reassembled on the chromosome during telophase. Moreover, the morpholino-mediated knockdown of HpCTCF resulted in mitotic arrest during the morula to blastula stage. Most of the arrested chromosomes were not phospholylated at serine 10 of histone H3, indicating that mitosis was arrested at the telophase by HpCTCF depletion. Furthermore, impaired sister chromatid segregation was observed using time-lapse imaging of HpCTCF-knockdown embryos. Thus, HpCTCF is essential for mitotic progression during the early development of sea urchins, especially during the telophase-to-interphase transition. However, the normal development of pluteus larvae in CRISPR-mediated HpCTCF-knockout embryos suggests that disruption of zygotic HpCTCF expression has little effect on embryonic and larval development.

Comprehensive Genome-Wide Analyses of Poplar R2R3-MYB Transcription Factors and Tissue-Specific Expression Patterns under Drought Stress

R2R3-type MYB transcription factors are implicated in drought stress, which is a primary factor limiting the growth and development of woody plants. The identification of R2R3-MYB genes in the Populus trichocarpa genome has been previously reported. Nevertheless, the diversity and complexity of the conserved domain of the MYB gene caused inconsistencies in these identification results. There is still a lack of drought-responsive expression patterns and functional studies of R2R3-MYB transcription factors in Populus species. In this study, we identified a total of 210 R2R3-MYB genes in the P. trichocarpa genome, of which 207 genes were unevenly distributed across all 19 chromosomes. These poplar R2R3-MYB genes were phylogenetically divided into 23 subgroups. Collinear analysis demonstrated that the poplar R2R3-MYB genes underwent rapid expansion and that whole-genome duplication events were a dominant factor in the process of rapid gene expansion. Subcellular localization assays indicated that poplar R2R3-MYB TFs mainly played a transcriptional regulatory role in the nucleus. Ten R2R3-MYB genes were cloned from P. deltoides × P. euramericana cv. Nanlin895, and their expression patterns were tissue-specific. A majority of the genes showed similar drought-responsive expression patterns in two out of three tissues. This study provides a valid cue for further functional characterization of drought-responsive R2R3-MYB genes in poplar and provides support for the development of new poplar genotypes with elevated drought tolerance.

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.

An ancestral Wnt-Brachyury feedback loop in axial patterning and recruitment of mesoderm-determining target genes

Transcription factors are crucial drivers of cellular differentiation during animal development and often share ancient evolutionary origins. The T-box transcription factor Brachyury plays a pivotal role as an early mesoderm determinant and neural repressor in vertebrates; yet, the ancestral function and key evolutionary transitions of the role of this transcription factor remain obscure. Here, we present a genome-wide target-gene screen using chromatin immunoprecipitation sequencing in the sea anemone Nematostella vectensis, an early branching non-bilaterian, and the sea urchin Strongylocentrotus purpuratus, a representative of the sister lineage of chordates. Our analysis reveals an ancestral gene regulatory feedback loop connecting Brachyury, FoxA and canonical Wnt signalling involved in axial patterning that predates the cnidarian-bilaterian split about 700 million years ago. Surprisingly, we also found that part of the gene regulatory network controlling the fate of neuromesodermal progenitors in vertebrates was already present in the common ancestor of cnidarians and bilaterians. However, while several endodermal and neuronal Brachyury target genes are ancestrally shared, hardly any of the key mesodermal downstream targets in vertebrates are found in the sea anemone or the sea urchin. Our study suggests that a limited number of target genes involved in mesoderm formation were newly acquired in the vertebrate lineage, leading to a dramatic shift in the function of this ancestral developmental regulator.

Gene Expression Analysis of the Stress Response to Lithium, Nickel, and Zinc in Paracentrotus lividus Embryos

Many anthropogenic pollutants such as metals are discharged into the marine environment through modern sources. Among these, lithium (Li), nickel (Ni), and zinc (Zn) can interfere with biological processes in many organisms when their concentration rises. These metals are toxic to sea urchin embryos, affecting their development. Indeed, animal/vegetal and dorso/ventral embryonic axes are differently perturbed: Li is a vegetalizing agent, Ni can disrupt dorso-ventral axis, Zn can be animalizing. To address the molecular response adopted by embryos to cope with these metals or involved in the gene networks regulating embryogenesis, and to detect new biomarkers for evaluating hazards in polluted environments in a well-known in vivo model, we applied a high-throughput screening approach to sea urchin embryos. After fertilization, Paracentrotus lividus embryos were exposed to Li, Ni, and Zn for 24/48 h. At both endpoints, RNAs were analyzed by NanoString nCounter technology. By in silico analyses, we selected a panel of 127 transcripts encoding for regulatory and structural proteins, ranked in categories: Apoptosis, Defense, Immune, Nervous, Development, and Biomineralization. The data analysis highlighted the dysregulation of many genes in a metal-dependent manner. A functional annotation analysis was performed by the KEEG Orthology database. This study provides a platform for research on metals biomarkers in sea urchins.

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.

The arm of the starfish: The far-reaching applications of Patiria miniata as a model system in evolutionary, developmental, and regenerative biology

Many species of echinoderms have long been considered model research organisms in biology. Historically, much of this research has focused on the embryology of sea urchins and the use of their extensive gene regulatory networks as a tool to understand how the genome controls cell state specification and patterning. The establishment of Patiria miniata, the bat sea star, as a research organism has allowed us to expand on the concepts explored with sea urchins, viewing these genetic networks through a comparative lens, gaining great insight into the evolutionary mechanisms that shape developmental diversity. Extensive molecular tools have been developed in P. miniata, designed to explore gene expression dynamics and build gene regulatory networks. Echinoderms also have a robust set of bioinformatic and computational resources, centered around echinobase.org, an extensive database containing multiomic, developmental, and experimental resources for researchers. In addition to comparative evolutionary development, P. miniata is a promising system in its own right for studying whole body regeneration, metamorphosis and body plan development, as well as marine disease.

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.

Developmental gene regulatory network connections predicted by machine learning from gene expression data alone

Gene regulatory network (GRN) inference can now take advantage of powerful machine learning algorithms to complement traditional experimental methods in building gene networks. However, the dynamical nature of embryonic development-representing the time-dependent interactions between thousands of transcription factors, signaling molecules, and effector genes-is one of the most challenging arenas for GRN prediction. In this work, we show that successful GRN predictions for a developmental network from gene expression data alone can be obtained with the Priors Enriched Absent Knowledge (PEAK) network inference algorithm. PEAK is a noise-robust method that models gene expression dynamics via ordinary differential equations and selects the best network based on information-theoretic criteria coupled with the machine learning algorithm Elastic Net. We test our GRN prediction methodology using two gene expression datasets for the purple sea urchin, Stronglyocentrotus purpuratus, and cross-check our results against existing GRN models that have been constructed and validated by over 30 years of experimental results. Our results find a remarkably high degree of sensitivity in identifying known gene interactions in the network (maximum 81.58%). We also generate novel predictions for interactions that have not yet been described, which provide a resource for researchers to use to further complete the sea urchin GRN. Published ChIPseq data and spatial co-expression analysis further support a subset of the top novel predictions. We conclude that GRN predictions that match known gene interactions can be produced using gene expression data alone from developmental time series experiments.

Coordinated Gene Expression and Chromatin Regulation during Hydra Head Regeneration

The cnidarian model organism Hydra has long been studied for its remarkable ability to regenerate its head, which is controlled by a head organizer located near the hypostome. The canonical Wnt pathway plays a central role in head organizer function during regeneration and during bud formation, which is the asexual mode of reproduction in Hydra. However, it is unclear how shared the developmental programs of head organizer genesis are in budding and regeneration. Time-series analysis of gene expression changes during head regeneration and budding revealed a set of 298 differentially expressed genes during the 48-h head regeneration and 72-h budding time courses. In order to understand the regulatory elements controlling Hydra head regeneration, we first identified 27,137 open-chromatin elements that are open in one or more sections of the organism body or regenerating tissue. We used histone modification ChIP-seq to identify 9,998 candidate proximal promoter and 3,018 candidate enhancer-like regions respectively. We show that a subset of these regulatory elements is dynamically remodeled during head regeneration and identify a set of transcription factor motifs that are enriched in the enhancer regions activated during head regeneration. Our results show that Hydra displays complex gene regulatory structures of developmentally dynamic enhancers, which suggests that the evolution of complex developmental enhancers predates the split of cnidarians and bilaterians.

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.

Inferring gene regulatory networks from single-cell RNA-seq temporal snapshot data requires higher-order moments

Single-cell RNA sequencing (scRNA-seq) has become ubiquitous in biology. Recently, there has been a push for using scRNA-seq snapshot data to infer the underlying gene regulatory networks (GRNs) steering cellular function. To date, this aspiration remains unrealized due to technical and computational challenges. In this work we focus on the latter, which is under-represented in the literature. We took a systemic approach by subdividing the GRN inference into three fundamental components: data pre-processing, feature extraction, and inference. We observed that the regulatory signature is captured in the statistical moments of scRNA-seq data and requires computationally intensive minimization solvers to extract it. Furthermore, current data pre-processing might not conserve these statistical moments. Although our moment-based approach is a didactic tool for understanding the different compartments of GRN inference, this line of thinking—finding computationally feasible multi-dimensional statistics of data—is imperative for designing GRN inference methods.

•

Single-cell RNA-seq temporal snapshot data for detecting regulation

•

Challenges in data pre-processing, feature extraction, and network inference for GRNs

•

Encoding of regulatory information in higher-order raw moments

•

Non-linear least-squares inference for temporal scRNA-seq snapshot data

Single-cell RNA sequencing (scRNA-seq) has become ubiquitous in biology. Recently, there has been a push for using scRNA-seq snapshot data to infer the underlying gene regulatory networks (GRNs) steering cellular function. A recent benchmark of 12 GRN methods demonstrated that the algorithms struggled to predict the ground-truth GRNs and speculated that the low performance was due to the insufficient resolution in the scRNA-seq data. Rather than proposing another method, this paper focuses on how to decompose a GRN problem into three subproblems (pre-processing, feature extraction, and inference), so that the gene regulatory information is preserved in each step. Subsequently, we discuss how to best approach each of the three subproblems.

FGF signalling plays similar roles in development and regeneration of the skeleton in the brittle star Amphiura filiformis

Regeneration as an adult developmental process is in many aspects similar to embryonic development. Although many studies point out similarities and differences, no large-scale, direct and functional comparative analyses between development and regeneration of a specific cell type or structure in one animal exist. Here, we use the brittle star Amphiura filiformis to characterise the role of the FGF signalling pathway during skeletal development in embryos and arm regeneration. In both processes, we find ligands expressed in ectodermal cells that flank underlying skeletal mesenchymal cells, which express the receptors. Perturbation of FGF signalling showed inhibited skeleton formation in both embryogenesis and regeneration, without affecting other key developmental processes. Differential transcriptome analysis finds mostly differentiation genes rather than transcription factors to be downregulated in both contexts. Moreover, comparative gene analysis allowed us to discover brittle star-specific differentiation genes. In conclusion, our results show that the FGF pathway is crucial for skeletogenesis in the brittle star, as in other deuterostomes, and provide evidence for the re-deployment of a developmental gene regulatory module during regeneration.

A preliminary gene expression analysis on Paracentrotus lividus embryos exposed to UVB, Cadmium and their combination

Paracentrotus lividus is a Mediterranean and Eastern Atlantic sea urchin species, very sensitive to chemical and physical environmental changes and widely used in eco-toxicological studies. Here, we applied a high throughput screening approach on P. lividus embryos exposed to UVB radiation (UV), Cadmium Chloride (Cd) and their combination (Cd/UV), to deeply characterize the molecular responses adopted by embryos to cope with these stressors. in vitro eco-toxicological assays were performed by exposing embryos to Cd (10-4 M) soon after fertilization, to UV (200 and 400J/m2) at early stage of development, while in co-exposure experiments, Cd-exposed embryos were irradiated with UV at 200 J/m2. By NanoString nCounter technology, custom-made probes were developed and hybridized on total RNA extracted from exposed embryos at 51h after fertilization. By in silico analyses, we selected and retrieved at the NCBI nucleotide database a panel of P. lividus transcripts encoding for many regulatory and structural proteins that we ranked in categories, i.e., Apoptosis, Biomineralization, Defense, Development, Immunity, Signaling and Transcription Factors. The analysis of 127 transcripts highlighted the dysregulation of many genes, some specifically activated to cope with stress agents, others involved in the complex molecular network of genes that regulate embryo development. We revealed the downregulation of Biomineralization and Development genes and the upregulation of Defensive genes in Cd and Cd/UV embryos. Our approach, using sea urchin embryo as an in vivomodel, contributes to advance our knowledge about cellular responses to UV, Cd and their combination.

Ultrastructural and molecular analysis of the origin and differentiation of cells mediating brittle star skeletal regeneration

BACKGROUND

Regeneration is the ability to re-grow body parts or tissues after trauma, and it is widespread across metazoans. Cells involved in regeneration can arise from a pool of undifferentiated proliferative cells or be recruited from pre-existing differentiated tissues. Both mechanisms have been described in different phyla; however, the cellular and molecular mechanisms employed by different animals to restore lost tissues as well as the source of cells involved in regeneration remain largely unknown. Echinoderms are a clade of deuterostome invertebrates that show striking larval and adult regenerative abilities in all extant classes. Here, we use the brittle star Amphiura filiformis to investigate the origin and differentiation of cells involved in skeletal regeneration using a combination of microscopy techniques and molecular markers.

RESULTS

Our ultrastructural analyses at different regenerative stages identify a population of morphologically undifferentiated cells which appear in close contact with the proliferating epithelium of the regenerating aboral coelomic cavity. These cells express skeletogenic marker genes, such as the transcription factor alx1 and the differentiation genes c-lectin and msp130L, and display a gradient of morphological differentiation from the aboral coelomic cavity towards the epidermis. Cells closer to the epidermis, which are in contact with developing spicules, have the morphology of mature skeletal cells (sclerocytes), and express several skeletogenic transcription factors and differentiation genes. Moreover, as regeneration progresses, sclerocytes show a different combinatorial expression of genes in various skeletal elements.

CONCLUSIONS

We hypothesize that sclerocyte precursors originate from the epithelium of the proliferating aboral coelomic cavity. As these cells migrate towards the epidermis, they differentiate and start secreting spicules. Moreover, our study shows that molecular and cellular processes involved in skeletal regeneration resemble those used during skeletal development, hinting at a possible conservation of developmental programmes during adult regeneration. Finally, we highlight that many genes involved in echinoderm skeletogenesis also play a role in vertebrate skeleton formation, suggesting a possible common origin of the deuterostome endoskeleton pathway.

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.

Regulation of dynamic pigment cell states at single-cell resolution

Cells bearing pigment have diverse roles and are often under strict evolutionary selection. Here, we explore the regulation of pigmented cells in the purple sea urchin Strongylocentrotus purpuratus, an emerging model for diverse pigment function. We took advantage of single cell RNA-seq (scRNAseq) technology and discovered that pigment cells in the embryo segregated into two distinct populations, a mitotic cluster and a post-mitotic cluster. Gcm is essential for expression of several genes important for pigment function, but is only transiently expressed in these cells. We discovered unique genes expressed by pigment cells and test their expression with double fluorescence in situ hybridization. These genes include new members of the fmo family that are expressed selectively in pigment cells of the embryonic and in the coelomic cells of the adult - both cell-types having immune functions. Overall, this study identifies nodes of molecular intersection ripe for change by selective evolutionary pressures.

Spatiotemporal control of gene expression boundaries using a feedforward loop

BACKGROUND

A feedforward loop (FFL) is commonly observed in several biological networks. The FFL network motif has been mostly studied with respect to variation of the input signal in time, with only a few studies of FFL activity in a spatially distributed system such as morphogen-mediated tissue patterning. However, most morphogen gradients also evolve in time.

RESULTS

We studied the spatiotemporal behavior of a coherent FFL in two contexts: (a) a generic, oscillating morphogen gradient and (b) the dorsal-ventral patterning of the early Drosophila embryo by a gradient of the NF-κB homolog dorsal with its early target Twist. In both models, we found features in the dynamics of the intermediate node-phase difference and noise filtering-that were largely independent of the parameterization of the models, and thus were functions of the structure of the FFL itself. In the dorsal gradient model, we also found that proper target gene expression was not possible without including the effect of maternal pioneer factor Zelda.

CONCLUSIONS

An FFL buffers fluctuation to changes in the morphogen signal ensuring stable gene expression boundaries.

Evolutionary Dynamics of the SKN-1 → MED → END-1,3 Regulatory Gene Cascade in Caenorhabditis Endoderm Specification

Gene regulatory networks and their evolution are important in the study of animal development. In the nematode, Caenorhabditis elegans, the endoderm (gut) is generated from a single embryonic precursor, E. Gut is specified by the maternal factor SKN-1, which activates the MED → END-1,3 → ELT-2,7 cascade of GATA transcription factors. In this work, genome sequences from over two dozen species within the Caenorhabditis genus are used to identify MED and END-1,3 orthologs. Predictions are validated by comparison of gene structure, protein conservation, and putative cis-regulatory sites. All three factors occur together, but only within the Elegans supergroup, suggesting they originated at its base. The MED factors are the most diverse and exhibit an unexpectedly extensive gene amplification. In contrast, the highly conserved END-1 orthologs are unique in nearly all species and share extended regions of conservation. The END-1,3 proteins share a region upstream of their zinc finger and an unusual amino-terminal poly-serine domain exhibiting high codon bias. Compared with END-1, the END-3 proteins are otherwise less conserved as a group and are typically found as paralogous duplicates. Hence, all three factors are under different evolutionary constraints. Promoter comparisons identify motifs that suggest the SKN-1, MED, and END factors function in a similar gut specification network across the Elegans supergroup that has been conserved for tens of millions of years. A model is proposed to account for the rapid origin of this essential kernel in the gut specification network, by the upstream intercalation of duplicate genes into a simpler ancestral network.

Classification-Based Inference of Dynamical Models of Gene Regulatory Networks

Cell-fate decisions during development are controlled by densely interconnected gene regulatory networks (GRNs) consisting of many genes. Inferring and predictively modeling these GRNs is crucial for understanding development and other physiological processes. Gene circuits, coupled differential equations that represent gene product synthesis with a switch-like function, provide a biologically realistic framework for modeling the time evolution of gene expression. However, their use has been limited to smaller networks due to the computational expense of inferring model parameters from gene expression data using global non-linear optimization. Here we show that the switch-like nature of gene regulation can be exploited to break the gene circuit inference problem into two simpler optimization problems that are amenable to computationally efficient supervised learning techniques. We present FIGR (Fast Inference of Gene Regulation), a novel classification-based inference approach to determining gene circuit parameters. We demonstrate FIGR’s effectiveness on synthetic data generated from random gene circuits of up to 50 genes as well as experimental data from the gap gene system of Drosophila melanogaster, a benchmark for inferring dynamical GRN models. FIGR is faster than global non-linear optimization by a factor of 600 and its computational complexity scales much better with GRN size. On a practical level, FIGR can accurately infer the biologically realistic gap gene network in under a minute on desktop-class hardware instead of requiring hours of parallel computing. We anticipate that FIGR would enable the inference of much larger biologically realistic GRNs than was possible before.

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.

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.

Encounters across networks: Windows into principles of genomic regulation

Gene regulatory networks account for the ability of the genome to program development in complex multi-cellular organisms. Such networks are based on principles of gene regulation by combinations of transcription factors that bind to specific cis-regulatory DNA sites to activate transcription. These cis-regulatory regions mediate logic processing at each network node, enabling progressive increases in organismal complexity with development. Gene regulatory network explanations of development have been shown to account for patterning and cell type diversification in fly and sea urchin embryonic systems, where networks are characterized by fast coupling between transcriptional inputs and changes in target gene transcription rates, and crucial cis-regulatory elements are concentrated relatively close to the protein coding sequences of the target genes, thus facilitating their identification. Stem cell-based development in post-embryonic mammalian systems also depends on gene networks, but differs from the fly and sea urchin systems. First, the number of regulatory elements per gene and the distances between regulatory elements and the genes they control are considerably larger, forcing searches via genome-wide transcription factor binding surveys rather than functional assays. Second, the intrinsic timing of network state transitions can be slowed considerably by the need to undo stem-cell chromatin configurations, which presumably add stability to stem-cell states but retard responses to transcription factor changes during differentiation. The dispersed, partially redundant cis-regulatory systems controlling gene expression and the slow state transition kinetics in these systems already reveal new insights and opportunities to extend understanding of the repertoire of gene networks and regulatory system logic.

Expression of exogenous mRNAs to study gene function in echinoderm embryos

With the completion of the genome sequencing projects, a new challenge for developmental biologists is to assign a function to the thousands of genes identified. Expression of exogenous mRNAs is a powerful, versatile and rapid technique that can be used to study gene function during development of the sea urchin. This chapter describes how this technique can be used to analyze gene function in echinoderm embryos, how it can be combined with cell transplantation to perform mosaic analysis and how it can be applied to identify downstream targets genes of transcription factors and signaling pathways. We describe specific examples of the use of overexpression of mRNA to analyze gene function, mention the benefits and current limitations of the technique and emphasize the importance of using different controls to assess the specificity of the effects observed. Finally, this chapter details the different steps, vectors and protocols for in vitro production of mRNA and phenotypic analysis.

Using ATAC-seq and RNA-seq to increase resolution in GRN connectivity

Echinoderms have some of the most complete reconstructed developmental gene regulatory networks (GRN) of any embryo, accounting for the formation of most embryo tissues and organs. Yet, many nodes (genes and regulators) and their regulatory interactions are still to be uncovered. Traditionally, knockdown/knockout experiments are performed to determine regulator-gene interactions, which are individually validated by cis-regulatory analysis. Differential RNA-seq, combined with perturbation analysis, allows for genome-wide reconstruction of a GRN around given regulators; however, this level of resolution cannot determine direct interactions. ChiP-chip or ChIP-seq is better equipped for determining, genome-wide, whether binding of a given transcription factor (TF) to cis-regulatory elements occurs. Antibodies for the TFs of interest must be available, and if not, this presents a limiting factor. ATAC-seq identifies regions of open chromatin, that are typically trimethylated at H3K4, H3K36 and H3K79 (Kouzarides, 2007), for a given time point, condition, or tissue. This technology combined with RNA-seq and perturbation analysis provides high resolution of the possible functional interactions occurring during development. Additionally, ATAC-seq is less expensive than ChIP-seq, requires less starting material, and provides a global view of regulatory regions. This chapter provides detailed steps to identify potential regulatory relationships between the nodes of a GRN, given a well assembled genome, annotated with gene models, and ATAC-seq data combined with RNA-seq and knockdown experiments.

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.

Methods for the experimental and computational analysis of gene regulatory networks in sea urchins

The discovery of gene regulatory networks (GRNs) has opened a gate to access the genomic mechanisms controlling development. GRNs are systems of transcriptional regulatory circuits that control the differential specification of cell fates during development by regulating gene expression. The experimental analysis of GRNs involves a collection of methods, each revealing aspects of the overall control process. This review provides an overview of experimental and computational methods that have been successfully applied for solving developmental GRNs in the sea urchin embryo. The key in this approach is to obtain experimental evidence for functional interactions between transcription factors and regulatory DNA. In the second part of this review, a more generally applicable strategy is discussed that shows a path from experimental evidence to annotation of regulatory linkages to the generation of GRN models.

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

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

Gene networks and the evolution of plant morphology

Elaboration of morphology depends on the precise orchestration of gene expression by key regulatory genes. The hierarchy and relationship among the participating genes is commonly known as gene regulatory network (GRN). Therefore, the evolution of morphology ultimately occurs by the rewiring of gene network structures or by the co-option of gene networks to novel domains. The availability of high-resolution expression data combined with powerful statistical tools have opened up new avenues to formulate and test hypotheses on how diverse gene networks influence trait development and diversity. Here we summarize recent studies based on both big-data and genetics approaches to understand the evolution of plant form and physiology. We also discuss recent genome-wide investigations on how studying open-chromatin regions may help study the evolution of gene expression patterns.

Germ layer-specific regulation of cell polarity and adhesion gives insight into the evolution of mesoderm

In triploblastic animals, Par-proteins regulate cell-polarity and adherens junctions of both ectodermal and endodermal epithelia. But, in embryos of the diploblastic cnidarian Nematostella vectensis, Par-proteins are degraded in all cells in the bifunctional gastrodermal epithelium. Using immunohistochemistry, CRISPR/Cas9 mutagenesis, and mRNA overexpression, we describe the functional association between Par-proteins, ß-catenin, and snail transcription factor genes in N. vectensis embryos. We demonstrate that the aPKC/Par complex regulates the localization of ß-catenin in the ectoderm by stabilizing its role in cell-adhesion, and that endomesodermal epithelial cells are organized by a different cell-adhesion system than overlying ectoderm. We also show that ectopic expression of snail genes, which are expressed in mesodermal derivatives in bilaterians, is sufficient to downregulate Par-proteins and translocate ß-catenin from the junctions to the cytoplasm in ectodermal cells. These data provide molecular insight into the evolution of epithelial structure and distinct cell behaviors in metazoan embryos.

Notch signaling in the division of germ layers in bilaterian embryos

Bilaterian embryos are triploblastic organisms which develop three complete germ layers (ectoderm, mesoderm, and endoderm). While the ectoderm develops mainly from the animal hemisphere, there is diversity in the location from where the endoderm and the mesoderm arise in relation to the animal-vegetal axis, ranging from endoderm being specified between the ectoderm and mesoderm in echinoderms, and the mesoderm being specified between the ectoderm and the endoderm in vertebrates. A common feature is that part of the mesoderm segregates from an ancient bipotential endomesodermal domain. The process of segregation is noisy during the initial steps but it is gradually refined. In this review, we discuss the role of the Notch pathway in the establishment and refinement of boundaries between germ layers in bilaterians, with special focus on its interaction with the Wnt/β-catenin pathway.

Regulatory states in the developmental control of gene expression

A growing body of evidence shows that gene expression in multicellular organisms is controlled by the combinatorial function of multiple transcription factors. This indicates that not the individual transcription factors or signaling molecules, but the combination of expressed regulatory molecules, the regulatory state, should be viewed as the functional unit in gene regulation. Here, I discuss the concept of the regulatory state and its proposed role in the genome-wide control of gene expression. Recent analyses of regulatory gene expression in sea urchin embryos have been instrumental for solving the genomic control of cell fate specification in this system. Some of the approaches that were used to determine the expression of regulatory states during sea urchin embryogenesis are reviewed. Significant developmental changes in regulatory state expression leading to the distinct specification of cell fates are regulated by gene regulatory network circuits. How these regulatory state transitions are encoded in the genome is illuminated using the sea urchin endoderm-mesoderms cell fate decision circuit as an example. These observations highlight the importance of considering developmental gene regulation, and the function of individual transcription factors, in the context of regulatory states.

Omics approaches to study gene regulatory networks for development in echinoderms

Gene regulatory networks (GRNs) describe the interactions for a developmental process at a given time and space. Historically, perturbation experiments represent one of the key methods for analyzing and reconstructing a GRN, and the GRN governing early development in the sea urchin embryo stands as one of the more deeply dissected so far. As technology progresses, so do the methods used to address different biological questions. Next-generation sequencing (NGS) has become a standard experimental technique for genome and transcriptome sequencing and studies of protein-DNA interactions and DNA accessibility. While several efforts have been made toward the integration of different omics approaches for the study of the regulatory genome in many animals, in a few cases, these are applied with the purpose of reconstructing and experimentally testing developmental GRNs. Here, we review emerging approaches integrating multiple NGS technologies for the prediction and validation of gene interactions within echinoderm GRNs. These approaches can be applied to both 'model' and 'non-model' organisms. Although a number of issues still need to be addressed, advances in NGS applications, such as assay for transposase-accessible chromatin sequencing, combined with the availability of embryos belonging to different species, all separated by various evolutionary distances and accessible to experimental regulatory biology, place echinoderms in an unprecedented position for the reconstruction and evolutionary comparison of developmental GRNs. We conclude that sequencing technologies and integrated omics approaches allow the examination of GRNs on a genome-wide scale only if biological perturbation and cis-regulatory analyses are experimentally accessible, as in the case of echinoderm embryos.

Loss of αB-crystallin function in zebrafish reveals critical roles in the development of the lens and stress resistance of the heart

Genetic mutations in the human small heat shock protein αB-crystallin have been implicated in autosomal cataracts and skeletal myopathies, including heart muscle diseases (cardiomyopathy). Although these mutations lead to modulation of their chaperone activity in vitro, the in vivo functions of αB-crystallin in the maintenance of both lens transparency and muscle integrity remain unclear. This lack of information has hindered a mechanistic understanding of these diseases. To better define the functional roles of αB-crystallin, we generated loss-of-function zebrafish mutant lines by utilizing the CRISPR/Cas9 system to specifically disrupt the two αB-crystallin genes, αBa and αBb We observed lens abnormalities in the mutant lines of both genes, and the penetrance of the lens phenotype was higher in αBa than αBb mutants. This finding is in contrast with the lack of a phenotype previously reported in αB-crystallin knock-out mice and suggests that the elevated chaperone activity of the two zebrafish orthologs is critical for lens development. Besides its key role in the lens, we uncovered another critical role for αB-crystallin in providing stress tolerance to the heart. The αB-crystallin mutants exhibited hypersusceptibility to develop pericardial edema when challenged by crowding stress or exposed to elevated cortisol stress, both of which activate glucocorticoid receptor signaling. Our work illuminates the involvement of αB-crystallin in stress tolerance of the heart presumably through the proteostasis network and reinforces the critical role of the chaperone activity of αB-crystallin in the maintenance of lens transparency.

Dynamic changes in the interchromosomal interaction of early histone gene loci during development of sea urchin

The nuclear positioning and chromatin dynamics of eukaryotic genes are closely related to the regulation of gene expression, but they have not been well examined during early development, which is accompanied by rapid cell cycle progression and dynamic changes in nuclear organization, such as nuclear size and chromatin constitution. In this study, we focused on the early development of the sea urchin Hemicentrotus pulcherrimus and performed three-dimensional fluorescence in situ hybridization of gene loci encoding early histones (one of the types of histone in sea urchin). There are two non-allelic early histone gene loci per sea urchin genome. We found that during the morula stage, when the early histone gene expression levels are at their maximum, interchromosomal interactions were often formed between the early histone gene loci on separate chromosomes and that the gene loci were directed to locate to more interior positions. Furthermore, these interactions were associated with the active transcription of the early histone genes. Thus, such dynamic interchromosomal interactions may contribute to the efficient synthesis of early histone mRNA during the morula stage of sea urchin development.

Fitting structure to function in gene regulatory networks

Cascades of transcriptional regulation are the common source of the forward drive in all developmental systems. Increases in complexity and specificity of gene expression at successive stages are based on the collaboration of varied combinations of transcription factors already expressed in the cells to turn on new genes, and the logical relationships between the transcription factors acting and becoming newly expressed from stage to stage are best visualized as gene regulatory networks. However, gene regulatory networks used in different developmental contexts underlie processes that actually operate through different sets of rules, which affect the kinetics, synchronicity, and logical properties of individual network nodes. Contrasting early embryonic development in flies and sea urchins with adult mammalian hematopoietic development from stem cells, major differences are seen in transcription factor dosage dependence, the silencing or damping impacts of repression, and the impact of cellular regulatory history on the parts of the genome that are accessible to transcription factor action in a given cell type. These different features not only affect the kinds of models that can illuminate developmental mechanisms in the respective biological systems, but also reflect the evolutionary needs of these biological systems to optimize different aspects of development.