Genome-wide profiling of DNA methylome and transcriptome reveals epigenetic regulation of Urechis unicinctus response to sulfide stress

Sulfide is a well-known environmental pollutant that can have detrimental effects on most organisms. However, few metazoans living in sulfide-rich environments have developed mechanisms to tolerate and adapt to sulfide stress. Epigenetic mechanisms, including DNA methylation, have been shown to play a vital role in environmental stress adaptation. Nevertheless, the precise function of DNA methylation in biological sulfide adaptation remains unclear. Urechis unicinctus, a benthic organism inhabiting sulfide-rich intertidal environments, is an ideal model organism for studying adaptation to sulfide environments. In this study, we conducted a comprehensive analysis of the DNA methylome and transcriptome of U. unicinctus after exposure to 50 μM sulfide. The results revealed dynamic changes in the DNA methylation (5-methylcytosine) landscape in response to sulfide stress, with U. unicinctus exhibiting elevated DNA methylation levels following stress exposure. Integrating differentially expressed genes (DEGs) and differentially methylated regions (DMRs), we identified a crucial role of gene body methylation in predicting gene expression. Furthermore, using a DNA methyltransferase inhibitor, we validated the involvement of DNA methylation in the sulfide stress response and the gene regulatory network influenced by DNA methylation. The results indicated that by modulating DNA methylation levels during sulfide stress, the expression of glutathione S-transferase, glutamyl aminopeptidase, and cytochrome c oxidase could be up-regulated, thereby facilitating the metabolism and detoxification of exogenous sulfides. Moreover, DNA methylation was found to regulate and enhance the oxidative phosphorylation pathway, including NADH dehydrogenase, isocitrate dehydrogenase, and ATP synthase. Additionally, DNA methylation influenced the regulation of Cytochrome P450 and macrophage migration inhibitory factor, both of which are closely associated with oxidative stress and stress resistance. Our findings not only emphasize the role of DNA methylation in sulfide adaptation but also provide novel insights into the potential mechanisms through which marine organisms adapt to environmental changes.

Temperature and sex shape reproductive barriers in a climate change hotspot

Climate change is altering species ranges and reproductive interactions in existing ranges, offering species new scope to mate and hybridize. The outcomes will depend on how environmental factors shape reproductive barriers across life stages, yet this is rarely assessed across the environments that species encounter in nature. We assess prezygotic and postzygotic barriers, and their dependence on temperature and parental sex, in species of a reef-building tubeworm (Galeolaria) from a fast-warming biodiversity hotspot in southern Australia. By replicating pure and reciprocal hybrid crosses across 5 temperatures spanning species' thermal ranges, we estimate thermal tolerance curves (defining niches) for crosses and reproductive isolation at each temperature. By also replicating crosses at 3 life stages, we partition the contributions of prezygotic barriers at fertilization, postzygotic barriers at embryogenesis, and postzygotic barriers at larval development to reproductive isolation. We show that barriers are weaker at fertilization and embryogenesis, but stronger and more temperature sensitive at larval development, as species diverge in thermal niche. Asymmetry of barriers between parental sexes, moreover, suggests a complex interplay between niche differentiation and maternal inheritance. Our findings point to a key role for temperature in reproductive isolation, but also challenges for predicting the fate of isolation in future climates.

Identification of the principal neuropeptide MIP and its action pathway in larval settlement of the echiuran worm Urechis unicinctus

BACKGROUND

Larval settlement and metamorphosis represent critical events in the life history of marine benthic animals. Myoinhibitory peptide (MIP) plays a pivotal role in larval settlement of marine invertebrates. However, the molecular mechanisms of MIP involved in this process are not well understood.

RESULTS

In this study, we evaluated the effects of thirteen MIP mature peptides on triggering the larval settlement of Urechis unicinctus (Xenopneusta, Urechidae), and determined that MIP2 was the principal neuropeptide. Transcriptomic analysis was employed to identify differentially expressed genes (DEGs) between the MIP2-treated larvae and normal early-segmentation larvae. Both cAMP and calcium signaling pathways were enriched in the DEGs of the MIP2-treated larvae, and two neuropeptide receptor genes (Spr, Fmrfar) were up-regulated in the MIP2-treated larvae. The activation of the SPR-cAMP pathway by MIP2 was experimentally validated in HEK293T cells. Furthermore, fourteen cilia-related genes, including Tctex1d2, Cfap45, Ift43, Ift74, Ift22, Cav1 and Mns1, etc. exhibited down-regulated expression in the MIP2-treated larvae. Whole-mount in situ hybridization identified two selected ciliary genes, Tctex1d2 and Cfap45, were specially expressed in circumoral ciliary cells of the early-segmentation larvae. Knocking down Tctex1d2 mRNA levels by in vivo RNA interference significantly increased the larval settlement rate.

CONCLUSION

Our findings suggest that MIP2 inhibits the function of the cilia-related genes, such as Tctex1d2, through the SPR-cAMP-PKA pathway, thereby inducing larval settlement in U. unicinctus. The study contributes important data to the understanding of neuropeptide regulation in larval settlement.

A whole-body transcriptome assembly of the annelid worm Hediste diversicolor

The Annelida phylum is composed of a myriad of species exhibiting key phenotypic adaptations. They occupy key ecological niches in a variety of marine, freshwater and terrestrial ecosystems. Importantly, the increment of omic resources is rapidly modifying the taxonomic landscape and knowledge of species belonging to this phylum. Here, we comprehensively characterised and annotated a transcriptome of the common ragworm, Hediste diversicolor (OF Müller). This species belongs to the family Nereididae and inhabits estuarine and lagoon areas on the Atlantic coasts of Europe and North America. Ecologically, H. diversicolor plays an important role in benthic food webs. Given its commercial value, H. diversicolor is a promising candidate for aquaculture development and production in farming facilities, under a circular economy framework. We used Illumina next-generation sequencing technology, to produce a total of 105 million (M) paired-end (PE) raw reads and generate the first whole-body transcriptome assembly of H. diversicolor species. This high-quality transcriptome contains 69,335 transcripts with an N50 transcript length of 2313 bp and achieved a BUSCO gene completeness of 97.7% and 96% in Eukaryota and Metazoa lineage-specific profile libraries. Our findings offer a valuable resource for multiple biological applications using this species.

A remarkable new deep-sea nereidid (Annelida: Nereididae) with gills

Nereidid polychaetes are well known from shallow marine habitats, but their diversity in the deep sea is poorly known. Here we describe an unusual new nereidid species found at methane seeps off the Pacific coast of Costa Rica. Specimens of Pectinereis strickrotti gen. nov., sp. nov. had been observed dating back to 2009 swimming just above the seafloor at ~1,000 m depth but were not successfully captured until 2018. Male epitokes were collected as well as a fragment of an infaunal female found in a pushcore sample. The specimens were all confirmed as the same species based on mitochondrial COI. Phylogenetic analyses, including one based on available whole mitochondrial genomes for nereidids, revealed no close relative, allowing for the placement of the new species in its own genus within the subfamily Nereidinae. This was supported by the unusual non-reproductive and epitokous morphology, including parapodial cirrostyles as pectinate gills, hooked aciculae, elfin-shoe-shaped ventral cirrophores, and elongate, fusiform dorsal ligules emerging sub-medially to enlarged cirrophores. Additionally, the gill-bearing subfamily Dendronereidinae, generally regarded as a junior synonym of Gymnonereidinae, is reviewed and it is here reinstated and as a monogeneric taxon.

Biomineralization in a cold environment: Insights from shield compositions and transcriptomics of polar sternaspids (Sternaspidae, Polychaeta)

The survival and physiological functions of polar marine organisms are impacted by global climate changes. Investigation of the adaptation mechanisms underlying biomineralization in polar organisms at low temperatures is important for understanding mineralized organismal sensitivity to climate change. Here, we performed electron probe analysis on the shields of Antarctic polychaete Sternaspis sendalli and Arctic polychaete Sternaspis buzhinskajae (Sternaspidae), and sequenced the transcriptomes of the tissues surrounding shields to examine biomineral characteristics and adaptive mechanisms in persistently cold environments. Compared to the temperate relative species, the relative abundance of iron, phosphorus, calcium, magnesium, nitrogen, sulfur and silicon in two polar sternaspid shields was similar to Sternaspis chinensis. However, the diversity and expression levels of biomineralization-related shell matrix proteins differed between the polar and temperate species, suggesting distinct molecular mechanisms underlying shield formation in cold environments. Tubulin and cyclophilin were upregulated compared to the temperate species. Furthermore, 42 positively selected genes were identified in Antarctic S. sendalli, with functions in cytoskeletal structure, DNA repair, immunity, transcription, translation, protein synthesis, and lipid metabolism. Highly expressed genes in both polar species were associated with cytoskeleton, macromolecular complexes and cellular component biosynthesis. Overall, this study reveals conserved elemental composition yet distinct biomineralization processes in the shields of polar sternaspids. The unique expression of biomineralization related genes and other cold-adaptation related genes provide molecular insights into biomineralization in cold marine environments.

The development of the adult nervous system in the annelid Owenia fusiformis

BACKGROUND

The evolutionary origins of animal nervous systems remain contentious because we still have a limited understanding of neural development in most major animal clades. Annelids - a species-rich group with centralised nervous systems - have played central roles in hypotheses about the origins of animal nervous systems. However, most studies have focused on adults of deeply nested species in the annelid tree. Recently, Owenia fusiformis has emerged as an informative species to reconstruct ancestral traits in Annelida, given its phylogenetic position within the sister clade to all remaining annelids.

METHODS

Combining immunohistochemistry of the conserved neuropeptides FVamide-lir, RYamide-lir, RGWamide-lir and MIP-lir with gene expression, we comprehensively characterise neural development from larva to adulthood in Owenia fusiformis.

RESULTS

The early larval nervous system comprises a neuropeptide-rich apical organ connected through peripheral nerves to a prototroch ring and the chaetal sac. There are seven sensory neurons in the prototroch. A bilobed brain forms below the apical organ and connects to the ventral nerve cord of the developing juvenile. During metamorphosis, the brain compresses, becoming ring-shaped, and the trunk nervous system develops several longitudinal cords and segmented lateral nerves.

CONCLUSIONS

Our findings reveal the formation and reorganisation of the nervous system during the life cycle of O. fusiformis, an early-branching annelid. Despite its apparent neuroanatomical simplicity, this species has a diverse peptidergic nervous system, exhibiting morphological similarities with other annelids, particularly at the larval stages. Our work supports the importance of neuropeptides in animal nervous systems and highlights how neuropeptides are differentially used throughout development.

A chromosome-level genome assembly of the Echiura Urechis unicinctus

Echiura is a distinctive family of unsegmented sausage-shaped marine worms whose phylogenetic relationship still needs strong evidence from the phylogenomic analysis. In this family, Urechis unicinctus is known for its high nutritional and medicinal value and adaptation to harsh intertidal conditions. Herein, we combined PacBio long-read, short-read Illumina and Hi-C sequencing, generating a high-quality chromosome-level genome assembly of U. unicinctus. The assembled genome spans ~1,138.6 Mb with a scaffold N50 of 68.3 Mb, of which 1,113.8 Mb (97.82%) were anchored into 17 pseudo-chromosomes. The BUSCO analysis demonstrated the completeness of the genome assembly and gene model prediction are 93.5% and 91.5%, respectively. A total of 482.1 Mb repetitive sequences, 21,524 protein-coding genes, 1,535 miRNAs, 3,431 tRNAs, 124 rRNAs, and 348 snRNAs were annotated. This study significantly improves the quality of U. unicinctus genome assembly, sets the footsteps for molecular breeding and further study in genome evolution, genetic and molecular biology of U. unicinctus.

Fast cycling culture of the annelid model Platynereis dumerilii

Platynereis dumerilii, a marine annelid, is a model animal that has gained popularity in various fields such as developmental biology, biological rhythms, nervous system organization and physiology, behaviour, reproductive biology, and epigenetic regulation. The transparency of P. dumerilii tissues at all developmental stages makes it easy to perform live microscopic imaging of all cell types. In addition, the slow-evolving genome of P. dumerilii and its phylogenetic position as a representative of the vast branch of Lophotrochozoans add to its evolutionary significance. Although P. dumerilii is amenable to transgenesis and CRISPR-Cas9 knockouts, its relatively long and indefinite life cycle, as well as its semelparous reproduction have been hindrances to its adoption as a reverse genetics model. To overcome this limitation, an adapted culturing method has been developed allowing much faster life cycling, with median reproductive age at 13-14 weeks instead of 25-35 weeks using the traditional protocol. A low worm density in boxes and a strictly controlled feeding regime are important factors for the rapid growth and health of the worms. This culture method has several advantages, such as being much more compact, not requiring air bubbling or an artificial moonlight regime for synchronized sexual maturation and necessitating only limited water change. A full protocol for worm care and handling is provided.

Disjunct habitat of cryptic Terebellides (Annelida, Trichobranchidae) species shows a phylogenetic link between polychaetes from the White and the North Seas

Understanding the distribution and biodiversity of marine species is crucial for developing effective conservation strategies and maintaining the health of global ecosystems. Advancements in molecular data utilization have significantly improved our understanding of biodiversity within the genus Terebellides. In this study, we conducted a phylogenetic analysis on polychaete samples from the Kandalaksha Bay, White Sea, revealing their affiliation with a putative undescribed species of the genus Terebellides found in two locations of the North Sea. Interestingly, this species was not detected in the Norwegian and Barents Seas, leading us to propose a disjunct distribution scenario for this Terebellides species. This unique distribution pattern might be attributed to the succession of polychaetes by new species, facilitated by the Gulf Stream and a climate change role in driving shifts in species' ranges and altering marine ecosystem dynamics.

Buried treasure in a public repository: Mining mitochondrial genes of 32 annelid species from sequence reads deposited in the Sequence Read Archive (SRA)

Background

The mitochondrial genomes (mitogenomes) of metazoans generally include the same set of protein-coding genes, which ensures the homology of mitochondrial genes between species. The mitochondrial genes are often used as reference data for species identification based on genetic data (DNA barcoding). The need for such reference data has been increasing due to the application of environmental DNA (eDNA) analysis for environmental assessments. Recently, the number of publicly available sequence reads obtained with next-generation sequencing (NGS) has been increasing in the public database (the NCBI Sequence Read Archive, SRA). Such freely available NGS reads would be promising sources for assembling mitochondrial protein-coding genes (mPCGs) of organisms whose mitochondrial genes are not available in GenBank. The present study aimed to assemble annelid mPCGs from raw data deposited in the SRA.

Methods

The recent progress in the classification of Annelida was briefly introduced. In the present study, the mPCGs of 32 annelid species of 19 families in clitellates and allies in Sedentaria (echiurans and polychaetes) were newly assembled from the reads deposited in the SRA. Assembly was performed with a recently published pipeline mitoRNA, which includes cycles of Bowtie2 mapping and Trinity assembly. Assembled mPCGs were deposited in GenBank as Third Party Data (TPA) data. A phylogenetic tree was reconstructed with maximum likelihood (ML) analysis, together with other mPCGs deposited in GenBank.

Results and Discussion

mPCG assembly was largely successful except for Travisia forbesii; only four genes were detected from the assembled contigs of the species probably due to the reads targeting its parasite. Most genes were largely successfully obtained, whereas atp8, nad2, and nad4l were only successful in 22-24 species. The high nucleotide substitution rates of these genes might be relevant to the failure in the assembly although nad6, which showed a similarly high substitution rate, was successfully assembled. Although the phylogenetic positions of several lineages were not resolved in the present study, the phylogenetic relationships of some polychaetes and leeches that were not inferred by transcriptomes were well resolved probably due to a more dense taxon sampling than previous phylogenetic analyses based on transcriptomes. Although NGS data are generally better sources for resolving phylogenetic relationships of both higher and lower classifications, there are ensuring needs for specific loci of the mitochondrial genes for analyses that do not require high resolutions, such as DNA barcoding, eDNA, and phylogenetic analysis among lower taxa. Assembly from publicly available NGS reads would help design specific primers for the mitochondrial gene sequences of species, whose mitochondrial genes are hard to amplify by Sanger sequencing using universal primers.

Morphological, histological and gene-expression analyses on stolonization in the Japanese Green Syllid, Megasyllis nipponica (Annelida, Syllidae)

Benthic annelids belonging to the family Syllidae (Annelida, Errantia, Phyllodocida) exhibit a unique reproduction mode called "schizogamy" or "stolonization", in which the posterior body part filled with gametes detaches from the original body, as a reproductive unit (stolon) that autonomously swims and spawns. In this study, morphological and histological observations on the developmental processes during stolonization were carried out in Megasyllis nipponica. Results suggest that the stolon formation started with maturation of gonads, followed by the formation of a head ganglion in the anteriormost segment of the developing stolon. Then, the detailed stolon-specific structures such as stolon eyes and notochaetae were formed. Furthermore, expression profiles of genes involved in the anterior-posterior identity (Hox genes), head determination, germ-line, and hormone regulation were compared between anterior and posterior body parts during the stolonization process. The results reveal that, in the posterior body part, genes for gonadal development were up-regulated, followed by hormone-related genes and head-determination genes. Unexpectedly, Hox genes known to identify body parts along the anterior-posterior axis showed no significant temporal expression changes. These findings suggest that during stolonization, gonad development induces the head formation of a stolon, without up-regulation of anterior Hox genes.

Transcriptomic landscape of posterior regeneration in the annelid Platynereis dumerilii

BACKGROUND

Restorative regeneration, the capacity to reform a lost body part following amputation or injury, is an important and still poorly understood process in animals. Annelids, or segmented worms, show amazing regenerative capabilities, and as such are a crucial group to investigate. Elucidating the molecular mechanisms that underpin regeneration in this major group remains a key goal. Among annelids, the nereididae Platynereis dumerilii (re)emerged recently as a front-line regeneration model. Following amputation of its posterior part, Platynereis worms can regenerate both differentiated tissues of their terminal part as well as a growth zone that contains putative stem cells. While this regeneration process follows specific and reproducible stages that have been well characterized, the transcriptomic landscape of these stages remains to be uncovered.

RESULTS

We generated a high-quality de novo Reference transcriptome for the annelid Platynereis dumerilii. We produced and analyzed three RNA-sequencing datasets, encompassing five stages of posterior regeneration, along with blastema stages and non-amputated tissues as controls. We included two of these regeneration RNA-seq datasets, as well as embryonic and tissue-specific datasets from the literature to produce a Reference transcriptome. We used this Reference transcriptome to perform in depth analyzes of RNA-seq data during the course of regeneration to reveal the important dynamics of the gene expression, process with thousands of genes differentially expressed between stages, as well as unique and specific gene expression at each regeneration stage. The study of these genes highlighted the importance of the nervous system at both early and late stages of regeneration, as well as the enrichment of RNA-binding proteins (RBPs) during almost the entire regeneration process.

CONCLUSIONS

In this study, we provided a high-quality de novo Reference transcriptome for the annelid Platynereis that is useful for investigating various developmental processes, including regeneration. Our extensive stage-specific transcriptional analysis during the course of posterior regeneration sheds light upon major molecular mechanisms and pathways, and will foster many specific studies in the future.

Updated phylogeny of Vestimentifera (Siboglinidae, Polychaeta, Annelida) based on mitochondrial genomes, with a new species

Siboglinid tubeworms are found at chemosynthetic environments worldwide and the Vestimentifera clade is particularly well known for their reliance on chemoautotrophic bacterial symbionts for nutrition. The mitochondrial genomes have been published for nine vestimentiferan species to date. This study provides new complete mitochondrial genomes for ten further Vestimentifera, including the first mitochondrial genomes sequenced for Alaysia spiralis, Arcovestia ivanovi, Lamellibrachia barhami, Lamellibrachia columna, Lamellibrachia donwalshi, and unnamed species of Alaysia and Oasisia. Phylogenetic analyses combining fifteen mitochondrial genes and the nuclear 18S rRNA gene recovered Lamellibrachia as sister to the remaining Vestimentifera and Riftia pachyptila as separate from the other vent-endemic taxa. Implications and auxiliary analyses regarding differing phylogenetic tree topologies, substitution saturation, ancestral state reconstruction, and divergence estimates are also discussed. Additionally, a new species of Alaysia is described from the Manus Basin.

First report of the mitogenome of the invasive reef-building polychaete Ficopomatus enigmaticus (Annelida: Serpulidae) and a cryptic lineage from the Japanese Archipelago

BACKGROUND

The mitochondrial genomes (mitogenomes) of the family Serpulidae are characterized by a high nucleotide sequence divergence and a significant number of gene order rearrangements compared with other families within the phylum Annelida. However, only two of 50 genera of serpulids have mitogenomes already sequenced. In this study, we report the first sequencing and assembly of the complete mitogenome of Ficopomatus, thus providing further knowledge on mitochondrial gene sequences of Serpulidae.

METHODS AND RESULTS

A mitogenome of the invasive reef-building polychaete Ficopomatus enigmaticus was amplified by long PCR and sequenced using the Illumina MiSeq System. It comprised 15,853 bp and consisted of 12 protein-coding genes (atp8 was not found), 23 tRNA, and two rRNA genes. The AT and GC skew values were infrequent when compared to annelid mitogenomes but similar to other serpulids sequenced to date (i.e., Spirobranchus and Hydroides). The mitochondrial gene order of F. enigmaticus was highly rearranged compared to other serpulids. To amplify 16S rRNA gene sequences, we developed a 16S rRNA primer set by modifying the universal primer set 16SarL/16SbrH. We detected the 16S rRNA sequence of F. enigmaticus deposited in GenBank erroneously characterized as of serpulid origin. We reported for the first time the presence of two lineages of F. enigmaticus in Japan, which have already been identified in California, Australia, and the Mediterranean.

CONCLUSIONS

The first mitochondrial genome of F. enigmaticus showed a unique gene order rearrangement, corroborating the remarkable diversity in the previously reported mitogenomes of other serpulid species. The presence of the two lineages of F. enigmaticus identified for the first time in Japan represents another case of cryptic invasion. The first 16S rRNA gene sequences of F. enigmaticus obtained in the present study can be used as reference sequences in future DNA metabarcoding studies.

Third-Generation Sequencing Reveals the Adaptive Role of the Epigenome in Three Deep-Sea Polychaetes

The roles of DNA methylation in invertebrates are poorly characterized, and critical data are missing for the phylum Annelida. We fill this knowledge gap by conducting the first genome-wide survey of DNA methylation in the deep-sea polychaetes dominant in deep-sea vents and seeps: Paraescarpia echinospica, Ridgeia piscesae, and Paralvinella palmiformis. DNA methylation calls were inferred from Oxford Nanopore sequencing after assembling high-quality genomes of these animals. The genomes of these worms encode all the key enzymes of the DNA methylation metabolism and possess a mosaic methylome similar to that of other invertebrates. Transcriptomic data of these polychaetes support the hypotheses that gene body methylation strengthens the expression of housekeeping genes and that promoter methylation acts as a silencing mechanism but not the hypothesis that DNA methylation suppresses the activity of transposable elements. The conserved epigenetic profiles of genes responsible for maintaining homeostasis under extreme hydrostatic pressure suggest DNA methylation plays an important adaptive role in these worms.

Annelid functional genomics reveal the origins of bilaterian life cycles

Indirect development with an intermediate larva exists in all major animal lineages1, which makes larvae central to most scenarios of animal evolution2-11. Yet how larvae evolved remains disputed. Here we show that temporal shifts (that is, heterochronies) in trunk formation underpin the diversification of larvae and bilaterian life cycles. We performed chromosome-scale genome sequencing in the annelid Owenia fusiformis with transcriptomic and epigenomic profiling during the life cycles of this and two other annelids. We found that trunk development is deferred to pre-metamorphic stages in the feeding larva of O. fusiformis but starts after gastrulation in the non-feeding larva with gradual metamorphosis of Capitella teleta and the direct developing embryo of Dimorphilus gyrociliatus. Accordingly, the embryos of O. fusiformis develop first into an enlarged anterior domain that forms larval tissues and the adult head12. Notably, this also occurs in the so-called 'head larvae' of other bilaterians13-17, with which the O. fusiformis larva shows extensive transcriptomic similarities. Together, our findings suggest that the temporal decoupling of head and trunk formation, as maximally observed in head larvae, facilitated larval evolution in Bilateria. This diverges from prevailing scenarios that propose either co-option9,10 or innovation11 of gene regulatory programmes to explain larva and adult origins.

Gotta Go Slow: Two Evolutionarily Distinct Annelids Retain a Common Hedgehog Pathway Composition, Outlining Its Pan-Bilaterian Core

Hedgehog signaling is one of the key regulators of morphogenesis, cell differentiation, and regeneration. While the Hh pathway is present in all bilaterians, it has mainly been studied in model animals such as Drosophila and vertebrates. Despite the conservatism of its core components, mechanisms of signal transduction and additional components vary in Ecdysozoa and Deuterostomia. Vertebrates have multiple copies of the pathway members, which complicates signaling implementation, whereas model ecdysozoans appear to have lost some components due to fast evolution rates. To shed light on the ancestral state of Hh signaling, models from the third clade, Spiralia, are needed. In our research, we analyzed the transcriptomes of two spiralian animals, errantial annelid Platynereis dumerilii (Nereididae) and sedentarian annelid Pygospio elegans (Spionidae). We found that both annelids express almost all Hh pathway components present in Drosophila and mouse. We performed a phylogenetic analysis of the core pathway components and built multiple sequence alignments of the additional key members. Our results imply that the Hh pathway compositions of both annelids share more similarities with vertebrates than with the fruit fly. Possessing an almost complete set of single-copy Hh pathway members, lophotrochozoan signaling composition may reflect the ancestral features of all three bilaterian branches.

First mitochondrial genomes of Chrysopetalidae (Annelida) from shallow-water and deep-sea chemosynthetic environments

Among Annelida, Chrysopetalidae is an ecologically and morphologically diverse group, which includes shallow-water, deep-sea, free-living, and symbiotic species. Here, the four first mitochondrial genomes of this group are presented and described. One of the free-living shallow-water species Chrysopetalum debile (Chrysopetalinae), one of the yet undescribed free-living deep-sea species Boudemos sp., and those of the two deep-sea bivalve endosymbionts Craseoschema thyasiricola and Iheyomytilidicola lauensis (Calamyzinae). An updated phylogeny of Chrysopetalidae is performed, which supports previous phylogenetic hypotheses within Chrysopetalinae and indicates a complex ecological evolution within Calamyzinae. Additionally, analyses of natural selection pressure in the four mitochondrial genomes and additional genes from the two shallow-water species Bhawania goodei and Arichlidon gathofi were performed. Relaxed selection pressure in the mitochondrion of deep-sea and symbiotic species was found, with many sites under selection identified in the COX3 gene of deep-sea species.

Genome-Wide Analyses of Heat Shock Protein Superfamily Provide New Insights on Adaptation to Sulfide-Rich Environments in Urechis unicinctus (Annelida, Echiura)

The intertidal zone is a transitional area of the land-sea continuum, in which physical and chemical properties vary during the tidal cycle and highly toxic sulfides are rich in sediments due to the dynamic regimes. As a typical species thriving in this habitat, Urechis unicinctus presents strong sulfide tolerance and is expected to be a model species for sulfide stress research. Heat shock proteins (HSPs) consist of a large group of highly conserved molecular chaperones, which play important roles in stress responses. In this study, we systematically analyzed the composition and expression of HSPs in U. unicinctus. A total of eighty-six HSP genes from seven families were identified, in which two families, including sHSP and HSP70, showed moderate expansion, and this variation may be related to the benthic habitat of the intertidal zone. Furthermore, expression analysis revealed that almost all the HSP genes in U. unicinctus were significantly induced under sulfide stress, suggesting that they may be involved in sulfide stress response. Weighted gene co-expression network analysis (WGCNA) showed that 12 HSPs, including 5 sHSP and 4 HSP70 family genes, were highly correlated with the sulfide stress response which was distributed in steelblue and green modules. Our data indicate that HSPs, especially sHSP and HSP70 families, may play significant roles in response to sulfide stress in U. unicinctus. This systematic analysis provides valuable information for further understanding of the function of the HSP gene family for sulfide adaptation in U. unicinctus and contributes a better understanding of the species adaptation strategies of marine benthos in the intertidal zone.

Variable heat shock response in Antarctic biofouling serpulid worms

The classical heat shock response (HSR) with up-regulation of hsp70 in response to warming is often absent in Antarctic marine species. Whilst in Antarctic fish, this is due to a mutation in the gene promoter region resulting in permanent constitutive expression of the inducible form of hsp70; there are further questions as to whether evolution to life below 0 °C has resulted in a generalised alteration to the HSR in Antarctic marine invertebrates. However, the number of species investigated to date is limited. In the first evaluation of the HSR in two spirorbid polychaetes Romanchella perrieri and Protolaeospira stalagmia, we show highly variable results of HSR induction depending on warming regimes. These animals were subjected to in situ warming (+ 1 °C and + 2 °C above ambient conditions) using heated settlement panels for 18 months, and then the HSR was tested in R. perrieri using acute and chronic temperature elevation trials. The classic HSR was not induced in response to acute thermal challenge in this species (2 h at 15 °C) and significant down-regulation of hsp90 occurred during chronic warming at 4 °C for 30 days. Analysis of heat shock protein (HSP) genes in a transcriptome study of P. stalagmia, which had been warmed in situ for 18 months, showed up-regulation of HSP70 and HSP90 family members, thus further emphasising the complexity of the response in Antarctic marine species. It is increasingly apparent that the Antarctic HSR has evolved in a species-specific manner to life in the cold.

The worm affair: fidelity and environmental adaptation in symbiont species that co-occur in vestimentiferan tubeworms

The symbioses between the vestimentiferan tubeworms and their chemosynthetic partners (Gammaproteobacteria, Chromatiales and Sedimenticolaceae) hallmark the success of these organisms in hydrothermal vent and hydrocarbon seep deep-sea habitats. The fidelity of these associations varies, as both the hosts and the symbionts can be loose in partner choice. Some tubeworms may host distinct symbiont phylotypes, which often co-occur in a single host individual. To better understand the genetic basis for the promiscuity of tubeworm symbioses, we assembled and investigated metagenome-assembled genomes of two symbiont phylotypes (species, based on the average nucleotide identity < 95%) in Lamellibrachia anaximandri, a vestimentiferan endemic to the Mediterranean Sea, in individuals collected from Palinuro hydrothermal vents (Italy) and hydrocarbon seeps (Eratosthenes seamount and Palmahim disturbance). Using comparative genomics, we show that mainly mobilome and genes involved in defence mechanisms distinguish the symbiont genotypes. While many central metabolic functions are conserved in the tubeworm symbionts, nitrate respiration (Nar, Nap and Nas proteins) is modular, yet this modularity is not linked to phylotype, but rather to geographic location, potentially implying adaptation to the local environment. Our results hint that variation in a single moonlighting protein may be responsible for the fidelity of these symbioses.

DNA methylation atlas and machinery in the developing and regenerating annelid Platynereis dumerilii

BACKGROUND

Methylation of cytosines in DNA (5mC methylation) is a major epigenetic modification that modulates gene expression and constitutes the basis for mechanisms regulating multiple aspects of embryonic development and cell reprogramming in vertebrates. In mammals, 5mC methylation of promoter regions is linked to transcriptional repression. Transcription regulation by 5mC methylation notably involves the nucleosome remodeling and deacetylase complex (NuRD complex) which bridges DNA methylation and histone modifications. However, less is known about regulatory mechanisms involving 5mC methylation and their function in non-vertebrate animals. In this paper, we study 5mC methylation in the marine annelid worm Platynereis dumerilii, an emerging evolutionary and developmental biology model capable of regenerating the posterior part of its body post-amputation.

RESULTS

Using in silico and experimental approaches, we show that P. dumerilii displays a high level of DNA methylation comparable to that of mammalian somatic cells. 5mC methylation in P. dumerilii is dynamic along the life cycle of the animal and markedly decreases at the transition between larval to post-larval stages. We identify a full repertoire of mainly single-copy genes encoding the machinery associated with 5mC methylation or members of the NuRD complex in P. dumerilii and show that this repertoire is close to the one inferred for the last common ancestor of bilaterians. These genes are dynamically expressed during P. dumerilii development and regeneration. Treatment with the DNA hypomethylating agent Decitabine impairs P. dumerilii larval development and regeneration and has long-term effects on post-regenerative growth.

CONCLUSIONS

Our data reveal high levels of 5mC methylation in the annelid P. dumerilii, highlighting that this feature is not specific to vertebrates in the bilaterian clade. Analysis of DNA methylation levels and machinery gene expression during development and regeneration, as well as the use of a chemical inhibitor of DNA methylation, suggest an involvement of 5mC methylation in P. dumerilii development and regeneration. We also present data indicating that P. dumerilii constitutes a promising model to study biological roles and mechanisms of DNA methylation in non-vertebrate bilaterians and to provide new knowledge about evolution of the functions of this key epigenetic modification in bilaterian animals.

RNA-Seq reveals divergent gene expression between larvae with contrasting trophic modes in the poecilogonous polychaete Boccardia wellingtonensis

The polychaete Boccardia wellingtonensis is a poecilogonous species that produces different larval types. Females may lay Type I capsules, in which only planktotrophic larvae are present, or Type III capsules that contain planktotrophic and adelphophagic larvae as well as nurse eggs. While planktotrophic larvae do not feed during encapsulation, adelphophagic larvae develop by feeding on nurse eggs and on other larvae inside the capsules and hatch at the juvenile stage. Previous works have not found differences in the morphology between the two larval types; thus, the factors explaining contrasting feeding abilities in larvae of this species are still unknown. In this paper, we use a transcriptomic approach to study the cellular and genetic mechanisms underlying the different larval trophic modes of B. wellingtonensis. By using approximately 624 million high-quality reads, we assemble the de novo transcriptome with 133,314 contigs, coding 32,390 putative proteins. We identify 5221 genes that are up-regulated in larval stages compared to their expression in adult individuals. The genetic expression profile differed between larval trophic modes, with genes involved in lipid metabolism and chaetogenesis over expressed in planktotrophic larvae. In contrast, up-regulated genes in adelphophagic larvae were associated with DNA replication and mRNA synthesis.

Baby makes three: Maternal, paternal, and zygotic genetic effects shape larval phenotypic evolution

The evolutionary potential of a population is shaped by the genetic architecture of its life-history traits. Early-life phenotypes are influenced by both maternal and offspring genotype, and efforts to understand life-history evolution therefore require consideration of the interactions between these separate but correlated genomes. We used a four-generation experimental pedigree to estimate the genetic architecture of early-life phenotypes in a species with dramatic variation in larval size and morphology. In the polychaete annelid Streblospio benedicti, females make either many small eggs that develop into complex larvae that feed in the plankton or few large eggs that develop into benthic juveniles without having to feed as larvae. By isolating the contributions of maternal, paternal, and zygotic genotype to larval traits, we determined that larval anatomical structures are governed by the offspring genotype at a small number of large-effect loci. Larval size is not shaped by the larva's own genotype but instead depends on loci that act in the mother, and at two genomic locations, by loci that act in the father. The overall phenotype of each larva thus depends on three separate genomes, and a population's response to selection on larval traits will reflect the interactions among them.

A novel transcription factor MRPS27 up-regulates the expression of sqr, a key gene of mitochondrial sulfide metabolism in echiuran worm Urechis unicinctus

Hydrogen sulfide is a natural, widely distributed, poisonous substance and sulfide: quinone oxidoreductase (SQR) is responsible for oxidizing hydrogen sulfide to less toxic sulfur compounds. The increase of SQR mRNA level is an important mechanism for organisms to adapt to hydrogen sulfide-rich environments. However, its transcriptional regulation mechanism is not very clear. In this study, a mitochondrial 28S ribosomal protein S27 (MRPS27), which has never been reported as a transcription factor, was screened by yeast one-hybrid experiment from the echiuran worm Urechis unicinctus, a benthic organism living in marine sediments. Western blotting indicated that UuMRPS27 contents increased significantly in the nuclear extract of hindgut under exposed to 150 μM sulfide. ChIP and EMSA assays demonstrated that UuMRPS27 did bind to the sqr proximal promoter, the key binding sequence was CTAGAG (+12 to +17 of the promoter) detected by DNase I footprinting assay as well as transient transfection experiments. Furthermore, UuMRPS27, as a transcription activator, exhibited the highest transcription activity compared with other reported sqr transcription factors. Our data revealed for the first time the role of MRPS27 acting as a transcription factor which expanded the understanding of sqr transcriptional regulation in sulfide metabolism mechanism.

Multilevel Selection on Offspring Size and the Maintenance of Variation

AbstractMultilevel selection on offspring size occurs when offspring fitness depends on both absolute size (hard selection) and size relative to neighbors (soft selection). We examined multilevel selection on egg size at two biological scales-within clutches and among clutches from different females-using an external fertilizing tube worm. We exposed clutches of eggs to two sperm environments (limiting and saturating) and measured their fertilization success. We then modeled environmental (sperm-dependent) differences in hard and soft selection on individual eggs as well as selection on clutch-level traits (means and variances). Hard and soft selection differed in strength and form depending on sperm availability-hard selection was consistently stabilizing; soft selection was directional and favored eggs relatively larger (sperm limitation) or smaller (sperm saturation) than the clutch mean. At the clutch level, selection on mean egg size was largely concave, while selection on within-clutch variance was weak but generally negative-although some correlational selection occurred between these two traits. Importantly, we found that the optimal clutch mean egg size differed for mothers and offspring, suggesting some antagonism between the levels of selection. We thus identify several pathways that may maintain offspring size variation: environmentally (sperm-) dependent soft selection, antagonistic multilevel selection, and correlational selection on clutch means and variances. Multilevel approaches are powerful but seldom-used tools for studies of offspring size, and we encourage their future use.

Ecosystem engineers drive differing microbial community composition in intertidal estuarine sediments

Intertidal systems are complex and dynamic environments with many interacting factors influencing biochemical characteristics and microbial communities. One key factor are the actions of resident fauna, many of which are regarded as ecosystem engineers because of their bioturbation, bioirrigation and sediment stabilising activities. The purpose of this investigation was to elucidate the evolutionary implications of the ecosystem engineering process by identifying, if any, aspects that act as selection pressures upon microbial communities. A mesocosm study was performed using the well characterised intertidal ecosystem engineers Corophium volutator, Hediste diversicolor, and microphytobenthos, in addition to manual turbation of sediments to compare effects of bioturbation, bioirrigation and stabilisation. A range of sediment functions and biogeochemical gradients were measured in conjunction with 16S rRNA sequencing and diatom taxonomy, with downstream bacterial metagenome function prediction, to identify selection pressures that incited change to microbial community composition and function. Bacterial communities were predominantly Proteobacteria, with the relative abundance of Bacteroidetes, Alphaproteobacteria and Verrucomicrobia being partially displaced by Deltaproteobacteria, Acidobacteria and Chloroflexi as dissolved oxygen concentration and redox potential decreased. Bacterial community composition was driven strongly by biogeochemistry; surface communities were affected by a combination of sediment functions and overlying water turbidity, and subsurface communities by biogeochemical gradients driven by sediment reworking. Diatom communities were dominated by Nitzschia laevis and Achnanthes sp., and assemblage composition was influenced by overlying water turbidity (manual or biogenic) rather than direct infaunal influences such as grazing.

Identifying conserved polychaete molecular markers of metal exposure: Comparative analyses using the Alitta virens (Annelida, Lophotrochozoa) transcriptome

Polychaetes are vital for evaluating the effects of toxic metals in marine systems, and sensitive molecular biomarkers should be integral to monitoring efforts. However, the few polychaete markers that exist are inconsistent, even within the same species, failing to identify gene expression changes in metal-exposed animals incurring clear metabolic costs. Comparing previously characterised polychaete metal-responsive genes with those of another carefully selected species could identify biomarkers applicable across polychaetes. The ragworm Alitta virens (Sars, 1835) is particularly suited for such comparisons due to its dominance of fully saline coastal areas, widespread distribution, large biomass, and its phylogenetic position relative to other polychaete 'omic' resources. A transcriptome atlas for A. virens was generated and an RNASeq-qPCR screening approach was used to characterise the response to chronic exposures of environmentally relevant concentrations of copper and zinc in controlled mesocosms. Genes presenting dramatic expression changes in A. virens were compared with known metal-responsive genes in other polychaetes to identify new possible biomarkers and assess those currently used. This revealed some current markers should probably be abandoned (e.g. Atox1), while others, such as GST-Omega, should be used with caution, as different polychaete species appear to upregulate distinct GST-Omega orthologues. In addition, the comparisons give some indication of genes that are induced by metal exposure across phylogenetically divergent polychaetes, including a suite of haemoglobin subunits and linker chains that could play conserved roles in metal-stress response. Although such newly identified markers need further characterisation, they offer alternatives to current markers that are plainly insufficient.

A Transcriptomic Approach to the Recruitment of Venom Proteins in a Marine Annelid

The growing number of known venomous marine invertebrates indicates that chemical warfare plays an important role in adapting to diversified ecological niches, even though it remains unclear how toxins fit into the evolutionary history of these animals. Our case study, the Polychaeta Eulalia sp., is an intertidal predator that secretes toxins. Whole-transcriptome sequencing revealed proteinaceous toxins secreted by cells in the proboscis and delivered by mucus. Toxins and accompanying enzymes promote permeabilization, coagulation impairment and the blocking of the neuromuscular activity of prey upon which the worm feeds by sucking pieces of live flesh. The main neurotoxins ("phyllotoxins") were found to be cysteine-rich proteins, a class of substances ubiquitous among venomous animals. Some toxins were phylogenetically related to Polychaeta, Mollusca or more ancient groups, such as Cnidaria. Some toxins may have evolved from non-toxin homologs that were recruited without the reduction in molecular mass and increased specificity of other invertebrate toxins. By analyzing the phylogeny of toxin mixtures, we show that Polychaeta is uniquely positioned in the evolution of animal venoms. Indeed, the phylogenetic models of mixed or individual toxins do not follow the expected eumetazoan tree-of-life and highlight that the recruitment of gene products for a role in venom systems is complex.

Detecting the effects of chronic metal exposure on benthic systems: Importance of biomarker and endpoint selection

Understanding metal toxicity to benthic systems is still an ecotoxicological priority and, although numerous biomarkers exist, a multi-biomarker and endpoint approach with sediment as the delivery matrix combined with life-history relevant exposure timescales is missing. Here we assess potential toxicity by measuring a suite of biomarkers and endpoints after exposing the ecologically important polychaete Alitta(Nereis)virens to sediment spiked with environmentally relevant concentrations of copper and zinc (and in combination) for 3, 6 and 9 months. We compared biomarker and endpoint sensitivity providing a guide to select the appropriate endpoints for the chosen time frame (exposure period) and concentration (relevant to Sediment Quality Guidelines) needed to identify effects for benthic polychaetes such as A. virens. Target bioavailable sediment and subsequent porewater concentrations reflect the global contamination range, whilst tissue concentrations, although elevated, were comparable with other polychaetes. Survival reduced as concentrations increased, but growth was not significantly different between treatments. Metabolic changes were restricted to significant reductions in protein after 9 months exposure across all copper concentrations, and reductions in lipid at high copper concentrations (3 months). Significant changes in feeding behaviour and increases in metallothionein-like protein concentration were limited to the medium and high copper and zinc concentrations, respectively, both after 6 months exposure. Despite data highlighting A. virens' metal tolerance, DNA damage and protein concentrations are the most sensitive biomarkers. Copper and zinc cause biomarker responses at concentrations routinely found in coastal sediments that are characterised as low contamination, suggesting a reappraisal of the current input sources (especially copper) is required.

Investigating cellular and molecular mechanisms of neurogenesis in Capitella teleta sheds light on the ancestor of Annelida

BACKGROUND

Diverse architectures of nervous systems (NSs) such as a plexus in cnidarians or a more centralized nervous system (CNS) in insects and vertebrates are present across Metazoa, but it is unclear what selection pressures drove evolution and diversification of NSs. One underlying aspect of this diversity lies in the cellular and molecular mechanisms driving neurogenesis, i.e. generation of neurons from neural precursor cells (NPCs). In cnidarians, vertebrates, and arthropods, homologs of SoxB and bHLH proneural genes control different steps of neurogenesis, suggesting that some neurogenic mechanisms may be conserved. However, data are lacking for spiralian taxa.

RESULTS

To that end, we characterized NPCs and their daughters at different stages of neurogenesis in the spiralian annelid Capitella teleta. We assessed cellular division patterns in the neuroectoderm using static and pulse-chase labeling with thymidine analogs (EdU and BrdU), which enabled identification of NPCs that underwent multiple rounds of division. Actively-dividing brain NPCs were found to be apically-localized, whereas actively-dividing NPCs for the ventral nerve cord (VNC) were found apically, basally, and closer to the ventral midline. We used lineage tracing to characterize the changing boundary of the trunk neuroectoderm. Finally, to start to generate a genetic hierarchy, we performed double-fluorescent in-situ hybridization (FISH) and single-FISH plus EdU labeling for neurogenic gene homologs. In the brain and VNC, Ct-soxB1 and Ct-neurogenin were expressed in a large proportion of apically-localized, EdU+ NPCs. In contrast, Ct-ash1 was expressed in a small subset of apically-localized, EdU+ NPCs and subsurface, EdU- cells, but not in Ct-neuroD+ or Ct-elav1+ cells, which also were subsurface.

CONCLUSIONS

Our data suggest a putative genetic hierarchy with Ct-soxB1 and Ct-neurogenin at the top, followed by Ct-ash1, then Ct-neuroD, and finally Ct-elav1. Comparison of our data with that from Platynereis dumerilii revealed expression of neurogenin homologs in proliferating NPCs in annelids, which appears different than the expression of vertebrate neurogenin homologs in cells that are exiting the cell cycle. Furthermore, differences between neurogenesis in the head versus trunk of C. teleta suggest that these two tissues may be independent developmental modules, possibly with differing evolutionary trajectories.

Molecular analysis of Spiophanes bombyx complex (Annelida: Spionidae) with description of a new species

Spiophanes bombyx (Claparède, 1870) from the Gulf of Naples, Tyrrhenian Sea, Italy, was the first described Spiophanes with fronto-lateral horns on the prostomium. It was also considered the only horned species occurring in European waters. Our sequence data of five gene fragments suggest the presence of two horned sibling Spiophanes species in northern Europe: S. cf. bombyx in the North and the Norwegian seas, and S. cf. convexus in Brittany, northern France, and Bay of Biscay, northern Spain. Spiophanes cf. bombyx worms are genetically close to a single examined specimen of S. bombyx from Venice Lagoon, Italy but their conspecificity should be verified by further study. Our sequence data show that horned Spiophanes from the North Pacific are genetically distant from horned European species, and that S. uschakowi Zachs, 1933, originally described from the Sea of Japan (East Sea) is a valid species. The data also suggest the presence of two horned sibling Spiophanes species in the North East Pacific: S. hakaiensis Radashevsky & Pankova, n. sp. distributed from Alaska south to about Point Conception, and S. norrisi Meißner & Blank, 2009, distributed from San Francisco Bay south to Baja California Sur, Mexico. Spiophanes from South America, morphologically similar to S. norrisi, are suggested to belong to a new species. Molecular data also suggest the presence of two sibling species among the worms from northern Europe identified by morphology as S. kroyeri Grube, 1860. Worms from the Barents Sea and northern part of the North Sea are tentatively referred to as S. cf. kroyeri; worms from the northern and central parts of the North Sea and from the Bay of Biscay, northern Spain, are tentatively referred to as S. cf. cirrata M. Sars in G.O. Sars, 1872. Sequence data also show that S. duplex from California is genetically different from morphologically similar worms from South America. The South American worms are referred to resurrected S. soederstroemi Hartman, 1953 which was originally described from off Rio Grande do Sul, Brazil, and then considered as a junior synonym of S. duplex. Analysis of divergence times of Spiophanes lineages suggested that the origin of the most recent common ancestor of horned Spiophanes with metameric nuchal organs was around 11.1 mya (95% HPD: 5.1–19.0 mya) and that the divergence of the North Atlantic and North Pacific lineages was around 7.9 mya (95% HPD: 4.1–13.3 mya). The North Atlantic lineage was estimated to have diverged 4.8 mya (95% HPD: 2.2–8.6 mya), resulting in the origin of S. cf. bombyx and S. cf. convexus. The North Pacific lineage was estimated to have diverged first by the isolation and speciation of S. norrisi 1.7 mya (95% HPD: 2.3–1.0 mya), and then by the isolation and speciation of S. uschakowi and S. hakaiensis n. sp. 1.3 mya (95% HPD: 2.0–0.7 mya). The estimates place the divergences soon after maximum glacial period in the North Pacific (2.4–3.0 mya).

First genetic assessment of brackish water polychaete Tylorrhynchus heterochaetus: mitochondrial COI sequences reveal strong genetic differentiation and population expansion in samples collected from southeast China and north Vietnam

Tylorrhynchus heterochaetus is a widespread benthic polychaete worm found in coastal brackish waters of the west Pacific. It has high ecological and economic value as a biomarker of water quality and as a high-quality feed in aquaculture and fisheries and is considered a delicacy in some areas of Asia. However, it has experienced a marked reduction in recent years due to overexploitation as well as changes in the environment and climate. Here, to comprehensively understand its genetic background and thus provide insights for better conservation and utilization of this species, we assessed the genetic variability and demographic history of T. heterochaetus individuals sampled from eight locations along the coasts of southeast China and north Vietnam based on mitochondrial cytochrome c oxidase I ( COI) sequences. We observed high haplotype diversity ( Hd), with an average of 0.926, but relatively low nucleotide diversity ( π), with a mean of 0.032 across all samples. A total of 94 polymorphic sites and 85 haplotypes were identified among 320 individuals. The pairwise genetic distances among haplotypes ranged from 0.001 to 0.067, with the high intraspecific divergence possibly reflecting geographic isolation and gene pool fragmentation. Significant genetic structures were revealed among the studied locations; specifically, the eight locations could be treated as six genetically different populations based on pairwise Φ ST results (0.026-0.951, P<0.01). A significant pattern of isolation-by-distance was detected between the genetic and geographic distances ( r=0.873, P=0.001). Three geographic lineages were defined based on phylogenetic tree and network analyses of COI haplotypes. AMOVA results indicated that genetic variations mainly occurred among the three lineages (89.96%). Tests of neutrality and mismatch distribution suggested that T. heterochaetus underwent recent population expansion. These results provide the first report on the genetic status of T. heterochaetus and will be valuable for the management of genetic resources and better understanding of the ecology and evolution in this species.

Intra-host symbiont diversity in eastern Pacific cold seep tubeworms identified by the 16S-V6 region, but undetected by the 16S-V4 region

Vestimentiferan tubeworms are key taxa in deep-sea chemosynthetic habitats worldwide. As adults they obtain their nutrition through their sulfide-oxidizing bacterial endosymbionts, which are acquired from the environment. Although horizontal transmission should favor infections by various symbiotic microbes, the current paradigm holds that every tubeworm harbors only one endosymbiotic 16S rRNA phylotype. Although previous studies based on traditional Sanger sequencing have questioned these findings, population level high-throughput analyses of the symbiont 16S diversity are still missing. To get further insights into the symbiont genetic variation and uncover hitherto hidden diversity we applied state-of-the-art 16S-V4 amplicon sequencing to populations of the co-occurring tubeworm species Lamellibrachia barhami and Escarpia spicata that were collected during E/V Nautilus and R/V Western Flyer cruises to cold seeps in the eastern Pacific Ocean. In agreement with earlier work our sequence data indicated that L. barhami and E. spicata share one monomorphic symbiont phylotype. However, complementary CARD-FISH analyses targeting the 16S-V6 region implied the existence of an additional phylotype in L. barhami. Our results suggest that the V4 region might not be sufficiently variable to investigate diversity in the intra-host symbiont population at least in the analyzed sample set. This is an important finding given that this region has become the standard molecular marker for high-throughput microbiome analyses. Further metagenomic research will be necessary to solve these issues and to uncover symbiont diversity that is hidden below the 16S rRNA level.

Genomic, transcriptomic, and proteomic insights into the symbiosis of deep-sea tubeworm holobionts

Deep-sea hydrothermal vents and methane seeps are often densely populated by animals that host chemosynthetic symbiotic bacteria, but the molecular mechanisms of such host-symbiont relationship remain largely unclear. We characterized the symbiont genome of the seep-living siboglinid Paraescarpia echinospica and compared seven siboglinid-symbiont genomes. Our comparative analyses indicate that seep-living siboglinid endosymbionts have more virulence traits for establishing infections and modulating host-bacterium interaction than the vent-dwelling species, and have a high potential to resist environmental hazards. Metatranscriptome and metaproteome analyses of the Paraescarpia holobiont reveal that the symbiont is highly versatile in its energy use and efficient in carbon fixation. There is close cooperation within the holobiont in production and supply of nutrients, and the symbiont may be able to obtain nutrients from host cells using virulence factors. Moreover, the symbiont is speculated to have evolved strategies to mediate host protective immunity, resulting in weak expression of host innate immunity genes in the trophosome. Overall, our results reveal the interdependence of the tubeworm holobiont through mutual nutrient supply, a pathogen-type regulatory mechanism, and host-symbiont cooperation in energy utilization and nutrient production, which is a key adaptation allowing the tubeworm to thrive in deep-sea chemosynthetic environments.

Genomic adaptations to chemosymbiosis in the deep-sea seep-dwelling tubeworm Lamellibrachia luymesi

BACKGROUND

Symbiotic relationships between microbes and their hosts are widespread and diverse, often providing protection or nutrients, and may be either obligate or facultative. However, the genetic mechanisms allowing organisms to maintain host-symbiont associations at the molecular level are still mostly unknown, and in the case of bacterial-animal associations, most genetic studies have focused on adaptations and mechanisms of the bacterial partner. The gutless tubeworms (Siboglinidae, Annelida) are obligate hosts of chemoautotrophic endosymbionts (except for Osedax which houses heterotrophic Oceanospirillales), which rely on the sulfide-oxidizing symbionts for nutrition and growth. Whereas several siboglinid endosymbiont genomes have been characterized, genomes of hosts and their adaptations to this symbiosis remain unexplored.

RESULTS

Here, we present and characterize adaptations of the cold seep-dwelling tubeworm Lamellibrachia luymesi, one of the longest-lived solitary invertebrates. We sequenced the worm's ~ 688-Mb haploid genome with an overall completeness of ~ 95% and discovered that L. luymesi lacks many genes essential in amino acid biosynthesis, obligating them to products provided by symbionts. Interestingly, the host is known to carry hydrogen sulfide to thiotrophic endosymbionts using hemoglobin. We also found an expansion of hemoglobin B1 genes, many of which possess a free cysteine residue which is hypothesized to function in sulfide binding. Contrary to previous analyses, the sulfide binding mediated by zinc ions is not conserved across tubeworms. Thus, the sulfide-binding mechanisms in sibgolinids need to be further explored, and B1 globins might play a more important role than previously thought. Our comparative analyses also suggest the Toll-like receptor pathway may be essential for tolerance/sensitivity to symbionts and pathogens. Several genes related to the worm's unique life history which are known to play important roles in apoptosis, cell proliferation, and aging were also identified. Last, molecular clock analyses based on phylogenomic data suggest modern siboglinid diversity originated in 267 mya (± 70 my) support previous hypotheses indicating a Late Mesozoic or Cenozoic origins of approximately 50-126 mya for vestimentiferans.

CONCLUSIONS

Here, we elucidate several specific adaptations along various molecular pathways that link phenome to genome to improve understanding of holobiont evolution. Our findings of adaptation in genomic mechanisms to reducing environments likely extend to other chemosynthetic symbiotic systems.

Evolution and function analysis of interleukin-17 gene from Pinctada fucata martensii

Interleukin-17 (IL-17) is a proinflammatory cytokine that plays an important role in immune responses. In this study, we identified 57 IL-17 genes from the genomes of six marine invertebrates, including Pinctada fucata martensii, Crassostrea gigas, Lottia gigantea, Capitella teleta, Mizuhopecten yessoensis, and Mytilus galloprovincialis. Phylogenetic analysis showed that all invertebrate IL-17 genes were clustered into one group, implying that invertebrate IL-17 evolved from one common ancestral gene. From the extron-intron analysis, we found many intronless IL-17 genes in mollusks, which may be caused by retroposition. Tissue and development transcriptomic analysis showed that the expression of PmIL-17 was tissue and developmental stage-specific. Moreover, we cloned the full length of the IL-17-2 gene from P. f. martensii (PmIL-17-2) and explored its function in the immune response. The full-length cDNA of PmIL-17-2 is 719 bp, containing an open reading frame of 564 bp, a 5' -untranslated region (UTR) of 31 bp, and a 3' -UTR of 124 bp with a 30 bp poly (A) tail. PmIL-17-2 had a strong response to lipopolysaccharide (LPS), indicating that the PmIL-17-2 participates in innate immune responses. In situ hybridization of hemocytes showed that PmIL-17-2 was mainly produced by granulosa cells, and the number of the stained granulosa increased after LPS stimulation. These results lay the foundation for the research of IL-17 family in marine invertebrates.

Transcriptome Analysis of Larval Segment Formation and Secondary Loss in the Echiuran Worm Urechis unicinctus

The larval segment formation and secondary loss in echiurans is a special phenomenon, which is considered to be one of the important characteristics in the evolutionary relationship between the Echiura and Annelida. To better understand the molecular mechanism of this phenomenon, we revealed the larval transcriptome profile of the echiuran worm Urechis unicinctus using RNA-Seq technology. Twelve cDNA libraries of U. unicinctus larvae, late-trochophore (LT), early-segmentation larva (ES), segmentation larva (SL), and worm-shaped larva (WL) were constructed. Totally 243,381 unigenes were assembled with an average length of 1125 bp and N50 of 1836 bp, and 149,488 unigenes (61.42%) were annotated. We obtained 70,517 differentially expressed genes (DEGs) by pairwise comparison of the larval transcriptome data at different developmental stages and clustered them into 20 gene expression profiles using STEM software. Based on the typical profiles during the larval segment formation and secondary loss, eight signaling pathways were enriched, and five of which, mTOR, PI3K-AKT, TGF-β, MAPK, and Dorso-ventral axis formation signaling pathway, were proposed for the first time to be involved in the segment formation. Furthermore, we identified 119 unigenes related to the segment formation of annelids, arthropods, and chordates, in which 101 genes were identified in Drosophila and annelids. The function of most segment polarity gene homologs (hedgehog, wingless, engrailed, etc.) was conserved in echiurans, annelids, and arthropods based on their expression profiles, while the gap and pair-rule gene homologs were not. Finally, we verified that strong positive signals of Hedgehog were indeed located on the boundary of larval segments using immunofluorescence. Data in this study provide molecular evidence for the understanding of larval segment development in echiurans and may serve as a blueprint for segmented ancestors in future research.

Species delimitation in Amblyosyllis (Annelida, Syllidae)

Amblyosyllis is a worldwide distributed group of annelids mainly found in coastal environments. It is well known among the polychaete specialists mostly because of its notable beauty, showing bright colourful patterns and outstanding long and coiled appendices. Amblyosyllis is a monophyletic genus easy to identify due to its distinct diagnostic features; however, the species and their boundaries are, in most cases, not well defined. Herein, we provide an extensive sample of Amblyosyllis material (115 specimens) from several world geographic areas. We have studied the morphological features of each specimen and photographed them alive. Two mitochondrial DNA markers (COI and 16S) and one nuclear gene fragment (28S, D1 region) were sequenced. We performed phylogenetic analyses based on each DNA partition, as well as the combined data sets, obtaining congruent results. Species delimitation methods such as distance analyses, statistical parsimony networks and multi-rate Poisson tree processes were also applied. The combined results obtained from different methodologies and data sets are used to differentiate between, at least, 19 lineages compatible with the separately evolving meta-populations species concept. Four of these lineages are identified as nominal species, including the type species of Amblyosyllis, A. rhombeata. For three other lineages previously synonymized names are recovered, and seven lineages are described as new species. All of these species are described and supported by appropriate iconography. We recognize several morphological characters useful to identify species of Amblyosyllis, which in some cases should also be combined with molecular methods for species delineation. The genetic divergence in the genus is high, contrary to the morphological homogeneity observed. Two species show a wide geographical distribution, while the rest have a more restricted distribution. There are several examples of species with overlapping distribution patterns.

Dose-dependent effects induced by cadmium in polychaete Perinereis aibuhitensis

Cadmium is a known metal contaminant in the Bohai Sea. In this study, the dose-dependent responses induced by Cd were characterized in marine polychaete Perinereis aibuhitensis using the endpoints, including activities of enzymes, expression levels of stress-responsive genes and metabolic responses. Both enzyme activities and gene expression levels exhibited the hormetic effects induced by Cd in P. aibuhitensis, as shown by the typical U-shaped or inverted U-shaped response profiles. The highest concentration (1280 µg/L) of Cd exposure induced obvious oxidative stresses. NMR-based metabolomics revealed that Cd induced both linear dose-dependent effects (69.13% of the total variation) and a relatively slight hormesis (5.54% of the total variation) in energy metabolism in P. aibuhitensis at metabolite level. In details, Cd exposures linearly reduced the consumption of amino acids and enhanced the consumption of glucose for energy supply, resulting in elevated contents of amino acids and depleted contents of glucose. Additionally, Cd treatments induced hormesis in the conversion of ATP hydrolysis to AMP. This work suggested that the hormetic effects should be considered in the ecological risk assessment for the environmental pollutants.

Distribution, diversity and functional dissociation of the mac genes in marine biofilms

Bacteria produce metamorphosis-associated contractile (MAC) structures to induce larval metamorphosis in Hydroides elegans. The distribution and diversity of mac gene homologs in marine environments are largely unexplored. In the present study mac genes were examined in marine environments by analyzing 101 biofilm and 91 seawater metagenomes. There were more mac genes in biofilms than in seawater, and substratum type, location, or sampling time did not affect the mac genes in biofilms. The mac gene clusters were highly diverse and often incomplete while the three MAC components co-occurred with other genes of different functions. Genomic analysis of four Pseudoalteromonas and two Streptomyces strains revealed the mac genes transfers among different microbial taxa. It is proposed that mac genes are more specific to biofilms; gene transfer among different microbial taxa has led to highly diverse mac gene clusters; and in most cases, the three MAC components function individually rather than forming a complex.

Chronic exposure to copper and zinc induces DNA damage in the polychaete Alitta virens and the implications for future toxicity of coastal sites

Copper and zinc are metals that have been traditionally thought of as past contamination legacies. However, their industrial use is still extensive and current applications (e.g. nanoparticles and antifouling paints) have become additional marine environment delivery routes. Determining a pollutant's genotoxicity is an ecotoxicological priority, but in marine benthic systems putative substances responsible for sediment genotoxicity have rarely been identified. Studies that use sediment as the delivery matrix combined with exposures over life-history relevant timescales are also missing for metals. Here we assess copper and zinc's genotoxicity by exposing the ecologically important polychaete Alitta virens to sediment spiked with environmentally relevant concentrations for 9 months. Target bioavailable sediment and subsequent porewater concentrations reflect the global contamination range for coasts, whilst tissue concentrations, although elevated, were comparable with other polychaetes. Survival generally reduced as concentrations increased, but monthly analyses show that growth was not significantly different between treatments. The differential treatment mortality may have enabled the surviving worms in the high concentration treatments to capture more food thus removing any concentration treatment effects for biomass. Using the alkaline comet assay we confirm that both metals via the sediment are genotoxic at concentrations routinely found in coastal regions and this is supported by elevated DNA damage in worms from field sites. However, combined with the growth data it also highlights the tolerance of A. virens to DNA damage. Finally, using long term (decadal) monitoring data we show stable or increasing sediment concentrations of these metals for many areas. This will potentially mean coastal sediment is a significant mutagenic hazard to the benthic community for decades to come. An urgent reappraisal of the current input sources for these 'old pollutants' is, therefore, required.

Morphological, cellular and molecular characterization of posterior regeneration in the marine annelid Platynereis dumerilii

Regeneration, the ability to restore body parts after an injury or an amputation, is a widespread but highly variable and complex phenomenon in animals. While having fascinated scientists for centuries, fundamental questions about the cellular basis of animal regeneration as well as its evolutionary history remain largely unanswered. Here, we present a study of regeneration of the marine annelid Platynereis dumerilii, an emerging comparative developmental biology model, which, like many other annelids, displays important regenerative abilities. When P. dumerilii worms are amputated, they are able to regenerate the posteriormost differentiated part of their body and a stem cell-rich growth zone that allows the production of new segments replacing the amputated ones. We show that posterior regeneration is a rapid process that follows a well reproducible path and timeline, going through specific stages that we thoroughly defined. Wound healing is achieved one day after amputation and a regeneration blastema forms one day later. At this time point, some tissue specification already occurs, and a functional posterior growth zone is re-established as early as three days after amputation. Regeneration timing is only influenced, in a minor manner, by worm size. Comparable regenerative abilities are found for amputations performed at different positions along the antero-posterior axis of the worm, except when amputation planes are very close to the pharynx. Regenerative abilities persist upon repeated amputations without important alterations of the process. We also show that intense cell proliferation occurs during regeneration and that cell divisions are required for regeneration to proceed normally. Finally, 5-ethynyl-2'-deoxyuridine (EdU) pulse and chase experiments suggest that blastemal cells mostly derive from the segment immediately abutting the amputation plane. The detailed characterization of P. dumerilii posterior body regeneration presented in this article provides the foundation for future mechanistic and comparative studies of regeneration in this species.

I-SceI Meganuclease-mediated transgenesis in the acorn worm, Saccoglossus kowalevskii

Hemichordates are a phylum of marine invertebrate deuterostomes that are closely related to chordates, and represent one of the most promising models to provide insights into early deuterostome evolution. The genome of the hemichordate, Saccoglossus kowalevskii, reveals an extensive set of non-coding elements conserved across all three deuterostome phyla. Functional characterization and cross-phyla comparisons of these putative regulatory elements will enable a better understanding of enhancer evolution, and subsequently how changes in gene regulation give rise to morphological innovation. Here, we describe an efficient method of transgenesis for the characterization of non-coding elements in S. kowalevskii. We first test the capacity of an I-SceI transgenesis system to drive ubiquitous or regionalized gene expression, and to label specific cell types. Finally, we identified a minimal promoter that can be used to test the capacity of putative enhancers in S. kowalevskii. This work demonstrates that this I-SceI transgenesis technique, when coupled with an understanding of chromatin accessibility, can be a powerful tool for studying how evolutionary changes in gene regulatory mechanisms contributed to the diversification of body plans in deuterostomes.

A novel transcription factor Rwdd1 and its SUMOylation inhibit the expression of sqr, a key gene of mitochondrial sulfide metabolism in Urechis unicinctus

Sulfide-quinone oxidoreductase (SQR) is a key enzyme of sulfide metabolism in metazoans, and responsible for oxidizing sulfide into thiosulfate and transmitting the generated electrons to the ubiquinone. It has been revealed that the sqr mRNA level increases significantly in echiuran worm Urechis unicinctus exposed to sulfide, and HSF1, NF1 and Sp1 have been verified to participate in its transcriptional regulation. In this study, we obtained 23 potential transcription factors interacting possibly with the proximal region (-391 to +50) of sqr promoter, and focused on the RWD domain-containing 1 (Rwdd1), a protein with the maximum number of clones in yeast one-hybrid (Y1H) screening, to investigate its transcriptional regulation to U. unincitus sqr. The ChIP and EMSA assays identified that the Rwdd1 can bind directly to the promoter (+18/+36) of U. unicinctus sqr. The point mutation and transient transfection experiments discovered that TACG was the key sequence of the DNA element bound by the Rwdd1. Furthermore, the U. unicinctus Rwdd1 (UuRwdd1) was identified to be a transcription repressor inhibiting the sqr promoter activity, and the SUMOylation of UuRwdd1 at the lysine of 90^th enhanced its inhibitory effect on sqr transcription further. Western blotting found Rwdd1 responded to sulfide in hindguts from U. unincitus, and the protein content showed a remarkable drop in hindgut nuclei in the early sulfide exposure, and then increased significantly both in the total protein and the nuclear protein extract. We suggested that the Rwdd1 is a novel transcription factor, and these data improve our understanding of the sqr transcriptional regulation and the mitochondrial sulfide metabolism.

Decoupled maternal and zygotic genetic effects shape the evolution of development

Evolutionary transitions from indirect to direct development involve changes in both maternal and zygotic genetic factors, with distinctive population-genetic implications, but empirical data on the genetics of such transitions are lacking. The polychaete Streblospio benedicti provides an opportunity to dissect a major transition in developmental mode using forward genetics. Females in this species produce either small eggs that develop into planktonic larvae or large eggs that develop into benthic juveniles. We identify large-effect loci that act maternally to influence larval size and independent, unlinked large-effect loci that act zygotically to affect discrete aspects of larval morphology. The likely fitness of zygotic alleles depends on their maternal background, creating a positive frequency-dependence that may homogenize local populations. Developmental and population genetics interact to shape larval evolution.

A Revised Spiralian Homeobox Gene Classification Incorporating New Polychaete Transcriptomes Reveals a Diverse TALE Class and a Divergent Hox Gene

The diversity of mechanisms and capacity for regeneration across the Metazoa present an intriguing challenge in evolutionary biology, impacting on the burgeoning field of regenerative medicine. Broad taxonomic sampling is essential to improve our understanding of regeneration, and studies outside of the traditional model organisms have proved extremely informative. Within the historically understudied Spiralia, the Annelida have an impressive variety of tractable regenerative systems. The biomeralizing, blastema-less regeneration of the head appendage (operculum) of the serpulid polychaete keelworm Spirobranchus (formerly Pomatoceros) lamarcki is one such system. To profile potential regulatory mechanisms, we classified the homeobox gene content of opercular regeneration transcriptomes. As a result of retrieving several difficult-to-classify homeobox sequences, we performed an extensive search and phylogenetic analysis of the TALE and PRD-class homeobox gene content of a broad selection of lophotrochozoan genomes. These analyses contribute to our increasing understanding of the diversity, taxonomic extent, rapid evolution, and radical flexibility of these recently discovered homeobox gene radiations. Our expansion and integration of previous nomenclature systems helps to clarify their cryptic orthology. We also describe an unusual divergent S. lamarcki Antp gene, a previously unclassified lophotrochozoan orphan gene family (Lopx), and a number of novel Nk class orphan genes. The expression and potential involvement of many of these lineage- and clade-restricted homeobox genes in S. lamarcki operculum regeneration provides an example of diversity in regenerative mechanisms, as well as significantly improving our understanding of homeobox gene evolution.

PdumBase: a transcriptome database and research tool for Platynereis dumerilii and early development of other metazoans

BACKGROUND

The marine polychaete annelid Platynereis dumerilii has recently emerged as a prominent organism for the study of development, evolution, stem cells, regeneration, marine ecology, chronobiology and neurobiology within metazoans. Its phylogenetic position within the spiralian/ lophotrochozoan clade, the comparatively high conservation of ancestral features in the Platynereis genome, and experimental access to any stage within its life cycle, make Platynereis an important model for elucidating the complex regulatory and functional molecular mechanisms governing early development, later organogenesis, and various features of its larval and adult life. High resolution RNA-seq gene expression data obtained from specific developmental stages can be used to dissect early developmental mechanisms. However, the potential for discovery of these mechanisms relies on tools to search, retrieve, and compare genome-wide information within Platynereis, and across other metazoan taxa.

RESULTS

To facilitate exploration and discovery by the broader scientific community, we have developed a web-based, searchable online research tool, PdumBase, featuring the first comprehensive transcriptome database for Platynereis dumerilii during early stages of development (2 h ~ 14 h). Our database also includes additional stages over the P. dumerilii life cycle and provides access to the expression data of 17,213 genes (31,806 transcripts) along with annotation information sourced from Swiss-Prot, Gene Ontology, KEGG pathways, Pfam domains, TmHMM, SingleP, and EggNOG orthology. Expression data for each gene includes the stage, the normalized FPKM, the raw read counts, and information that can be leveraged for statistical analyses of differential gene expression and the construction of genome-wide co-expression networks. In addition, PdumBase offers early stage transcriptome expression data from five further species as a valuable resource for investigators interested in comparing early development in different organisms. To understand conservation of Platynereis gene models and to validate gene annotation, most Platynereis gene models include a comprehensive phylogenetic analysis across 18 species representing diverse metazoan taxa.

CONCLUSIONS

PdumBase represents the first online resource for the early developmental transcriptome of Platynereis dumerilii. It serves as a research platform for discovery and exploration of gene expression during early stages, throughout the Platynereis life cycle, and enables comparison to other model organisms. PdumBase is freely available at http://pdumbase.gdcb.iastate.edu .

Adhesive gland transcriptomics uncovers a diversity of genes involved in glue formation in marine tube-building polychaetes

Tube-building sabellariid polychaetes are hermatypic organisms capable of forming vast reefs in highly turbulent marine habitats. Sabellariid worms assemble their tube by gluing together siliceous and calcareous clastic particles using a polyelectrolytic biocement. Here, we performed transcriptomic analyses to investigate the genes that are differentially expressed in the parathorax region, which contains the adhesive gland and tissues, from the rest of the body. We found a large number of candidate genes to be involved in the composition and formation of biocement in two species: Sabellaria alveolata and Phragmatopoma caudata. Our results indicate that the glue is likely to be composed by a large diversity of cement-related proteins, including Poly(S), GY-rich, H-repeat and miscellaneous categories. However, sequences divergence and differences in expression profiles between S. alveolata and P. caudata of cement-related proteins may reflect adaptation to the type of substratum used to build their tube, and/or to their habitat (temperate vs tropical, amplitude of pH, salinity …). Related to the L-DOPA metabolic pathways and linked with the genes that were differentially expressed in the parathorax region, we found that tyrosinase and peroxidase gene families may have undergone independent expansion in the two Sabellariidae species investigated. Our data also reinforce the importance of protein modifications in cement formation. Altogether these new genomic resources help to identify novel transcripts encoding for cement-related proteins, but also important enzymes putatively involved in the chemistry of the adhesion process, such as kinases, and may correspond to new targets to develop biomimetic approaches.

STATEMENTS OF SIGNIFICANCE

The diversity of bioadhesives elaborated by marine invertebrates is a tremendous source of inspiration to develop biomimetic approaches for biomedical and technical applications. Recent studies on the adhesion system of mussel, barnacle and sea star had highlighted the usefulness of high-throughput RNA sequencing in accelerating the development of biomimetic adhesives. Adhesion in sandcastle worms, which involves catechol and phosphate chemistries, polyelectrolyte complexes, supramolecular architectures, and a coacervation process, is a useful model to develop multipurpose wet adhesives. Using transcriptomic tools, we have explored the diversity of genes encoding for structural and catalytic proteins involved in cement formation of two sandcastle worm species, Sabellaria alveolata and Phragmatopoma caudata. The important genomic resource generated should help to design novel "blue" adhesives.