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

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.

Expression and Transcriptional Targets of TGFβ-RII in Paracentrotus lividus Larval Skeletogenesis

Organisms from the five kingdoms of life use minerals to harden their tissues and make teeth, shells and skeletons, in the process of biomineralization. The sea urchin larval skeleton is an excellent system to study the biological regulation of biomineralization and its evolution. The gene regulatory network (GRN) that controls sea urchin skeletogenesis is known in great details and shows similarity to the GRN that controls vertebrates' vascularization while it is quite distinct from the GRN that drives vertebrates' bone formation. Yet, transforming growth factor beta (TGF-β) signaling regulates both sea urchin and vertebrates' skeletogenesis. Here, we study the upstream regulation and identify transcriptional targets of TGF-β in the Mediterranean Sea urchin species, Paracentrotus lividus. TGF-βRII is transiently active in the skeletogenic cells downstream of vascular endothelial growth factor (VEGF) signaling, in P. lividus. Continuous perturbation of TGF-βRII activity significantly impairs skeletal elongation and the expression of key skeletogenic genes. Perturbation of TGF-βRII after skeletal initiation leads to a delay in skeletal elongation and minor changes in gene expression. TGF-β targets are distinct from its transcriptional targets during vertebrates' bone formation, suggesting that the role of TGF-β in biomineralization in these two phyla results from convergent evolution.

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.

PVC pellet leachates affect adult immune system and embryonic development but not reproductive capacity in the sea urchin Paracentrotus lividus

Microplastic pollution is a major concern of our age, eliciting a range of effects on organisms including during embryonic development. Plastic preproduction pellets stunt the development of sea urchins through the leaching of teratogenic compounds. However, the effect of these leachates on adult sea urchins and their fertility is unknown. Here we investigate the effect of PVC leachates on the capacity to produce normal embryos, and demonstrate that adults kept in contaminated water still produce viable offspring. However, we observe a cumulative negative effect by continued exposure to highly polluted water: adult animals had lower counts and disturbed morphological profiles of immune cells, were under increased oxidative stress, and produced embryos less tolerant of contaminated environments. Our findings suggest that even in highly polluted areas, sea urchins are fertile, but that sublethal effects seen in the adults may lead to transgenerational effects that reduce developmental robustness of the embryos.

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.

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.

Polymer aging affects the bioavailability of microplastics-associated contaminants in sea urchin embryos

Microplastics (MPs) in the marine environment undergo complex weathering factors that can affect their ability to interact with different coexisting environmental contaminants (termed here co-contaminants). In this study, the influence of artificially aging using UV on the sorption of a complex mixture of co-contaminants onto MPs was investigated in order to provide meaningful hypotheses on their individual and combined toxicities on sea urchin embryos. A mixture of artificially aged MPs (PS particles and PA microfibers) combined with 2,2',4,4'-tetrabromodiphenyl ether (BDE-47), or Cd or Cu, both alone and in a mix, were used to expose embryos of Paracentrotus lividus. The effects of polymer aging on co-contaminants bioavailability were assessed by measuring changes in the transcriptional profile of genes involved in oxidative-stress response and skeletogenic and endo-mesodermal specification. Changes in the sorption ability of MPs to co-contaminants in the aqueous phase highlighted that aging did not affect the sorption of BDE-47 and Cd on MPs, although a certain influence on Cu sorption was found. Despite no morphological effects in embryos at the gastrula stage after MPs/contaminants combinatorial exposure emerged, the greatest influence of the aging process was mainly found for combined exposures which included BDE-47. Finally, the exposure to multiple contaminants generated transcriptional profiles poorly related to those activated by single contaminant, at times suggesting a mixture-dependent different aging influence. These results open new scenarios on the controversial role of vector of co-contaminants for MPs, especially when complex and different types of mixtures were considered.

Recent reconfiguration of an ancient developmental gene regulatory network in Heliocidaris sea urchins

Changes in developmental gene regulatory networks (dGRNs) underlie much of the diversity of life, but the evolutionary mechanisms that operate on regulatory interactions remain poorly understood. Closely related species with extreme phenotypic divergence provide a valuable window into the genetic and molecular basis for changes in dGRNs and their relationship to adaptive changes in organismal traits. Here we analyse genomes, epigenomes and transcriptomes during early development in two Heliocidaris sea urchin species that exhibit highly divergent life histories and in an outgroup species. Positive selection and chromatin accessibility modifications within putative regulatory elements are enriched on the branch leading to the derived life history, particularly near dGRN genes. Single-cell transcriptomes reveal a dramatic delay in cell fate specification in the derived state, which also has far fewer open chromatin regions, especially near conserved cell fate specification genes. Experimentally perturbing key transcription factors reveals profound evolutionary changes to early embryonic patterning events, disrupting regulatory interactions previously conserved for ~225 million years. These results demonstrate that natural selection can rapidly reshape developmental gene expression on a broad scale when selective regimes abruptly change. More broadly, even highly conserved dGRNs and patterning mechanisms in the early embryo remain evolvable under appropriate ecological circumstances.

Notch Signalling Under Maternal-to-Zygotic Transition

The development of all animal embryos is initially directed by the gene products supplied by their mothers. With the progression of embryogenesis, the embryo's genome is activated to command subsequent developments. This transition, which has been studied in many model animals, is referred to as the Maternal-to-Zygotic Transition (MZT). In many organisms, including flies, nematodes, and sea urchins, genes involved in Notch signaling are extensively influenced by the MZT. This signaling pathway is highly conserved across metazoans; moreover, it regulates various developmental processes. Notch signaling defects are commonly associated with various human diseases. The maternal contribution of its factors was first discovered in flies. Subsequently, several genes were identified from mutant embryos with a phenotype similar to Notch mutants only upon the removal of the maternal contributions. Studies on these maternal genes have revealed various novel steps in the cascade of Notch signal transduction. Among these genes, pecanex and almondex have been functionally characterized in recent studies. Therefore, in this review, we will focus on the roles of these two maternal genes in Notch signaling and discuss future research directions on its maternal function.

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 skeletogenic transcription factors form specific regulatory states in various skeletogenic sub‐populations.

Late VEGF signaling controls the regulatory states at the tips of the post‐oral and anterolateral skeletal rods.

VEGF signaling controls the expression of the transcription factors, MyoD1 and Pitx1.

Pitx1 activity is required for normal skeletal elongation and for the expression of some of VEGF target genes.

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.

Dorsal-ventral axis formation in sea urchin embryos

Most sea urchin species produce planktonic feeding larvae with distinct dorsal-ventral polarity. Such morphological indicators of polarity arise after gastrulation, when several morphogenesis and cell differentiation events occur differentially along the dorsal-ventral axis. For instance, the gut bends toward the ventral side where the mouth will form, skeletogenesis occurs initially near the ventral side with the forming skeleton extending dorsally, and pigment cells differentiate and embed in the dorsal ectoderm. The patterning mechanisms and gene regulatory networks underlying these events have been extensively studied. Two opposing TGF-β signaling pathways, Nodal and BMP, play key roles in all three germ layers to respectively pattern the sea urchin ventral and dorsal sides. In this chapter, I describe our current understanding of sea urchin dorsal-ventral patterning mechanisms. Additionally, differences in the patterning mechanisms observed in lecithotrophic sea urchins (nonfeeding larvae) and in cidaroid sea urchins are also discussed, along with evolutionary insights gained from comparative analyses.

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.

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

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

ATAC-Seq for Assaying Chromatin Accessibility Protocol Using Echinoderm Embryos

Cis-regulatory elements (CREs) and transcription factors (TFs) associated with them determine temporal and spatial domains of gene expression. Therefore, identification of these CREs and TFs is crucial to elucidating transcriptional programs across taxa. With chromatin accessibility facilitating transcription factor access to DNA, the identification of regions of open chromatin sheds light both on the function of the regulatory elements and their evolution, thus allowing the recognition of potential CREs. Buenrostro and colleagues have developed a novel method for exploring chromatin accessibility: assay for transposase-accessible chromatin with high-throughput sequencing (ATAC-seq), which can be used for the purpose of identifying putative CREs. This method was shown to have considerable advantages when compared to traditional methods such as sequence conservation analyses or functional assays. Here we present the adaptation of the ATAC-seq method to echinoderm species and discuss how it can be used for CRE discovery.

Maternal control of early patterning in sea urchin embryos

Sea urchin development has been studied extensively for more than a century and considered regulative since the first experimental evidence. Further investigations have repeatedly supported this standpoint by revealing the presence of inductive mechanisms that alter cell fate decisions at early cleavage stages and flexibility of development in response to environmental conditions. Some features indicate that sea urchin development is not completely regulative, but actually includes determinative events. In 16-cell embryos, mesomeres and macromeres represent multipotency, while the cell fate of most vegetal micromeres is restricted. It is known that the mature sea urchin eggs are polarized by the asymmetrical distribution of some maternal mRNAs and proteins. Spatially-distributed maternal factors are necessary for the orientation of the primary animal-vegetal axis, which is established by both maternal and zygotic mechanisms later in development. The secondary dorsal-ventral axis is conditionally specified later in development. Dorsal-ventral polarity is very liable during the early cleavages, though more recent data argue that its direction may be oriented by maternal asymmetry. In this review, we focus on the role of maternal factors in initial embryonic patterning during the first cleavages of sea urchin embryos before activation of the embryonic genome.

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

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

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.

Trans-Species Polymorphism in Mitochondrial Genome of Camarodont Sea Urchins

Mitochondrial (mt) genomes of the sea urchins Strongylocentrotus intermedius and Mesocentrotus nudus demonstrate the identical patterns of intraspecific length variability of the ND6 gene, consisting of 489 bp (S variant) and 498 bp (L variant), respectively. For both species, the ND6 length difference is due to the 488A>G substitution, which changes the stop codon TAG in S variant for a tryptophan codon TGG in L variant and elongates the corresponding ND6 protein by three additional amino acids, Trp-Leu-Trp. The phylogenetic analysis based on mt genomes of sea urchins and related echinoderm groups from GenBank has shown the S and L ND6 variants as shared among the camarodont sea urchins; the rest of the echinoderms demonstrate the S variant only. The data suggest that the ND6 488A>G substitution can be the first example of the trans-species polymorphism in sea urchins, persisting at least since the time of the Odontophora diversification at the Eocene/Oligocene boundary (approximately 34 million years ago), which was characterized by an abrupt climate change and significant global ocean cooling. Alternative hypotheses, including the convergent RNA editing and/or codon reassignment, are not supported by direct comparisons of the ND6 gene sequences with the corresponding transcripts using the basic local alignment search tool (BLAST) of full sea urchin transcriptomes.

The maternal-to-zygotic transition revisited

The development of animal embryos is initially directed by maternal gene products. Then, during the maternal-to-zygotic transition (MZT), developmental control is handed to the zygotic genome. Extensive research in both vertebrate and invertebrate model organisms has revealed that the MZT can be subdivided into two phases, during which very different modes of gene regulation are implemented: initially, regulation is exclusively post-transcriptional and post-translational, following which gradual activation of the zygotic genome leads to predominance of transcriptional regulation. These changes in the gene expression program of embryos are precisely controlled and highly interconnected. Here, we review current understanding of the mechanisms that underlie handover of developmental control during the MZT.

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.

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.

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

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

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

Across the animal kingdom, embryos of closely related species show high morphological similarity despite genetic and environmental distances. Deciphering the molecular mechanisms that underlie morphological conservation and those that support embryonic adaptation are keys to understand developmental robustness and evolution. Comparative studies of developmental gene regulatory networks can track the genetic changes that lead to evolutionary novelties. However, these studies are limited to a relatively small set of genes and demand extensive experimental efforts. An alternative approach enabled by next-generation sequencing, is to compare the expression kinetic of large sets of genes between different species. The advantages of these comparisons are that they can be done relatively easily, for any species and they provide information of all expressed genes. The challenge in these experiments is to compare the kinetic profiles of thousands of genes between species that develop in different rates. Here we review recent comparative studies that tackled the challenges of accurate staging and large-scale analyses using different computational approaches. These studies reveal how correct temporal scaling exposes the striking conservation of developmental gene expression between morphologically similar species. Different clustering approaches are used to address various comparative questions and identify the conservation and divergence of large gene sets. We discuss the unexpected contribution of housekeeping genes to the interspecies correlations and how this contribution distorts the hourglass pattern generated by developmental genes. Overall, we demonstrate how comparative studies of gene expression kinetics can provide novel insights into the developmental constraints and plasticity that shape animal body plans.

The developmental transcriptomes of two sea biscuit species with differing larval types

BACKGROUND

Larval developmental patterns are extremely varied both between and within phyla, however the genetic mechanisms leading to this diversification are poorly understood. We assembled and compared the developmental transcriptomes for two sea biscuit species (Echinodermata: Echinoidea) with differing patterns of larval development, to provide a resource for investigating the evolution of alternate life cycles. One species (Clypeaster subdepressus) develops via an obligately feeding larva which metamorphoses 3-4 weeks after fertilization; the other (Clypeaster rosaceus) develops via a rare, intermediate larval type-facultative feeding- and can develop through metamorphosis entirely based on egg provisioning in under one week.

RESULTS

Overall, the two transcriptomes are highly similar, containing largely orthologous contigs with similar functional annotation. However, we found distinct differences in gene expression patterns between the two species. Larvae from C. rosaceus, the facultative planktotroph, turned genes on at earlier stages and had less differentiation in gene expression between larval stages, whereas, C. subdepressus showed a higher degree of stage-specific gene expression.

CONCLUSION

This study is the first genetic analysis of a species with facultatively feeding larvae. Our results are consistent with known developmental differences between the larval types and raise the question of whether earlier onset of developmental genes is a key step in the evolution of a reduced larval period. By publishing a transcriptome for this rare, intermediate, larval type, this study adds developmental breadth to the current genetic resources, which will provide a valuable tool for future research on echinoderm development as well as studies on the evolution of development in general.

Translatome analysis at the egg-to-embryo transition in sea urchin

Early embryogenesis relies on the translational regulation of maternally stored mRNAs. In sea urchin, fertilization triggers a dramatic rise in translation activity, necessary for the onset of cell division. Here, the full spectrum of the mRNAs translated upon fertilization was investigated by polysome profiling and sequencing. The translatome of the early sea urchin embryo gave a complete picture of the polysomal recruitment dynamics following fertilization. Our results indicate that only a subset of maternal mRNAs were selectively recruited onto polysomes, with over-represented functional categories in the translated set. The increase in translation upon fertilization depends on the formation of translation initiation complexes following mTOR pathway activation. Surprisingly, mTOR pathway inhibition differentially affected polysomal recruitment of the newly translated mRNAs, which thus appeared either mTOR-dependent or mTOR-independent. Therefore, our data argue for an alternative to the classical cap-dependent model of translation in early development. The identification of the mRNAs translated following fertilization helped assign translational activation events to specific mRNAs. This translatome is the first step to a comprehensive analysis of the molecular mechanisms governing translation upon fertilization and the translational regulatory networks that control the egg-to-embryo transition as well as the early steps of embryogenesis.

Paleogenomics of echinoids reveals an ancient origin for the double-negative specification of micromeres in sea urchins

Establishing a timeline for the evolution of novelties is a common, unifying goal at the intersection of evolutionary and developmental biology. Analyses of gene regulatory networks (GRNs) provide the ability to understand the underlying genetic and developmental mechanisms responsible for the origin of morphological structures both in the development of an individual and across entire evolutionary lineages. Accurately dating GRN novelties, thereby establishing a timeline for GRN evolution, is necessary to answer questions about the rate at which GRNs and their subcircuits evolve, and to tie their evolution to paleoenvironmental and paleoecological changes. Paleogenomics unites the fossil record and all aspects of deep time, with modern genomics and developmental biology to understand the evolution of genomes in evolutionary time. Recent work on the regulatory genomic basis of development in cidaroid echinoids, sand dollars, heart urchins, and other nonmodel echinoderms provides an ideal dataset with which to explore GRN evolution in a comparative framework. Using divergence time estimation and ancestral state reconstructions, we have determined the age of the double-negative gate (DNG), the subcircuit which specifies micromeres and skeletogenic cells in Strongylocentrotus purpuratus We have determined that the DNG has likely been used for euechinoid echinoid micromere specification since at least the Late Triassic. The innovation of the DNG thus predates the burst of post-Paleozoic echinoid morphological diversification that began in the Early Jurassic. Paleogenomics has wide applicability for the integration of deep time and molecular developmental data, and has wide utility in rigorously establishing timelines for GRN evolution.

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.

Asymmetric distribution of hypoxia-inducible factor α regulates dorsoventral axis establishment in the early sea urchin embryo

Hypoxia signaling is an ancient pathway by which animals can respond to low oxygen. Malfunction of this pathway disturbs hypoxic acclimation and can result in various diseases, including cancers. The role of hypoxia signaling in early embryogenesis remains unclear. Here, we show that in the blastula of the sea urchin Strongylocentrotus purpuratus, hypoxia-inducible factor α (HIFα), the downstream transcription factor of the hypoxia pathway, is localized and transcriptionally active on the future dorsal side. This asymmetric distribution is attributable to its oxygen-sensing ability. Manipulations of the HIFα level entrained the dorsoventral axis, as the side with the higher level of HIFα tends to develop into the dorsal side. Gene expression analyses revealed that HIFα restricts the expression of nodal to the ventral side and activates several genes encoding transcription factors on the dorsal side. We also observed that intrinsic hypoxic signals in the early embryos formed a gradient, which was disrupted under hypoxic conditions. Our results reveal an unprecedented role of the hypoxia pathway in animal development.

Echinoderm development and evolution in the post-genomic era

The highly recognizable animals within the phylum Echinodermata encompass an enormous disparity of adult and larval body plans. The extensive knowledge of sea urchin development has culminated in the description of the exquisitely detailed gene regulatory network (GRN) that governs the specification of various embryonic territories. This information provides a unique opportunity for comparative studies in other echinoderm taxa to understand the evolution and developmental mechanisms underlying body plan change. This review focuses on recent work that has utilized new genomic resources and systems-level experiments to address questions of evolutionary developmental biology. In particular, we synthesize the results of several recent studies from various echinoderm classes that have explored the development and evolution of the larval skeleton, which is a major feature that distinguishes the two predominant larval subtypes within the Phylum. We specifically examine the ways in which GRNs can evolve, either through cis regulatory and/or protein-level changes in transcription factors. We also examine recent work comparing evolution across shorter time scales that occur within and between species of sea urchin, and highlight the kinds of questions that can be addressed by these comparisons. The advent of new genomic and transcriptomic datasets in additional species from all classes of echinoderm will continue to empower the use of these taxa for evolutionary developmental studies.

Divergence of ectodermal and mesodermal gene regulatory network linkages in early development of sea urchins

Developmental gene regulatory networks (GRNs) are assemblages of gene regulatory interactions that direct ontogeny of animal body plans. Studies of GRNs operating in the early development of euechinoid sea urchins have revealed that little appreciable change has occurred since their divergence ∼90 million years ago (mya). These observations suggest that strong conservation of GRN architecture was maintained in early development of the sea urchin lineage. Testing whether this holds for all sea urchins necessitates comparative analyses of echinoid taxa that diverged deeper in geological time. Recent studies highlighted extensive divergence of skeletogenic mesoderm specification in the sister clade of euechinoids, the cidaroids, suggesting that comparative analyses of cidaroid GRN architecture may confer a greater understanding of the evolutionary dynamics of developmental GRNs. Here I report spatiotemporal patterning of 55 regulatory genes and perturbation analyses of key regulatory genes involved in euechinoid oral-aboral patterning of nonskeletogenic mesodermal and ectodermal domains in early development of the cidaroid Eucidaris tribuloides These results indicate that developmental GRNs directing mesodermal and ectodermal specification have undergone marked alterations since the divergence of cidaroids and euechinoids. Notably, statistical and clustering analyses of echinoid temporal gene expression datasets indicate that regulation of mesodermal genes has diverged more markedly than regulation of ectodermal genes. Although research on indirect-developing euechinoid sea urchins suggests strong conservation of GRN circuitry during early embryogenesis, this study indicates that since the divergence of cidaroids and euechinoids, developmental GRNs have undergone significant, cell type-biased alterations.

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.

Comparative Developmental Transcriptomics Reveals Rewiring of a Highly Conserved Gene Regulatory Network during a Major Life History Switch in the Sea Urchin Genus Heliocidaris

The ecologically significant shift in developmental strategy from planktotrophic (feeding) to lecithotrophic (nonfeeding) development in the sea urchin genus Heliocidaris is one of the most comprehensively studied life history transitions in any animal. Although the evolution of lecithotrophy involved substantial changes to larval development and morphology, it is not known to what extent changes in gene expression underlie the developmental differences between species, nor do we understand how these changes evolved within the context of the well-defined gene regulatory network (GRN) underlying sea urchin development. To address these questions, we used RNA-seq to measure expression dynamics across development in three species: the lecithotroph Heliocidaris erythrogramma, the closely related planktotroph H. tuberculata, and an outgroup planktotroph Lytechinus variegatus. Using well-established statistical methods, we developed a novel framework for identifying, quantifying, and polarizing evolutionary changes in gene expression profiles across the transcriptome and within the GRN. We found that major changes in gene expression profiles were more numerous during the evolution of lecithotrophy than during the persistence of planktotrophy, and that genes with derived expression profiles in the lecithotroph displayed specific characteristics as a group that are consistent with the dramatically altered developmental program in this species. Compared to the transcriptome, changes in gene expression profiles within the GRN were even more pronounced in the lecithotroph. We found evidence for conservation and likely divergence of particular GRN regulatory interactions in the lecithotroph, as well as significant changes in the expression of genes with known roles in larval skeletogenesis. We further use coexpression analysis to identify genes of unknown function that may contribute to both conserved and derived developmental traits between species. Collectively, our results indicate that distinct evolutionary processes operate on gene expression during periods of life history conservation and periods of life history divergence, and that this contrast is even more pronounced within the GRN than across the transcriptome as a whole.

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.