Chromosomes of Japanese starfishes

Population genomic analysis reveals a polygenic sex determination system in Apostichopus japonicus

The sea cucumber Apostichopus japonicus, a key species in Chinese aquaculture, plays a significant evolutionary role within the Echinodermata phylum. However, the sex determination mechanism in this species remains poorly understood. Here, we conducted extensive sex surveys and sampling of eight wild populations, investigating the sex-related SNPs and insertion or deletions (indels) through bulk segregation analysis (BSA) and genome-wide association study (GWAS) analysis. Our findings suggest that A. japonicus employs a polygenic sex determination (PSD) system, with solute carrier family 8 (SLC8A) being the candidate gene for sex determination, encoding sodium-calcium exchanger (NCX1). The analysis of normalized sequencing depth reveals variations across chromosomes 6, 13, 14, 16, and 18, supporting the PSD system. We also identified 541.656 kb of male-specific sequences and screened five markers (C77185, C98086, C64977, C125, and C876) for molecular sex identification. Overall, this study provides new insights into A. japonicus sex determination, highlighting a complex multi-gene mechanism rather than a simple XX/XY system.

A reference genome for ecological restoration of the sunflower sea star, Pycnopodia helianthoides

Wildlife diseases, such as the sea star wasting (SSW) epizootic that outbroke in the mid-2010s, appear to be associated with acute and/or chronic abiotic environmental change; dissociating the effects of different drivers can be difficult. The sunflower sea star, Pycnopodia helianthoides, was the species most severely impacted during the SSW outbreak, which overlapped with periods of anomalous atmospheric and oceanographic conditions, and there is not yet a consensus on the cause(s). Genomic data may reveal underlying molecular signatures that implicate a subset of factors and, thus, clarify past events while also setting the scene for effective restoration efforts. To advance this goal, we used Pacific Biosciences HiFi long sequencing reads and Dovetail Omni-C proximity reads to generate a highly contiguous genome assembly that was then annotated using RNA-seq-informed gene prediction. The genome assembly is 484 Mb long, with contig N50 of 1.9 Mb, scaffold N50 of 21.8 Mb, BUSCO completeness score of 96.1%, and 22 major scaffolds consistent with prior evidence that sea star genomes comprise 22 autosomes. These statistics generally fall between those of other recently assembled chromosome-scale assemblies for two species in the distantly related asteroid genus Pisaster. These novel genomic resources for P. helianthoides will underwrite population genomic, comparative genomic, and phylogenomic analyses-as well as their integration across scales-of SSW and environmental stressors.

Chromosome-level genome assembly of the northern Pacific seastar Asterias amurensis

Asterias amurensis has attracted widespread concern because of its population outbreaks, which has impacted fisheries and aquaculture, as well as disrupting local ecosystems. A high-quality reference genome is necessary to better investigate mechanisms of outbreak and adaptive changes. Combining PacBio HiFi and Hi-C sequencing data, we generated a chromosome-level A. amurensis genome with a size of 491.53 Mb. The contig N50 and scaffold N50 were 8.05 and 23.75 Mb, respectively. The result of BUSCO analysis revealed a completeness score of 98.85%. A total of 16,531 protein-coding genes were predicted in the genome, of which 94.63% were functionally annotated. The high-quality genome assembly resulting from this study will provide a valuable genetic resource for future research on the mechanism of population outbreaks and invasion ecology.

The first high-quality genome assembly and annotation of Patiria pectinifera

The blue bat star, a highly adaptive species in the East Sea of Korea, has displayed remarkable success in adapting to recent climate change. The genetic mechanisms behind this success were not well-understood, prompting our report on the first chromosome-level assembly of the Patiria genus. We assembled the genome using Nanopore and Illumina sequences, yielding a total length of 615 Mb and a scaffold N50 of 24,204,423 bp. Hi-C analysis allowed us to anchor the scaffold sequences onto 22 pseudochromosomes. K-mer based analysis revealed 5.16% heterozygosity rate of the genome, higher than any previously reported echinoderm species. Our transposable element analysis exposed a substantial number of genome-wide retrotransposons and DNA transposons. These results offer valuable resources for understanding the evolutionary mechanisms behind P. pectinifera's successful adaptation in fluctuating environments.

A chromosome-level genome assembly of a deep-sea starfish (Zoroaster cf. ophiactis)

Understanding of adaptation and evolution of organisms in the deep sea requires more genomic resources. Zoroaster cf. ophiactis is a sea star in the family Zoroasteridae occurring exclusively in the deep sea. In this study, a chromosome-level genome assembly for Z. cf. ophiactis was generated by combining Nanopore long-read, Illumina short-read, and Hi-C sequencing data. The final assembly was 1,002.0 Mb in length, with a contig N50 of 376 Kb and a scaffold N50 of 40.4 Mb, and included 22 pseudo-chromosomes, covering 92.3% of the assembly. Completeness analysis evaluated with BUSCO revealed that 95.91% of the metazoan conserved genes were complete. Additionally, 39,426 protein-coding genes were annotated for this assembly. This chromosome-level genome assembly represents the first high-quality genome for the deep-sea Asteroidea, and will provide a valuable resource for future studies on evolution and adaptation of deep-sea echinoderms.

A chromosome-level reference genome for the giant pink sea star, Pisaster brevispinus, a species severely impacted by wasting

Efforts to protect the ecologically and economically significant California Current Ecosystem from global change will greatly benefit from data about patterns of local adaptation and population connectivity. To facilitate that work, we present a reference-quality genome for the giant pink sea star, Pisaster brevispinus, a species of ecological importance along the Pacific west coast of North America that has been heavily impacted by environmental change and disease. We used Pacific Biosciences HiFi long sequencing reads and Dovetail Omni-C proximity reads to generate a highly contiguous genome assembly of 550 Mb in length. The assembly contains 127 scaffolds with a contig N50 of 4.6 Mb and a scaffold N50 of 21.4 Mb; the BUSCO completeness score is 98.70%. The P. brevispinus genome assembly is comparable to the genome of the congener species P. ochraceus in size and completeness. Both Pisaster assemblies are consistent with previously published karyotyping results showing sea star genomes are organized into 22 autosomes. The reference genome for P. brevispinus is an important first step toward the goal of producing a comprehensive, population genomics view of ecological and evolutionary processes along the California coast. This resource will help scientists, managers, and policy makers in their task of understanding and protecting critical coastal regions from the impacts of global change.

Characterization of Coelomic Fluid Cell Types in the Starfish Marthasterias glacialis Using a Flow Cytometry/Imaging Combined Approach

Coelomocytes is the generic name for a collection of cellular morphotypes, present in many coelomate animals, and highly variable among echinoderm classes. The roles attributed to the major types of these free circulating cells present in the coelomic fluid of echinoderms include immune response, phagocytic digestion and clotting. Our main aim in this study was to characterize coelomocytes found in the coelomic fluid of Marthasterias glacialis (class Asteroidea) by using a combination of flow cytometry (FC), imaging flow cytometry (IFC) and fluorescence plus transmission electron microscopy (TEM). Two coelomocyte populations (P1 and P2) identified through flow cytometry were subsequently studied in terms of abundance, morphology, ultrastructure, cell viability and cell cycle profiles. Ultrastructurally, P2 diploid cells were present as two main morphotypes, similar to phagocytes and vertebrate thrombocytes, whereas the smaller P1 cellular population was characterized by low mitotic activity, a relatively undifferentiated cytotype and a high nucleus/cytoplasm ratio. In the present study we could not rule out possible similarities between haploid P1 cells and stem-cell types in other animals. Additionally, we report the presence of two other morphotypes in P2 that could only be detected by fluorescence microscopy, as well as a morphotype revealed via combined microscopy/FC. This integrative experimental workflow combined cells physical separation with different microscopic image capture technologies, enabling us to better tackle the characterization of the heterogeneous composition of coelomocytes populations.

An initial comparative genomic autopsy of wasting disease in sea stars

Beginning in 2013, sea stars throughout the Eastern North Pacific were decimated by wasting disease, also known as "asteroid idiopathic wasting syndrome" (AIWS) due to its elusive aetiology. The geographic extent and taxonomic scale of AIWS meant events leading up to the outbreak were heterogeneous, multifaceted, and oftentimes unobserved; progression from morbidity to death was rapid, leaving few tell-tale symptoms. Here, we take a forensic genomic approach to discover candidate genes that may help explain sea star wasting syndrome. We report the first genome and annotation for Pisaster ochraceus, along with differential gene expression (DGE) analyses in four size classes, three tissue types, and in symptomatic and asymptomatic individuals. We integrate nucleotide polymorphisms associated with survivors of the wasting disease outbreak, DGE associated with temperature treatments in P. ochraceus, and DGE associated with wasting in another asteroid Pycnopodia helianthoides. In P. ochraceus, we found DGE across all tissues, among size classes, and between asymptomatic and symptomatic individuals; the strongest wasting-associated DGE signal was in pyloric caecum. We also found previously identified outlier loci co-occur with differentially expressed genes. In cross-species comparisons of symptomatic and asymptomatic individuals, consistent responses distinguish genes associated with invertebrate innate immunity and chemical defence, consistent with context-dependent stress responses, defensive apoptosis, and tissue degradation. Our analyses thus highlight genomic constituents that may link suspected environmental drivers (elevated temperature) with intrinsic differences among individuals (age/size, alleles associated with susceptibility) that elicit organismal responses (e.g., coelomocyte proliferation) and manifest as sea star wasting mass mortality.

Current Status of Echinoderm Genome Analysis - What do we Know?

Echinoderms have long served as model organisms for a variety of biological research, especially in the field of developmental biology. Although the genome of the purple sea urchin Strongylocentrotus purpuratus has been sequenced, it is the only echinoderm whose whole genome sequence has been reported. Nevertheless, data is rapidly accumulating on the chromosomes and genomic sequences of all five classes of echinoderms, including the mitochondrial genomes and Hox genes. This blossoming new data will be essential for estimating the phylogenetic relationships among echinoderms, and also to examine the underlying mechanisms by which the diverse morphologies of echinoderms have arisen.