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.

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.

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.

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.

Single-cell RNA sequencing of the Strongylocentrotus purpuratus larva reveals the blueprint of major cell types and nervous system of a non-chordate deuterostome.

Identifying the molecular fingerprint of organismal cell types is key for understanding their function and evolution. Here, we use single-cell RNA sequencing (scRNA-seq) to survey the cell types of the sea urchin early pluteus larva, representing an important developmental transition from non-feeding to feeding larva. We identify 21 distinct cell clusters, representing cells of the digestive, skeletal, immune, and nervous systems. Further subclustering of these reveal a highly detailed portrait of cell diversity across the larva, including the identification of neuronal cell types. We then validate important gene regulatory networks driving sea urchin development and reveal new domains of activity within the larval body. Focusing on neurons that co-express Pdx-1 and Brn1/2/4, we identify an unprecedented number of genes shared by this population of neurons in sea urchin and vertebrate endocrine pancreatic cells. Using differential expression results from Pdx-1 knockdown experiments, we show that Pdx1 is necessary for the acquisition of the neuronal identity of these cells. We hypothesize that a network similar to the one orchestrated by Pdx1 in the sea urchin neurons was active in an ancestral cell type and then inherited by neuronal and pancreatic developmental lineages in sea urchins and vertebrates.

microRNA-31 regulates skeletogenesis by direct suppression of Eve and Wnt1.

microRNAs (miRNAs) play a critical role in a variety of biological processes, including embryogenesis and the physiological functions of cells. Evolutionarily conserved microRNA-31 (miR-31) has been found to be involved in cancer, bone formation, and lymphatic development. We previously discovered that, in the sea urchin, miR-31 knockdown (KD) embryos have shortened dorsoventral connecting rods, mispatterned skeletogenic primary mesenchyme cells (PMCs) and shifted and expanded Vegf3 expression domain. Vegf3 itself does not contain miR-31 binding sites; however, we identified its upstream regulators Eve and Wnt1 to be directly suppressed by miR-31. Removal of miR-31's suppression of Eve and Wnt1 resulted in skeletal and PMC patterning defects, similar to miR-31 KD phenotypes. Additionally, removal of miR-31's suppression of Eve and Wnt1 results in an expansion and anterior shift in expression of Veg1 ectodermal genes, including Vegf3 in the blastulae. This indicates that miR-31 indirectly regulates Vegf3 expression through directly suppressing Eve and Wnt1. Furthermore, removing miR-31 suppression of Eve is sufficient to cause skeletogenic defects, revealing a novel regulatory role of Eve in skeletogenesis and PMC patterning. Overall, this study provides a proposed molecular mechanism of miR-31's regulation of skeletogenesis and PMC patterning through its cross-regulation of a Wnt signaling ligand and a transcription factor of the endodermal and ectodermal gene regulatory network.

CRISPR-Cas9 editing of non-coding genomic loci as a means of controlling gene expression in the sea urchin.

We seek to manipulate gene function here through CRISPR-Cas9 editing of cis-regulatory sequences, rather than the more typical mutation of coding regions. This approach would minimize secondary effects of cellular responses to nonsense mediated decay pathways or to mutant protein products by premature stops. This strategy also allows for reducing gene activity in cases where a complete gene knockout would result in lethality, and it can be applied to the rapid identification of key regulatory sites essential for gene expression. We tested this strategy here with genes of known function as a proof of concept, and then applied it to examine the upstream genomic region of the germline gene Nanos2 in the sea urchin, Strongylocentrotus purpuratus. We first used CRISPR-Cas9 to target established genomic cis-regulatory regions of the skeletogenic cell transcription factor, Alx1, and the TGF-ß signaling ligand, Nodal, which produce obvious developmental defects when altered in sea urchin embryos. Importantly, mutation of cis-activator sites (Alx1) and cis-repressor sites (Nodal) result in the predicted decreased and increased transcriptional output, respectively. Upon identification of efficient gRNAs by genomic mutations, we then used the same validated gRNAs to target a deadCas9-VP64 transcriptional activator to increase Nodal transcription directly. Finally, we paired these new methodologies with a more traditional, GFP reporter construct approach to further our understanding of the transcriptional regulation of Nanos2, a key gene required for germ cell identity in S. purpuratus. With a series of reporter assays, upstream Cas9-promoter targeted mutagenesis, coupled with qPCR and in situ RNA hybridization, we concluded that the promoter of Nanos2 drives strong mRNA expression in the sea urchin embryo, indicating that its primordial germ cell (PGC)-specific restriction may rely instead on post-transcriptional regulation. Overall, we present a proof-of-principle tool-kit of Cas9-mediated manipulations of promoter regions that should be applicable in most cells and embryos for which CRISPR-Cas9 is employed.

Regulation of dynamic pigment cell states at single-cell resolution.

Cells bearing pigment have diverse roles and are often under strict evolutionary selection. Here, we explore the regulation of pigmented cells in the purple sea urchin Strongylocentrotus purpuratus, an emerging model for diverse pigment function. We took advantage of single cell RNA-seq (scRNAseq) technology and discovered that pigment cells in the embryo segregated into two distinct populations, a mitotic cluster and a post-mitotic cluster. Gcm is essential for expression of several genes important for pigment function, but is only transiently expressed in these cells. We discovered unique genes expressed by pigment cells and test their expression with double fluorescence in situ hybridization. These genes include new members of the fmo family that are expressed selectively in pigment cells of the embryonic and in the coelomic cells of the adult - both cell-types having immune functions. Overall, this study identifies nodes of molecular intersection ripe for change by selective evolutionary pressures.

Unique Genomic and Phenotypic Responses to Extreme and Variable pH Conditions in Purple Urchin Larvae.

Environmental variation experienced by a species across space and time can promote the maintenance of genetic diversity that may be adaptive in future global change conditions. Selection experiments have shown that purple sea urchin, Strongylocentrotus purpuratus, populations have adaptive genetic variation for surviving pH conditions at the "edge" (pH 7.5) of conditions experienced in nature. However, little is known about whether populations have genetic variation for surviving low-pH events beyond those currently experienced in nature or how variation in pH conditions affects organismal and genetic responses. Here, we quantified survival, growth, and allele frequency shifts in experimentally selected developing purple sea urchin larvae in static and variable conditions at three pH levels: pH 8.1 (control), pH 7.5 (edge-of-range), and pH 7.0 (extreme). Variable treatments recovered body size relative to static treatments, but resulted in higher mortality, suggesting a potential tradeoff between survival and growth under pH stress. However, within each pH level, allele frequency changes were overlapping between static and variable conditions, suggesting a shared genetic basis underlying survival to mean pH regardless of variability. In contrast, genetic responses to pH 7.5 (edge) versus pH 7.0 (extreme) conditions were distinct, indicating a unique genetic basis of survival. In addition, loci under selection were more likely to be in exonic regions than regulatory, indicating that selection targeted protein-coding variation. Loci under selection in variable pH 7.5 conditions, more similar to conditions periodically experienced in nature, performed functions related to lipid biosynthesis and metabolism, while loci under selection in static pH 7.0 conditions performed functions related to transmembrane and mitochondrial processes. While these results are promising in that purple sea urchin populations possess genetic variation for surviving extreme pH conditions not currently experienced in nature, they caution that increased acidification does not result in a linear response but elicits unique physiological stresses and survival mechanisms.