Conservation of cis-Regulatory Syntax Underlying Deuterostome Gastrulation.

Throughout embryonic development, the shaping of the functional and morphological characteristics of embryos is orchestrated by an intricate interaction between transcription factors and cis-regulatory elements. In this study, we conducted a comprehensive analysis of deuterostome cis-regulatory landscapes during gastrulation, focusing on four paradigmatic species: the echinoderm Strongylocentrotus purpuratus, the cephalochordate Branchiostoma lanceolatum, the urochordate Ciona intestinalis, and the vertebrate Danio rerio. Our approach involved comparative computational analysis of ATAC-seq datasets to explore the genome-wide blueprint of conserved transcription factor binding motifs underlying gastrulation. We identified a core set of conserved DNA binding motifs associated with 62 known transcription factors, indicating the remarkable conservation of the gastrulation regulatory landscape across deuterostomes. Our findings offer valuable insights into the evolutionary molecular dynamics of embryonic development, shedding light on conserved regulatory subprograms and providing a comprehensive perspective on the conservation and divergence of gene regulation underlying the gastrulation process.

From molecules to morphology: How food supply influences the larvae of sea urchins across all levels of biological organization.

An important goal of many studies in molecular ecology is to utilize molecular tools to elucidate how critical traits like metabolism and growth are affected by environmental stressors and how organisms offset these stresses by adaptive molecular-level responses. Stress from food deprivation may be critical for early developmental stages that require a continued supply of substrates for energy metabolism and growth if development is to be completed. In a 'From the Cover' article in this issue of Molecular Ecology, Li et al. (2023) examined the effects of withholding food (unicellular algae) on 10 traits of larvae of the purple sea urchin (Strongylocentrotus purpuratus), ranging from the molecular level (gene expression) to morphology. Overall, this study sheds new light on the plasticity of larval development and the tight linkages that exist among traits as they respond to changes in food availability. Importantly, shifts in the sources of food utilized under different dietary treatments show the plasticity of these larvae to alter reliance on endogenous energy stores and dissolved organic matter (DOM) as algae deprivation continues. The effects of global change on the amounts and phenology of productivity in the seas make this type of integrated, multi-level analysis an important tool for predicting the future states of marine ecosystems.

Recombinant SpTransformer proteins are functionally diverse for binding and phagocytosis by three subtypes of sea urchin phagocytes.

Introduction: The California purple sea urchin, Strongylocentrotus purpuratus, relies solely on an innate immune system to combat the many pathogens in the marine environment. One aspect of their molecular defenses is the SpTransformer (SpTrf) gene family that is upregulated in response to immune challenge. The gene sequences are highly variable both within and among animals and likely encode thousands of SpTrf isoforms within the sea urchin population. The native SpTrf proteins bind foreign targets and augment phagocytosis of a marine Vibrio. A recombinant (r)SpTrf-E1-Ec protein produced by E. coli also binds Vibrio but does not augment phagocytosis. Methods: To address the question of whether other rSpTrf isoforms function as opsonins and augment phagocytosis, six rSpTrf proteins were expressed in insect cells. Results: The rSpTrf proteins are larger than expected, are glycosylated, and one dimerized irreversibly. Each rSpTrf protein cross-linked to inert magnetic beads (rSpTrf::beads) results in different levels of surface binding and phagocytosis by phagocytes. Initial analysis shows that significantly more rSpTrf::beads associate with cells compared to control BSA::beads. Binding specificity was verified by pre-incubating the rSpTrf::beads with antibodies, which reduces the association with phagocytes. The different rSpTrf::beads show significant differences for cell surface binding and phagocytosis by phagocytes. Furthermore, there are differences among the three distinct types of phagocytes that show specific vs. constitutive binding and phagocytosis. Conclusion: These findings illustrate the complexity and effectiveness of the sea urchin innate immune system driven by the natSpTrf proteins and the phagocyte cell populations that act to neutralize a wide range of foreign pathogens.

Local Genomic Instability of the SpTransformer Gene Family in the Purple Sea Urchin Inferred from BAC Insert Deletions.

The SpTransformer (SpTrf) gene family in the purple sea urchin, Strongylocentrotus purpuratus, encodes immune response proteins. The genes are clustered, surrounded by short tandem repeats, and some are present in genomic segmental duplications. The genes share regions of sequence and include repeats in the coding exon. This complex structure is consistent with putative local genomic instability. Instability of the SpTrf gene cluster was tested by 10 days of growth of Escherichia coli harboring bacterial artificial chromosome (BAC) clones of sea urchin genomic DNA with inserts containing SpTrf genes. After the growth period, the BAC DNA inserts were analyzed for size and SpTrf gene content. Clones with multiple SpTrf genes showed a variety of deletions, including loss of one, most, or all genes from the cluster. Alternatively, a BAC insert with a single SpTrf gene was stable. BAC insert instability is consistent with variations in the gene family composition among sea urchins, the types of SpTrf genes in the family, and a reduction in the gene copy number in single coelomocytes. Based on the sequence variability among SpTrf genes within and among sea urchins, local genomic instability of the family may be important for driving sequence diversity in this gene family that would be of benefit to sea urchins in their arms race with marine microbes.

Spotting disease disrupts the microbiome of infected purple sea urchins, Strongylocentrotus purpuratus.

BACKGROUND: Spotting disease infects a variety of sea urchin species across many different marine locations. The disease is characterized by discrete lesions on the body surface composed of discolored necrotic tissue that cause the loss of all surface appendages within the lesioned area. A similar, but separate disease of sea urchins called bald sea urchin disease (BSUD) has overlapping symptoms with spotting disease, resulting in confusions in distinguishing the two diseases. Previous studies have focus on identifying the underlying causative agent of spotting disease, which has resulted in the identification of a wide array of pathogenic bacteria that vary based on location and sea urchin species. Our aim was to investigate the spotting disease infection by characterizing the microbiomes of the animal surface and various tissues. RESULTS: We collected samples of the global body surface, the lesion surface, lesioned and non-lesioned body wall, and coelomic fluid, in addition to samples from healthy sea urchins. 16S rRNA gene was amplified and sequenced from the genomic DNA. Results show that the lesions are composed mainly of Cyclobacteriaceae, Cryomorphaceae, and a few other taxa, and that the microbial composition of lesions is the same for all infected sea urchins. Spotting disease also alters the microbial composition of the non-lesioned body wall and coelomic fluid of infected sea urchins. In our closed aquarium systems, sea urchins contracted spotting disease and BSUD separately and therefore direct comparisons could be made between the microbiomes from diseased and healthy sea urchins. CONCLUSION: Results show that spotting disease and BSUD are separate diseases with distinct symptoms and distinct microbial compositions.

Whole-Genome Sequencing Reveals That Regulatory and Low Pleiotropy Variants Underlie Local Adaptation to Environmental Variability in Purple Sea Urchins.

AbstractOrganisms experience environments that vary across both space and time. Such environmental heterogeneity shapes standing genetic variation and may influence species' capacity to adapt to rapid environmental change. However, we know little about the kind of genetic variation that is involved in local adaptation to environmental variability. To address this gap, we sequenced the whole genomes of 140 purple sea urchins (Strongylocentrotus purpuratus) from seven populations that vary in their degree of pH variability. Despite no evidence of global population structure, we found a suite of single-nucleotide polymorphisms (SNPs) tightly correlated with local pH variability (outlier SNPs), which were overrepresented in regions putatively involved in gene regulation (long noncoding RNA and enhancers), supporting the idea that variation in regulatory regions is important for local adaptation to variability. In addition, outliers in genes were found to be (i) enriched for biomineralization and ion homeostasis functions related to low pH response, (ii) less central to the protein-protein interaction network, and (iii) underrepresented among genes highly expressed during early development. Taken together, these results suggest that loci that underlie local adaptation to pH variability in purple sea urchins fall in regions with potentially low pleiotropic effects (based on analyses involving regulatory regions, network centrality, and expression time) involved in low pH response (based on functional enrichment).

Widespread priming of transcriptional regulatory elements by incipient accessibility or RNA polymerase II pause in early embryos of the sea urchin Strongylocentrotus purpuratus.

Transcriptional regulatory elements (TREs) are the primary nodes that control developmental gene regulatory networks. In embryo stages, larvae, and adult differentiated red spherule cells of the sea urchin Strongylocentrotus purpuratus, transcriptionally engaged TREs are detected by Precision Run-On Sequencing (PRO-seq), which maps genome-wide at base pair resolution the location of paused or elongating RNA polymerase II (Pol II). In parallel, TRE accessibility is estimated by the Assay for Transposase-Accessible Chromatin using Sequencing (ATAC-seq). Our analysis identifies surprisingly early and widespread TRE accessibility in 4-cell cleavage embryos that is not necessarily followed by concurrent or subsequent transcription. TRE transcriptional differences identified by PRO-seq provide more contrast among embryonic stages than ATAC-seq accessibility differences, in agreement with the apparent excess of accessible but inactive TREs during embryogenesis. Global TRE accessibility reaches a maximum around the 20-hour late blastula stage, which coincides with the consolidation of major embryo regionalizations and peak histone variant H2A.Z expression. A transcriptional potency model based on labile nucleosome TRE occupancy driven by DNA sequences and the prevalence of histone variants is proposed in order to explain the basal accessibility of transcriptionally inactive TREs during embryogenesis. However, our results would not reconcile well with labile nucleosome models based on simple A/T sequence enrichment. In addition, a large number of distal TREs become transcriptionally disengaged during developmental progression, in support of an early Pol II paused model for developmental gene regulation that eventually resolves in transcriptional activation or silencing. Thus, developmental potency in early embryos may be facilitated by incipient accessibility and transcriptional pause at TREs.

Thyroid hormone membrane receptor binding and transcriptional regulation in the sea urchin Strongylocentrotus purpuratus.

Thyroid hormones (THs) are small amino acid derived signaling molecules with broad physiological and developmental functions in animals. Specifically, their function in metamorphic development, ion regulation, angiogenesis and many others have been studied in detail in mammals and some other vertebrates. Despite extensive reports showing pharmacological responses of invertebrate species to THs, little is known about TH signaling mechanisms outside of vertebrates. Previous work in sea urchins suggests that non-genomic mechanisms are activated by TH ligands. Here we show that several THs bind to sea urchin (Strongylocentrotus purpuratus) cell membrane extracts and are displaced by ligands of RGD-binding integrins. A transcriptional analysis across sea urchin developmental stages shows activation of genomic and non-genomic pathways in response to TH exposure, suggesting that both pathways are activated by THs in sea urchin embryos and larvae. We also provide evidence associating TH regulation of gene expression with TH response elements in the genome. In ontogeny, we found more differentially expressed genes in older larvae compared to gastrula stages. In contrast to gastrula stages, the acceleration of skeletogenesis by thyroxine in older larvae is not fully inhibited by competitive ligands or inhibitors of the integrin membrane receptor pathway, suggesting that THs likely activate multiple pathways. Our data confirms a signaling function of THs in sea urchin development and suggests that both genomic and non-genomic mechanisms play a role, with genomic signaling being more prominent during later stages of larval development.

A case of hermaphroditism in the gonochoristic sea urchin, Strongylocentrotus purpuratus, reveals key mechanisms of sex determination†.

Sea urchins are usually gonochoristic, with all of their five gonads either testes or ovaries. Here, we report an unusual case of hermaphroditism in the purple sea urchin, Strongylocentrotus purpuratus. The hermaphrodite is self-fertile, and one of the gonads is an ovotestis; it is largely an ovary with a small segment containing fully mature sperm. Molecular analysis demonstrated that each gonad producedviable gametes, and we identified for the first time a somatic sex-specific marker in this phylum: Doublesex and mab-3 related transcription factor 1 (DMRT1). This finding also enabled us to analyze the somatic tissues of the hermaphrodite, and we found that the oral tissues (including gut) were out of register with the aboral tissues (including tube feet) enabling a genetic lineage analysis. Results from this study support a genetic basis of sex determination in sea urchins, the viability of hermaphroditism, and distinguish gonad determination from somatic tissue organization in the adult.

Sunflower sea star predation on urchins can facilitate kelp forest recovery.

The recent collapse of predatory sunflower sea stars (Pycnopodia helianthoides) owing to sea star wasting disease (SSWD) is hypothesized to have contributed to proliferation of sea urchin barrens and losses of kelp forests on the North American west coast. We used experiments and a model to test whether restored Pycnopodia populations may help recover kelp forests through their consumption of nutritionally poor purple sea urchins (Strongylocentrotus purpuratus) typical of barrens. Pycnopodia consumed 0.68 S. purpuratus d-1, and our model and sensitivity analysis shows that the magnitude of recent Pycnopodia declines is consistent with urchin proliferation after modest sea urchin recruitment, and even small Pycnopodia recoveries could generally lead to lower densities of sea urchins that are consistent with kelp-urchin coexistence. Pycnopodia seem unable to chemically distinguish starved from fed urchins and indeed have higher predation rates on starved urchins owing to shorter handling times. These results highlight the importance of Pycnopodia in regulating purple sea urchin populations and maintaining healthy kelp forests through top-down control. The recovery of this important predator to densities commonly found prior to SSWD, whether through natural means or human-assisted reintroductions, may therefore be a key step in kelp forest restoration at ecologically significant scales.

Molecular and Cellular Characterization of the TH Pathway in the Sea Urchin Strongylocentrotus purpuratus.

Thyroid Hormones (THs) are a class of signaling molecules produced by coupling iodine with tyrosine residues. In vertebrates, extensive data support their important role in a variety of processes such as metabolism, development and metamorphosis. On the other hand, in invertebrates, the synthesis and role of the THs have been, so far, poorly investigated, thus limiting our understanding of the function and evolution of this important animal signaling pathway. In sea urchins, for example, while several studies focused on the availability and function of external sources of iodotyrosines, preliminary evidence suggests that an endogenous TH pathway might be in place. Here, integrating available literature with an in silico analysis, various homologous genes of the vertebrate TH molecular toolkit have been identified in the sea urchin Strongylocentrotus purpuratus. They include genes involved in the synthesis (Sp-Pxdn), metabolism (Sp-Dios), transport (Sp-Ttrl, Sp-Mct7/8/10) and response (Sp-Thr, Sp-Rxr and Sp-Integrin αP) to thyroid hormones. To understand the cell type(s) involved in TH synthesis and/or response, we studied the spatial expression of the TH toolkit during urchin development. Exploiting single-cell transcriptomics data in conjunction with in situ hybridization and immunohistochemistry, we identified cell types that are potentially producing or responding to THs in the sea urchin. Finally, growing sea urchin embryos until the larva stage with and without a source of inorganic iodine, we provided evidence that iodine organification is important for larval skeleton growth.

Identification of the genes encoding candidate septate junction components expressed during early development of the sea urchin, Strongylocentrotus purpuratus, and evidence of a role for Mesh in the formation of the gut barrier.

Septate junctions (SJs) evolved as cell-cell junctions that regulate the paracellular barrier and integrity of epithelia in invertebrates. Multiple morphological variants of SJs exist specific to different epithelia and/or phyla but the biological significance of varied SJ morphology is unclear because the knowledge of the SJ associated proteins and their functions in non-insect invertebrates remains largely unknown. Here we report cell-specific expression of nine candidate SJ genes in the early life stages of the sea urchin Strongylocentrotus purpuratus. By use of in situ RNA hybridization and single cell RNA-seq we found that the expression of selected genes encoding putatively SJ associated transmembrane and cytoplasmic scaffold molecules was dynamically regulated during epithelial development in the embryos and larvae with different epithelia expressing different cohorts of SJ genes. We focused a functional analysis on SpMesh, a homolog of the Drosophila smooth SJ component Mesh, which was highly enriched in the endodermal epithelium of the mid- and hindgut. Functional perturbation of SpMesh by both CRISPR/Cas9 mutagenesis and vivo morpholino-mediated knockdown shows that loss of SpMesh does not disrupt the formation of the gut epithelium during gastrulation. However, loss of SpMesh resulted in a severely reduced gut-paracellular barrier as quantitated by increased permeability to 3-5 â€‹kDa FITC-dextran. Together, these studies provide a first look at the molecular SJ physiology during the development of a marine organism and suggest a shared role for Mesh-homologous proteins in forming an intestinal barrier in invertebrates. Results have implications for consideration of the traits underlying species-specific sensitivity of marine larvae to climate driven ocean change.

Functional characterization of the MyD88 homologs in Strongylocentrotus purpuratus.

Toll-like receptor signaling is an evolutionarily conserved pathway to induce the expression of immune mediators in response to encounters with pathogens. MyD88 is a central adapter connecting the intracellular domain of the receptors to downstream kinases. Here, we conducted a comprehensive assessment of the MyD88 family in an echinoderm, Strongylocentrotus purpuratus. Of five SpMyD88s only two closely related proteins, SpMyD88A and SpMyD88B, are functional in mammalian cell lines as their overexpression facilitates the activation of the downstream transcription factor NF-κB. This requires the presence of the endogenous mammalian MyD88s, and domain swapping indicated that the death domains of S. purpuratus MyD88 are unable to efficiently connect to the respective domains of the vertebrate IRAK kinases. This suggests that the interaction surfaces between the signaling mediators in this conserved signaling pathway are not as conserved as previously thought but were likely shaped and evolved by pathogenic selection over evolutionary timescales.

Coelomocyte populations in the sea urchin, Strongylocentrotus purpuratus, undergo dynamic changes in response to immune challenge.

The sea urchin, Strongylocentrotus purpuratus has seven described populations of distinct coelomocytes in the coelomic fluid that are defined by morphology, size, and for some types, by known functions. Of these subtypes, the large phagocytes are thought to be key to the sea urchin cellular innate immune response. The concentration of total coelomocytes in the coelomic fluid increases in response to pathogen challenge. However, there is no quantitative analysis of how the respective coelomocyte populations change over time in response to immune challenge. Accordingly, coelomocytes collected from immunoquiescent, healthy sea urchins were evaluated by flow cytometry for responses to injury and to challenge with either heat-killed Vibrio diazotrophicus, zymosan A, or artificial coelomic fluid, which served as the vehicle control. Responses to the initial injury of coelomic fluid collection or to injection of V. diazotrophicus show significant increases in the concentration of large phagocytes, small phagocytes, and red spherule cells after one day. Responses to zymosan A show decreases in the concentration of large phagocytes and increases in the concentration of small phagocytes. In contrast, responses to injections of vehicle result in decreased concentration of large phagocytes. When these changes in coelomocytes are evaluated based on proportions rather than concentration, the respective coelomocyte proportions are generally maintained in response to injection with V. diazotrophicus and vehicle. However, this is not observed in response to zymosan A and this lack of correspondence between proportions and concentrations may be an outcome of clearing these large particles by the large phagocytes. Variations in coelomocyte populations are also noted for individual sea urchins evaluated at different times for their responses to immune challenge compared to the vehicle. Together, these results demonstrate that the cell populations in sea urchin immune cell populations undergo dynamic changes in vivo in response to distinct immune stimuli and to injury and that these changes are driven by the responses of the large phagocyte populations.

Genetic variation underlies plastic responses to global change drivers in the purple sea urchin, Strongylocentrotus purpuratus.

Phenotypic plasticity and adaptive evolution enable population persistence in response to global change. However, there are few experiments that test how these processes interact within and across generations, especially in marine species with broad distributions experiencing spatially and temporally variable temperature and pCO2. We employed a quantitative genetics experiment with the purple sea urchin, Strongylocentrotus purpuratus, to decompose family-level variation in transgenerational and developmental plastic responses to ecologically relevant temperature and pCO2. Adults were conditioned to controlled non-upwelling (high temperature, low pCO2) or upwelling (low temperature, high pCO2) conditions. Embryos were reared in either the same conditions as their parents or the crossed environment, and morphological aspects of larval body size were quantified. We find evidence of family-level phenotypic plasticity in response to different developmental environments. Among developmental environments, there was substantial additive genetic variance for one body size metric when larvae developed under upwelling conditions, although this differed based on parental environment. Furthermore, cross-environment correlations indicate significant variance for genotype-by-environment interactive effects. Therefore, genetic variation for plasticity is evident in early stages of S. purpuratus, emphasizing the importance of adaptive evolution and phenotypic plasticity in organismal responses to global change.

Architecture and evolution of the cis-regulatory system of the echinoderm kirrelL gene.

The gene regulatory network (GRN) that underlies echinoderm skeletogenesis is a prominent model of GRN architecture and evolution. KirrelL is an essential downstream effector gene in this network and encodes an Ig-superfamily protein required for the fusion of skeletogenic cells and the formation of the skeleton. In this study, we dissected the transcriptional control region of the kirrelL gene of the purple sea urchin, Strongylocentrotus purpuratus. Using plasmid- and bacterial artificial chromosome-based transgenic reporter assays, we identified key cis-regulatory elements (CREs) and transcription factor inputs that regulate Sp-kirrelL, including direct, positive inputs from two key transcription factors in the skeletogenic GRN, Alx1 and Ets1. We next identified kirrelL cis-regulatory regions from seven other echinoderm species that together represent all classes within the phylum. By introducing these heterologous regulatory regions into developing sea urchin embryos we provide evidence of their remarkable conservation across ~500 million years of evolution. We dissected in detail the kirrelL regulatory region of the sea star, Patiria miniata, and demonstrated that it also receives direct inputs from Alx1 and Ets1. Our findings identify kirrelL as a component of the ancestral echinoderm skeletogenic GRN. They support the view that GRN subcircuits, including specific transcription factor-CRE interactions, can remain stable over vast periods of evolutionary history. Lastly, our analysis of kirrelL establishes direct linkages between a developmental GRN and an effector gene that controls a key morphogenetic cell behavior, cell-cell fusion, providing a paradigm for extending the explanatory power of GRNs.

Four plastic additives reduce larval growth and survival in the sea urchin Strongylocentrotus purpuratus.

Plastic additives are utilized during the production of plastic to modify the attributes and stability of the polymer. As oceanic plastic waste degrades, these additives can leach, and are harmful to global marine ecosystems. Despite the high abundance of additives leached into the marine environment, little is known about their direct impact on marine zooplankton. Here we test for impacts of four plastic additives, UV-327, Irganox 1010, DEHP, and methylparaben, all commonly used in plastic manufacturing, on purple sea urchin (Strongylocentrotus purpuratus) larval growth and survival in a serial dose response for 4 days. Methylparaben, UV-327, and Irganox 1010 significantly reduced larval body length by about 5% for at least one dose. In contrast, all compounds reduced larval survival by 20-70% with strongest effects at intermediate rather than high doses. Our results highlight that plastic additives should be tested for their effects on marine organisms.

A flow cytometry based approach to identify distinct coelomocyte subsets of the purple sea urchin, Strongylocentrotus purpuratus.

The sea urchin, Strongylocentrotus purpuratus, possesses at least seven distinguishable cell populations in the coelomic fluid, which vary in morphology, size, and function. Of these, the large phagocytes, small phagocytes, and red spherule cells are thought to be key to the echinoid immune response. Because there are currently no effective and rapid means of evaluating sea urchin coelomocytes, we developed a flow cytometry based approach to identify these subsets from unseparated, unstained, live cells. In particular our gating strategy distinguishes between the large phagocytes, small phagocytes, red spherule cells, and a mixed population of vibratile cells and colorless spherule cells. This flow cytometry based analysis increases the speed and improves the reliability of coelomocyte analysis compared to differential cell counts by microscopy.

Sequence Diversity, Locus Structure, and Evolutionary History of the SpTransformer Genes in the Sea Urchin Genome.

The generation of large immune gene families is often driven by evolutionary pressure exerted on host genomes by their pathogens, which has been described as the immunological arms race. The SpTransformer (SpTrf) gene family from the California purple sea urchin, Strongylocentrotus purpuratus, is upregulated upon immune challenge and encodes the SpTrf proteins that interact with pathogens during an immune response. Native SpTrf proteins bind both bacteria and yeast, and augment phagocytosis of a marine Vibrio, while a recombinant SpTrf protein (rSpTrf-E1) binds a subset of pathogens and a range of pathogen associated molecular patterns. In the sequenced sea urchin genome, there are four SpTrf gene clusters for a total of 17 genes. Here, we report an in-depth analysis of these genes to understand the sequence complexities of this family, its genomic structure, and to derive a putative evolutionary history for the formation of the gene clusters. We report a detailed characterization of gene structure including the intron type and UTRs with conserved transcriptional start sites, the start codon and multiple stop codons, and locations of polyadenylation signals. Phylogenetic and percent mismatch analyses of the genes and the intergenic regions allowed us to predict the last common ancestral SpTrf gene and a theoretical evolutionary history of the gene family. The appearance of the gene clusters from the theoretical ancestral gene may have been driven by multiple duplication and deletion events of regions containing SpTrf genes. Duplications and ectopic insertion events, indels, and point mutations in the exons likely resulted in the extant genes and family structure. This theoretical evolutionary history is consistent with the involvement of these genes in the arms race in responses to pathogens and suggests that the diversification of these genes and their encoded proteins have been selected for based on the survival benefits of pathogen binding and host protection.

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

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