Developmental gene regulatory networks in sea urchins and what we can learn from them

Stable germline transgenesis using the Minos Tc1/mariner element in the sea urchin Lytechinus pictus

Stable transgenesis is a transformative tool in model organism biology. Although the sea urchin is one of the oldest animal models in cell and developmental biology, studies in this animal have largely relied on transient manipulation of wild animals, without a strategy for stable transgenesis. Here, we build on recent progress to develop a more genetically tractable sea urchin species, Lytechinus pictus, and establish a robust transgene integration method. Three commonly used transposons (Minos, Tol2 and piggyBac) were tested for non-autonomous transposition, using plasmids containing a polyubiquitin promoter upstream of a H2B-mCerulean nuclear marker. Minos was the only transposable element that resulted in significant expression beyond metamorphosis. F0 animals were raised to sexual maturity, and spawned to determine germline integration and transgene inheritance frequency, and to characterize expression patterns of the transgene in F1 progeny. The results demonstrate transgene transmission through the germline, the first example of a germline transgenic sea urchin and, indeed, of any echinoderm. This milestone paves the way for the generation of diverse transgenic resources that will dramatically enhance the utility, reproducibility and efficiency of sea urchin research.

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

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

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

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

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

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

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

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

Comparative Analysis of Bivalve and Sea Urchin Genetics and Development: Investigating the Dichotomy in Bilateria

This comprehensive review presents a comparative analysis of early embryogenesis in Protostomia and Deuterostomia, the first of which exhibit a mosaic pattern of development, where cells are fated deterministically, while Deuterostomia display a regulatory pattern of development, where the fate of cells is indeterminate. Despite these fundamental differences, there are common transcriptional mechanisms that underline their evolutionary linkages, particularly in the field of functional genomics. By elucidating both conserved and unique regulatory strategies, this review provides essential insights into the comparative embryology and developmental dynamics of these groups. The objective of this review is to clarify the shared and distinctive characteristics of transcriptional regulatory mechanisms. This will contribute to the extensive areas of functional genomics, evolutionary biology and developmental biology, and possibly lay the foundation for future research and discussion on this seminal topic.

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.

Early expression onset of tissue-specific effector genes during the specification process in sea urchin embryos

BACKGROUND

In the course of animal developmental processes, various tissues are differentiated through complex interactions within the gene regulatory network. As a general concept, differentiation has been considered to be the endpoint of specification processes. Previous works followed this view and provided a genetic control scheme of differentiation in sea urchin embryos: early specification genes generate distinct regulatory territories in an embryo to express a small set of differentiation driver genes; these genes eventually stimulate the expression of tissue-specific effector genes, which provide biological identity to differentiated cells, in each region. However, some tissue-specific effector genes begin to be expressed in parallel with the expression onset of early specification genes, raising questions about the simplistic regulatory scheme of tissue-specific effector gene expression and the current concept of differentiation itself.

RESULTS

Here, we examined the dynamics of effector gene expression patterns during sea urchin embryogenesis. Our transcriptome-based analysis indicated that many tissue-specific effector genes begin to be expressed and accumulated along with the advancing specification GRN in the distinct cell lineages of embryos. Moreover, we found that the expression of some of the tissue-specific effector genes commences before cell lineage segregation occurs.

CONCLUSIONS

Based on this finding, we propose that the expression onset of tissue-specific effector genes is controlled more dynamically than suggested in the previously proposed simplistic regulation scheme. Thus, we suggest that differentiation should be conceptualized as a seamless process of accumulation of effector expression along with the advancing specification GRN. This pattern of effector gene expression may have interesting implications for the evolution of novel cell types.

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

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

Pattern of Repetitive Element Transcription Segregate Cell Lineages during the Embryogenesis of Sea Urchin Strongylocentrotus purpuratus

Repetitive elements (REs) occupy a significant part of eukaryotic genomes and are shown to play diverse roles in genome regulation. During embryogenesis of the sea urchin, a large number of REs are expressed, but the role of these elements in the regulation of biological processes remains unknown. The aim of this study was to identify the RE expression at different stages of embryogenesis. REs occupied 44% of genomic DNA of Strongylocentrotus purpuratus. The most prevalent among these elements were the unknown elements-in total, they contributed 78.5% of REs (35% in total genome occupancy). It was revealed that the transcription pattern of genes and REs changes significantly during gastrulation. Using the de novo transcriptome assembly, we showed that the expression of RE is independent of its copy number in the genome. We also identified copies that are expressed. Only active RE copies were used for mapping and quantification of RE expression in the single-cell RNA sequencing data. REs expression was observed in all cell lineages and they were detected as population markers. Moreover, the primary mesenchyme cell (PMC) line had the greatest diversity of REs among the markers. Our data suggest a role for RE in the organization of developmental domains during the sea urchin embryogenesis at the single-cell resolution level.

A nomenclature for echinoderm genes

Echinoderm embryos and larvae are prominent experimental model systems for studying developmental mechanisms. High-quality, assembled, annotated genome sequences are now available for several echinoderm species, including representatives from most classes. The increased availability of these data necessitates the development of a nomenclature that assigns universally interpretable gene symbols to echinoderm genes to facilitate cross-species comparisons of gene functions, both within echinoderms and across other phyla. This paper describes the implementation of an improved set of echinoderm gene nomenclature guidelines that both communicates meaningful orthology information in protein-coding gene symbols and names and establishes continuity with nomenclatures developed for major vertebrate model organisms, including humans. Differences between the echinoderm gene nomenclature guidelines and vertebrate guidelines are examined and explained. This nomenclature incorporates novel solutions to allow for several types of orthologous relationships, including the single echinoderm genes with multiple vertebrate co-orthologs that result from whole-genome-duplication events. The current version of the Echinoderm Gene Nomenclature Guidelines can be found at https://www.echinobase.org/gene/static/geneNomenclature.jsp

Database URL https://www.echinobase.org/

Global patterns of enhancer activity during sea urchin embryogenesis assessed by eRNA profiling

We used capped analysis of gene expression with sequencing (CAGE-seq) to profile eRNA expression and enhancer activity during embryogenesis of a model echinoderm: the sea urchin, Strongylocentrotus purpuratus. We identified more than 18,000 enhancers that were active in mature oocytes and developing embryos and documented a burst of enhancer activation during cleavage and early blastula stages. We found that a large fraction (73.8%) of all enhancers active during the first 48 h of embryogenesis were hyperaccessible no later than the 128-cell stage and possibly even earlier. Most enhancers were located near gene bodies, and temporal patterns of eRNA expression tended to parallel those of nearby genes. Furthermore, enhancers near lineage-specific genes contained signatures of inputs from developmental gene regulatory networks deployed in those lineages. A large fraction (60%) of sea urchin enhancers previously shown to be active in transgenic reporter assays was associated with eRNA expression. Moreover, a large fraction (50%) of a representative subset of enhancers identified by eRNA profiling drove tissue-specific gene expression in isolation when tested by reporter assays. Our findings provide an atlas of developmental enhancers in a model sea urchin and support the utility of eRNA profiling as a tool for enhancer discovery and regulatory biology. The data generated in this study are available at Echinobase, the public database of information related to echinoderm genomics.

Sea Urchin as a Universal Model for Studies of Gene Networks

The purple sea urchin Strongylocentrotus purpuratus has been used for over 150 years as a model organism in developmental biology. Using this model species, scientists have been able to describe, in detail, the mechanisms of cell cycle control and cell adhesion, fertilization, calcium signaling, cell differentiation, and death. Massive parallel sequencing of the sea urchin genome enabled the deciphering of the main components of gene regulatory networks during the activation of embryonic signaling pathways. This knowledge helped to extrapolate aberrations in somatic cells that may lead to diseases, including cancer in humans. Furthermore, since many, if not all, developmental signaling pathways were shown to be controlled by non-coding RNAs (ncRNAs), the sea urchin organism represents an attractive experimental model. In this review, we discuss the main discoveries in the genetics, genomics, and transcriptomics of sea urchins during embryogenesis with the main focus on the role of ncRNAs. This information may be useful for comparative studies between different organisms, and may help identify new regulatory networks controlled by ncRNAs.

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.

Embryo microinjection of the lecithotrophic sea urchin Heliocidaris erythrogramma

Microinjection is a common embryological technique used for many types of experiments, including lineage tracing, manipulating gene expression, or genome editing. Injectable reagents include mRNA overexpression, mis-expression, or dominant-negative experiments to examine a gene of interest, a morpholino antisense oligo to prevent translation of an mRNA or spliceoform of interest and CRISPR-Cas9 reagents. Thus, the technique is broadly useful for basic embryological studies, constructing gene regulatory networks, and directly testing hypotheses about cis-regulatory and coding sequence changes underlying the evolution of development. However, the methods for microinjection in typical planktotrophic marine invertebrates may not work well in the highly modified eggs and embryos of lecithotrophic species. This protocol is optimized for the lecithotrophic sea urchin Heliocidaris erythrogramma.

Gene expression profiling suggests differences in molecular mechanisms of fin elongation between cichlid species

Comparative analyses of gene regulation inform about the molecular basis of phenotypic trait evolution. Here, we address a fin shape phenotype that evolved multiple times independently across teleost fish, including several species within the family Cichlidae. In a previous study, we proposed a gene regulatory network (GRN) involved in the formation and regeneration of conspicuous filamentous elongations adorning the unpaired fins of the Neolamprologus brichardi. Here, we tested the members of this network in the blockhead cichlid, Steatocranus casuarius, which displays conspicuously elongated dorsal and moderately elongated anal fins. Our study provided evidence for differences in the anatomy of fin elongation and suggested gene regulatory divergence between the two cichlid species. Only a subset of the 20 genes tested in S. casuarius showed the qPCR expression patterns predicted from the GRN identified in N. brichardi, and several of the gene-by-gene expression correlations differed between the two cichlid species. In comparison to N. brichardi, gene expression patterns in S. casuarius were in better (but not full) agreement with gene regulatory interactions inferred in zebrafish. Within S. casuarius, the dorsoventral asymmetry in ornament expression was accompanied by differences in gene expression patterns, including potential regulatory differentiation, between the anal and dorsal fin.

Encounters across networks: Windows into principles of genomic regulation

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

From genome to anatomy: The architecture and evolution of the skeletogenic gene regulatory network of sea urchins and other echinoderms

The skeletogenic gene regulatory network (GRN) of sea urchins and other echinoderms is one of the most intensively studied transcriptional networks in any developing organism. As such, it serves as a preeminent model of GRN architecture and evolution. This review summarizes our current understanding of this developmental network. We describe in detail the most comprehensive model of the skeletogenic GRN, one developed for the euechinoid sea urchin Strongylocentrotus purpuratus, including its initial deployment by maternal inputs, its elaboration and stabilization through regulatory gene interactions, and its control of downstream effector genes that directly drive skeletal morphogenesis. We highlight recent comparative studies that have leveraged the euechinoid GRN model to examine the evolution of skeletogenic programs in diverse echinoderms, studies that have revealed both conserved and divergent features of skeletogenesis within the phylum. Last, we summarize the major insights that have emerged from analysis of the structure and evolution of the echinoderm skeletogenic GRN and identify key, unresolved questions as a guide for future work.

Perturbation of Developmental Regulatory Gene Expression by a G-Quadruplex DNA Inducer in the Sea Urchin Embryo

The G-quadruplex (G4) is a four-stranded DNA structure identified in vivo in guanine-rich regions located in the promoter of a number of genes. Intriguing evidence suggested that small molecules acting as G4-targeting ligands could potentially regulate multiple cellular processes via either stabilizing or disruptive effects on G4 motifs. Research in this field aims to prove the direct role of G4 ligands and/or structures on a specific biological process in a complex living organism. In this study, we evaluate in vivo the effects of a nickel(II)-salnaphen-like complex, named Nisaln, a potent G4 binder and stabilizer, during embryogenesis of the sea urchin embryo. We describe developmental defects inflicted by Nisaln and correlate them with variation in the expression of several regulatory genes. It is worth mentioning that we show that Nisaln binds a G4 structure in the promoter of hbox12-a, a gene lying at the top of the developmental regulatory hierarchy, inducing overexpression of this gene.

Cdc42 controls primary mesenchyme cell morphogenesis in the sea urchin embryo

In the sea urchin embryo, gastrulation is characterized by the ingression and directed cell migration of primary mesenchyme cells (PMCs), as well as the primary invagination and convergent extension of the endomesoderm. Like all cell shape changes, individual and collective cell motility is orchestrated by Rho family GTPases and their modulation of the actomyosin cytoskeleton. And while endomesoderm specification has been intensively studied in echinoids, much less is known about the proximate regulators driving cell motility. Toward these ends, we employed anti-sense morpholinos, mutant alleles and pharmacological inhibitors to assess the role of Cdc42 during sea urchin gastrulation. While inhibition of Cdc42 expression or activity had only mild effects on PMC ingression, PMC migration, alignment and skeletogenesis were disrupted in the absence of Cdc42, as well as elongation of the archenteron. PMC migration and patterning of the larval skeleton relies on the extension of filopodia, and Cdc42 was required for filopodia in vivo as well as in cultured PMCs. Lastly, filopodial extension required both Arp2/3 and formin actin-nucleating factors, supporting models of filopodial nucleation observed in other systems. Together, these results suggest that Cdc42 plays essential roles during PMC cell motility and organogenesis.

Two genetic changes in cis-regulatory elements caused evolution of petal spot position in Clarkia

A major premise in evolutionary developmental biology is that regulatory changes, often involving cis-regulatory elements, are responsible for much morphological evolution. This premise is supported by recent investigations of animal development, but information is just beginning to accumulate regarding whether it also applies to the evolution of plant morphology^1-4. Here, we identify the genetic differences between species in the genus Clarkia that are responsible for evolutionary change in an ecologically important element of floral colour patterns: spot position. The evolutionary shift in spot position was due to two simple genetic changes that resulted in the appearance of a transcription factor binding site mutation in the R2R3 Myb gene that changes spot formation. These genetic changes caused R2R3 Myb to be activated by a different transcription factor that is expressed in a different position in the petal. These results suggest that the regulatory rewiring paradigm is as applicable to plants as it is to animals, and support the hypothesis that cis-regulatory changes may often play a role in plant morphological evolution.

New biomarkers of post-settlement growth in the sea urchin Strongylocentrotus purpuratus

Some sea urchins, including the purple sea urchin Strongylocentrotus purpuratus, have been successfully used in aquaculture, but their slow growth and late reproduction are challenging to overcome when developing efficient aquaculture production techniques. S. purpuratus develops via an indirect life history that is characterized by a drastic settlement process at the end of a larval period that lasts for several weeks. During this transition, the bilateral larva is transformed into a pentaradial juvenile, which will start feeding and growing in the benthic habitat. Due to predation and other ecological factors, settlement is typically associated with high mortality rates in juvenile populations. Additionally, juveniles require several days to develop a functional mouth and digestive system. During this perimetamorphic period, juveniles use up larval resources until they are capable to digest adult food. Mechanisms underlying the onset of juvenile feeding and metabolism have implications for the recruitment of natural populations as well as aquaculture and are relatively poorly understood in S. purpuratus. The insulin/insulin-like growth factor signalling (IIS)/Target of Rapamycin (TOR) pathway (IIS/TOR) is well conserved among animal phyla and regulates physiological and developmental functions, such as growth, reproduction, aging and nutritional status. We analyzed the expression of FoxO, TOR, and ILPs in post-settlement juveniles in conjunction with their early growth trajectories. We also tested how pre-settlement starvation affected post-settlement expression of IIS. We found that FoxO provides a useful molecular marker in early juveniles as its expression is strongly correlated with juvenile growth. We also found that pre-settlement starvation affects juvenile growth trajectories as well as IIS. Our findings provide preliminary insights into the mechanisms underlying post-settlement growth and metabolism in S. purpuratus. They also have important implications for sea urchin aquaculture, as they show that pre-settlement nutrient environment significantly affects both early growth trajectories and gene expression. This information can be used to develop new biomarkers for juvenile health in sea urchin population ecology and aquaculture aquaculture.

Diversification of spatiotemporal expression and copy number variation of the echinoid hbox12/pmar1/micro1 multigene family

Changes occurring during evolution in the cis-regulatory landscapes of individual members of multigene families might impart diversification in their spatiotemporal expression and function. The archetypal member of the echinoid hbox12/pmar1/micro1 family is hbox12-a, a homeobox-containing gene expressed exclusively by dorsal blastomeres, where it governs the dorsal/ventral gene regulatory network during embryogenesis of the sea urchin Paracentrotus lividus. Here we describe the inventory of the hbox12/pmar1/micro1 genes in P. lividus, highlighting that gene copy number variation occurs across individual sea urchins of the same species. We show that the various hbox12/pmar1/micro1 genes group into three subfamilies according to their spatiotemporal expression, which ranges from broad transcription throughout development to transient expression in either the animal hemisphere or micromeres of the early embryo. Interestingly, the promoter regions of those genes showing comparable expression patterns are highly similar, while differing from those of the other subfamilies. Strikingly, phylogenetic analysis suggests that the hbox12/pmar1/micro1 genes are species-specific, exhibiting extensive divergence in their noncoding, but not in their coding, sequences across three distinct sea urchin species. In spite of this, two micromere-specific genes of P. lividus possess a TCF/LEF-binding motif in a similar position, and their transcription relies on Wnt/β-catenin signaling, similar to the pmar1 and micro1 genes, which in other sea urchin species are involved in micromere specification. Altogether, our findings suggest that the hbox12/pmar1/micro1 gene family evolved rather rapidly, generating paralogs whose cis-regulatory sequences diverged following multiple rounds of duplication from a common ancestor.

A gene-centered C. elegans protein-DNA interaction network provides a framework for functional predictions

Transcription factors (TFs) play a central role in controlling spatiotemporal gene expression and the response to environmental cues. A comprehensive understanding of gene regulation requires integrating physical protein–DNA interactions (PDIs) with TF regulatory activity, expression patterns, and phenotypic data. Although great progress has been made in mapping PDIs using chromatin immunoprecipitation, these studies have only characterized ~10% of TFs in any metazoan species. The nematode C. elegans has been widely used to study gene regulation due to its compact genome with short regulatory sequences. Here, we delineated the largest gene‐centered metazoan PDI network to date by examining interactions between 90% of C. elegans TFs and 15% of gene promoters. We used this network as a backbone to predict TF binding sites for 77 TFs, two‐thirds of which are novel, as well as integrate gene expression, protein–protein interaction, and phenotypic data to predict regulatory and biological functions for multiple genes and TFs.