Single-cell RNA-sequencing analysis of early sea star development.

Echinoderms represent a broad phylum with many tractable features to test evolutionary changes and constraints. Here, we present a single-cell RNA-sequencing analysis of early development in the sea star Patiria miniata, to complement the recent analysis of two sea urchin species. We identified 20 cell states across six developmental stages from 8 hpf to mid-gastrula stage, using the analysis of 25,703 cells. The clusters were assigned cell states based on known marker gene expression and by in situ RNA hybridization. We found that early (morula, 8-14 hpf) and late (blastula-to-mid-gastrula) cell states are transcriptionally distinct. Cells surrounding the blastopore undergo rapid cell state changes that include endomesoderm diversification. Of particular import to understanding germ cell specification is that we never see Nodal pathway members within Nanos/Vasa-positive cells in the region known to give rise to the primordial germ cells (PGCs). The results from this work contrast the results of PGC specification in the sea urchin, and the dataset presented here enables deeper comparative studies in tractable developmental models for testing a variety of developmental mechanisms.

Gene regulatory divergence amongst echinoderms underlies appearance of pigment cells in sea urchin development.

Larvae of the sea urchin, Strongylocentrotus purpuratus, have pigmented migratory cells implicated in immune defense and gut patterning. The transcription factor SpGcm activates the expression of many pigment cell-specific genes, including those involved in pigment biosynthesis (SpPks1 and SpFmo3) and immune related genes (e.g. SpMif5). Despite the importance of this cell type in sea urchins, pigmented cells are absent in larvae of the sea star, Patiria miniata. In this study, we tested the premises that sea stars lack genes to synthesize echinochrome pigment, that the genes are present but are not expressed in the larvae, or rather that the homologous gene expression does not contribute to echinochrome synthesis. Our results show that orthologs of sea urchin pigment cell-specific genes (PmPks1, PmFmo3-1 and PmMifL1-2) are present in the sea star genome and expressed in the larvae. Although no cell lineage homologous to migratory sea urchin pigment cells is present, dynamic gene activation accomplishes a similar spatial and temporal expression profile. The mechanisms regulating the expression of these genes, though, is highly divergent. In sea stars, PmGcm lacks the central role in pigment gene expression since it is not expressed in PmPks1 and PmFmo3-1-positive cells, and knockdown of Gcm does not abrogate pigment gene expression. Pigment genes are instead expressed in the coelomic mesoderm early in development before later being expressed in the ectoderm. These findings were supported by in situ RNA hybridization and comparative scRNA-seq analyses. We conclude that simply the coexpression of Pks1 and Fmo3 orthologs in cells of the sea star is not sufficient to underlie the emergence of the larval pigment cell in the sea urchin.

A single-cell RNA-seq analysis of Brachyury-expressing cell clusters suggests a morphogenesis-associated signal center of oral ectoderm in sea urchin embryos.

Brachyury is a T-box family transcription factor and plays pivotal roles in morphogenesis. In sea urchin embryos, Brachyury is expressed in the invaginating endoderm, and in the oral ectoderm of the invaginating mouth opening. The oral ectoderm is hypothesized to serve as a signaling center for oral (ventral)-aboral (dorsal) axis formation and to function as a ventral organizer. Our previous results of a single-cell RNA-seq (scRNA-seq) atlas of early Strongylocentrotus purpuratus embryos categorized the constituent cells into 22 clusters, in which the endoderm consists of three clusters and the oral ectoderm four clusters (Foster et al., 2020). Here we examined which clusters of cells expressed Brachyury in relation to the morphogenesis and the identity of the ventral organizer. Our results showed that cells of all three endoderm clusters expressed Brachyury in blastulae. Based on expression profiles of genes involved in the gene regulatory networks (GRNs) of sea urchin embryos, the three clusters are distinguishable, two likely derived from the Veg2 tier and one from the Veg1 tier. On the other hand, of the four oral-ectoderm clusters, cells of two clusters expressed Brachyury at the gastrula stage and genes that are responsible for the ventral organizer at the late blastula stage, but the other two clusters did not. At a single-cell level, most cells of the two oral-ectoderm clusters expressed organizer-related genes, nearly a half of which coincidently expressed Brachyury. This suggests that the ventral organizer contains Brachyury-positive cells which invaginate to form the stomodeum. This scRNA-seq study therefore highlights significant roles of Brachyury-expressing cells in body-plan formation of early sea urchin embryos, though cellular and molecular mechanisms for how Brachyury functions in these processes remain to be elucidated in future studies.

A single cell RNA sequencing resource for early sea urchin development.

Identifying cell states during development from their mRNA profiles provides insight into their gene regulatory network. Here, we leverage the sea urchin embryo for its well-established gene regulatory network to interrogate the embryo using single cell RNA sequencing. We tested eight developmental stages in Strongylocentrotus purpuratus, from the eight-cell stage to late in gastrulation. We used these datasets to parse out 22 major cell states of the embryo, focusing on key transition stages for cell type specification of each germ layer. Subclustering of these major embryonic domains revealed over 50 cell states with distinct transcript profiles. Furthermore, we identified the transcript profile of two cell states expressing germ cell factors, one we conclude represents the primordial germ cells and the other state is transiently present during gastrulation. We hypothesize that these cells of the Veg2 tier of the early embryo represent a lineage that converts to the germ line when the primordial germ cells are deleted. This broad resource will hopefully enable the community to identify other cell states and genes of interest to expose the underpinning of developmental mechanisms.

Conservation and contrast in cell states of echinoderm ovaries.

Echinoderms produce functional gametes throughout their lifespan, in some cases exceeding 200 years. The histology and ultrastructure of echinoderm ovaries has been described but how these ovaries function and maintain the production of high-quality gametes remains a mystery. Here, we present the first single cell RNA sequencing data sets of mature ovaries from two sea urchin species (Strongylocentrotus purpuratus [Sp] and Lytechinus variegatus [Lv]), and one sea star species (Patiria miniata [Pm]). We find 14 cell states in the Sp ovary, 16 cell states in the Lv ovary and 13 cell states in the ovary of the sea star. This resource is essential to understand the structure and functional biology of the ovary in echinoderms, and better informs decisions in the utilization of in situ RNA hybridization probes selective for various cell types. We link key genes with cell clusters in validation of this approach. This resource also aids in the identification of the stem cells for prolonged and continuous gamete production, is a foundation for testing changes in the annual reproductive cycle, and is essential for understanding the evolution of reproduction of this important phylum.

Single cell RNA-seq in the sea urchin embryo show marked cell-type specificity in the Delta/Notch pathway.

Sea urchin embryos are excellent for in vivo functional studies because of their transparency and tractability in manipulation. They are also favorites for pharmacological approaches since they develop in an aquatic environment and addition of test substances is straightforward. A concern in many pharmacological tests though is the potential for pleiotropic effects that confound the conclusions drawn from the results. Precise cellular interpretations are often not feasible because the impact of the perturbant is not known. Here we use single-cell mRNA (messenger RNA) sequencing as a metric of cell types in the embryo and to determine the selectivity of two commonly used inhibitors, one each for the Wnt and the Delta-Notch pathways, on these nascent cell types. We identified 11 distinct cell types based on mRNA profiling, and that the cell lineages affected by Wnt and Delta/Notch inhibition were distinct from each other. These data support specificity and distinct effects of these signaling pathways in the embryo and illuminate how these conserved pathways selectively regulate cell lineages at a single cell level. Overall, we conclude that single cell RNA-seq analysis in this embryo is revealing of the cell types present during development, of the changes in the gene regulatory network resulting from inhibition of various signaling pathways, and of the selectivity of these pathways in influencing developmental trajectories.

Elements of divergence in germline determination in closely related species.

Evolutionary transitions are particularly important in development of the germ line, cells which directly impact sexual reproduction. Differences in the primordial germ cells (PGCs) of two sea urchin species were examined here by stage-matched, integrated, single cell RNA-seq (scRNA-seq) datasets. Even though both species rely on inherited mechanisms to specify their germ line, this analysis revealed a variety of differences in germline gene expression, including a broader expression of the germline factor Nanos2 (Nan2) in Lytechinus variegatus (Lv) compared to Strongylocentrotus purpuratus (Sp). In Sp, Nan2 mRNA expression is highly restricted to the PGCs by a lability element in its 3'UTR, which is lacking in the mRNA of Lv-Nan2, thus explaining the difference. We discovered that the Lv-Nan2 3'UTR instead leads to its specific translation in the PGCs. The results emphasize that regulatory mechanisms resulting in germline specification rely greatly on post-transcriptional restrictions of key gene products.

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.

Identifying gene expression from single cells to single genes.

Gene regulatory networks reveal how transcription factors contribute to a dynamic cascade of cellular information processing. Recent advances in technologies have enhanced the toolkit for testing GRN mechanisms and connections. Here we emphasize three approaches that we have found important for interrogating transcriptional mechanisms in echinoderms: single cell mRNA sequencing (drop-seq), nascent RNA detection and identification, and chromatin immunoprecipitation (ChIP). We present these applications in order since it is a logical experimental protocol. With preliminary information from bulk mRNA transcriptome analysis and differential gene expression studies (DE-seq), one may need to test in what specific cells important genes may be expressed and to use single cell sequencing to define such links. Nascent RNA analysis with the Click-iT chemistry allows the investigator to deduce when the RNA was transcribed, not just identify its presence, and ChIP allows the investigator to study direct interactions of putative transcriptional regulators with the gene promoter of interest. This flow of thinking, and the technologies to support it, is presented here for echinoderms. While many of the procedures are general and applicable to many organisms and cell types, we emphasize unique aspects of the protocols for consideration in using echinoderm embryos, larvae, and adult tissues.