Complete mitochondrial genome sequences for Crown-of-thorns starfish Acanthaster planci and Acanthaster brevispinus

Complete Mitochondrial Genomes and Phylogenetic Analysis of Genus Henricia (Asteroidea: Spinulosida: Echinasteridae)

The genus Henricia is known to have intraspecific morphological variations, making species identification difficult. Therefore, molecular phylogeny analysis based on genetic characteristics is valuable for species identification. We present complete mitochondrial genomic sequences of Henricia longispina aleutica, H. reniossa, and H. sanguinolenta for the first time in this study. This study will make a significant contribution to our understanding of Henricia species and its relationships within the class Asteroidea. Lengths of mitochondrial genomes of the three species are 16,217, 16,223, and 16,194 bp, respectively, with a circular form. These genomes contained 13 protein-coding genes, two ribosomal RNA genes, 22 transfer RNA genes, and a D-loop. The gene order and direction aligned with other asteroid species. Phylogenetic relationship analysis showed that our Henricia species were in a monophyletic clade with other Henricia species and in a large clade with species (Echinaster brasiliensis) from the same family. These findings provide valuable insight into understanding the phylogenetic relationships of species in the genus Henricia.

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.

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).

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.

Mitogenomes Reveal Alternative Initiation Codons and Lineage-Specific Gene Order Conservation in Echinoderms

The mitochondrial genetic code is much more varied than the standard genetic code. The invertebrate mitochondrial code, for instance, comprises six initiation codons, including five alternative start codons. However, only two initiation codons are known in the echinoderm and flatworm mitochondrial code, the canonical ATG and alternative GTG. Here, we analyzed 23 Asteroidea mitogenomes, including ten newly sequenced species and unambiguously identified at least two other start codons, ATT and ATC, both of which also initiate translation of mitochondrial genes in other invertebrates. These findings underscore the diversity of the genetic code and expand upon the suite of initiation codons among echinoderms to avoid erroneous annotations. Our analyses have also uncovered the remarkable conservation of gene order among asteroids, echinoids, and holothuroids, with only an interchange between two gene positions in asteroids over ∼500 Ma of echinoderm evolution.

An Investigation into the Genetic History of Japanese Populations of Three Starfish, Acanthaster planci, Linckia laevigata, and Asterias amurensis, Based on Complete Mitochondrial DNA Sequences

Crown-of-thorns starfish, Acanthaster planci (COTS), are common in coral reefs of Indo-Pacific Ocean. Since they are highly fecund predators of corals, periodic outbreaks of COTS cause substantial loss of healthy coral reefs. Using complete mitochondrial DNA sequences, we here examined how COTS outbreaks in the Ryukyu Archipelago, Japan are reflected by the profile of their population genetics. Population genetics of the blue starfish, Linckia laevigata, which lives in the Ryukyu Archipelago, but not break out and the northern Pacific sea star, Asterias amurensis, which lives in colder seawater around the main Islands of Japan, were also examined as controls. Our results showed that As. amurensis has at least two local populations that diverged approximately 4.7 million years ago (MYA), and no genetic exchanges have occurred between the populations since then. Linckia laevigata shows two major populations in the Ryukyu Archipelago that likely diverged ∼6.8 MYA. The two populations, each comprised of individuals collected from coast of the Okinawa Island and those from the Ishigaki Island, suggest the presence of two cryptic species in the Ryukyu Archipelago. On the other hand, population genetics of COTS showed a profile quite different from those of Asterias and Linckia. At least five lineages of COTS have arisen since their divergence ∼0.7 MYA, and each of the lineages is present at the Okinawa Island, Miyako Island, and Ishigaki Island. These results suggest that COTS have experienced repeated genetic bottlenecks that may be associated with or caused by repeated outbreaks.

An initial comparative genomic autopsy of wasting disease in sea stars

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

Nonspecific expression of fertilization genes in the crown-of-thorns Acanthaster cf. solaris: Unexpected evidence of hermaphroditism in a coral reef predator

The characterization of gene expression in gametes has advanced our understanding of the molecular basis for ecological variation in reproductive success and the evolution of reproductive isolation. These advances are especially significant for ecologically important keystone predators such as the coral-eating crown-of-thorns sea stars (COTS, Acanthaster) which are the most influential predator species in Indo-Pacific coral reef ecosystems and the focus of intensive management efforts. We used RNA-seq and transcriptome assemblies to characterize the expression of genes in mature COTS gonads. We described the sequence and domain organization of eight genes with sex-specific expression and well known functions in fertilization in other echinoderms. We found unexpected expression of genes in one ovary transcriptome that are characteristic of males and sperm, including genes that encode the sperm-specific guanylate cyclase receptor for an egg pheromone, and the sperm acrosomal protein bindin. In a reassembly of previously published RNA-seq data from COTS testes, we found a complementary pattern: strong expression of four genes that are otherwise well known to encode egg-specific fertilization proteins, including the egg receptor for bindin (EBR1) and the acrosome reaction-inducing substance in the egg coat (ARIS1, ARIS2, ARIS3). We also found histological evidence of both eggs and sperm developing in the same gonad in several COTS individuals from a parallel study. These results suggest the occurrence of hermaphrodites, and the potential for reproductive assurance via self-fertilization. Our findings have implications for management of COTS populations, especially in consideration of the large size and massive fecundity of these sea stars.

Complete mitochondrial genome of the sea star Archaster typicus (Asteroidea: Archasteridae)

The complete mitochondrial genome of the widespread and common Indo-Pacific sea star Archaster typicus has been sequenced in this study. The mitogenome is 16,230 base pairs (bp) in length, with 13 protein coding genes (PCGs), 22 tRNAs and 2 rRNAs. Gene order of its PCGs and rRNAs matches those of nine other asteroid taxa included for comparison in this study, and it has a similar nucleotide composition of 33.08% A, 26.38% T, 25.53% C and 15.01% G nucleotides. Phylogenetic analyses place A. typicus as the sister group to Acanthaster spp., consistent with previous inferences.

The first complete mitochondrial genome of the Mariana Trench Freyastera benthophila (Asteroidea: Brisingida: Brisingidae) allows insights into the deep-sea adaptive evolution of Brisingida

Starfish (phylum Echinodermata) are ecologically important and diverse members of marine ecosystems in all of the world's oceans, from the shallow water to the hadal zone. The deep sea is recognized as an extremely harsh environment on earth. In this study, we present the mitochondrial genome sequence of Mariana Trench starfish Freyastera benthophila, and this study is the first to explore in detail the mitochondrial genome of a deep-sea member of the order Brisingida. Similar to other starfish, it contained 13 protein-coding genes, two ribosomal RNA genes, and 22 transfer RNA genes (duplication of two tRNAs: trnL and trnS). Twenty-two of these genes are encoded on the positive strand, while the other 15 are encoded on the negative strand. The gene arrangement was identical to those of sequenced starfish. Phylogenetic analysis showed the deep-sea Brisingida as a sister taxon to the traditional members of the Asteriidae. Positive selection analysis indicated that five residues (8 N and 16 I in atp8, 47 D and 196 V in nad2, 599 N in nad5) were positively selected sites with high posterior probabilities. Compared these features with shallow sea starfish, we predict that variation specifically in atp8, nad2, and nad5 may play an important role in F. benthophila's adaptation to deep-sea environment.

Complete mitochondrial genome analysis of Distolasterias nipon (Echinodermata, Asteroidea, Forcipulatida)

In this study, we determined the complete mitochondrial genome sequences of Distolasterias nipon belonging to the class Asteroidea in the phylum Echinodermata. The complete mitogenome of D. nipon was 16,193 bp in length and consisted of 13 protein-coding genes (PCGs), 22 tRNAs, and 2 rRNAs. The orders of PCGs and rRNAs were identical to those of the eight recorded mitogenomes of asteroids. The phylogenetic tree determined by the maximum likelihood method revealed that D. nipon was clearly grouped with Aphelasterias japonica and Asterias amurensis, which belong to the family Asteriidae.

Conservation of mitochondrial genome arrangements in brittle stars (Echinodermata, Ophiuroidea)

Brittle stars are conspicuous members of benthic ecosystems, fill many ecological niches and are the most speciose of all classes of echinoderms. With high levels of biodiversity, elucidating the evolutionary history of this group is important. Understanding of higher-level relationships within Ophiuroidea has been aided by multilocus nuclear data and DNA barcoding. However, the degree of consistency between mitochondrial and nuclear data within ophiuroids remains unclear and deserves further assessment. In this study, 17 mitochondrial genomes spanning the taxonomic breadth of Ophiuroidea were utilized to explore evolutionary relationships through maximum likelihood analyses, Bayesian inference and comparative assessment of gene order. Our phylogenetic analyses, based on both nucleotide and amino acid residues, support recent findings based on multilocus nuclear data and morphology, in that the brittle star clades Ophintegrida and Euryophiurida were recovered as monophyletic with the latter comprising Euyalida, Ophiuridae and Ophiopyrgidae. Only three different arrangements of the 13 protein coding and 2 ribosomal RNA genes were observed. As expected, tRNA genes were more likely to have undergone rearrangement but the order of all 37 genes was found to be conserved in all sampled Euryalida and Ophiuridae. Both Euryalida and the clade comprised of Ophiuridae and Ophiopyrgidae, each had their own conserved rearrangement of protein coding genes and ribosomal genes, after divergence from their last common ancestor. Euryalida has a rearrangement of the two ribosomal RNA genes, rrnS and rrnL, in contrast to Ophiuridae and Ophiopyrgidae, which had an inversion of the genes nad1, nad2, and cob relative to Ophintegrida. Further, our data support the gene order found in all sampled Euryalida as the most likely ancestral order for all Ophiuroidea.

Complexity of Yolk Proteins and Their Dynamics in the Sea Star Patiria miniata

Oviparous animals store yolk proteins within the developing oocyte. These proteins are used in gametogenesis and as a nutritional source for embryogenesis. Vitellogenin and the major yolk protein are two of the most important yolk proteins among diverse species of invertebrates and vertebrates. Among the echinoderms, members of the subphyla Echinozoa (sea urchins and sea cucumbers) express the major yolk protein (MYP) but not vitellogenin (Vtg), while an initial report has documented that two Asterozoa (sea stars) express a vitellogenin. Our results show that sea stars contain two vitellogenins, Vtg 1 and Vtg 2, and MYP. In Patiria miniata, these genes are differentially expressed in the somatic and germ cells of the ovary: Vtg 1 is enriched in the somatic cells of the ovary but not in the oocytes, and Vtg 2 accumulates in both oocytes and somatic cells; MYP is not robustly present in either. Remarkably, Vtg 2 and MYP mRNA reappear in larvae; Vtg 2 is detected within cells of the ectoderm, and MYP accumulates in the coelomic pouches, the intestine, and the posterior enterocoel (PE), the site of germ line formation in this animal. Additionally, the Vtg 2 protein is present in oocytes, follicle cells, and developing embryos, but becomes undetectable following gastrulation. These results help elucidate the mechanisms involved in yolk dynamics, and provide molecular information that allows for greater understanding of the evolution of these important gene products.

Phylogeography of the crown-of-thorns starfish in the Indian Ocean

Background

Understanding the limits and population dynamics of closely related sibling species in the marine realm is particularly relevant in organisms that require management. The crown-of-thorns starfish Acanthaster planci, recently shown to be a species complex of at least four closely related species, is a coral predator infamous for its outbreaks that have devastated reefs throughout much of its Indo-Pacific distribution.

Methodology/Principal Findings

In this first Indian Ocean-wide genetic study of a marine organism we investigated the genetic structure and inferred the paleohistory of the two Indian Ocean sister-species of Acanthaster planci using mitochondrial DNA sequence analyses. We suggest that the first of two main diversification events led to the formation of a Southern and Northern Indian Ocean sister-species in the late Pliocene-early Pleistocene. The second led to the formation of two internal clades within each species around the onset of the last interglacial. The subsequent demographic history of the two lineages strongly differed, the Southern Indian Ocean sister-species showing a signature of recent population expansion and hardly any regional structure, whereas the Northern Indian Ocean sister-species apparently maintained a constant size with highly differentiated regional groupings that were asymmetrically connected by gene flow.

Conclusions/Significance

Past and present surface circulation patterns in conjunction with ocean primary productivity were identified as the processes most likely to have shaped the genetic structure between and within the two Indian Ocean lineages. This knowledge will help to understand the biological or ecological differences of the two sibling species and therefore aid in developing strategies to manage population outbreaks of this coral predator in the Indian Ocean.

The complete mitochondrial genome of Galba pervia (Gastropoda: Mollusca), an intermediate host snail of Fasciola spp

Complete mitochondrial (mt) genomes and the gene rearrangements are increasingly used as molecular markers for investigating phylogenetic relationships. Contributing to the complete mt genomes of Gastropoda, especially Pulmonata, we determined the mt genome of the freshwater snail Galba pervia, which is an important intermediate host for Fasciola spp. in China. The complete mt genome of G. pervia is 13,768 bp in length. Its genome is circular, and consists of 37 genes, including 13 genes for proteins, 2 genes for rRNA, 22 genes for tRNA. The mt gene order of G. pervia showed novel arrangement (tRNA-His, tRNA-Gly and tRNA-Tyr change positions and directions) when compared with mt genomes of Pulmonata species sequenced to date, indicating divergence among different species within the Pulmonata. A total of 3655 amino acids were deduced to encode 13 protein genes. The most frequently used amino acid is Leu (15.05%), followed by Phe (11.24%), Ser (10.76%) and IIe (8.346%). Phylogenetic analyses using the concatenated amino acid sequences of the 13 protein-coding genes, with three different computational algorithms (maximum parsimony, maximum likelihood and Bayesian analysis), all revealed that the families Lymnaeidae and Planorbidae are closely related two snail families, consistent with previous classifications based on morphological and molecular studies. The complete mt genome sequence of G. pervia showed a novel gene arrangement and it represents the first sequenced high quality mt genome of the family Lymnaeidae. These novel mtDNA data provide additional genetic markers for studying the epidemiology, population genetics and phylogeographics of freshwater snails, as well as for understanding interplay between the intermediate snail hosts and the intra-mollusca stages of Fasciola spp..

There's no place like home: crown-of-thorns outbreaks in the central pacific are regionally derived and independent events

One of the most significant biological disturbances on a tropical coral reef is a population outbreak of the fecund, corallivorous crown-of-thorns sea star, Acanthaster planci. Although the factors that trigger an initial outbreak may vary, successive outbreaks within and across regions are assumed to spread via the planktonic larvae released from a primary outbreak. This secondary outbreak hypothesis is predominantly based on the high dispersal potential of A. planci and the assertion that outbreak populations (a rogue subset of the larger population) are genetically more similar to each other than they are to low-density non-outbreak populations. Here we use molecular techniques to evaluate the spatial scale at which A. planci outbreaks can propagate via larval dispersal in the central Pacific Ocean by inferring the location and severity of gene flow restrictions from the analysis of mtDNA control region sequence (656 specimens, 17 non-outbreak and six outbreak locations, six archipelagos, and three regions). Substantial regional, archipelagic, and subarchipelagic-scale genetic structuring of A. planci populations indicate that larvae rarely realize their dispersal potential and outbreaks in the central Pacific do not spread across the expanses of open ocean. On a finer scale, genetic partitioning was detected within two of three islands with multiple sampling sites. The finest spatial structure was detected at Pearl & Hermes Atoll, between the lagoon and forereef habitats (<10 km). Despite using a genetic marker capable of revealing subtle partitioning, we found no evidence that outbreaks were a rogue genetic subset of a greater population. Overall, outbreaks that occur at similar times across population partitions are genetically independent and likely due to nutrient inputs and similar climatic and ecological conditions that conspire to fuel plankton blooms.

Identification and localization of the bilitranslocase homologue in white grape berries (Vitis vinifera L.) during ripening

A homologue of the mammalian bilirubin transporter bilitranslocase (BTL) (TCDB 2.A.65.1.1), able to perform an apparent secondary active transport of flavonoids, has previously been found in carnation petals and red grape berries. In the present work, a BTL homologue was also shown in white berries from Vitis vinifera L. cv. Tocai/Friulano, using anti-sequence antibodies specific for rat liver BTL. This transporter, similarly to what found in red grape, was localized in the first layers of the epidermal tissue and in the vascular bundle cells of the mesocarp. In addition, a strong immunochemical reaction was detected in the placental tissue and particularly in peripheral integuments of the seed. The protein was expressed during the last maturation stages in both skin and pulp tissues and exhibited an apparent molecular mass of c. 31 kDa. Furthermore, the transport activity of such a carrier, measured as bromosulphophthalein (BSP) uptake, was detected in berry pulp microsomes, where it was inhibited by specific anti-BTL antibodies. The BTL homologue activity exhibited higher values, for both K(m) and V(max), than those found in the red cultivar. Moreover, two non-pigmented flavonoids, such as quercetin (a flavonol) and eriodictyol (a flavanone), inhibited the uptake of BSP in an uncompetitive manner. Such results strengthen the hypothesis that this BTL homologue acts as a carrier involved also in the membrane transport of colourless flavonoids and demonstrate the presence of such a carrier in different organs and tissues.