The crown-of-thorns starfish genome as a guide for biocontrol of this coral reef pest

The brittle star genome illuminates the genetic basis of animal appendage regeneration

Species within nearly all extant animal lineages are capable of regenerating body parts. However, it remains unclear whether the gene expression programme controlling regeneration is evolutionarily conserved. Brittle stars are a species-rich class of echinoderms with outstanding regenerative abilities, but investigations into the genetic bases of regeneration in this group have been hindered by the limited genomic resources. Here we report a chromosome-scale genome assembly for the brittle star Amphiura filiformis. We show that the brittle star genome is the most rearranged among echinoderms sequenced so far, featuring a reorganized Hox cluster reminiscent of the rearrangements observed in sea urchins. In addition, we performed an extensive profiling of gene expression during brittle star adult arm regeneration and identified sequential waves of gene expression governing wound healing, proliferation and differentiation. We conducted comparative transcriptomic analyses with other invertebrate and vertebrate models for appendage regeneration and uncovered hundreds of genes with conserved expression dynamics, particularly during the proliferative phase of regeneration. Our findings emphasize the crucial importance of echinoderms to detect long-range expression conservation between vertebrates and classical invertebrate regeneration model systems.

A new strategy based on a cascade amplification strategy biosensor for on-site eDNA detection and outbreak warning of crown-of-thorns starfish

Population outbreaks of the crown-of-thorns starfish (COTS) seriously threaten the sustainability of coral reef ecosystems. However, traditional ecological monitoring techniques cannot provide early warning before the outbreaks, thus preventing timely intervention. Therefore, there is an urgent need for a more accurate and faster technology to predict the outbreaks of COTS. In this work, we developed an electrochemical biosensor based on a programmed catalytic hairpin assembly (CHA) and hybridization chain reaction (HCR) cyclic amplification strategy for sensitive and selective detection of COTS environmental DNA (eDNA) in water bodies. This biosensor exhibited excellent electrochemical characteristics, including a low limit of detection (LOD = 18.4 fM), low limit of quantification (LOQ = 41.1 fM), and wide linear range (50 fM - 10 nM). The biosensing technology successfully allowed the detection of COTS eDNA in the aquarium environment, and the results also demonstrated a significant correlation between eDNA concentration and COTS number (r = 0.990; P < 0.001). The reliability and accuracy of the biosensor results have been further validated through comparison with digital droplet PCR (ddPCR). Moreover, the applicability and accuracy of the biosensor were reconfirmed in field tests at the COTS outbreak site in the South China Sea, which has shown potential application in dynamically monitoring the larvae before the COTS outbreak. Therefore, this efficient electrochemical biosensing technology offers a new solution for on-site monitoring and early warning of the COTS outbreak.

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.

Chemosensory behaviour of juvenile crown-of-thorns sea star (Acanthaster sp.), attraction to algal and coral food and avoidance of adult conspecifics

Intraspecific and habitat-mediated responses to chemical cues play key roles in structuring populations of marine species. We investigated the behaviour of herbivorous-stage juvenile crown-of-thorns sea stars (COTS; Acanthaster sp.) in flow-through choice chambers to determine if chemical cues from their habitat influence movement and their transition to become coral predators. Juveniles at the diet transition stage were exposed to cues from their nursery habitat (coral rubble-crustose coralline algae (CCA)), live coral and adult COTS to determine if waterborne cues influence movement. In response to CCA and coral as sole cues, juveniles moved towards the cue source and when these cues were presented in combination, they exhibited a preference for coral. Juveniles moved away from adult COTS cues. Exposure to food cues (coral, CCA) in the presence of adult cues resulted in variable responses. Our results suggest a feedback mechanism whereby juvenile behaviour is mediated by adult chemical cues. Cues from the adult population may deter juveniles from the switch to corallivory. As outbreaks wane, juveniles released from competition may serve as a proximate source of outbreaks, supporting the juveniles-in-waiting hypothesis. The accumulation of juveniles within the reef infrastructure is an underappreciated potential source of COTS outbreaks that devastate coral reefs.

Protecting Great Barrier Reef resilience through effective management of crown-of-thorns starfish outbreaks

Resilience-based management is essential to protect ecosystems in the Anthropocene. Unlike large-scale climate threats to Great Barrier Reef (GBR) corals, outbreaks of coral-eating crown-of-thorns starfish (COTS; Acanthaster cf. solaris) can be directly managed through targeted culling. Here, we evaluate the outcomes of a decade of strategic COTS management in suppressing outbreaks and protecting corals during the 4th COTS outbreak wave at reef and regional scales (sectors). We compare COTS density and coral cover dynamics during the 3rd and 4th outbreak waves. During the 4th outbreak wave, sectors that received limited to no culling had sustained COTS outbreaks causing significant coral losses. In contrast, in sectors that received timely and sufficient cull effort, coral cover increased substantially, and outbreaks were suppressed with COTS densities up to six-fold lower than in the 3rd outbreak wave. In the Townsville sector for example, despite exposure to comparable disturbance regimes during the 4th outbreak wave, effective outbreak suppression coincided with relative increases in sector-wide coral cover (44%), versus significant coral cover declines (37%) during the 3rd outbreak wave. Importantly, these estimated increases span entire sectors, not just reefs with active COTS control. Outbreaking reefs with higher levels of culling had net increases in coral cover, while the rate of coral loss was more than halved on reefs with lower levels of cull effort. Our results also indicate that outbreak wave progression to adjoining sectors has been delayed, probably via suppression of COTS larval supply. Our findings provide compelling evidence that proactive, targeted, and sustained COTS management can effectively suppress COTS outbreaks and deliver coral growth and recovery benefits at reef and sector-wide scales. The clear coral protection outcomes demonstrate the value of targeted manual culling as both a scalable intervention to mitigate COTS outbreaks, and a potent resilience-based management tool to "buy time" for coral reefs, protecting reef ecosystem functions and biodiversity as the climate changes.

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.

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.

Induction of larval settlement in crown-of-thorns starfish is not mediated by conspecific cues

Population irruptions of crown-of-thorns starfish (COTS; Acanthaster spp.) remain a major cause of coral reef degradation throughout the Pacific and Indian Oceans and are inherently modulated by larval settlement and recruitment success. Gregarious larval settlement, as exhibited by many other ecologically important marine invertebrates, can catalyse population growth and replenishment. However, whether conspecific cues induce or influence the settlement of COTS larvae remains a critical information gap. This experimental study examined the induction of COTS settlement in response to a range of conspecific cues associated with early- and late-stage herbivorous juveniles, corallivorous juveniles and adults. Competent COTS larvae were generally not induced to settle by the presence of conspecifics or cues associated with conspecifics, while the settlement success of COTS in the presence of coralline algae was not inhibited or enhanced by adding conspecific conditioned seawater. Rather than being reinforced by gregarious settlement, the recruitment of COTS populations appears dependent on associative settlement cues (i.e., coralline algae and/or associated microbial communities) signalling suitable benthic habitat.

Crown-of-thorns starfish spines secrete defence proteins

Background

The crown-of-thorns starfish (COTS; Acanthaster species) is a slow-moving corallivore protected by an extensive array of long, sharp toxic spines. Envenomation can result in nausea, numbness, vomiting, joint aches and sometimes paralysis. Small molecule saponins and the plancitoxin proteins have been implicated in COTS toxicity.

Methods

Brine shrimp lethality assays were used to confirm the secretion of spine toxin biomolecules. Histological analysis, followed by spine-derived proteomics helped to explain the source and identity of proteins, while quantitative RNA-sequencing and phylogeny confirmed target gene expression and relative conservation, respectively.

Results

We demonstrate the lethality of COTS spine secreted biomolecules on brine shrimp, including significant toxicity using aboral spine semi-purifications of >10 kDa (p > 0.05, 9.82 µg/ml), supporting the presence of secreted proteins as toxins. Ultrastructure observations of the COTS aboral spine showed the presence of pores that could facilitate the distribution of secreted proteins. Subsequent purification and mass spectrometry analysis of spine-derived proteins identified numerous secretory proteins, including plancitoxins, as well as those with relatively high gene expression in spines, including phospholipase A2, protease inhibitor 16-like protein, ependymin-related proteins and those uncharacterized. Some secretory proteins (e.g., vitellogenin and deleted in malignant brain tumor protein 1) were not highly expressed in spine tissue, yet the spine may serve as a storage or release site. This study contributes to our understanding of the COTS through functional, ultrastructural and proteomic analysis of aboral spines.

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.

Genome-wide identification and comparative analysis of Dmrt genes in echinoderms

The Dmrt (Doublesex-mab3-related transcription factor) gene family is a class of crucial transcription factors characterized by one or several conserved DM (Doublesex/Mab-3) domains. Dmrt family genes can participate in various physiological developmental processes, especially in sex determination/differentiation. Echinoderms are extremely important research objects in various fields, such as sex determination/differentiation and neuroscience. However, to date, the genome-wide characterization and analysis of Dmrt genes in echinoderms have not been investigated. In this study, the identification and analysis of Dmrt genes in 11 representative echinoderms were performed using bioinformatics methods. A total of 43 Dmrt genes have been found in the studied echinoderms, and the number of Dmrt genes in different species ranges from 2 to 5. The phylogenetic tree showed that all Dmrt genes from echinoderms can be subdivided into 5 classes, the Dmrt2-like class, Dmrt3-like class, Dmrt4/5-like class, Dsx-like class, and a novel Dmrt (starfish-specific) class. Furthermore, selective pressure assessment suggested that the Dmrt genes underwent purifying selection pressure. In general, this study provides a molecular basis for echinoderm Dmrt genes and may serve as a reference for in-depth phylogenomics.

The Holothuria leucospilota genome elucidates sacrificial organ expulsion and bioadhesive trap enriched with amyloid-patterned proteins

Some tropical sea cucumbers of the family Holothuriidae can efficiently repel or even fatally ensnare predators by sacrificially ejecting a bioadhesive matrix termed the Cuvierian organ (CO), so named by the French zoologist Georges Cuvier who first described it in 1831. Still, the precise mechanisms for how adhesiveness genetically arose in CO and how sea cucumbers perceive and transduce danger signals for CO expulsion during defense have remained unclear. Here, we report the first high-quality, chromosome-level genome assembly of Holothuria leucospilota, an ecologically significant sea cucumber with prototypical CO. The H. leucospilota genome reveals characteristic long-repeat signatures in CO-specific outer-layer proteins, analogous to fibrous proteins of disparate species origins, including spider spidroin and silkworm fibroin. Intriguingly, several CO-specific proteins occur with amyloid-like patterns featuring extensive intramolecular cross-β structures readily stainable by amyloid indicator dyes. Distinct proteins within the CO connective tissue and outer surface cooperate to give the expelled matrix its apparent tenacity and adhesiveness, respectively. Genomic evidence offers further hints that H. leucospilota directly transduces predator-induced mechanical pressure onto the CO surface through mediation by transient receptor potential channels, which culminates in acetylcholine-triggered CO expulsion in part or in entirety. Evolutionarily, innovative events in two distinct regions of the H. leucospilota genome have apparently spurred CO's differentiation from the respiratory tree to a lethal defensive organ against predators.

Molecular evolution of the ependymin-related gene epdl2 in African weakly electric fish

Gene duplication and subsequent molecular evolution can give rise to taxon-specific gene specializations. In previous work, we found evidence that African weakly electric fish (Mormyridae) may have as many as three copies of the epdl2 gene, and the expression of two epdl2 genes is correlated with electric signal divergence. Epdl2 belongs to the ependymin-related family (EPDR), a functionally diverse family of secretory glycoproteins. In this study, we first describe vertebrate EPDR evolution and then present a detailed evolutionary history of epdl2 in Mormyridae with emphasis on the speciose genus Paramormyrops. Using Sanger sequencing, we confirm three apparently functional epdl2 genes in Paramormyrops kingsleyae. Next, we developed a nanopore-based amplicon sequencing strategy and bioinformatics pipeline to obtain and classify full-length epdl2 gene sequences (N = 34) across Mormyridae. Our phylogenetic analysis proposes three or four epdl2 paralogs dating from early Paramormyrops evolution. Finally, we conducted selection tests which detected positive selection around the duplication events and identified ten sites likely targeted by selection in the resulting paralogs. These sites' locations in our modeled 3D protein structure involve four sites in ligand binding and six sites in homodimer formation. Together, these findings strongly imply an evolutionary mechanism whereby epdl2 genes underwent selection-driven functional specialization after tandem duplications in the rapidly speciating Paramormyrops. Considering previous evidence, we propose that epdl2 may contribute to electric signal diversification in mormyrids, an important aspect of species recognition during mating.

Structure and proteomic analysis of the crown-of-thorns starfish (Acanthaster sp.) radial nerve cord

The nervous system of the Asteroidea (starfish or seastar) consists of radial nerve cords (RNCs) that interconnect with a ring nerve. Despite its relative simplicity, it facilitates the movement of multiple arms and numerous tube feet, as well as regeneration of damaged limbs. Here, we investigated the RNC ultrastructure and its molecular components within the of Pacific crown-of-thorns starfish (COTS; Acanthaster sp.), a well-known coral predator that in high-density outbreaks has major ecological impacts on coral reefs. We describe the presence of an array of unique small bulbous bulbs (40-100 μm diameter) that project from the ectoneural region of the adult RNC. Each comprise large secretory-like cells and prominent cilia. In contrast, juvenile COTS and its congener Acanthaster brevispinus lack these features, both of which are non-corallivorous. Proteomic analysis of the RNC (and isolated neural bulbs) provides the first comprehensive echinoderm protein database for neural tissue, including numerous secreted proteins associated with signalling, transport and defence. The neural bulbs contained several neuropeptides (e.g., bombyxin-type, starfish myorelaxant peptide, secretogranin 7B2-like, Ap15a-like, and ApNp35) and Deleted in Malignant Brain Tumor 1-like proteins. In summary, this study provides a new insight into the novel traits of COTS, a major pest on coral reefs, and a proteomics resource that can be used to develop (bio)control strategies and understand molecular mechanisms of regeneration.

Biochemical metabolomic profiling of the Crown-of-Thorns Starfish (Acanthaster): New insight into its biology for improved pest management

The Crown-of-Thorns Starfish (COTS), Acanthaster species, is a voracious coral predator that destroys coral reefs when in outbreak status. The baseline metabolite and lipid biomolecules of 10 COTS tissues, including eggs from gravid females, were investigated in this study to provide insight into their biology and identify avenues for control. Targeted and untargeted metabolite- and lipidomics-based mass spectrometry approaches were used to obtain tissue-specific metabolite and lipid profiles. Across all COTS tissues, 410 metabolites and 367 lipids were identified. Most abundant were amino acids and peptides (18.7%) and wax esters (17%). There were 262 metabolites and 192 lipids identified in COTS eggs. Wax esters were more abundant in the eggs (30%) followed by triacylglycerols (TG), amino acids, and peptides. The diversity of asterosaponins in eggs (34) was higher than in tissues (2). Several asterosaponins known to modulate sperm acrosome reaction were putatively identified, including glycoside B, asterosaponin-4 (Co-Aris III), and regularoside B (asterosaponin A). The saponins saponin A, thornasteroside A, hillaside B, and non-saponins dictyol J and axinellamine B which have been shown to possess defensive properties, were found in abundance in gonads, skin, and radial nerve tissues. Inosine and 2-hexyldecanoic acid are the most abundant metabolites in tissues and eggs. As a secondary metabolite of purine degradation, inosine plays an important role in purine biosynthesis, while 2-hexyldecanoic acid is known to suppress side-chain crystallization during the synthesis of amphiphilic macromolecules (i.e., phospholipids). These significant spatial changes in metabolite, lipid, and asterosaponin profiles enabled unique insights into key biological tissue-specific processes that could be manipulated to better control COTS populations. Our findings highlight COTS as a novel source of molecules with therapeutic and cosmetic properties (ceramides, sphingolipids, carnosine, and inosine). These outcomes will be highly relevant for the development of strategies for COTS management including chemotaxis-based biocontrol and exploitation of COTS carcasses for the extraction of commercial molecules.

A Proteomic Analysis for the Red Seaweed Asparagopsis taxiformis

The red seaweed Asparagopsis taxiformis is a promising ruminant feed additive with anti-methanogenic properties that could contribute to global climate change solutions. Genomics has provided a strong foundation for in-depth molecular investigations, including proteomics. Here, we investigated the proteome of A. taxiformis (Lineage 6) in both sporophyte and gametophyte stages, using soluble and insoluble extraction methods. We identified 741 unique non-redundant proteins using a genome-derived database and 2007 using a transcriptome-derived database, which included numerous proteins predicted to be of fungal origin. We further investigated the genome-derived proteins to focus on seaweed-specific proteins. Ontology analysis indicated a relatively large proportion of ion-binding proteins (i.e., iron, zinc, manganese, potassium and copper), which may play a role in seaweed heavy metal tolerance. In addition, we identified 58 stress-related proteins (e.g., heat shock and vanadium-dependent haloperoxidases) and 44 photosynthesis-related proteins (e.g., phycobilisomes, photosystem I, photosystem II and ATPase), which were in general more abundantly identified from female gametophytes. Forty proteins were predicted to be secreted, including ten rhodophyte collagen-alpha-like proteins (RCAPs), which displayed overall high gene expression levels. These findings provide a comprehensive overview of expressed proteins in A. taxiformis, highlighting the potential for targeted protein extraction and functional characterisation for future biodiscovery.

The genome of a hadal sea cucumber reveals novel adaptive strategies to deep-sea environments

How organisms cope with coldness and high pressure in the hadal zone remains poorly understood. Here, we sequenced and assembled the genome of hadal sea cucumber Paelopatides sp. Yap with high quality and explored its potential mechanisms for deep-sea adaptation. First, the expansion of ACOX1 for rate-limiting enzyme in the DHA synthesis pathway, increased DHA content in the phospholipid bilayer, and positive selection of EPT1 may maintain cell membrane fluidity. Second, three genes for translation initiation factors and two for ribosomal proteins underwent expansion, and three ribosomal protein genes were positively selected, which may ameliorate the protein synthesis inhibition or ribosome dissociation in the hadal zone. Third, expansion and positive selection of genes associated with stalled replication fork recovery and DNA repair suggest improvements in DNA protection. This is the first genome sequence of a hadal invertebrate. Our results provide insights into the genetic adaptations used by invertebrate in deep oceans.

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.

Genome Survey Sequencing and Genetic Background Characterization of Ilex chinensis Sims (Aquifoliaceae) Based on Next-Generation Sequencing

Ilex chinensis Sims. is an evergreen arbor species with high ornamental and medicinal value that is widely distributed in China. However, there is a lack of molecular and genomic data for this plant, which severely restricts the development of its relevant research. To obtain the whole reference genome, we first conducted a genome survey of I. chinensis by next-generation sequencing (NGS) to perform de novo whole-genome sequencing. As a result, our estimates using k-mer and flow cytometric analysis suggested the genome size of I. chinensis to be around 618-655 Mb, with the GC content, heterozygous rate, and repeat sequence rate of 37.52%, 1.1%, and 38%, respectively. A total of 334,649 microsatellite motifs were detected from the I. chinensis genome data, which will provide basic molecular markers for germplasm characterization, genetic diversity, and QTL mapping studies for I. chinensis. In summary, the I. chinensis genome is complex with high heterozygosity and few repeated sequences. Overall, this is the first report on the genome features of I. chinensis, and the information may lay a strong groundwork for future whole-genome sequencing and molecular breeding studies of this species.

The Impact of Phenotypic Heterogeneity on Chemotactic Self-Organisation

The capacity to aggregate through chemosensitive movement forms a paradigm of self-organisation, with examples spanning cellular and animal systems. A basic mechanism assumes a phenotypically homogeneous population that secretes its own attractant, with the well known system introduced more than five decades ago by Keller and Segel proving resolutely popular in modelling studies. The typical assumption of population phenotypic homogeneity, however, often lies at odds with the heterogeneity of natural systems, where populations may comprise distinct phenotypes that vary according to their chemotactic ability, attractant secretion, etc. To initiate an understanding into how this diversity can impact on autoaggregation, we propose a simple extension to the classical Keller and Segel model, in which the population is divided into two distinct phenotypes: those performing chemotaxis and those producing attractant. Using a combination of linear stability analysis and numerical simulations, we demonstrate that switching between these phenotypic states alters the capacity of a population to self-aggregate. Further, we show that switching based on the local environment (population density or chemoattractant level) leads to diverse patterning and provides a route through which a population can effectively curb the size and density of an aggregate. We discuss the results in the context of real world examples of chemotactic aggregation, as well as theoretical aspects of the model such as global existence and blow-up of solutions.

Twinkle twinkle brittle star: the draft genome of Ophioderma brevispinum (Echinodermata: Ophiuroidea) as a resource for regeneration research

Background

Echinoderms are established models in experimental and developmental biology, however genomic resources are still lacking for many species. Here, we present the draft genome of Ophioderma brevispinum, an emerging model organism in the field of regenerative biology. This new genomic resource provides a reference for experimental studies of regenerative mechanisms.

Results

We report a de novo nuclear genome assembly for the brittle star O. brevispinum and annotation facilitated by the transcriptome assembly. The final assembly is 2.68 Gb in length and contains 146,703 predicted protein-coding gene models. We also report a mitochondrial genome for this species, which is 15,831 bp in length, and contains 13 protein-coding, 22 tRNAs, and 2 rRNAs genes, respectively. In addition, 29 genes of the Notch signaling pathway are identified to illustrate the practical utility of the assembly for studies of regeneration.

Conclusions

The sequenced and annotated genome of O. brevispinum presented here provides the first such resource for an ophiuroid model species. Considering the remarkable regenerative capacity of this species, this genome will be an essential resource in future research efforts on molecular mechanisms regulating regeneration.

Supplementary Information

The online version contains supplementary material available at (10.1186/s12864-022-08750-y).

New Three-Finger Protein from Starfish Asteria rubens Shares Structure and Pharmacology with Human Brain Neuromodulator Lynx2

Three-finger proteins (TFPs) are small proteins with characteristic three-finger β-structural fold stabilized by the system of conserved disulfide bonds. These proteins have been found in organisms from different taxonomic groups and perform various important regulatory functions or act as components of snake venoms. Recently, four TFPs (Lystars 1–4) with unknown function were identified in the coelomic fluid proteome of starfish A. rubens. Here we analyzed the genomes of A. rubens and A. planci starfishes and predicted additional five and six proteins containing three-finger domains, respectively. One of them, named Lystar5, is expressed in A. rubens coelomocytes and has sequence homology to the human brain neuromodulator Lynx2. The three-finger structure of Lystar5 close to the structure of Lynx2 was confirmed by NMR. Similar to Lynx2, Lystar5 negatively modulated α4β2 nicotinic acetylcholine receptors (nAChRs) expressed in X. laevis oocytes. Incubation with Lystar5 decreased the expression of acetylcholine esterase and α4 and α7 nAChR subunits in the hippocampal neurons. In summary, for the first time we reported modulator of the cholinergic system in starfish.

A Review of Asteroid Biology in the Context of Sea Star Wasting: Possible Causes and Consequences

AbstractSea star wasting-marked in a variety of sea star species as varying degrees of skin lesions followed by disintegration-recently caused one of the largest marine die-offs ever recorded on the west coast of North America, killing billions of sea stars. Despite the important ramifications this mortality had for coastal benthic ecosystems, such as increased abundance of prey, little is known about the causes of the disease or the mechanisms of its progression. Although there have been studies indicating a range of causal mechanisms, including viruses and environmental effects, the broad spatial and depth range of affected populations leaves many questions remaining about either infectious or non-infectious mechanisms. Wasting appears to start with degradation of mutable connective tissue in the body wall, leading to disintegration of the epidermis. Here, we briefly review basic sea star biology in the context of sea star wasting and present our current knowledge and hypotheses related to the symptoms, the microbiome, the viruses, and the associated environmental stressors. We also highlight throughout the article knowledge gaps and the data needed to better understand sea star wasting mechanistically, its causes, and potential management.

Biosynthesis of saponin defensive compounds in sea cucumbers

Soft-bodied slow-moving sea creatures such as sea stars and sea cucumbers lack an adaptive immune system and have instead evolved the ability to make specialized protective chemicals (glycosylated steroids and triterpenes) as part of their innate immune system. This raises the intriguing question of how these biosynthetic pathways have evolved. Sea star saponins are steroidal, while those of the sea cucumber are triterpenoid. Sterol biosynthesis in animals involves cyclization of 2,3-oxidosqualene to lanosterol by the oxidosqualene cyclase (OSC) enzyme lanosterol synthase (LSS). Here we show that sea cucumbers lack LSS and instead have two divergent OSCs that produce triterpene saponins and that are likely to have evolved from an ancestral LSS by gene duplication and neofunctionalization. We further show that sea cucumbers make alternate sterols that confer protection against self-poisoning by their own saponins. Collectively, these events have enabled sea cucumbers to evolve the ability to produce saponins and saponin-resistant sterols concomitantly.

Sea stars and sea cucumbers biosynthesize protective glycosylated steroids and triterpenes via divergent oxidosqualene cyclases (OSCs) that produce these distinct saponins in different species as well as in different tissues of a single species.

Proteomic changes associated with predator-induced morphological defenses in oysters

Inducible prey defenses occur when organisms undergo plastic changes in phenotype to reduce predation risk. When predation pressure varies persistently over space or time, such as when predator and prey co-occur over only part of their biogeographic ranges, prey populations can become locally adapted in their inducible defenses. In California estuaries, native Olympia oyster (Ostrea lurida) populations have evolved disparate phenotypic responses to an invasive predator, the Atlantic oyster drill (Urosalpinx cinerea). In this study, oysters from an estuary with drills, and oysters from an estuary without drills, were reared for two generations in a laboratory common garden, and subsequently exposed to cues from Atlantic drills. Comparative proteomics was then used to investigate molecular mechanisms underlying conserved and divergent aspects of their inducible defenses. Both populations developed smaller, thicker, and harder shells after drill exposure, and these changes in shell phenotype were associated with up-regulation of calcium transport proteins that could influence biomineralization. Inducible defenses evolve in part because defended phenotypes incur fitness costs when predation risk is low. Immune proteins were down-regulated by both oyster populations after exposure to drills, implying a trade-off between biomineralization and immune function. Following drill exposure, oysters from the population that co-occurs with drills grew smaller shells than oysters inhabiting the estuary not yet invaded by the predator. Variation in the response to drills between populations was associated with isoform-specific protein expression. This trend suggests that a stronger inducible defense response evolved in oysters that co-occur with drills through modification of an existing mechanism.

Metazoans and Intrinsic Apoptosis: An Evolutionary Analysis of the Bcl-2 Family

The B-cell lymphoma-2 (Bcl-2) family is a group of genes regulating intrinsic apoptosis, a process controlling events such as development, homeostasis and the innate and adaptive immune responses in metazoans. In higher organisms, Bcl-2 proteins coordinate intrinsic apoptosis through their regulation of the integrity of the mitochondrial outer membrane; this function appears to have originated in the basal metazoans. Bcl-2 genes predate the cnidarian-bilaterian split and have been identified in porifera, placozoans and cnidarians but not ctenophores and some nematodes. The Bcl-2 family is composed of two groups of proteins, one with an α-helical Bcl-2 fold that has been identified in porifera, placozoans, cnidarians, and almost all higher bilaterians. The second group of proteins, the BH3-only group, has little sequence conservation and less well-defined structures and is found in cnidarians and most bilaterians, but not porifera or placozoans. Here we examine the evolutionary relationships between Bcl-2 proteins. We show that the structures of the Bcl-2-fold proteins are highly conserved over evolutionary time. Some metazoans such as the urochordate Oikopleura dioica have lost all Bcl-2 family members. This gene loss indicates that Bcl-2 regulated apoptosis is not an absolute requirement in metazoans, a finding mirrored in recent gene deletion studies in mice. Sequence analysis suggests that at least some Bcl-2 proteins lack the ability to bind BH3-only antagonists and therefore potentially have other non-apoptotic functions. By examining the foundations of the Bcl-2 regulated apoptosis, functional relationships may be clarified that allow us to understand the role of specific Bcl-2 proteins in evolution and disease.

The genome of an apodid holothuroid (Chiridota heheva) provides insights into its adaptation to a deep-sea reducing environment

Cold seeps and hydrothermal vents are deep-sea reducing environments that are characterized by lacking oxygen and photosynthesis-derived nutrients. Most animals acquire nutrition in cold seeps or hydrothermal vents by maintaining epi- or endosymbiotic relationship with chemoautotrophic microorganisms. Although several seep- and vent-dwelling animals hosting symbiotic microbes have been well-studied, the genomic basis of adaptation to deep-sea reducing environment in nonsymbiotic animals is still lacking. Here, we report a high-quality genome of Chiridota heheva Pawson & Vance, 2004, which thrives by extracting organic components from sediment detritus and suspended material, as a reference for nonsymbiotic animal's adaptation to deep-sea reducing environments. The expansion of the aerolysin-like protein family in C. heheva compared with other echinoderms might be involved in the disintegration of microbes during digestion. Moreover, several hypoxia-related genes (Pyruvate Kinase M2, PKM2; Phospholysine Phosphohistidine Inorganic Pyrophosphate Phosphatase, LHPP; Poly(A)-specific Ribonuclease Subunit PAN2, PAN2; and Ribosomal RNA Processing 9, RRP9) were subject to positive selection in the genome of C. heheva, which contributes to their adaptation to hypoxic environments.

Crown of thorns starfish life-history traits contribute to outbreaks, a continuing concern for coral reefs

Crown of thorns starfish (COTS, Acanthaster sp.) are notorious for their destructive consumption of coral that decimates tropical reefs, an attribute unique among tropical marine invertebrates. Their populations can rapidly increase from 0–1 COTS ha⁻¹ to more than 10–1000 COTS ha⁻¹ in short order causing a drastic change to benthic communities and reducing the functional and species diversity of coral reef ecosystems. Population outbreaks were first identified to be a significant threat to coral reefs in the 1960s. Since then, they have become one of the leading causes of coral loss along with coral bleaching. Decades of research and significant investment in Australia and elsewhere, particularly Japan, have been directed towards identifying, understanding, and managing the potential causes of outbreaks and designing population control methods. Despite this, the drivers of outbreaks remain elusive. What is becoming increasingly clear is that the success of COTS is tied to their inherent biological traits, especially in early life. Survival of larval and juvenile COTS is likely to be enhanced by their dietary flexibility and resilience to variable food conditions as well as their phenotypically plastic growth dynamics, all magnified by the extreme reproductive potential of COTS. These traits enable COTS to capitalise on anthropogenic disturbances to reef systems as well as endure less favourable conditions.

Somatostatin-type and allatostatin-C-type neuropeptides are paralogous and have opposing myoregulatory roles in an echinoderm

Somatostatin (SS) and allatostatin-C (ASTC) are related neuropeptide hormones that act as inhibitory regulators of physiological processes in chordates (e.g., humans) and protostome invertebrates (e.g., insects), respectively. We have discovered that echinoderms (e.g., starfish) uniquely have both SS-type and ASTC-type neuropeptides, which act as inhibitory and excitatory regulators of muscle activity, respectively. Our findings suggest that SS-type and ASTC-type neuropeptides evolved by duplication of a common ancestral encoding gene. Then, one of the neuropeptides was lost in protostomes and chordates, probably because of their functional redundancy as inhibitory regulators. Conversely, the unique retention of both neuropeptide types in echinoderms may be explained by evolution of an excitatory role for ASTC-type neuropeptides mediated by yet-to-be-determined signaling mechanisms.

Somatostatin (SS) and allatostatin-C (ASTC) are inhibitory neuropeptides in chordates and protostomes, respectively, which hitherto were identified as orthologs. However, echinoderms have two SS/ASTC-type neuropeptides (SS1 and SS2), and here, our analysis of sequence data indicates that SS1 is an ortholog of ASTC and SS2 is an ortholog of SS. The occurrence of both SS-type and ASTC-type neuropeptides in echinoderms provides a unique context to compare their physiological roles. Investigation of the expression and actions of the ASTC-type neuropeptide ArSS1 in the starfish Asterias rubens revealed that it causes muscle contraction (myoexcitation), contrasting with myoinhibitory effects of the SS-type neuropeptide ArSS2. Our findings suggest that SS-type and ASTC-type neuropeptides are paralogous and originated by gene duplication in a common ancestor of the Bilateria, with only one type being retained in chordates (SS) and protostomes (ASTC) but with both types being retained in echinoderms. Loss of ASTC-type and SS-type neuropeptides in chordates and protostomes, respectively, may have been due to their functional redundancy as inhibitory regulators of physiological processes. Conversely, the retention of both neuropeptide types in echinoderms may be a consequence of the evolution of a myoexcitatory role for ASTC-type neuropeptides mediated by as yet unknown signaling mechanisms.

Knowledge Gaps in the Biology, Ecology, and Management of the Pacific Crown-of-Thorns Sea Star Acanthaster sp. on Australia's Great Barrier Reef

AbstractCrown-of-thorns sea stars (Acanthaster sp.) are among the most studied coral reef organisms, owing to their propensity to undergo major population irruptions, which contribute to significant coral loss and reef degradation throughout the Indo-Pacific. However, there are still important knowledge gaps pertaining to the biology, ecology, and management of Acanthaster sp. Renewed efforts to advance understanding and management of Pacific crown-of-thorns sea stars (Acanthaster sp.) on Australia's Great Barrier Reef require explicit consideration of relevant and tractable knowledge gaps. Drawing on established horizon scanning methodologies, this study identified contemporary knowledge gaps by asking active and/or established crown-of-thorns sea star researchers to pose critical research questions that they believe should be addressed to improve the understanding and management of crown-of-thorns sea stars on the Great Barrier Reef. A total of 38 participants proposed 246 independent research questions, organized into 7 themes: feeding ecology, demography, distribution and abundance, predation, settlement, management, and environmental change. Questions were further assigned to 48 specific topics nested within the 7 themes. During this process, redundant questions were removed, which reduced the total number of distinct research questions to 172. Research questions posed were mostly related to themes of demography (46 questions) and management (48 questions). The dominant topics, meanwhile, were the incidence of population irruptions (16 questions), feeding ecology of larval sea stars (15 questions), effects of elevated water temperature on crown-of-thorns sea stars (13 questions), and predation on juveniles (12 questions). While the breadth of questions suggests that there is considerable research needed to improve understanding and management of crown-of-thorns sea stars on the Great Barrier Reef, the predominance of certain themes and topics suggests a major focus for new research while also providing a roadmap to guide future research efforts.

Cloning and Selfing Affect Population Genetic Variation in Simulations of Outcrossing, Sexual Sea Stars

AbstractMany sea stars are well known for facultative or obligate asexual reproduction in both the adult and larval life-cycle stages. Some species and lineages are also capable of facultative or obligate hermaphroditic reproduction with self-fertilization. However, models of population genetic variation and empirical analyses of genetic data typically assume only sexual reproduction and outcrossing. A recent reanalysis of previously published empirical data (microsatellite genotypes) from two studies of one of the most well-known sea star species (the crown-of-thorns sea star; Acanthaster sp.) concluded that cloning and self-fertilization in that species are rare and contribute little to patterns of population genetic variation. Here we reconsider that conclusion by simulating the contribution of cloning and selfing to genetic variation in a series of models of sea star demography. Simulated variation in two simple models (analogous to previous analyses of empirical data) was consistent with high rates of cloning or selfing or both. More realistic scenarios that characterize population flux in sea stars of ecological significance, including outbreaks of crown-of-thorns sea stars that devastate coral reefs, invasions by Asterias amurensis, and epizootics of sea star wasting disease that kill Pisaster ochraceus, also showed significant but smaller effects of cloning and selfing on variation within subpopulations and differentiation between subpopulations. Future models or analyses of genetic variation in similar study systems might benefit from simulation modeling to characterize possible contributions of cloning or selfing to genetic variation in population samples or to understand the limits on inferring the effects of cloning or selfing in nature.

Echidnas of the Sea: The Defensive Behavior of Juvenile and Adult Crown-of-Thorns Sea Stars

AbstractCrown-of-thorns sea stars are one of the most ecologically important tropical marine invertebrates, with boom-bust population dynamics that influence the community structure of coral reefs. Although predation is likely to influence the development of population outbreaks, little is known about the defensive behavior of crown-of-thorns sea stars. Righting behavior after being overturned, a key defensive response in echinoderms, was investigated for the newly settled herbivorous juvenile, the corallivorous juvenile, and adult stages of crown-of-thorns sea stars. The average righting time of the newly settled juveniles (0.3-1.0-mm diameter) was 2.74 minutes. For the coral-eating juveniles (15-55-mm diameter), the righting time (mean = 6.24 min) was faster in larger juveniles, and the mean righting time of the adults was 6.28 minutes. During righting and in response to being lifted off of the substrate, the juveniles and adults exhibited an arm curling response, during which their arms closed over their oral side, often forming a spine ball, a feature not known for other asteroids. The righting and curling responses of the corallivorous juveniles were influenced by the presence of a natural enemy, a coral guard crab, which caused the juveniles to spend more time with their arms curled. These behaviors indicate that crown-of-thorns sea stars use their spines to protect the soft tissue of their oral side. The highly defended morphology and behavioral adaptations of crown-of-thorns sea stars are likely to have evolved as antipredator mechanisms. This points to the potential importance of predators in regulating their populations, which may have decreased in recent times due to fishing, a factor that may contribute to outbreaks.

The impact of rheotaxis and flow on the aggregation of organisms

Dispersed populations often need to organize into groups. Chemical attractants provide one means of directing individuals into an aggregate, but whether these structures emerge can depend on various factors, such as there being a sufficiently large population or the response to the attractant being sufficiently sensitive. In an aquatic environment, fluid flow may heavily impact on population distribution and many aquatic organisms adopt a rheotaxis response when exposed to a current, orienting and swimming according to the flow field. Consequently, flow-induced transport could be substantially different for the population members and any aggregating signal they secrete. With the aim of investigating how flows and rheotaxis responses impact on an aquatic population's ability to form and maintain an aggregated profile, we develop and analyse a mathematical model that incorporates these factors. Through a systematic analysis into the effect of introducing rheotactic behaviour under various forms of environmental flow, we demonstrate that each of flow and rheotaxis can act beneficially or detrimentally on the ability to form and maintain a cluster. Synthesizing these findings, we test a hypothesis that density-dependent rheotaxis may be optimal for group formation and maintenance, in which individuals increase their rheotactic effort as they approach an aggregated state.

A new subfamily of ionotropic glutamate receptors unique to the echinoderms with putative sensory role

Chemosensation is a critical signalling process in animals and especially important in sea cucumbers, a group of ecologically and economically important marine echinoderms (class Holothuroidea), which lack audio and visual organs and rely on chemical sensing for survival, feeding and reproduction. The ionotropic receptors are a recently identified family of chemosensory receptors in insects and other protostomes, related to the ionotropic glutamate receptor family (iGluR), a large family of membrane receptors in metazoan. Here we characterize the echinoderm iGluR subunits and consider their possible role in chemical communication in sea cucumbers. Sequence similarity searches revealed that sea cucumbers have in general a higher number of iGluR subunits when compared to other echinoderms. Phylogenetic analysis and sequence comparisons revealed GluH as a specific iGluR subfamily present in all echinoderms. Homologues of the vertebrate GluA (aka α-amino-3-hydroxy-5-methyl-4-isoxazolepropionic acid, AMPA), GluK (aka kainate) and GluD (aka delta) were also identified. The GluN (aka N-methyl-d-aspartate, NMDA) as well as the invertebrate deuterostome subfamily GluF (aka phi) are absent in echinoderms. The echinoderm GluH subfamily shares conserved structural protein organization with vertebrate iGluRs and the ligand binding domain (LBD) is the most conserved region; genome analysis indicates evolution via lineage and species-specific tandem gene duplications. GluH genes (named Grih) are the most highly expressed iGluRs subunit genes in tissues in the sea cucumber Holothuria arguinesis, with Griha1, Griha2 and Griha5 exclusively expressed in tentacles, making them candidates to have a chemosensory role in this species. The multiple GluH subunits may provide alternative receptor assembly combinations, thus expanding the functional possibilities and widening the range of compounds detected during aggregation and spawning in echinoderms.

A highly-contiguous genome assembly of the Eurasian spruce bark beetle, Ips typographus, provides insight into a major forest pest

Conifer-feeding bark beetles are important herbivores and decomposers in forest ecosystems. These species complete their life cycle in nutritionally poor substrates and some can kill enormous numbers of trees during population outbreaks. The Eurasian spruce bark beetle (Ips typographus) can destroy >100 million m3 of spruce in a single year. We report a 236.8 Mb I. typographus genome assembly using PacBio long-read sequencing. The final phased assembly has a contig N50 of 6.65 Mb in 272 contigs and is predicted to contain 23,923 protein-coding genes. We reveal expanded gene families associated with plant cell wall degradation, including pectinases, aspartyl proteases, and glycosyl hydrolases. This genome sequence from the genus Ips provides timely resources to address questions about the evolutionary biology of the true weevils (Curculionidae), one of the most species-rich animal families. In forests of today, increasingly stressed by global warming, this draft genome may assist in developing pest control strategies to mitigate outbreaks.

The Development and Neuronal Complexity of Bipinnaria Larvae of the Sea Star Asterias rubens

Free-swimming planktonic larvae are a key stage in the development of many marine phyla, and studies of these organisms have contributed to our understanding of major genetic and evolutionary processes. Although transitory, these larvae often attain a remarkable degree of tissue complexity, with well-defined musculature and nervous systems. Among the best studied are larvae belonging to the phylum Echinodermata, but with work largely focused on the pluteus larvae of sea urchins (class Echinoidea). The greatest diversity of larval strategies among echinoderms is found in the class Asteroidea (sea stars), organisms that are rapidly emerging as experimental systems for genetic and developmental studies. However, the bipinnaria larvae of sea stars have only been studied in detail in a small number of species and although they have been relatively well described neuro-anatomically, they are poorly understood neurochemically. Here, we have analyzed embryonic development and bipinnaria larval anatomy in the common North Atlantic sea star Asterias rubens, using a variety of staining methods in combination with confocal microscopy. Importantly, the chemical complexity of the nervous system of bipinnaria larvae was revealed through use of a diverse set of antibodies, with identification of at least three centers of differing neurochemical signature within the previously described nervous system: the anterior apical organ, oral region, and ciliary bands. Furthermore, the anatomy of the musculature and sites of cell division in bipinnaria larvae was analyzed. Comparisons of developmental progression and molecular anatomy across the Echinodermata provided a basis for hypotheses on the shared evolutionary and developmental processes that have shaped this group of animals. We conclude that bipinnaria larvae appear to be remarkably conserved across ∼200 million years of evolutionary time and may represent a strong evolutionary and/or developmental constraint on species utilizing this larval strategy.

Elucidation of the speciation history of three sister species of crown-of-thorns starfish (Acanthaster spp.) based on genomic analysis

The crown-of-thorns starfish (COTS) is a coral predator that is widely distributed in Indo-Pacific Oceans. A previous phylogenetic study using partial mitochondrial sequences suggested that COTS had diverged into four distinct species, but a nuclear genome-based analysis to confirm this was not conducted. To address this, COTS species nuclear genome sequences were analysed here, sequencing Northern Indian Ocean (NIO) and Red Sea (RS) species genomes for the first time, followed by a comparative analysis with the Pacific Ocean (PO) species. Phylogenetic analysis and ADMIXTURE analysis revealed clear divergences between the three COTS species. Furthermore, within the PO species, the phylogenetic position of the Hawaiian sample was further away from the other Pacific-derived samples than expected based on the mitochondrial data, suggesting that it may be a PO subspecies. The pairwise sequentially Markovian coalescent model showed that the trajectories of the population size diverged by region during the Mid-Pleistocene transition when the sea-level was dramatically decreased, strongly suggesting that the three COTS species experienced allopatric speciation. Analysis of the orthologues indicated that there were remarkable genes with species-specific positive selection in the genomes of the PO and RS species, which suggested that there may be local adaptations in the COTS species.

TIAMMAt: Leveraging biodiversity to revise protein domain models, evidence from innate immunity

Sequence annotation is fundamental for studying the evolution of protein families, particularly when working with nonmodel species. Given the rapid, ever-increasing number of species receiving high-quality genome sequencing, accurate domain modeling that is representative of species diversity is crucial for understanding protein family sequence evolution and their inferred function(s). Here, we describe a bioinformatic tool called Taxon-Informed Adjustment of Markov Model Attributes (TIAMMAt) which revises domain profile hidden Markov models (HMMs) by incorporating homologous domain sequences from underrepresented and nonmodel species. Using innate immunity pathways as a case study, we show that revising profile HMM parameters to directly account for variation in homologs among underrepresented species provides valuable insight into the evolution of protein families. Following adjustment by TIAMMAt, domain profile HMMs exhibit changes in their per-site amino acid state emission probabilities and insertion/deletion probabilities while maintaining the overall structure of the consensus sequence. Our results show that domain revision can heavily impact evolutionary interpretations for some families (i.e., NLR’s NACHT domain), whereas impact on other domains (e.g., rel homology domain and interferon regulatory factor domains) is minimal due to high levels of sequence conservation across the sampled phylogenetic depth (i.e., Metazoa). Importantly, TIAMMAt revises target domain models to reflect homologous sequence variation using the taxonomic distribution under consideration by the user. TIAMMAt’s flexibility to revise any subset of the Pfam database using a user-defined taxonomic pool will make it a valuable tool for future protein evolution studies, particularly when incorporating (or focusing) on nonmodel species.

Facilitating population genomics of non-model organisms through optimized experimental design for reduced representation sequencing

BACKGROUND

Genome-wide data are invaluable to characterize differentiation and adaptation of natural populations. Reduced representation sequencing (RRS) subsamples a genome repeatedly across many individuals. However, RRS requires careful optimization and fine-tuning to deliver high marker density while being cost-efficient. The number of genomic fragments created through restriction enzyme digestion and the sequencing library setup must match to achieve sufficient sequencing coverage per locus. Here, we present a workflow based on published information and computational and experimental procedures to investigate and streamline the applicability of RRS.

RESULTS

In an iterative process genome size estimates, restriction enzymes and size selection windows were tested and scaled in six classes of Antarctic animals (Ostracoda, Malacostraca, Bivalvia, Asteroidea, Actinopterygii, Aves). Achieving high marker density would be expensive in amphipods, the malacostracan target taxon, due to the large genome size. We propose alternative approaches such as mitogenome or target capture sequencing for this group. Pilot libraries were sequenced for all other target taxa. Ostracods, bivalves, sea stars, and fish showed overall good coverage and marker numbers for downstream population genomic analyses. In contrast, the bird test library produced low coverage and few polymorphic loci, likely due to degraded DNA.

CONCLUSIONS

Prior testing and optimization are important to identify which groups are amenable for RRS and where alternative methods may currently offer better cost-benefit ratios. The steps outlined here are easy to follow for other non-model taxa with little genomic resources, thus stimulating efficient resource use for the many pressing research questions in molecular ecology.

A nomenclature for echinoderm genes

Echinoderm embryos and larvae are prominent experimental model systems for studying developmental mechanisms. High-quality, assembled, annotated genome sequences are now available for several echinoderm species, including representatives from most classes. The increased availability of these data necessitates the development of a nomenclature that assigns universally interpretable gene symbols to echinoderm genes to facilitate cross-species comparisons of gene functions, both within echinoderms and across other phyla. This paper describes the implementation of an improved set of echinoderm gene nomenclature guidelines that both communicates meaningful orthology information in protein-coding gene symbols and names and establishes continuity with nomenclatures developed for major vertebrate model organisms, including humans. Differences between the echinoderm gene nomenclature guidelines and vertebrate guidelines are examined and explained. This nomenclature incorporates novel solutions to allow for several types of orthologous relationships, including the single echinoderm genes with multiple vertebrate co-orthologs that result from whole-genome-duplication events. The current version of the Echinoderm Gene Nomenclature Guidelines can be found at https://www.echinobase.org/gene/static/geneNomenclature.jsp

Database URL https://www.echinobase.org/

Ambulacrarian insulin-related peptides and their putative receptors suggest how insulin and similar peptides may have evolved from insulin-like growth factor

Background

Some insulin/IGF-related peptides (irps) stimulate a receptor tyrosine kinase (RTK) that transfers the extracellular hormonal signal into an intracellular response. Other irps, such as relaxin, do not use an RTK, but a G-protein coupled receptor (GPCR). This is unusual since evolutionarily related hormones typically either use the same or paralogous receptors. In arthropods three different irps, i.e. arthropod IGF, gonadulin and Drosophila insulin-like peptide 7 (dilp7), likely evolved from a gene triplication, as in several species genes encoding these three peptides are located next to one another on the same chromosomal fragment. These arthropod irps have homologs in vertebrates, suggesting that the initial gene triplication was perhaps already present in the last common ancestor of deuterostomes and protostomes. It would be interesting to know whether this is indeed so and how insulin might be related to this trio of irps.

Methodology

Genes encoding irps as well as their putative receptors were identified in genomes and transcriptomes from echinoderms and hemichordates.

Results

A similar triplet of genes coding for irps also occurs in some ambulacrarians. Two of these are orthologs of arthropod IGF and dilp7 and the third is likely a gonadulin ortholog. In echinoderms, two novel irps emerged, gonad stimulating substance (GSS) and multinsulin, likely from gene duplications of the IGF and dilp7-like genes respectively. The structures of GSS diverged considerably from IGF, which would suggest they use different receptors from IGF, but no novel irp receptors evolved. If IGF and GSS use different receptors, and the evolution of GSS from a gene duplication of IGF is not associated with the appearance of a novel receptor, while irps are known to use two different types of receptors, the ancestor of GSS and IGF might have acted on both types of receptors while one or both of its descendants act on only one. There are three ambulacrarian GPCRs that have amino acid sequences suggestive of being irp GPCRs, two of these are orthologs of the gonadulin and dilp7 receptors. This suggests that the third might be an IGF receptor, and that by deduction, GSS only acts on the RTK. The evolution of GSS from IGF may represent a pattern, where IGF gene duplications lead to novel genes coding for shorter peptides that activate an RTK. It is likely this is how insulin and the insect neuroendocrine irps evolved independently from IGF.

Conclusion

The local gene triplication described from arthropods that yielded three genes encoding irps was already present in the last common ancestor of protostomes and deuterostomes. It seems plausible that irps, such as those produced by neuroendocrine cells in the brain of insects and echinoderm GSS evolved independently from IGF and, thus, are not true orthologs, but the result of convergent evolution.

The Use of Larval Sea Stars and Sea Urchins in the Discovery of Shared Mechanisms of Metazoan Whole-Body Regeneration

The ability to regenerate is scattered among the metazoan tree of life. Further still, regenerative capacity varies widely within these specific organisms. Numerous organisms, all with different regenerative capabilities, have been studied at length and key similarities and disparities in how regeneration occurs have been identified. In order to get a better grasp on understanding regeneration as a whole, we must search for new models that are capable of extensive regeneration, as well as those that have been under sampled in the literature. As invertebrate deuterostomes, echinoderms fit both of these requirements. Multiple members regenerate various tissue types at all life stages, including examples of whole-body regeneration. Interrogations in two highly studied echinoderms, the sea urchin and the sea star, have provided knowledge of tissue and whole-body regeneration at various life stages. Work has begun to examine regeneration in echinoderm larvae, a potential new system for understanding regenerative mechanisms in a basal deuterostome. Here, we review the ways these two animals’ larvae have been utilized as a model of regeneration.

Chemical defense in developmental stages and adult of the sea star Echinaster (Othilia) brasiliensis

To date, evidence regarding the performance of secondary metabolites from larval stages of sea stars as an anti-predation defense relates only to a few species/specimens from a few geographic ranges. Unfortunately, this hinders a comprehensive global understanding of this inter-specific predator-prey interaction. Here, we present laboratory experimental evidence of chemical defense action in the early developmental stages and adults of the sea star Echinaster (Othilia) brasiliensis from Brazil against sympatric and allopatric invertebrate consumers. Blastulae, early and late brachiolarias of E. (O.) brasiliensis were not consumed by the sympatric and allopatric crabs Mithraculus forceps. Blastulae were also avoided by the sympatric and allopatric individuals of the anemone Anemonia sargassensis, but not the larval stages. Extracts from embryos (blastula) and brachiolarias of E. (O.) brasiliensis from one sampled population (João Fernandes beach) significantly inhibited the consumption by sympatric M. forceps, but not by allopatric crabs and A. sargassensi anemone. In this same site, extracts from adults E. (O.) brasiliensis significantly inhibited the consumption by sympatric and allopatric specimens of the crab in a range of concentrations. Whereas equivalent extract concentrations of E. (O.) brasiliensis from other population (Itaipu beach)inhibited the predation by allopatric M. forceps, while sympatric individuals of this crab avoided the only the higher level tested. Then, early stages and adult specimens of E. (O.) brasiliensis can be chemically defended against consumers, but this action is quite variable, depending on the type (anemone or crab) and the origin of the consumer (sympatric or allopatric).

Development and Interrogation of a Transcriptomic Resource for the Giant Triton Snail (Charonia tritonis)

Gastropod molluscs are among the most abundant species that inhabit coral reef ecosystems. Many are specialist predators, along with the giant triton snail Charonia tritonis (Linnaeus, 1758) whose diet consists of Acanthaster planci (crown-of-thorns starfish), a corallivore known to consume enormous quantities of reef-building coral. C. tritonis are considered vulnerable due to overexploitation, and a decline in their populations is believed to have contributed to recurring A. planci population outbreaks. Aquaculture is considered one approach that could help restore natural populations of C. tritonis and mitigate coral loss; however, numerous questions remain unanswered regarding their life cycle, including the molecular factors that regulate their reproduction and development. In this study, we have established a reference C. tritonis transcriptome derived from developmental stages (embryo and veliger) and adult tissues. This was used to identify genes associated with cell signalling, such as neuropeptides and G protein-coupled receptors (GPCRs), involved in endocrine and olfactory signalling. A comparison of developmental stages showed that several neuropeptide precursors are exclusively expressed in post-hatch veligers and functional analysis found that FFamide stimulated a significant (20.3%) increase in larval heart rate. GPCRs unique to veligers, and a diversity of rhodopsin-like GPCRs located within adult cephalic tentacles, all represent candidate olfactory receptors. In addition, the cytochrome P450 superfamily, which participates in the biosynthesis and degradation of steroid hormones and lipids, was also found to be expanded with at least 91 genes annotated, mostly in gill tissue. These findings further progress our understanding of C. tritonis with possible application in developing aquaculture methods.

The protein and volatile components of trail mucus in the Common Garden Snail, Cornu aspersum

The Common or Brown Garden Snail, Cornu aspersum, is an invasive land snail that has successfully colonized a diverse range of global environments. Like other invasive land snails, it is a significant pest of a variety of agricultural crops, including citrus, grapes and canola. Cornu aspersum secretes a mucus trail when mobile that facilitates locomotion. The involvement of the trail in conspecific chemical communication has also been postulated. Our study found that anterior tentacle contact with conspecific mucus elicited a significant increase in heart rate from 46.9 to 51 beats per minute. In order to gain a better understanding of the constituents of the trail mucus and the role it may play in snail communication, the protein and volatile components of mucus trails were investigated. Using two different protein extraction methods, mass spectrometry analysis yielded 175 different proteins, 29 of which had no significant similarity to any entries in the non-redundant protein sequence database. Of the mucus proteins, 22 contain features consistent with secreted proteins, including a perlucin-like protein. The eight most abundant volatiles detected using gas chromatography were recorded (including propanoic acid and limonene) and their potential role as putative pheromones are discussed. In summary, this study has provided an avenue for further research pertaining to the role of trail mucus in snail communication and provides a useful repository for land snail trail mucus components. This may be utilized for further research regarding snail attraction and dispersal, which may be applied in the fields of agriculture, ecology and human health.

Divalent metal transporter-related protein restricts animals to marine habitats

Utilization and regulation of metals from seawater by marine organisms are important physiological processes. To better understand metal regulation, we searched the crown-of-thorns starfish genome for the divalent metal transporter (DMT) gene, a membrane protein responsible for uptake of divalent cations. We found two DMT-like sequences. One is an ortholog of vertebrate DMT, but the other is an unknown protein, which we named DMT-related protein (DMTRP). Functional analysis using a yeast expression system demonstrated that DMT transports various metals, like known DMTs, but DMTRP does not. In contrast, DMTRP reduced the intracellular concentration of some metals, especially zinc, suggesting its involvement in negative regulation of metal uptake. Phylogenetic distribution of the DMTRP gene in various metazoans, including sponges, protostomes, and deuterostomes, indicates that it originated early in metazoan evolution. However, the DMTRP gene is only retained in marine species, and its loss seems to have occurred independently in ecdysozoan and vertebrate lineages from which major freshwater and land animals appeared. DMTRP may be an evolutionary and ecological limitation, restricting organisms that possess it to marine habitats, whereas its loss may have allowed other organisms to invade freshwater and terrestrial habitats.

Mieko Sassa et al. report a novel divalent metal transporter protein (DMTRP) in the crown-of-thorns starfish genome and determine that all organisms with a DMTRP gene are located in marine habitats. They also show in a functional yeast system that DMTRP can prevent uptake of certain metals, bringing insight into the evolution of metal regulation for marine organisms.

Quantitative study of the behavior of two broadcast spawners, the sea urchins Strongylocentrotus intermedius and Mesocentrotus nudus, during mass spawning events in situ

Background

The spatial distribution of spawners and temporal parameters of spawning in motile invertebrates with external fertilization might influence reproductive success. However, to date, data on the prespawning and spawning behaviors of broadcast spawners in the field have been scarce and mostly qualitative. The present study was intended to clarify the behavioral adaptations of two sea urchin species, Strongylocentrotus intermedius and Mesocentrotus nudus, using quantitative analysis of their behavior during mass spawning events under natural conditions.

Methods

We analyzed in situ video recordings of sea urchin behavior obtained during six spawning seasons (2014–2019). The total number of specimens of each sea urchin species and the numbers of spawning males and females were counted. Quantitative parameters of sea urchin spawning (numbers of gamete batches, release duration of one gamete batch, time intervals between gamete batches and total duration of spawning) and movement (step length of spawners and nonspawners before and during spawning and changes in distances between males/nonspawners and females) were determined.

Results

For each species, 12 mass spawning events were recorded in which 10 or more individuals participated. The temporal dynamics of the numbers of males and females participating in mass spawning were well synchronized in both species; however, males began to spawn earlier and ended their spawning later than females. In both species, the most significant intersex difference was the longer spawning duration in males due to the longer pause between gamete batches. The total duration of gamete release did not differ significantly between sexes. The average duration of sperm release during mass spawning events was longer than solitary male spawning. Males and females showed significant increases in the locomotion rate 35 min before the start of spawning and continued to actively move during spawning. An increase in movement rate before spawning in males and females was induced by environmental factor(s). Nonspawners of both species showed increased locomotion activity but in the presence of spawning neighbors and less prominently than spawners. On a vertical surface, both echinoids moved strictly upward. On flat surfaces, males, females and nonspawners of both echinoids became closer during spawning.

Discussion

We showed that two sea urchin species with planktotrophic larvae display similar behavioral adaptations aimed at enhancing reproductive success. The high sensitivity of sea urchins, primarily males, to some environmental factors, most likely phytoplankton, may be considered a large-scale adaptation promoting the development of mass spawning events. The longer spawning duration in males and increased movement activity before and during spawning in both sexes may be considered small-scale adaptations promoting approach of males and females and enhancing the chances of egg fertilization.