ROCK and the actomyosin network control biomineral growth and morphology during sea urchin skeletogenesis.

Biomineralization had apparently evolved independently in different phyla, using distinct minerals, organic scaffolds, and gene regulatory networks (GRNs). However, diverse eukaryotes from unicellular organisms, through echinoderms to vertebrates, use the actomyosin network during biomineralization. Specifically, the actomyosin remodeling protein, Rho-associated coiled-coil kinase (ROCK) regulates cell differentiation and gene expression in vertebrates' biomineralizing cells, yet, little is known on ROCK's role in invertebrates' biomineralization. Here, we reveal that ROCK controls the formation, growth, and morphology of the calcite spicules in the sea urchin larva. ROCK expression is elevated in the sea urchin skeletogenic cells downstream of the Vascular Endothelial Growth Factor (VEGF) signaling. ROCK inhibition leads to skeletal loss and disrupts skeletogenic gene expression. ROCK inhibition after spicule formation reduces the spicule elongation rate and induces ectopic spicule branching. Similar skeletogenic phenotypes are observed when ROCK is inhibited in a skeletogenic cell culture, indicating that these phenotypes are due to ROCK activity specifically in the skeletogenic cells. Reduced skeletal growth and enhanced branching are also observed under direct perturbations of the actomyosin network. We propose that ROCK and the actomyosin machinery were employed independently, downstream of distinct GRNs, to regulate biomineral growth and morphology in Eukaryotes.

Analysis of the P. lividus sea urchin genome highlights contrasting trends of genomic and regulatory evolution in deuterostomes.

Sea urchins are emblematic models in developmental biology and display several characteristics that set them apart from other deuterostomes. To uncover the genomic cues that may underlie these specificities, we generated a chromosome-scale genome assembly for the sea urchin Paracentrotus lividus and an extensive gene expression and epigenetic profiles of its embryonic development. We found that, unlike vertebrates, sea urchins retained ancestral chromosomal linkages but underwent very fast intrachromosomal gene order mixing. We identified a burst of gene duplication in the echinoid lineage and showed that some of these expanded genes have been recruited in novel structures (water vascular system, Aristotle's lantern, and skeletogenic micromere lineage). Finally, we identified gene-regulatory modules conserved between sea urchins and chordates. Our results suggest that gene-regulatory networks controlling development can be conserved despite extensive gene order rearrangement.

Distinct regulatory states control the elongation of individual skeletal rods in the sea urchin embryo.

BACKGROUND: Understanding how gene regulatory networks (GRNs) control developmental progression is a key to the mechanistic understanding of morphogenesis. The sea urchin larval skeletogenesis provides an excellent platform to tackle this question. In the early stages of sea urchin skeletogenesis, skeletogenic genes are uniformly expressed in the skeletogenic lineage. Yet, during skeletal elongation, skeletogenic genes are expressed in distinct spatial sub-domains. The regulation of differential gene expression during late skeletogenesis is not well understood. RESULTS: Here we reveal the dynamic expression of the skeletogenic regulatory genes that define a specific regulatory state for each pair of skeletal rods, in the sea urchin Paracentrotus lividus. The vascular endothelial growth factor (VEGF) signaling, essential for skeleton formation, specifically controls the migration of cells that form the postoral and distal anterolateral skeletogenic rods. VEGF signaling also controls the expression of regulatory genes in cells at the tips of the postoral rods, including the transcription factors Pitx1 and MyoD1. Pitx1 activity is required for normal skeletal elongation and for the expression of some of VEGF target genes. CONCLUSIONS: Our study illuminates the fine-tuning of the regulatory system during the transition from early to late skeletogenesis that gives rise to rod-specific regulatory states.

The Evolution of Biomineralization through the Co-Option of Organic Scaffold Forming Networks.

Biomineralization is the process in which organisms use minerals to generate hard structures like teeth, skeletons and shells. Biomineralization is proposed to have evolved independently in different phyla through the co-option of pre-existing developmental programs. Comparing the gene regulatory networks (GRNs) that drive biomineralization in different species could illuminate the molecular evolution of biomineralization. Skeletogenesis in the sea urchin embryo was extensively studied and the underlying GRN shows high conservation within echinoderms, larval and adult skeletogenesis. The organic scaffold in which the calcite skeletal elements form in echinoderms is a tubular compartment generated by the syncytial skeletogenic cells. This is strictly different than the organic cartilaginous scaffold that vertebrates mineralize with hydroxyapatite to make their bones. Here I compare the GRNs that drive biomineralization and tubulogenesis in echinoderms and in vertebrates. The GRN that drives skeletogenesis in the sea urchin embryo shows little similarity to the GRN that drives bone formation and high resemblance to the GRN that drives vertebrates' vascular tubulogenesis. On the other hand, vertebrates' bone-GRNs show high similarity to the GRNs that operate in the cells that generate the cartilage-like tissues of basal chordate and invertebrates that do not produce mineralized tissue. These comparisons suggest that biomineralization in deuterostomes evolved through the phylum specific co-option of GRNs that control distinct organic scaffolds to mineralization.

The biological regulation of sea urchin larval skeletogenesis - From genes to biomineralized tissue.

Biomineralization is the process in which soft organic tissues use minerals to produce shells, skeletons and teeth for various functions such as protection and physical support. The ability of the cells to control the time and place of crystal nucleation as well as crystal orientation and stiffness is far beyond the state-of-the art of human technologies. Thus, understanding the biological control of biomineralization will promote our understanding of embryo development as well as provide novel approaches for material engineering. Sea urchin larval skeletogenesis offers an excellent platform for functional analyses of both the molecular control system and mineral uptake and deposition. Here we describe the current understanding of the genetic, molecular and cellular processes that underlie sea urchin larval skeletogenesis. We portray the regulatory genes that define the specification of the skeletogenic cells and drive the various morphogenetic processes that occur in the skeletogenic lineage, including: epithelial to mesenchymal transition, cell migration, spicule cavity formation and mineral deposition into the spicule cavity. We describe recent characterizations of the size, motion and mineral concentration of the calcium-bearing vesicles in the skeletogenic cells. We review the distinct specification states within the skeletogenic lineage that drive localized skeletal growth at the tips of the spicules. Finally, we discuss the surprising similarity between the regulatory network and cellular processes that drive sea urchin skeletogenesis and those that control vertebrate vascularization. Overall, we illustrate the novel insights on the biological regulation and evolution of biomineralization, gained from studies of the sea urchin larval skeletogenesis.

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

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

Calcium-vesicles perform active diffusion in the sea urchin embryo during larval biomineralization.

Biomineralization is the process by which organisms use minerals to harden their tissues and provide them with physical support. Biomineralizing cells concentrate the mineral in vesicles that they secret into a dedicated compartment where crystallization occurs. The dynamics of vesicle motion and the molecular mechanisms that control it, are not well understood. Sea urchin larval skeletogenesis provides an excellent platform for investigating the kinetics of mineral-bearing vesicles. Here we used lattice light-sheet microscopy to study the three-dimensional (3D) dynamics of calcium-bearing vesicles in the cells of normal sea urchin embryos and of embryos where skeletogenesis is blocked through the inhibition of Vascular Endothelial Growth Factor Receptor (VEGFR). We developed computational tools for displaying 3D-volumetric movies and for automatically quantifying vesicle dynamics. Our findings imply that calcium vesicles perform an active diffusion motion in both, calcifying (skeletogenic) and non-calcifying (ectodermal) cells of the embryo. The diffusion coefficient and vesicle speed are larger in the mesenchymal skeletogenic cells compared to the epithelial ectodermal cells. These differences are possibly due to the distinct mechanical properties of the two tissues, demonstrated by the enhanced f-actin accumulation and myosinII activity in the ectodermal cells compared to the skeletogenic cells. Vesicle motion is not directed toward the biomineralization compartment, but the vesicles slow down when they approach it, and probably bind for mineral deposition. VEGFR inhibition leads to an increase of vesicle volume but hardly changes vesicle kinetics and doesn't affect f-actin accumulation and myosinII activity. Thus, calcium vesicles perform an active diffusion motion in the cells of the sea urchin embryo, with diffusion length and speed that inversely correlate with the strength of the actomyosin network. Overall, our studies provide an unprecedented view of calcium vesicle 3D-dynamics and point toward cytoskeleton remodeling as an important effector of the motion of mineral-bearing vesicles.

VEGF signaling activates the matrix metalloproteinases, MmpL7 and MmpL5 at the sites of active skeletal growth and MmpL7 regulates skeletal elongation.

Organisms can uptake minerals, shape them in different forms and generate teeth, skeletons or shells that support and protect them. Mineral uptake, trafficking and nucleation are tightly regulated by the biomineralizing cells through networks of specialized proteins. Specifically, matrix metalloproteases (MMPs) digest various extracellular substrates and allow for mineralization in the vertebrates' teeth and bones, but little is known about their role in invertebrates' systems. The sea urchin embryo provides an excellent invertebrate model for genetic and molecular studies of biomineralization. MMP inhibition prevents the growth of the calcite spicules of the sea urchin larval skeleton, however, the molecular mechanisms and genes that underlie this response are not well understood. Here we study the spatial expression and regulation of two membrane type MMPs that were found to be occluded in the sea urchin spicules, Pl-MmpL7 and Pl-MmpL5, and investigate the function of Pl-MmpL7 in skeletogenesis. The inhibition of MMPs does not change the volume of the calcium vesicles in the skeletogenic cells. The expression of Pl-MmpL7 and Pl-MmpL5 is regulated by the Vascular Endothelial Growth Factor (VEGF) signaling, from the time of skeleton initiation and on. The expression of these genes is localized to the subsets of skeletogenic cells where active spicule growth occurs throughout skeletogenesis. Downregulation of Pl-MmpL7 expression delays the growth of the skeletal rods and in some cases, strongly perturbs skeletal shape. The localized expression of Pl-MmpL7 and Pl-MmpL5 to the active growth zone and the effect of Pl-MmpL7 perturbations on skeletal growth, suggest that these genes are essential for normal spicule elongation in the sea urchin embryo.

Developmental transcriptomes of the sea star, Patiria miniata, illuminate how gene expression changes with evolutionary distance.

Understanding how changes in developmental gene expression alter morphogenesis is a fundamental problem in development and evolution. A promising approach to address this problem is to compare the developmental transcriptomes between related species. The echinoderm phylum consists of several model species that have significantly contributed to the understanding of gene regulation and evolution. Particularly, the regulatory networks of the sea star, Patiria miniata (P. miniata), have been extensively studied, however developmental transcriptomes for this species were lacking. Here we generated developmental transcriptomes of P. miniata and compared these with those of two sea urchins species. We demonstrate that the conservation of gene expression depends on gene function, cell type and evolutionary distance. With increasing evolutionary distance the interspecies correlations in gene expression decreases. The reduction is more severe in the correlations between morphologically equivalent stages (diagonal elements) than in the correlation between morphologically distinct stages (off-diagonal elements). This could reflect a decrease in the morphological constraints compared to other constraints that shape gene expression at large evolutionary divergence. Within this trend, the interspecies correlations of developmental control genes maintain their diagonality at large evolutionary distance, and peak at the onset of gastrulation, supporting the hourglass model of phylotypic stage conservation.

Possible cooption of a VEGF-driven tubulogenesis program for biomineralization in echinoderms.

Biomineralization is the process by which living organisms use minerals to form hard structures that protect and support them. Biomineralization is believed to have evolved rapidly and independently in different phyla utilizing preexisting components. The mechanistic understanding of the regulatory networks that drive biomineralization and their evolution is far from clear. Sea urchin skeletogenesis is an excellent model system for studying both gene regulation and mineral uptake and deposition. The sea urchin calcite spicules are formed within a tubular cavity generated by the skeletogenic cells controlled by vascular endothelial growth factor (VEGF) signaling. The VEGF pathway is essential for biomineralization in echinoderms, while in many other phyla, across metazoans, it controls tubulogenesis and vascularization. Despite the critical role of VEGF signaling in sea urchin spiculogenesis, the downstream program it activates was largely unknown. Here we study the cellular and molecular machinery activated by the VEGF pathway during sea urchin spiculogenesis and reveal multiple parallels to the regulation of vertebrate vascularization. Human VEGF rescues sea urchin VEGF knockdown, vesicle deposition into an internal cavity plays a significant role in both systems, and sea urchin VEGF signaling activates hundreds of genes, including biomineralization and interestingly, vascularization genes. Moreover, five upstream transcription factors and three signaling genes that drive spiculogenesis are homologous to vertebrate factors that control vascularization. Overall, our findings suggest that sea urchin spiculogenesis and vertebrate vascularization diverged from a common ancestral tubulogenesis program, broadly adapted for vascularization and specifically coopted for biomineralization in the echinoderm phylum.

Corrigendum: Comparative Studies of Gene Expression Kinetics: Methodologies and Insights on Development and Evolution.

[This corrects the article DOI: 10.3389/fgene.2018.00339.].

The network remains.

Eric Davidson was a legend both in his science and his personality. He inspired and challenged a new generation of developmental biologists and I was lucky to be one of them. He changed the way we think about biological interactions by synthesizing a large scale, almost incomprehensible set of data into a causal model of a gene regulatory network. While his death leaves a big hole in our lives, his contribution to the conceptualization of regulatory biology will inspire developmental and evolutionary biologists for decades to come.

Parallel embryonic transcriptional programs evolve under distinct constraints and may enable morphological conservation amidst adaptation.

Embryonic development evolves by balancing stringent morphological constraints with genetic and environmental variation. The design principle that allows developmental transcriptional programs to conserve embryonic morphology while adapting to environmental changes is still not fully understood. To address this fundamental challenge, we compare developmental transcriptomes of two sea urchin species, Paracentrotus lividus and Strongylocentrotus purpuratus, that shared a common ancestor about 40 million years ago and are geographically distant yet show similar morphology. We find that both developmental and housekeeping genes show highly dynamic and strongly conserved temporal expression patterns. The expression of other gene sets, including homeostasis and response genes, show divergent expression which could result from either evolutionary drift or adaptation to local environmental conditions. The interspecies correlations of developmental gene expressions are highest between morphologically similar developmental time points whereas the interspecies correlations of housekeeping gene expression are high between all the late zygotic time points. Relatedly, the position of the phylotypic stage varies between these two groups of genes: developmental gene expression shows highest conservation at mid-developmental stage, in agreement with the hourglass model while the conservation of housekeeping genes keeps increasing with developmental time. When all genes are combined, the relationship between conservation of gene expression and morphological similarity is partially masked by housekeeping genes and genes with diverged expression. Our study illustrates various transcriptional programs that coexist in the developing embryo and evolve under different constraints. Apparently, morphological constraints underlie the conservation of developmental gene expression while embryonic fitness requires the conservation of housekeeping gene expression and the species-specific adjustments of homeostasis gene expression. The distinct evolutionary forces acting on these transcriptional programs enable the conservation of similar body plans while allowing adaption.

Regulatory heterochronies and loose temporal scaling between sea star and sea urchin regulatory circuits.

It has long been argued that heterochrony, a change in relative timing of a developmental process, is a major source of evolutionary innovation. Heterochronic changes of regulatory gene activation could be the underlying molecular mechanism driving heterochronic changes through evolution. Here, we compare the temporal expression profiles of key regulatory circuits between sea urchin and sea star, representative of two classes of Echinoderms that shared a common ancestor about 500 million years ago. The morphologies of the sea urchin and sea star embryos are largely comparable, yet, differences in certain mesodermal cell types and ectodermal patterning result in distinct larval body plans. We generated high resolution temporal profiles of 17 mesodermally-, endodermally- and ectodermally-expressed regulatory genes in the sea star, Patiria miniata, and compared these to their orthologs in the Mediterranean sea urchin, Paracentrotus lividus. We found that the maternal to zygotic transition is delayed in the sea star compared to the sea urchin, in agreement with the longer cleavage stage in the sea star. Interestingly, the order of gene activation shows the highest variation in the relatively diverged mesodermal circuit, while the correlations of expression dynamics are the highest in the strongly conserved endodermal circuit. We detected loose scaling of the developmental rates of these species and observed interspecies heterochronies within all studied regulatory circuits. Thus, after 500 million years of parallel evolution, mild heterochronies between the species are frequently observed and the tight temporal scaling observed for closely related species no longer holds.

Mature maternal mRNAs are longer than zygotic ones and have complex degradation kinetics in sea urchin.

Early in embryogenesis, maternally deposited transcripts are degraded and new zygotic transcripts are generated during the maternal to zygotic transition. Recent works have shown that early zygotic transcripts are short compared to maternal transcripts, in zebrafish and Drosophila species. The reduced zygotic transcript length was attributed to the short cell cycle in these organisms that prevents the transcription of long primary transcripts (intron delay). Here we study the length of maternal mRNAs and their degradation kinetics in two sea urchin species to further the understanding of maternal gene usage and processing. Early zygotic primary transcripts and mRNAs are shorter than maternal ones in the sea urchin, Strongylocentrotus purpuratus. Yet, while primary transcripts length increases when cell cycle lengthens, typical for intron delay, the relatively short length of zygotic mRNAs is consistent. The enhanced mRNA length is due to significantly longer maternal open reading frames and 3'UTRs compared to the zygotic lengths, a ratio that does not change with developmental time. This implies unique usage of both coding sequences and regulatory information in the maternal stage compared to the zygotic stages. We extracted the half-lifetimes due to maternal and zygotic degradation mechanisms from high-density time course of a set of maternal mRNAs in Paracentrotus lividus. The degradation rates due to maternal and zygotic degradation mechanisms are not correlated, indicating that these mechanisms are independent and relay on different regulatory information. Our studies illuminate specific structural and kinetic properties of sea urchin maternal mRNAs that might be broadly shared by other organisms.

Robustness and Accuracy in Sea Urchin Developmental Gene Regulatory Networks.

Developmental gene regulatory networks robustly control the timely activation of regulatory and differentiation genes. The structure of these networks underlies their capacity to buffer intrinsic and extrinsic noise and maintain embryonic morphology. Here I illustrate how the use of specific architectures by the sea urchin developmental regulatory networks enables the robust control of cell fate decisions. The Wnt-ßcatenin signaling pathway patterns the primary embryonic axis while the BMP signaling pathway patterns the secondary embryonic axis in the sea urchin embryo and across bilateria. Interestingly, in the sea urchin in both cases, the signaling pathway that defines the axis controls directly the expression of a set of downstream regulatory genes. I propose that this direct activation of a set of regulatory genes enables a uniform regulatory response and a clear cut cell fate decision in the endoderm and in the dorsal ectoderm. The specification of the mesodermal pigment cell lineage is activated by Delta signaling that initiates a triple positive feedback loop that locks down the pigment specification state. I propose that the use of compound positive feedback circuitry provides the endodermal cells enough time to turn off mesodermal genes and ensures correct mesoderm vs. endoderm fate decision. Thus, I argue that understanding the control properties of repeatedly used regulatory architectures illuminates their role in embryogenesis and provides possible explanations to their resistance to evolutionary change.

Quantitative developmental transcriptomes of the Mediterranean sea urchin Paracentrotus lividus.

Embryonic development progresses through the timely activation of thousands of differentially activated genes. Quantitative developmental transcriptomes provide the means to relate global patterns of differentially expressed genes to the emerging body plans they generate. The sea urchin is one of the classic model systems for embryogenesis and the models of its developmental gene regulatory networks are of the most comprehensive of their kind. Thus, the sea urchin embryo is an excellent system for studies of its global developmental transcriptional profiles. Here we produced quantitative developmental transcriptomes of the sea urchin Paracentrotus lividus (P. lividus) at seven developmental stages from the fertilized egg to prism stage. We generated de-novo reference transcriptome and identified 29,817 genes that are expressed at this time period. We annotated and quantified gene expression at the different developmental stages and confirmed the reliability of the expression profiles by QPCR measurement of a subset of genes. The progression of embryo development is reflected in the observed global expression patterns and in our principle component analysis. Our study illuminates the rich patterns of gene expression that participate in sea urchin embryogenesis and provide an essential resource for further studies of the dynamic expression of P. lividus genes.

Comparative Study of Regulatory Circuits in Two Sea Urchin Species Reveals Tight Control of Timing and High Conservation of Expression Dynamics.

Accurate temporal control of gene expression is essential for normal development and must be robust to natural genetic and environmental variation. Studying gene expression variation within and between related species can delineate the level of expression variability that development can tolerate. Here we exploit the comprehensive model of sea urchin gene regulatory networks and generate high-density expression profiles of key regulatory genes of the Mediterranean sea urchin, Paracentrotus lividus (Pl). The high resolution of our studies reveals highly reproducible gene initiation times that have lower variation than those of maximal mRNA levels between different individuals of the same species. This observation supports a threshold behavior of gene activation that is less sensitive to input concentrations. We then compare Mediterranean sea urchin gene expression profiles to those of its Pacific Ocean relative, Strongylocentrotus purpuratus (Sp). These species shared a common ancestor about 40 million years ago and show highly similar embryonic morphologies. Our comparative analyses of five regulatory circuits operating in different embryonic territories reveal a high conservation of the temporal order of gene activation but also some cases of divergence. A linear ratio of 1.3-fold between gene initiation times in Pl and Sp is partially explained by scaling of the developmental rates with temperature. Scaling the developmental rates according to the estimated Sp-Pl ratio and normalizing the expression levels reveals a striking conservation of relative dynamics of gene expression between the species. Overall, our findings demonstrate the ability of biological developmental systems to tightly control the timing of gene activation and relative dynamics and overcome expression noise induced by genetic variation and growth conditions.

Gene regulatory control in the sea urchin aboral ectoderm: spatial initiation, signaling inputs, and cell fate lockdown.

The regulation of oral-aboral ectoderm specification in the sea urchin embryo has been extensively studied in recent years. The oral-aboral polarity is initially imposed downstream of a redox gradient induced by asymmetric maternal distribution of mitochondria. Two TGF-ß signaling pathways, Nodal and BMP, are then respectively utilized in the generation of oral and aboral regulatory states. However, a causal understanding of the regulation of aboral ectoderm specification has been lacking. In this work control of aboral ectoderm regulatory state specification was revealed by combining detailed regulatory gene expression studies, perturbation and cis-regulatory analyses. Our analysis illuminates a dynamic system where different factors dominate at different developmental times. We found that the initial activation of aboral genes depends directly on the redox sensitive transcription factor, hypoxia inducible factor 1α (HIF-1α). Two BMP ligands, BMP2/4 and BMP5/8, then significantly enhance aboral regulatory gene transcription. Ultimately, encoded feedback wiring lockdown the aboral ectoderm regulatory state. Our study elucidates the different regulatory mechanisms that sequentially dominate the spatial localization of aboral regulatory states.

Perturbation analysis analyzed--athematical modeling of intact and perturbed gene regulatory circuits for animal development.

Gene regulatory networks for animal development are the underlying mechanisms controlling cell fate specification and differentiation. The architecture of gene regulatory circuits determines their information processing properties and their developmental function. It is a major task to derive realistic network models from exceedingly advanced high throughput experimental data. Here we use mathematical modeling to study the dynamics of gene regulatory circuits to advance the ability to infer regulatory connections and logic function from experimental data. This study is guided by experimental methodologies that are commonly used to study gene regulatory networks that control cell fate specification. We study the effect of a perturbation of an input on the level of its downstream genes and compare between the cis-regulatory execution of OR and AND logics. Circuits that initiate gene activation and circuits that lock on the expression of genes are analyzed. The model improves our ability to analyze experimental data and construct from it the network topology. The model also illuminates information processing properties of gene regulatory circuits for animal development.