Seasonal tissue-specific gene expression in wild crown-of-thorns starfish reveals reproductive and stress-related transcriptional systems.

Animals are influenced by the season, yet we know little about the changes that occur in most species throughout the year. This is particularly true in tropical marine animals that experience relatively small annual temperature and daylight changes. Like many coral reef inhabitants, the crown-of-thorns starfish (COTS), well known as a notorious consumer of corals and destroyer of coral reefs, reproduces exclusively in the summer. By comparing gene expression in 7 somatic tissues procured from wild COTS sampled on the Great Barrier Reef, we identified more than 2,000 protein-coding genes that change significantly between summer and winter. COTS genes that appear to mediate conspecific communication, including both signalling factors released into the surrounding sea water and cell surface receptors, are up-regulated in external secretory and sensory tissues in the summer, often in a sex-specific manner. Sexually dimorphic gene expression appears to be underpinned by sex- and season-specific transcription factors (TFs) and gene regulatory programs. There are over 100 TFs that are seasonally expressed, 87% of which are significantly up-regulated in the summer. Six nuclear receptors are up-regulated in all tissues in the summer, suggesting that systemic seasonal changes are hormonally controlled, as in vertebrates. Unexpectedly, there is a suite of stress-related chaperone proteins and TFs, including HIFa, ATF3, C/EBP, CREB, and NF-κB, that are uniquely and widely co-expressed in gravid females. The up-regulation of these stress proteins in the summer suggests the demands of oogenesis in this highly fecund starfish affects protein stability and turnover in somatic cells. Together, these circannual changes in gene expression provide novel insights into seasonal changes in this coral reef pest and have the potential to identify vulnerabilities for targeted biocontrol.

A sea star is only a head.

Where are the front and back ends in a sea star? Formery et al. recently tackled this long-standing mystery using state-of-the-art molecular tools, leading them to suggest that a sea star may be constructed from components that, in other animals, would constitute only the head.

High germline mutation rates, but not extreme population outbreaks, influence genetic diversity in a keystone coral predator.

Lewontin's paradox, the observation that levels of genetic diversity (π) do not scale linearly with census population size (Nc) variation, is an evolutionary conundrum. The most extreme mismatches between π and Nc are found for highly abundant marine invertebrates. Yet, the influences of new mutations on π relative to extrinsic processes such as Nc fluctuations are unknown. Here, we provide the first germline mutation rate (µ) estimate for a marine invertebrate in corallivorous crown-of-thorns sea stars (Acanthaster cf. solaris). We use high-coverage whole-genome sequencing of 14 parent-offspring trios alongside empirical estimates of Nc in Australia's Great Barrier Reef to jointly examine the determinants of π in populations undergoing extreme Nc fluctuations. The A. cf. solaris mean µ was 9.13 x 10-09 mutations per-site per-generation (95% CI: 6.51 x 10-09 to 1.18 x 10-08), exceeding estimates for other invertebrates and showing greater concordance with vertebrate mutation rates. Lower-than-expected Ne (~70,000-180,000) and low Ne/Nc values (0.0047-0.048) indicated weak influences of population outbreaks on long-term π. Our findings are consistent with elevated µ evolving in response to reduced Ne and generation time length, with important implications for explaining high mutational loads and the determinants of genetic diversity in marine invertebrate taxa.

Development and Testing of Species-Specific Primers for Detecting the Presence of the Northern Pacific Sea Star (Asterias amurensis) from Environmental DNA.

The starfish Asterias amurensis, a well-known predator of molluscan species in intertidal ecosystems, has caused substantial ecological and economic losses in North China such as offshore Qingdao. Effective monitoring and prevention measures are urged to minimize its negative impacts. Compared with traditional biomonitoring methods, environmental DNA technology has emerged as a powerful and cost-efficient tool for inferring species' presence and abundance. In this study, we developed a pair of species-specific primers (i.e., Ast-F and Ast-R) for the A. amurensis mitochondrial COI gene and tested its utility in amplifying and quantifying the DNA fragments from environmental samples under both laboratory and field conditions. The results of controlled water tank experiments demonstrated that the amount of eDNA released by A. amurensis was positively related to its biomass; after the removal of the starfish, the eDNA degraded significantly in 24 h and remained detectable for 8 days. The number of eDNA copies enriched tended to increase with smaller pore size of filter membrane and larger volume of filtered water. For field tests, we confirmed the validation of our approach in six locations in Qingdao by filtering 1000 ml water per sample with a 0.45-µm pore size filtration. All the amplification products generated a single and bright band via gel electrophoresis, and the quantitative PCR results unveiled significant differences in eDNA copies. This study provided an eDNA-based approach for investigating the distribution and biomass of A. amurensis, which may help to formulate early warning and management strategies in coastal Qingdao and other regions.

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.

Systematics of deep-sea starfish order Brisingida (Echinodermata: Asteroidea), with a revised classification and assessments of morphological characters.

Brisingida Fisher 1928 is one of the seven currently recognised starfish orders, and one of the least known taxa as being exclusive deep-sea inhabitants. Modern deep-sea expeditions revealed their common occurrences in various deep-sea settings including seamounts, basins and hydrothermal vent peripheral, underlining the necessity of clarifying their global diversity and phylogeny. In this study, we present a comprehensive molecular phylogeny of Brisingida which encompasses the highest taxonomic diversity to date. DNA sequences (COI, 16S, 12S and 28S) were obtained from 225 specimens collected in the global ocean, identified as 58 species spanning 15 of the 17 extant genera. Phylogenetic relationship was inferred using both maximum likelihood and Bayesian inference methods, revealing polyphyletic families and genera and indicating nonnegligible bias in prior morphology-based systematics. Based on the new phylogeny, a novel classification of the order, consisting of 5 families and 17 genera, is proposed. Families Odinellidae, Brisingasteridae and Novodiniidae (sensu Clark and Mah, 2001) were resurrected to encompass the genera Odinella, Brisingaster and Novodinia. Brisingidae and Freyellidae were revised to include 11 and 3 genera, respectively. A new genus and species, two new subgenera and seven new combinations are described and a key to each genus and family is provided. Transformations of morphological traits were evaluated under the present phylogenetic hypothesis. A series of paedomorphic characters were found in many genera and species, which led to a high degree of homoplasy across phylogenetically distant groups. Our results provide new insights in the phylogeny and ontogeny of the order, and highlight the necessity to evaluate character convergence under sound phylogenetic hypothesis.

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.

New hypotheses of cell type diversity and novelty from orthology-driven comparative single cell and nuclei transcriptomics in echinoderms.

Cell types are the building blocks of metazoan biodiversity and offer a powerful perspective for inferring evolutionary phenomena. With the development of single-cell transcriptomic techniques, new definitions of cell types are emerging. This allows a conceptual reassessment of traditional definitions of novel cell types and their evolution. Research in echinoderms, particularly sea star and sea urchin embryos has contributed significantly to understanding the evolution of novel cell types, through the examination of skeletogenic mesenchyme and pigment cells, which are found in sea urchin larvae, but not sea star larvae. This paper outlines the development of a gene expression atlas for the bat sea star, Patiria miniata, using single nuclear RNA sequencing (snRNA-seq) of embryonic stages. The atlas revealed 23 cell clusters covering all expected cell types from the endoderm, mesoderm, and ectoderm germ layers. In particular, four distinct neural clusters, an immune-like cluster, and distinct right and left coelom clusters were revealed as distinct cell states. A comparison with Strongylocentrotus purpuratus embryo single-cell transcriptomes was performed using 1:1 orthologs to anchor and then compare gene expression patterns. The equivalent of S. purpuratus piwil3+ Cells were not detected in P. miniata, while the Left Coelom of P. miniata has no equivalent cell cluster in S. purpuratus. These differences may reflect changes in developmental timing between these species. While considered novel morphologically, the Pigment Cells of S. purpuratus map to clusters containing Immune-like Mesenchyme and Neural cells of P. miniata, while the Skeletogenic Mesenchyme of S. purpuratus are revealed as orthologous to the Right Coelom cluster of P. miniata. These results suggest a new interpretation of the evolution of these well-studied cell types and a reflection on the definition of novel cell types.

Detection of TALEN-mediated genome cleavage during the early embryonic stage of the starfish Patiria pectinifera.

BACKGROUND: Echinoderms have long been utilized as experimental materials to study the genetic control of developmental processes and their evolution. Among echinoderms, the molecular study of starfish embryos has received considerable attention across research topics such as gene regulatory network evolution and larval regeneration. Recently, experimental techniques to manipulate gene functions have been gradually established in starfish as the feasibility of genome editing methods was reported. However, it is still unclear when these techniques cause genome cleavage during the development of starfish, which is critical to understand the timeframe and applicability of the experiment during early development of starfish. RESULTS: We herein reported that gene functions can be analyzed by the genome editing method TALEN in early embryos, such as the blastula of the starfish Patiria pectinifera. We injected the mRNA of TALEN targeting rar, which was previously constructed, into eggs of P. pectinifera and examined the efficiency of genome cleavage through developmental stages from 6 to 48 hours post fertilization. CONCLUSION: The results will be key knowledge not only when designing TALEN-based experiments but also when assessing the results.

Coelomic fluid of asteroid echinoderms: Current knowledge and future perspectives on its utility for disease and mortality investigations.

Coelomic fluid surrounds the internal organs of asteroid echinoderms (asteroids, otherwise known as sea stars or starfish) and plays an essential role in the immune system, as well as in the transport of respiratory gases, nutrients, waste products, and reproductive mediators. Due to its importance in physiology and accessibility for nonlethal diagnostic sampling, coelomic fluid of asteroids provides an excellent sample matrix for health evaluations and can be particularly useful in disease and mortality investigations. This is especially important in light of recent increases in the number of affected individuals and species, larger geographic scope, and increased observed frequency of sea star wasting events compared with historic accounts of wasting. This review summarizes the current knowledge about coelomocytes, the effector cell of the asteroid immune system; coelomic fluid electrolytes, osmolality, acid-base status and respiratory gases, and microbiota; and genomic, transcriptomic, and proteomic investigations of coelomic fluid. The utility of coelomic fluid analysis for assessing stressor responses, diseases, and mortality investigations is considered with knowledge gaps and future directions identified. This complex body fluid provides an exciting opportunity to increase our understanding of this unique and ecologically important group of animals.

Molecular mechanisms of tubulogenesis revealed in the sea star hydro-vascular organ.

A fundamental goal in the organogenesis field is to understand how cells organize into tubular shapes. Toward this aim, we have established the hydro-vascular organ in the sea star Patiria miniata as a model for tubulogenesis. In this animal, bilateral tubes grow out from the tip of the developing gut, and precisely extend to specific sites in the larva. This growth involves cell migration coupled with mitosis in distinct zones. Cell proliferation requires FGF signaling, whereas the three-dimensional orientation of the organ depends on Wnt signaling. Specification and maintenance of tube cell fate requires Delta/Notch signaling. Moreover, we identify target genes of the FGF pathway that contribute to tube morphology, revealing molecular mechanisms for tube outgrowth. Finally, we report that FGF activates the Six1/2 transcription factor, which serves as an evolutionarily ancient regulator of branching morphogenesis. This study uncovers distinct mechanisms of tubulogenesis in vivo and we propose that cellular dynamics in the sea star hydro-vascular organ represents a key comparison for understanding the evolution of vertebrate organs.

Captivity induces a sweeping and sustained genomic response in a starfish.

Marine animals in the wild are often difficult to access, so they are studied in captivity. However, the implicit assumption that physiological processes of animals in artificial environments are not different from those in the wild has rarely been tested. Here, we investigate the extent to which an animal is impacted by captivity by comparing global gene expression in wild and captive crown-of-thorns starfish (COTS). In a preliminary analysis, we compared transcriptomes of three external tissues obtained from multiple wild COTS with a single captive COTS maintained in aquaria for at least 1 week. On average, an astonishingly large 24% of the coding sequences in the genome were differentially expressed. This led us to conduct a replicated experiment to test more comprehensively the impact of captivity on gene expression. Specifically, a comparison of 13 wild with 8 captive COTS coelomocyte transcriptomes revealed significant differences in the expression of 20% of coding sequences. Coelomocyte transcriptomes in captive COTS remain different from those in wild COTS for more than 30 days and show no indication of reverting back to a wild state (i.e. no evidence of acclimation). Genes upregulated in captivity include those involved in oxidative stress and energy metabolism, whereas genes downregulated are involved in cell signalling. These changes in gene expression indicate that being translocated and maintained in captivity has a marked impact on the physiology and health of these echinoderms. This study suggests that caution should be exercised when extrapolating results from captive aquatic invertebrates to their wild counterparts.

Ultrasensitive and on-site eDNA detection for the monitoring of crown-of-thorns starfish densities at the pre-outbreak stage using an electrochemical biosensor.

The coral reef crisis has significantly intensified over the last decades, mainly due to severe outbreaks of crown-of-thorns starfish (COTS). Current ecological monitoring has failed to detect COTS densities at the pre-outbreak stage, thus preventing early intervention. In this work, we developed an effective electrochemical biosensor modified by a MoO2/C nanomaterial, as well as a specific DNA probe that could detect trace COTS environmental DNA (eDNA) at a lower detection limit (LOD = 0.147 ng/µL) with excellent specificity. The reliability and accuracy of the biosensor were validated against the standard methods by an ultramicro spectrophotometer and droplet digital PCR (p > 0.05). The biosensor was then utilized for the on-site analysis of seawater samples from SYM-LD and SY sites in the South China Sea. For the SYM-LD site suffering an outbreak, the COTS eDNA concentrations were 0.33 ng/µL (1 m, depth) and 0.26 ng/µL (10 m, depth), respectively. According to the ecological survey, the COTS density was 500 ind/hm2 at the SYM-LD site, verifying the accuracy of our measurements. At the SY site, COTS eDNA was also detected at 0.19 ng/µL, but COTS was not found by the traditional survey. Hence, larvae were possibly present in this region. Therefore, this electrochemical biosensor could be used to monitor COTS populations at the pre-outbreak stages, and potentially serve as a revolutionary early warning method. We will continue to improve this method for picomolar or even femtomolar detection of COTS eDNA.

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.

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.

Sex-specific expression of pheromones and other signals in gravid starfish.

BACKGROUND: Many echinoderms form seasonal aggregations prior to spawning. In some fecund species, a spawning event can lead to population outbreaks with detrimental ecosystem impacts. For instance, outbreaks of crown-of-thorns starfish (COTS), a corallivore, can destroy coral reefs. Here, we examine the gene expression in gravid male and female COTS prior to spawning in the wild, to identify genome-encoded factors that may regulate aggregation and spawning. This study is informed by a previously identified exoproteome that attracts conspecifics. To capture the natural gene expression profiles, we isolated RNAs from gravid female and male COTS immediately after they were removed from the Great Barrier Reef.  RESULTS: Sexually dimorphic gene expression is present in all seven somatic tissues and organs that we surveyed and in the gonads. Approximately 40% of the exoproteome transcripts are differentially expressed between sexes. Males uniquely upregulate an additional 68 secreted factors in their testes. A suite of neuropeptides in sensory organs, coelomocytes and gonads is differentially expressed between sexes, including the relaxin-like gonad-stimulating peptide and gonadotropin-releasing hormones. Female sensory tentacles-chemosensory organs at the distal tips of the starfish arms-uniquely upregulate diverse receptors and signalling molecules, including chemosensory G-protein-coupled receptors and several neuropeptides, including kisspeptin, SALMFamide and orexin. CONCLUSIONS: Analysis of 103 tissue/organ transcriptomes from 13 wild COTS has revealed genes that are consistently differentially expressed between gravid females and males and that all tissues surveyed are sexually dimorphic at the molecular level. This finding is consistent with female and male COTS using sex-specific pheromones to regulate reproductive aggregations and synchronised spawning events. These pheromones appear to be received primarily by the sensory tentacles, which express a range of receptors and signalling molecules in a sex-specific manner. Furthermore, coelomocytes and gonads differentially express signalling and regulatory factors that control gametogenesis and spawning in other echinoderms.

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.

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.

Developing an effective marine eDNA monitoring: eDNA detection at pre-outbreak densities of corallivorous seastar (Acanthaster cf. solaris).

Outbreaks of the corallivorous Crown-of-Thorns Seastar (CoTS) Acanthaster cf. solaris contribute significantly to coral reef loss. Control of outbreaks is hampered because standard monitoring techniques do not detect outbreaks at early (low density) stages, thus preventing early intervention. We previously demonstrated that eDNA monitoring can detect CoTS at intermediate densities. Here, we test whether detection probability can be improved by (i) targeted site selection or collection at specific times and (ii) moving from an average eDNA copy number approach (based on the limit of quantification) to a presence/absence approach (based on the limit of detection). Using a dataset collected over three years and multiple reef sites, we demonstrated that adding water residence age, sea surface level and temperature into generalized linear models explained low amounts of variance of eDNA copy numbers. Site specific CoTS density, by contrast, was a significant predictor for eDNA copy numbers. Bayesian multi-scale occupancy modelling of the presence/absence data demonstrated that the probability of sample capture (θ) on most reefs with intermediate or high CoTS densities was >0.8. Thus, confirming CoTS presence on these reefs would only require 2-3 samples. Sample capture decreased with decreasing CoTS density. Collecting ten filters was sufficient to reliably (based on the lower 95 % Credibility Interval) detect CoTS below nominal outbreak levels (3 Ind. ha-1). Copy number-based estimates may be more relevant to quantify CoTS at higher densities. Although water residence age did contribute little to our models, sites with higher residence times may serve as sentinel sites accumulating eDNA. The approach based on presence or absence of eDNA facilitates eDNA monitoring to detect CoTS densities below outbreak thresholds and we continue to further develop this method for quantification.

Chromosome-level genome assembly of Plazaster borealis sheds light on the morphogenesis of multiarmed starfish and its regenerative capacity.

BACKGROUND: Plazaster borealis has a unique morphology, displaying multiple arms with a clear distinction between disk and arms, rather than displaying pentaradial symmetry, a remarkable characteristic of echinoderms. Herein we report the first chromosome-level reference genome of P. borealis and an essential tool to further investigate the basis of the divergent morphology. FINDINGS: In total, 57.76 Gb of a long read and 70.83 Gb of short-read data were generated to assemble a de novo 561-Mb reference genome of P. borealis, and Hi-C sequencing data (57.47 Gb) were used for scaffolding into 22 chromosomal scaffolds comprising 92.38% of the genome. The genome completeness estimated by BUSCO was 98.0% using the metazoan set, indicating a high-quality assembly. Through the comparative genome analysis, we identified evolutionary accelerated genes known to be involved in morphogenesis and regeneration, suggesting their potential role in shaping body pattern and capacity of regeneration. CONCLUSION: This first chromosome-level genome assembly of P. borealis provides fundamental insights into echinoderm biology, as well as the genomic mechanism underlying its unique morphology and regeneration.