New insights from a high-resolution look at gastrulation in the sea urchin, Lytechinus variegatus

A model for boundary-driven tissue morphogenesis

Tissue deformations during morphogenesis can be active, driven by internal processes, or passive, resulting from stresses applied at their boundaries. Here, we introduce the Drosophila hindgut primordium as a model for studying boundary-driven tissue morphogenesis. We characterize its deformations and show that its complex shape changes can be a passive consequence of the deformations of the active regions of the embryo that surround it. First, we find an intermediate characteristic triangular shape in the 3D deformations of the hindgut. We construct a minimal model of the hindgut primordium as an elastic ring deformed by active midgut invagination and germ band extension on an ellipsoidal surface, which robustly captures the symmetry-breaking into this triangular shape. We then quantify the 3D kinematics of the tissue by a set of contours and discover that the hindgut deforms in two stages: an initial translation on the curved embryo surface followed by a rapid breaking of shape symmetry. We extend our model to show that the contour kinematics in both stages are consistent with our passive picture. Our results suggest that the role of in-plane deformations during hindgut morphogenesis is to translate the tissue to a region with anisotropic embryonic curvature and show that uniform boundary conditions are sufficient to generate the observed nonuniform shape change. Our work thus provides a possible explanation for the various characteristic shapes of blastopore-equivalents in different organisms and a framework for the mechanical emergence of global morphologies in complex developmental systems.

Modeling Epithelial Morphogenesis and Cell Rearrangement during Zebrafish Epiboly: Tissue Deformation, Cell-Cell Coupling, and the Mechanical Response to Stress

Morphogenesis in early development involves complex and extreme deformations in response to intra- and intercellular forces. Zebrafish epiboly, the spreading of the blastoderm to cover and engulf the large yolk cell, is a key early event that sets the stage for gastrulation, but is poorly understood. The enveloping layer (EVL), the thin squamous outer epithelium of the blastoderm, bound at its edges to the yolk cell by tight junctions, plays a central role. The marginal EVL cells receive through their tight junctions, forces generated in the yolk cell, and propagate those forces through the rest of the EVL; and the resulting EVL expansion is required for epiboly of the deep cells within the blastoderm. However, we don't know how these exogenous forces are generated and regulated; several alternative force generation mechanisms have been proposed. In order to model their relative contributions, we need a mechanical model of the EVL capable of responding to such forces and undergoing the drastic deformation of epiboly. The expanding EVL more than doubles its surface area and experiences significant shear as it deforms from a thin cap at the pole to become a complete sphere, necessarily requiring extensive internal rearrangement. We constructed an agent-based model of the EVL and its response to exogenous forces using the center-based simulation framework, Tissue Forge. Our model captures the large viscoelastic deformation of the EVL by cell rearrangement, and incorporates algorithmic strategies to accommodate these dynamic changes while maintaining tissue cohesion. Features observed in living embryos, such as the straightening of the initially ragged leading edge, also emerge in the model. We identified two key components required for realistic epiboly in the model: first, a mechanism to enable tissue remodeling by cell rearrangement without tearing the tissue, and second, a regulatory negative feedback on the forces driving EVL expansion, to synchronize the advancement of the EVL margin. We discuss the implications of these findings for the behavior of living EVL and the mechanisms that drive epiboly.

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 single-cell transcriptome analysis (scRNA-seq) developmental time courses from H. erythrogramma and Lytechinus variegatus (representing the derived and ancestral states, respectively) reveal numerous evolutionary changes in embryonic patterning. The earliest cell fate specification events and the primary signaling center are co-localized in the ancestral developmental gene regulatory network; 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.

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.

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.

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.

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.

Live imaging of echinoderm embryos to illuminate evo-devo

Echinoderm embryos have been model systems for cell and developmental biology for over 150 years, in good part because of their optical clarity. Discoveries that shaped our understanding of fertilization, cell division and cell differentiation were only possible because of the transparency of sea urchin eggs and embryos, which allowed direct observations of intracellular structures. More recently, live imaging of sea urchin embryos, coupled with fluorescence microscopy, has proven pivotal to uncovering mechanisms of epithelial to mesenchymal transition, cell migration and gastrulation. However, live imaging has mainly been performed on sea urchin embryos, while echinoderms include numerous experimentally tractable species that present interesting variation in key aspects of morphogenesis, including differences in embryo compaction and mechanisms of blastula formation. The study of such variation would allow us not only to understand how tissues are formed in echinoderms, but also to identify which changes in cell shape, cell-matrix and cell-cell contact formation are more likely to result in evolution of new embryonic shapes. Here we argue that adapting live imaging techniques to more echinoderm species will be fundamental to exploit such an evolutionary approach to the study of morphogenesis, as it will allow measuring differences in dynamic cellular behaviors - such as changes in cell shape and cell adhesion - between species. We briefly review existing methods for live imaging of echinoderm embryos and describe in detail how we adapted those methods to allow long-term live imaging of several species, namely the sea urchin Lytechinus pictus and the sea stars Patiria miniata and Patiriella regularis. We outline procedures to successfully label, mount and image early embryos for 10-16 h, from cleavage stages to early blastula. We show that data obtained with these methods allows 3D segmentation and tracking of individual cells over time, the first step to analyze how cell shape and cell contact differ among species. The methods presented here can be easily adopted by most cell and developmental biology laboratories and adapted to successfully image early embryos of additional species, therefore broadening our understanding of the evolution of morphogenesis.

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.

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.

Transcriptomic analysis of sea star development through metamorphosis to the highly derived pentameral body plan with a focus on neural transcription factors

The Echinodermata is characterized by a secondarily evolved pentameral body plan. While the evolutionary origin of this body plan has been the subject of debate, the molecular mechanisms underlying its development are poorly understood. We assembled a de novo developmental transcriptome from the embryo through metamorphosis in the sea star Parvulastra exigua. We use the asteroid model as it represents the basal-type echinoderm body architecture. Global variation in gene expression distinguished the gastrula profile and showed that metamorphic and juvenile stages were more similar to each other than to the pre-metamorphic stages, pointing to the marked changes that occur during metamorphosis. Differential expression and gene ontology (GO) analyses revealed dynamic changes in gene expression throughout development and the transition to pentamery. Many GO terms enriched during late metamorphosis were related to neurogenesis and signalling. Neural transcription factor genes exhibited clusters with distinct expression patterns. A suite of these genes was up-regulated during metamorphosis (e.g. Pax6, Eya, Hey, NeuroD, FoxD, Mbx, and Otp). In situ hybridization showed expression of neural genes in the CNS and sensory structures. Our results provide a foundation to understand the metamorphic transition in echinoderms and the genes involved in development and evolution of pentamery.

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.

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.

A wnt2 ortholog in the sea urchin Paracentrotus lividus

Members of the wnt gene family encode secreted glycoproteins that mediate critical intercellular communications in metazoans. Large-scale genome and transcriptome analyses have shown that this family is composed of 13 distinct subfamilies. These analyses have further established that the number of wnt genes per subfamily varies significantly between metazoan phyla, highlighting that gene duplication and gene loss events have shaped the complements of wnt genes during evolution. In sea urchins, for example, previous work reported the absence of representatives of both the WNT2 and WNT11 subfamilies in two different species, Paracentrotus lividus and Strongylocentrotus purpuratus. Recently, however, we identified a gene encoding a WNT2 ortholog in P. lividus and, based on that finding, we also reanalyzed the genome of S. purpuratus. Yet, we found no evidence of a bona fide wnt2 gene in S. purpuratus. Furthermore, we established that the P. lividus wnt2 gene is selectively expressed in vegetal tissues during embryogenesis, in a pattern that is similar, although not identical, to that of other P. lividus wnt genes. Taken together, this study amends previous work on the P. lividus wnt complement and reveals an unexpected variation in the number of wnt genes between closely related sea urchin species.

Cell rearrangement induced by filopodial tension accounts for the late phase of convergent extension in the sea urchin archenteron

George Oster was a pioneer in using mechanical models to interrogate morphogenesis in animal embryos. Convergent extension is a particularly important morphogenetic process to which George Oster gave significant attention. Late elongation of the sea urchin archenteron is a classic example of convergent extension in a monolayered tube, which has been proposed to be driven by extrinsic axial tension due to the activity of secondary mesenchyme cells. Using a vertex-based mechanical model, we show that key features of archenteron elongation can be accounted for by passive cell rearrangement due to applied tension. The model mimics the cell elongation and the Poisson effect (necking) that occur in actual archenterons. We also show that, as predicted by the model, ablation of secondary mesenchyme cells late in archenteron elongation does not result in extensive elastic recoil. Moreover, blocking the addition of cells to the base of the archenteron late in archenteron elongation leads to excessive cell rearrangement consistent with tension-induced rearrangement of a smaller cohort of cells. Our mechanical simulation suggests that responsive rearrangement can account for key features of archenteron elongation and provides a useful starting point for designing future experiments to examine the mechanical properties of the archenteron.

Theoretical tool bridging cell polarities with development of robust morphologies

Despite continual renewal and damages, a multicellular organism is able to maintain its complex morphology. How is this stability compatible with the complexity and diversity of living forms? Looking for answers at protein level may be limiting as diverging protein sequences can result in similar morphologies. Inspired by the progressive role of apical-basal and planar cell polarity in development, we propose that stability, complexity, and diversity are emergent properties in populations of proliferating polarized cells. We support our hypothesis by a theoretical approach, developed to effectively capture both types of polar cell adhesions. When applied to specific cases of development – gastrulation and the origins of folds and tubes – our theoretical tool suggests experimentally testable predictions pointing to the strength of polar adhesion, restricted directions of cell polarities, and the rate of cell proliferation to be major determinants of morphological diversity and stability.

Cells have the power to organise themselves to form complex and stable structures, whether it is to create a fully shaped baby from a single egg, or to allow adult salamanders to grow a new limb after losing a leg. This ability has been scrutinised at many different levels. For example, researchers have looked at the chemical messages exchanged by cells, or they have recorded the different shapes an embryo goes through during development. However, it is still difficult to reconcile the information from these approaches into a description that makes sense at multiple scales.

When an embryo develops, sheets of cells fold and unfold to create complex 3D shapes, like the tubes that make our lungs. Moulding sheets into tubes relies on interactions between cells that are not the same in all directions. In fact, two types of asymmetry (or polarity) guide these interactions. Apical-basal polarity runs across a sheet of cells, which means that the top surface of the sheet differs from the bottom. Planar cell polarity runs along the sheet and distinguishes one end from the other. For instance, apical-basal polarity marks the inner and outer surfaces of our skin, while planar cell polarity controls the direction in which our hair grows.

Nissen et al. set out to investigate how these polarities help cells in an embryo organise themselves to form complicated folds and tubes. To do this, simple mathematical representations of both apical-basal and planar cell polarities were designed. The representations were then combined to create computer simulations of groups of cells as these divide and interact with each other.

Simulations of ‘cells’ with only apical-basal polarity were able to generate different shapes in the ‘tissues’ produced, including many found in living organisms. External conditions, such as how cells were arranged to start with, determined the resulting shape. With both apical-basal and planar cell polarities, the simulations reproduced an important change that occurs during early development. They also replicated how the tubes that transport nutrients and oxygen form.

These results show that simple properties of individual cells, such as polarities, can produce different shapes in developing tissues and organs, without the need for a complicated overarching program. Abnormal changes in cell polarity are also associated with diseases such as cancer. The mathematical model developed by Nissen et al. could therefore be a useful tool to study these events.