A single cell RNA sequencing resource for early sea urchin development

mRNA splicing variants of the transcription factor Blimp1 differentially regulate germline genes in echinoderms

Germ cell specification is an essential step in sexually reproducing animals. Echinoderms possess diverse representatives of the main mechanisms that result in this cell fate determination. Sea urchins use an inherited mechanism, whereas sea stars rely on the ancestral, induced mechanism. Blimp1 (B lymphocyte-induced maturation protein-1) is a transcriptional regulator reported in mice to function in the induction of germline cells. Here, we identify the dynamic function of Blimp1 during development in a comparative approach using the purple sea urchin, Strongylocentrotus purpuratus (inherited germline) and the batstar, Patiria miniata (induced germline). We found that Blimp1 is important for germ cell specification in both species and that multiple Blimp1 isoforms result from differential mRNA splicing in each animal. Each isoform of Blimp1 functions in distinct expression of germline determinants, including Vasa and Nanos. These results show that Blimp1 is a conserved and key regulator for germ cell specification, but divergent in function as a result of post-transcriptional modification. Overall, we conclude that Blimp1 is an intersectional node in diverse germline specification strategies and supports the concept that differential mRNA splicing is an essential mechanism in germ cell formation.

The Cell States of Sea Urchin During Metamorphosis Revealed by Single-Cell RNA Sequencing

Metamorphosis is a key process in the life history of sea urchin Heliocidaris crassispina. However, the understanding of its molecular mechanisms is still lacking, especially the basic cell biology pre-metamorphosis and post-metamorphosis. Therefore, we employed single-cell RNA sequencing to delineate the cellular states of larvae and juveniles of H. crassispina. Our investigation revealed that the cell composition in sea urchins comprises six primary populations, encompassing nerve cells, skeletogenic cells, immune cells, digestive cells, germ cells, and muscle cells. Subsequently, we identified subpopulations within these cells. Our findings indicated that the larval peripheral nerves were discarded during metamorphosis. A decrease in the number of spicules was observed during this process. Additionally, we examined the differences between larval and adult pigment cells. Meanwhile, cellulase is highlighted as an essential factor for the development of competent juveniles. In summary, this study not only serves as a valuable resource for future research on sea urchins but also deepens our understanding of the intricate metamorphosis process.

An ancient apical patterning system sets the position of the forebrain in chordates

The evolutionary origin of the vertebrate brain remains a major subject of debate, as its development from a dorsal tubular neuroepithelium is unique to chordates. To shed light on the evolutionary emergence of the vertebrate brain, we compared anterior neuroectoderm development across deuterostome species, using available single-cell datasets from sea urchin, amphioxus, and zebrafish embryos. We identified a conserved gene co-expression module, comparable to the anterior gene regulatory network (aGRN) controlling apical organ development in ambulacrarians, and spatially mapped it by multiplexed in situ hybridization to the developing retina and hypothalamus of chordates. Using functional approaches, we show Wnt signaling regulating this co-expression module in amphioxus, like the aGRN in echinoderms, and that its overactivation suppresses forebrain identity. This suggests a previously undescribed role for Wnt signaling in amphioxus in determining the position of the forebrain. We propose this Wnt-regulated gene co-expression module as a possible mechanism by which the brain set antero-dorsally early in chordate evolution.

Benchmarking cross-species single-cell RNA-seq data integration methods: towards a cell type tree of life

Cross-species single-cell RNA-seq data hold immense potential for unraveling cell type evolution and transferring knowledge between well-explored and less-studied species. However, challenges arise from interspecific genetic variation, batch effects stemming from experimental discrepancies and inherent individual biological differences. Here, we benchmarked nine data-integration methods across 20 species, encompassing 4.7 million cells, spanning eight phyla and the entire animal taxonomic hierarchy. Our evaluation reveals notable differences between the methods in removing batch effects and preserving biological variance across taxonomic distances. Methods that effectively leverage gene sequence information capture underlying biological variances, while generative model-based approaches excel in batch effect removal. SATURN demonstrates robust performance across diverse taxonomic levels, from cross-genus to cross-phylum, emphasizing its versatility. SAMap excels in integrating species beyond the cross-family level, especially for atlas-level cross-species integration, while scGen shines within or below the cross-class hierarchy. As a result, our analysis offers recommendations and guidelines for selecting suitable integration methods, enhancing cross-species single-cell RNA-seq analyses and advancing algorithm development.

Contrasting the development of larval and adult body plans during the evolution of biphasic lifecycles in sea urchins

Biphasic lifecycles are widespread among animals, but little is known about how the developmental transition between larvae and adults is regulated. Sea urchins are a unique system for studying this phenomenon because of the stark differences between their bilateral larval and pentaradial adult body plans. Here, we use single-cell RNA sequencing to analyze the development of Heliocidaris erythrogramma (He), a sea urchin species with an accelerated, non-feeding mode of larval development. The sequencing time course extends from embryogenesis to roughly a day before the onset of metamorphosis in He larvae, which is a period that has not been covered by previous datasets. We find that the non-feeding developmental strategy of He is associated with several changes in the specification of larval cell types compared to sea urchins with feeding larvae, such as the loss of a larva-specific skeletal cell population. Furthermore, the development of the larval and adult body plans in sea urchins may utilize largely different sets of regulatory genes. These findings lay the groundwork for extending existing developmental gene regulatory networks to cover additional stages of biphasic lifecycles.

Identification and expression analysis of Sox family genes in echinoderms

The Sox gene family, a collection of transcription factors widely distributed throughout the animal kingdom, plays a crucial role in numerous developmental processes. Echinoderms occupy a pivotal position in many research fields, such as neuroscience, sex determination and differentiation, and embryonic development. However, to date, no comprehensive study has been conducted to characterize and analyze Sox genes in echinoderms. In the present study, the evolution and expression of Sox family genes across 11 echinoderms were analyzed using bioinformatics methods. The results revealed a total of 70 Sox genes, with counts ranging from 5 to 8 across different echinoderms. Phylogenetic analysis revealed that the identified Sox genes could be categorized into seven distinct classes: the SoxB1 class, SoxB2 class, SoxC class, SoxD class, SoxE class, SoxF class and SoxH class. Notably, the SoxB1, SoxB2, and SoxF genes were ubiquitously present in all the echinoderms studied, which suggests that these genes may be conserved in echinoderms. The spatiotemporal expression patterns observed for Sox genes in the three echinoderms indicated that various Sox members perform distinct functional roles. Notably, SoxB1 is likely involved in echinoderm ovary development, while SoxH may play a crucial role in testis development in starfish and sea cucumber. In general, the present investigation provides a molecular foundation for exploring the Sox gene in echinoderms, providing a valuable resource for future phylogenetic and genomic studies.

Echinobase: a resource to support the echinoderm research community

Echinobase (www.echinobase.org) is a model organism knowledgebase serving as a resource for the community that studies echinoderms, a phylum of marine invertebrates that includes sea urchins and sea stars. Echinoderms have been important experimental models for over 100 years and continue to make important contributions to environmental, evolutionary, and developmental studies, including research on developmental gene regulatory networks. As a centralized resource, Echinobase hosts genomes and collects functional genomic data, reagents, literature, and other information for the community. This third-generation site is based on the Xenbase knowledgebase design and utilizes gene-centric pages to minimize the time and effort required to access genomic information. Summary gene pages display gene symbols and names, functional data, links to the JBrowse genome browser, and orthology to other organisms and reagents, and tabs from the Summary gene page contain more detailed information concerning mRNAs, proteins, diseases, and protein-protein interactions. The gene pages also display 1:1 orthologs between the fully supported species Strongylocentrotus purpuratus (purple sea urchin), Lytechinus variegatus (green sea urchin), Patiria miniata (bat star), and Acanthaster planci (crown-of-thorns sea star). JBrowse tracks are available for visualization of functional genomic data from both fully supported species and the partially supported species Anneissia japonica (feather star), Asterias rubens (sugar star), and L. pictus (painted sea urchin). Echinobase serves a vital role by providing researchers with annotated genomes including orthology, functional genomic data aligned to the genomes, and curated reagents and data. The Echinoderm Anatomical Ontology provides a framework for standardizing developmental data across the phylum, and knowledgebase content is formatted to be findable, accessible, interoperable, and reusable by the research community.

Single-cell transcriptomics reveals evolutionary reconfiguration of embryonic cell fate specification in the sea urchin Heliocidaris erythrogramma

Altered regulatory interactions during development likely underlie a large fraction of phenotypic diversity within and between species, yet identifying specific evolutionary changes remains challenging. Analysis of single-cell developmental transcriptomes from multiple species provides a powerful framework for unbiased identification of evolutionary changes in developmental mechanisms. Here, we leverage a "natural experiment" in developmental evolution in sea urchins, where a major life history switch recently evolved in the lineage leading to Heliocidaris erythrogramma, precipitating extensive changes in early development. Comparative analyses of scRNA-seq developmental time courses from H. erythrogramma and Lytechinus variegatus (representing the derived and ancestral states respectively) reveals numerous evolutionary changes in embryonic patterning. The earliest cell fate specification events, and the primary signaling center are co-localized in the ancestral dGRN but remarkably, in H. erythrogramma they are spatially and temporally separate. Fate specification and differentiation are delayed in most embryonic cell lineages, although in some cases, these processes are conserved or even accelerated. Comparative analysis of regulator-target gene co-expression is consistent with many specific interactions being preserved but delayed in H. erythrogramma, while some otherwise widely conserved interactions have likely been lost. Finally, specific patterning events are directly correlated with evolutionary changes in larval morphology, suggesting that they are directly tied to the life history shift. Together, these findings demonstrate that comparative scRNA-seq developmental time courses can reveal a diverse set of evolutionary changes in embryonic patterning and provide an efficient way to identify likely candidate regulatory interactions for subsequent experimental validation.

scRNA-seq analysis of cells comprising the amphioxus notochord

Cephalochordates occupy a key phylogenetic position for deciphering the origin and evolution of chordates, since they diverged earlier than urochordates and vertebrates. The notochord is the most prominent feature of chordates. The amphioxus notochord features coin-shaped cells bearing myofibrils. Notochord-derived hedgehog signaling contributes to patterning of the dorsal nerve cord, as in vertebrates. However, properties of constituent notochord cells remain unknown at the single-cell level. We examined these properties using Iso-seq analysis, single-cell RNA-seq analysis, and in situ hybridization (ISH). Gene expression profiles broadly categorize notochordal cells into myofibrillar cells and non-myofibrillar cells. Myofibrillar cells occupy most of the central portion of the notochord, and some cells extend the notochordal horn to both sides of the ventral nerve cord. Some notochord myofibrillar genes are not expressed in myotomes, suggesting an occurrence of myofibrillar genes that are preferentially expressed in notochord. On the other hand, non-myofibrillar cells contain dorsal, lateral, and ventral Müller cells, and all three express both hedgehog and Brachyury. This was confirmed by ISH, although expression of hedgehog in ventral Müller cells was minimal. In addition, dorsal Müller cells express neural transmission-related genes, suggesting an interaction with nerve cord. Lateral Müller cells express hedgehog and other signaling-related genes, suggesting an interaction with myotomes positioned lateral to the notochord. Ventral Müller cells also expressed genes for FGF- and EGF-related signaling, which may be associated with development of endoderm, ventral to the notochord. Lateral Müller cells were intermediate between dorsal/ventral Müller cells. Since vertebrate notochord contributes to patterning and differentiation of ectoderm (nerve cord), mesoderm (somite), and endoderm, this investigation provides evidence that an ancestral or original form of vertebrate notochord is present in extant cephalochordates.

Marine Invertebrates One Cell at A Time: Insights from Single-Cell Analysis

Over the past decade, single-cell RNA-sequencing (scRNA-seq) has made it possible to study the cellular diversity of a broad range of organisms. Technological advances in single-cell isolation and sequencing have expanded rapidly, allowing the transcriptomic profile of individual cells to be captured. As a result, there has been an explosion of cell type atlases created for many different marine invertebrate species from across the tree of life. Our focus in this review is to synthesize current literature on marine invertebrate scRNA-seq. Specifically, we provide perspectives on key insights from scRNA-seq studies, including descriptive studies of cell type composition, how cells respond in dynamic processes such as development and regeneration, and the evolution of new cell types. Despite these tremendous advances, there also lie several challenges ahead. We discuss the important considerations that are essential when making comparisons between experiments, or between datasets from different species. Finally, we address the future of single-cell analyses in marine invertebrates, including combining scRNA-seq data with other 'omics methods to get a fuller understanding of cellular complexities. The full diversity of cell types across marine invertebrates remains unknown and understanding this diversity and evolution will provide rich areas for future study.

Solute carrier (SLC) expression reveals skeletogenic cell diversity

Within the developing embryo is a microcosm of cell type diversity. Single cell RNA-sequencing (scRNA-seq) is used to reveal cell types, typically by grouping cells according to their gene regulatory states. However, both across and within these regulatory states are additional layers of cellular diversity represented by the differential expression of genes that govern cell function. Here, we analyzed scRNA-seq data representing the late gastrula stage of Strongylocentrotus purpuratus (purple sea urchin) to understand the patterning of transporters belonging to the ABC and SLC families. These transporters handle diverse substrates from amino acids to signaling molecules, nutrients and xenobiotics. Using transporter-based clustering, we identified unique transporter patterns that are both shared across cell lineages, as well as those that were unique to known cell types. We further explored three patterns of transporter expression in mesodermal cells including secondary mesenchyme cells (pigment cells and blastocoelar cells) and skeletogenic cells (primary mesenchyme cells). The results revealed the enrichment of SMTs potentially involved in nutrient absorption (SLC5A9, SLC7A11, SLC35F3, and SLC52A3) and skeletogenesis (SLC9A3, SLC13A2/3/5, and SLC39A13) in pigment cells and blastocoelar cells respectively. The results indicated that the strategy of clustering by cellular activity can be useful for discovering cellular populations that would otherwise remain obscured.

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.

The crucial role of CTCF in mitotic progression during early development of sea urchin

CCCTC-binding factor (CTCF), an insulator protein with 11 zinc fingers, is enriched at the boundaries of topologically associated domains (TADs) in eukaryotic genomes. In this study, we isolated and analyzed the cDNAs encoding HpCTCF, the CTCF homolog in the sea urchin Hemicentrotus pulcherrimus, to investigate its expression patterns and functions during the early development of sea urchin. HpCTCF contains nine zinc fingers corresponding to fingers 2-10 of the vertebrate CTCF. Expression pattern analysis revealed that HpCTCF mRNA was detected at all developmental stages and in the entire embryo. Upon expressing the HpCTCF-GFP fusion protein in early embryos, we observed its uniform distribution within interphase nuclei. However, during mitosis, it disappeared from the chromosomes and subsequently reassembled on the chromosome during telophase. Moreover, the morpholino-mediated knockdown of HpCTCF resulted in mitotic arrest during the morula to blastula stage. Most of the arrested chromosomes were not phospholylated at serine 10 of histone H3, indicating that mitosis was arrested at the telophase by HpCTCF depletion. Furthermore, impaired sister chromatid segregation was observed using time-lapse imaging of HpCTCF-knockdown embryos. Thus, HpCTCF is essential for mitotic progression during the early development of sea urchins, especially during the telophase-to-interphase transition. However, the normal development of pluteus larvae in CRISPR-mediated HpCTCF-knockout embryos suggests that disruption of zygotic HpCTCF expression has little effect on embryonic and larval development.

Single-Cell Transcriptomic Analysis Reveals the Molecular Profile of Go-Opsin Photoreceptor Cells in Sea Urchin Larvae

The ability to perceive and respond to light stimuli is fundamental not only for spatial vision but also to many other light-mediated interactions with the environment. In animals, light perception is performed by specific cells known as photoreceptors and, at molecular level, by a group of GPCRs known as opsins. Sea urchin larvae possess a group of photoreceptor cells (PRCs) deploying a Go-Opsin (Opsin3.2) which have been shown to share transcription factors and morphology with PRCs of the ciliary type, raising new questions related to how this sea urchin larva PRC is specified and whether it shares a common ancestor with ciliary PRCs or it if evolved independently through convergent evolution. To answer these questions, we combined immunohistochemistry and fluorescent in situ hybridization to investigate how the Opsin3.2 PRCs develop in the sea urchin Strongylocentrotus purpuratus larva. Subsequently, we applied single-cell transcriptomics to investigate the molecular signature of the Sp-Opsin3.2-expressing cells and show that they deploy an ancient regulatory program responsible for photoreceptors specification. Finally, we also discuss the possible functions of the Opsin3.2-positive cells based on their molecular fingerprint, and we suggest that they are involved in a variety of signaling pathways, including those entailing the thyrotropin-releasing hormone.

Tracing the history of cell types

A study of sea urchin and sea star larvae paves the way for understanding how cell types evolve and give rise to novel morphologies.

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.

Evolution of homology: From archetype towards a holistic concept of cell type

The concept of homology lies in the heart of comparative biological science. The distinction between homology as structure and analogy as function has shaped the evolutionary paradigm for a century and formed the axis of comparative anatomy and embryology, which accept the identity of structure as a ground measure of relatedness. The advent of single-cell genomics overturned the classical view of cell homology by establishing a backbone regulatory identity of cell types, the basic biological units bridging the molecular and phenotypic dimensions, to reveal that the cell is the most flexible unit of living matter and that many approaches of classical biology need to be revised to understand evolution and diversity at the cellular level. The emerging theory of cell types explicitly decouples cell identity from phenotype, essentially allowing for the divergence of evolutionarily related morphotypes beyond recognition, as well as it decouples ontogenetic cell lineage from cell-type phylogeny, whereby explicating that cell types can share common descent regardless of their structure, function or developmental origin. The article succinctly summarizes current progress and opinion in this field and formulates a more generalistic view of biological cell types as avatars, transient or terminal cell states deployed in a continuum of states by the developmental programme of one and the same omnipotent cell, capable of changing or combining identities with distinct evolutionary histories or inventing ad hoc identities that never existed in evolution or development. It highlights how the new logic grounded in the regulatory nature of cell identity transforms the concepts of cell homology and phenotypic stability, suggesting that cellular evolution is inherently and massively network-like, with one-to-one homologies being rather uncommon and restricted to shallower levels of the animal tree of life.

A single-cell RNA-seq analysis of early larval cell-types of the starfish, Patiria pectinifera: Insights into evolution of the chordate body plan

Ambulacrarians (echinoderms and hemichordates) are a sister group to chordates; thus, their larval cell-types may provide clues about evolution of chordate body plans. Although most genic information accumulated to date pertains to sea urchin embryogenesis, starfish embryogenesis represents a more ancestral mode than that of sea urchins. We performed single-cell RNA-seq analysis of cell-types from gastrulae and bipinnarial larvae of the starfish, Patiria pectinifera, and categorized them into 22 clusters, each of which is composed of cells with specific, shared profiles of development-relevant gene expression. Oral and aboral ectoderm, apical plate, hindgut or archenteron, midgut or intestine, pharynx, endomesoderm, stomodeum, and mesenchyme of the gastrulae, and neurons, ciliary bands, enterocoel and muscle of larvae were characterized by expression profiles of at least two relevant transcription factor genes and signaling molecular genes. Expression of Hox2, Hox7, Hox9/10, and Hox11/13b was detected in cells of clusters that form the larval enterocoel. By comparing homologous gene expression profiles in chordate embryos, we discuss and propose how the chordate body plan evolved from a deuterostome ancestor, from which the echinoderm body plan also evolved.

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.

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.

Different sea urchin RAG-like genes were domesticated to carry out different functions

The closely linked recombination activating genes (RAG1 and RAG2) in vertebrates encode the core of the RAG recombinase that mediates the V(D)J recombination of the immunoglobulin and T-cell receptor genes. RAG1 and RAG2 homologues (RAG1L and RAG2L) are present in multiple invertebrate phyla, including mollusks, nemerteans, cnidarians, and sea urchins. However, the function of the invertebrates' RAGL proteins is yet unknown. The sea urchins contain multiple RAGL genes that presumably originated in a common ancestral transposon. In this study, we demonstrated that two different RAG1L genes in the sea urchin Paracentrutus lividus (PlRAG1La and PlRAG1Lb) lost their mobility and, along with PlRAG2L, were fully domesticated to carry out different functions. We found that the examined echinoid RAGL homologues have distinct expression profiles in early developmental stages and in adult tissues. Moreover, the predicted structure of the proteins suggests that while PlRAG1La could maintain its endonuclease activity and create a heterotetramer with PlRAG2L, the PlRAG1Lb adopted a different function that does not include an interaction with DNA nor a collaboration with PlRAG2L. By characterizing the different RAG homologues in the echinoid lineage, we hope to increase the knowledge about the evolution of these genes and shed light on their domestication processes.

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.

Single-Cell RNA Sequencing Atlas From a Bivalve Larva Enhances Classical Cell Lineage Studies

Ciliated trochophore-type larvae are widespread among protostome animals with spiral cleavage. The respective phyla are often united into the superclade Spiralia or Lophotrochozoa that includes, for example, mollusks, annelids, and platyhelminths. Mollusks (bivalves, gastropods, cephalopods, polyplacophorans, and their kin) in particular are known for their morphological innovations and lineage-specific plasticity of homologous characters (e.g., radula, shell, foot, neuromuscular systems), raising questions concerning the cell types and the molecular toolkit that underlie this variation. Here, we report on the gene expression profile of individual cells of the trochophore larva of the invasive freshwater bivalve Dreissena rostriformis as inferred from single cell RNA sequencing. We generated transcriptomes of 632 individual cells and identified seven transcriptionally distinct cell populations. Developmental trajectory analyses identify cell populations that, for example, share an ectodermal origin such as the nervous system, the shell field, and the prototroch. To annotate these cell populations, we examined ontology terms from the gene sets that characterize each individual cluster. These were compared to gene expression data previously reported from other lophotrochozoans. Genes expected to be specific to certain tissues, such as Hox1 (in the shell field), Caveolin (in prototrochal cells), or FoxJ (in other cillia-bearing cells) provide evidence that the recovered cell populations contribute to various distinct tissues and organs known from morphological studies. This dataset provides the first molecular atlas of gene expression underlying bivalve organogenesis and generates an important framework for future comparative studies into cell and tissue type development in Mollusca and Metazoa as a whole.

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.

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.

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.

Developmental single-cell transcriptomics in the Lytechinus variegatus sea urchin embryo

Using scRNA-seq coupled with computational approaches, we studied transcriptional changes in cell states of sea urchin embryos during development to the larval stage. Eighteen closely spaced time points were taken during the first 24 h of development of Lytechinus variegatus (Lv). Developmental trajectories were constructed using Waddington-OT, a computational approach to 'stitch' together developmental time points. Skeletogenic and primordial germ cell trajectories diverged early in cleavage. Ectodermal progenitors were distinct from other lineages by the 6th cleavage, although a small percentage of ectoderm cells briefly co-expressed endoderm markers that indicated an early ecto-endoderm cell state, likely in cells originating from the equatorial region of the egg. Endomesoderm cells also originated at the 6th cleavage and this state persisted for more than two cleavages, then diverged into distinct endoderm and mesoderm fates asynchronously, with some cells retaining an intermediate specification status until gastrulation. Seventy-nine out of 80 genes (99%) examined, and included in published developmental gene regulatory networks (dGRNs), are present in the Lv-scRNA-seq dataset and are expressed in the correct lineages in which the dGRN circuits operate.

A Preliminary Single-Cell RNA-Seq Analysis of Embryonic Cells That Express Brachyury in the Amphioxus, Branchiostoma japonicum

Among chordate taxa, the cephalochordates diverged earlier than urochordates and vertebrates; thus, they retain unique, primitive developmental features. In particular, the amphioxus notochord has muscle-like properties, a feature not seen in urochordates or vertebrates. Amphioxus contains two Brachyury genes, Bra1 and Bra2. Bra2 is reportedly expressed in the blastopore, notochord, somites, and tail bud, in contrast to a low level of Bra1 expression only in notochord. To distinguish the expression profiles of the two Brachyury genes at the single-cell level, we carried out single-cell RNA-seq (scRNA-seq) analysis using the amphioxus, Branchiostoma japonicum. This scRNA-seq analysis classified B. japonicum embryonic cells into 15 clusters at developmental stages from midgastrula to early swimming larva. Brachyury was expressed in cells of clusters 4, 5, 8, and 9. We first confirmed that cluster 8 comprises cells that form somites since this cluster specifically expresses four myogenic factor genes. Cluster 9 contains a larger number of cells with high levels of Bra2 expression and a smaller number of cells with Bra1 expression. Simultaneous expression in cluster 9 of tool-kit genes, including FoxA, Goosecoid, and hedgehog, showed that this cluster comprises cells that form the notochord. Expression of Bra2, but not Bra1, in cells of clusters 4 and 5 at the gastrula stage together with expression of Wnt1 and Caudal indicates that clusters 4 and 5 comprise cells of the blastopore, which contiguously form the tail bud. In addition, Hox1, Hox3, and Hox4 were highly expressed in Bra2-expressing clusters 4, 5, 8, and 9 in a temporally coordinated manner, suggesting roles of anterior Hox genes in specification of mesodermal organs, including somites, notochord, and tail bud. This scRNA-seq analysis therefore highlights differences between the two Brachyury genes in relation to embryonic regions in which they are expressed and their levels of expression. Bra2 is the ancestral Brachyury in amphioxus, since expression in the blastopore is shared with other deuterostomes. On the other hand, Bra1 is a duplicate copy and likely evolved a supplementary function in notochord and somite formation in the Branchiostoma lineage.

New techniques for creating parthenogenetic larvae of the sea urchin Lytechinus pictus for gene expression studies

BACKGROUND

Sea urchins are model organisms for studying the spatial-temporal control of gene activity during development. The Southern California species, Lytechinus pictus, has a sequenced genome and can be raised in the laboratory from egg to egg in 4 to 5 months.

RESULTS

Here, we present new techniques for generating parthenogenetic larvae of this species and include a gallery of photomicrographs of morphologically abnormal larvae that could be used for transcriptomic analysis.

CONCLUSIONS

Comparison of gene expression in parthenogenotes to larvae produced by fertilization could provide novel insights into gene expression controls contributed by sperm in this important model organism. Knowledge gained from transcriptomics of sea urchin parthenogenotes could contribute to parthenogenetic studies of mammalian embryos.

Evolutionary Cell Type Mapping with Single-Cell Genomics

A fundamental characteristic of animal multicellularity is the spatial coexistence of functionally specialized cell types that are all encoded by a single genome sequence. Cell type transcriptional programs are deployed and maintained by regulatory mechanisms that control the asymmetric, differential access to genomic information in each cell. This genome regulation ultimately results in specific cellular phenotypes. However, the emergence, diversity, and evolutionary dynamics of animal cell types remain almost completely unexplored beyond a few species. Single-cell genomics is emerging as a powerful tool to build comprehensive catalogs of cell types and their associated gene regulatory programs in non-traditional model species. We review the current state of sampling efforts across the animal tree of life and challenges ahead for the comparative study of cell type programs. We also discuss how the phylogenetic integration of cell atlases can lead to the development of models of cell type evolution and a phylogenetic taxonomy of cells.

Resolving Transcriptional States and Predicting Lineages in the Annelid Capitella teleta Using Single-Cell RNAseq

Evolution and diversification of cell types has contributed to animal evolution. However, gene regulatory mechanisms underlying cell fate acquisition during development remains largely uncharacterized in spiralians. Here we use a whole-organism, single-cell transcriptomic approach to map larval cell types in the annelid Capitella teleta at 24- and 48-h post gastrulation (stages 4 and 5). We identified eight unique cell clusters (undifferentiated precursors, ectoderm, muscle, ciliary-band, gut, neurons, neurosecretory cells, and protonephridia), thus helping to identify uncharacterized molecular signatures such as previously unknown neurosecretory cell markers in C. teleta. Analysis of coregulatory programs in individual clusters revealed gene interactions that can be used for comparisons of cell types across taxa. We examined the neural and neurosecretory clusters more deeply and characterized a differentiation trajectory starting from dividing precursors to neurons using Monocle3 and velocyto. Pseudotime analysis along this trajectory identified temporally-distinct cell states undergoing progressive gene expression changes over time. Our data revealed two potentially distinct neural differentiation trajectories including an early trajectory for brain neurosecretory cells. This work provides a valuable resource for future functional investigations to better understanding neurogenesis and the transitions from neural precursors to neurons in an annelid.

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.