Adaptation and evolution of deep-sea scale worms (Annelida: Polynoidae): insights from transcriptome comparison with a shallow-water species

Genetic adaptations of marine invertebrates to hydrothermal vent habitats

Hydrothermal vents are unique habitats like an oases of life compared with typical deep-sea, soft-sediment environments. Most animals that live in these habitats are invertebrates, and they have adapted to extreme vent environments that include high temperatures, hypoxia, high sulfide, high metal concentration, and darkness. The advent of next-generation sequencing technology, especially the coming of the new era of omics, allowed more studies to focus on the molecular adaptation of these invertebrates to vent habitats. Many genes linked to hydrothermal adaptation have been studied. We summarize the findings related to these genetic adaptations and discuss which new techniques can facilitate studies in the future.

Phylogenetics of Lepidonotopodini (Macellicephalinae, Polynoidae, Annelida) and Comparative Mitogenomics of Shallow-Water vs. Deep-Sea Scaleworms (Aphroditiformia)

Within Polynoidae, a diverse aphroditiform family, the subfamily Macellicephalinae comprises anchialine cave-dwelling and deep-sea scaleworms. In this study, Lepidonotopodinae is synonymized with Macellicephalinae, and the tribe Lepidonotopodini is applied to a well-supported clade inhabiting deep-sea chemosynthetic-based ecosystems. Newly sequenced "genome skimming" data for 30 deep-sea polynoids and the comparatively shallow living Eulagisca gigantea is used to bioinformatically assemble their mitogenomes. When analyzed with existing scaleworm mitogenomes, deep-sea scaleworms exhibit increased gene order rearrangement events compared to shallow-water relatives. Additionally, comparative analyses of shallow-water vs. deep-sea polynoid substitution rates in mitochondrial protein-coding genes show an overall relaxed purifying selection and a positive selection of several amino acid sites in deep-sea species, indicating that polynoid mitogenomes have undergone selective pressure to evolve metabolic adaptations suited to deep-sea environments. Furthermore, the inclusion of skimming data for already known Lepidonotopodini species allowed for an increased coverage of DNA data and a representation of the taxa necessary to create a more robust phylogeny using 18 genes, as opposed to the six genes previously used. The phylogenetic results support the erection of Cladopolynoe gen. nov., Mamiwata gen. nov., Photinopolynoe gen. nov., Stratigos gen. nov., and Themis gen. nov., and emended diagnoses for Branchinotogluma, Branchipolynoe, Lepidonotopodium, and Levensteiniella.

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.

Genomic Analysis of a Scale Worm Provides Insights into Its Adaptation to Deep-Sea Hydrothermal Vents

Deep-sea polynoid scale worms endemic to hydrothermal vents have evolved an adaptive strategy to the chronically hypoxic environment, but its underlying molecular mechanisms remain elusive. Here, we assembled a chromosome-scale genome of the vent-endemic scale worm Branchipolynoe longqiensis (the first annotated genome in the subclass Errantia) and annotated two shallow-water polynoid genomes, aiming to elucidate the adaptive mechanisms. We present a genome-wide molecular phylogeny of Annelida which calls for extensive taxonomy revision by including more genomes from key lineages. The B. longqiensis genome with a genome size of 1.86 Gb and 18 pseudochromosomes is larger than the genomes of two shallow-water polynoids, possibly due to the expansion of various transposable elements (TEs) and transposons. We revealed two interchromosomal rearrangements in B. longqiensis when compared with the two shallow-water polynoid genomes. The intron elongation and interchromosomal rearrangement can influence a number of biological processes, such as vesicle transport, microtubules, and transcription factors. Furthermore, the expansion of cytoskeleton-related gene families may favor the cell structure maintenance of B. longqiensis in the deep ocean. The expansion of synaptic vesicle exocytosis genes has possibly contributed to the unique complex structure of the nerve system in B. longqiensis. Finally, we uncovered an expansion of single-domain hemoglobin and a unique formation of tetra-domain hemoglobin via tandem duplications, which may be related to the adaptation to a hypoxic environment.

Adaptation and evolution of the sea anemone Alvinactis sp. to deep-sea hydrothermal vents: A comparison using transcriptomes

Sea anemones are diverse and ecologically successful members of Anthozoa. They are often found in intertidal and shallow waters, although a few of them inhabit harsher living conditions, such as deep-sea hydrothermal vents. Here, we sequenced the transcriptome of the vent sea anemone Alvinactis sp., which was collected from Edmond vent along the central Indian Ocean ridge at a depth of 3275 m, to explore the molecular mechanisms related to adaptation to vents. Compared with another deep-sea anemone (Paraphelliactis xishaensis) and five shallow water sea anemones, a total of 117 positively selected genes and 46 significantly expanded gene families were found in Alvinactis sp. specifically that may be related to its vent-specific aspect of adaptation. In addition, 127 positively selected genes and 23 significantly expanded gene families that were found in both Alvinactis sp. and P. xishaensis. Among these, vent-specific adaptations of Alvinactis sp. may involve genetic alterations in peroxisome, ubiquitin-mediated protein degradation, oxidative phosphorylation, and endocytosis, and its deep-sea adaptation may involve changes in genetic information processing. Differentially expressed genes between Alvinactis sp. and the deep-sea anemone P. xishaensis were enriched in a variety of pathways related to adaptation, such as energy metabolism, genetic information processing, endocytosis, and peroxisomes. Overall, we provided the first transcriptome of sea anemones that inhabit vents, which enriches our knowledge of deep-sea hydrothermal vent adaptation and the diversity of sea anemones.

Effects of scale worm parasitism on interactions between the symbiotic gill microbiome and gene regulation in deep sea mussel hosts

Diverse adaptations to the challenging deep sea environment are expected to be found across all deep sea organisms. Scale worms Branchipolynoe pettiboneae are believed to adapt to the deep sea environment by parasitizing deep sea mussels; this biotic interaction is one of most known in the deep sea chemosynthetic ecosystem. However, the mechanisms underlying the effects of scale worm parasitism on hosts are unclear. Previous studies have revealed that the microbiota plays an important role in host adaptability. Here, we compared gill-microbiota, gene expression and host-microorganism interactions in a group of deep sea mussels (Gigantidas haimaensis) parasitized by scale worm (PA group) and a no parasitic control group (NPA group). The symbiotic microorganism diversity of the PA group significantly decreased than NPA group, while the relative abundance of chemoautotrophic symbiotic bacteria that provide the host with organic carbon compounds significantly increased in PA. Interestingly, RNA-seq revealed that G. haimaensis hosts responded to B. pettiboneaei parasitism through significant upregulation of protein and lipid anabolism related genes, and that this parasitism may enhance host mussel nutrient anabolism but inhibit the host’s ability to absorb nutrients, thus potentially helping the parasite obtain nutrients from the host. In an integrated analysis of the interactions between changes in the microbiota and host gene dysregulation, we found an agreement between the microbiota and transcriptomic responses to B. pettiboneaei parasitism. Together, our findings provide new insights into the effects of parasite scale worms on changes in symbiotic bacteria and gene expression in deep sea mussel hosts. We explored the potential role of host-microorganism interactions between scale worms and deep sea mussels, and revealed the mechanisms through which scale worm parasitism affects hosts in deep sea chemosynthetic ecosystem.

Full-Length Transcriptome Comparison Provides Novel Insights into the Molecular Basis of Adaptation to Different Ecological Niches of the Deep-Sea Hydrothermal Vent in Alvinocaridid Shrimps

The deep-sea hydrothermal vent ecosystem is one of the extreme chemoautotrophic environments. Shinkaicaris leurokolos Kikuchi and Hashimoto, 2000, and Alvinocaris longirostris Kikuchi and Ohta, 1995, are typically co-distributed and closely related alvinocaridid shrimps in hydrothermal vent areas with different ecological niches, providing an excellent model for studying the adaptive evolution mechanism of animals in the extreme deep-sea hydrothermal vent environment. The shrimp S. leurokolos lives in close proximity to the chimney vent discharging high-temperature fluid, while A. longirostris inhabits the peripheral areas of hydrothermal vents. In this study, full-length transcriptomes of S. leurokolos and A. longirostris were generated using a combination of single-molecule real-time (SMRT) and Illumina RNA-seq technology. Expression analyses of the transcriptomes showed that among the top 30% of highly expressed genes of each species, more genes related to sulfide and heavy metal metabolism (sulfide: quinone oxidoreductase, SQR; persulfide dioxygenase, ETHE1; thiosulfate sulfurtransferase, TST, and ferritin, FRI) were specifically highly expressed in S. leurokolos, while genes involved in maintaining epibiotic bacteria or pathogen resistance (beta-1,3-glucan-binding protein, BGBP; endochitinase, CHIT; acidic mammalian chitinase, CHIA, and anti-lipopolysaccharide factors, ALPS) were highly expressed in A. longirostris. Gene family expansion analysis revealed that genes related to anti-oxidant metabolism (cytosolic manganese superoxide dismutase, SODM; glutathione S-transferase, GST, and glutathione peroxidase, GPX) and heat stress (heat shock cognate 70 kDa protein, HSP70 and heat shock 70 kDa protein cognate 4, HSP7D) underwent significant expansion in S. leurokolos, while CHIA and CHIT involved in pathogen resistance significantly expanded in A. longirostris. Finally, 66 positively selected genes (PSGs) were identified in the vent shrimp S. leurokolos. Most of the PSGs were involved in DNA repair, antioxidation, immune defense, and heat stress response, suggesting their function in the adaptive evolution of species inhabiting the extreme vent microhabitat. This study provides abundant genetic resources for deep-sea invertebrates, and is expected to lay the foundation for deep decipherment of the adaptive evolution mechanism of shrimps in a deep-sea chemosynthetic ecosystem based on further whole-genome comparison.

Phylogenomics resolves ambiguous relationships within Aciculata (Errantia, Annelida)

Aciculata (Eunicida + Phyllodocida) is among the largest clades of annelids, comprising almost half of the known diversity of all marine annelids. Despite the group's large size and biological importance, most phylogenomic studies on Annelida to date have had a limited sampling of this clade. The phylogenetic placement of many clades within Phyllodocida in particular has remained poorly understood. To resolve the relationships within Aciculata we conducted a large-scale phylogenomic analysis based on 24 transcriptomes (13 new), chosen to represent many family-ranked taxa that have never been included in a broad phylogenomic study. Our sampling also includes several enigmatic taxa with challenging phylogenetic placement, such as Histriobdella, Struwela, Lacydonia, Pilargis and the holopelagic worms Lopadorrhynchus, Travisiopsis and Tomopteris. Our robust phylogeny allows us to name and place some of these problematic clades and has significant implications on the systematics of the group. Within Eunicida we reinstate the names Eunicoidea and Oenonoidea. Within Phyllodocida we delineate Phyllodociformia, Glyceriformia, Nereidiformia, Nephtyiformia and Aphroditiformia. Phyllodociformia now includes: Lacydonia, Typhloscolecidae, Lopadorrhynchidae and Phyllodocidae. Nephtyiformia includes Nephtyidae and Pilargidae. We also broaden the delineation of Glyceriformia to include Sphaerodoridae, Tomopteridae and Glyceroidea (Glyceridae + Goniadidae). Furthermore, our study demonstrates and explores how conflicting, yet highly supported topologies can result from confounding signals in gene trees.

Muscular adaptations in swimming scale worms (Polynoidae, Annelida)

Annelids are predominantly found along with the seafloor, but over time have colonized a vast diversity of habitats, such as the water column, where different modes of locomotion are necessary. Yet, little is known about their potential muscular adaptation to the continuous swimming behaviour required in the water column. The musculature and motility were examined for five scale worm species of Polynoidae (Aphroditiformia, Annelida) found in shallow waters, deep sea or caves and which exhibit crawling, occasional swimming or continuous swimming, respectively. Their parapodial musculature was reconstructed using microCT and computational three-dimensional analyses, and the muscular functions were interpreted from video recordings of their locomotion. Since most benthic scale worms are able to swim for short distances using body and parapodial muscle movements, suitable musculature for swimming is already present. Our results indicate that rather than rearrangements or addition of muscles, a shift to a pelagic lifestyle is mainly accompanied by structural loss of muscle bundles and density, as well as elongation of extrinsic dorsal and ventral parapodial muscles. Our study documents clear differences in locomotion and musculature among closely related annelids with different lifestyles as well as points to myoanatomical adaptations for accessing the water column.

Novel perspectives of environmental proteomics

The connection among genome expression, proteome alteration, metabolism regulation and phenotype change under environmental stresses is very vague. It is a tough task for the traditional research approaches to reveal the related scientific mechanisms of the above connection at molecular and systematic levels. Proteomics approach is an insightful tool for revealing the biological functions, metabolic networks and functional protein interaction networks of cells and organisms under stresses at the systematic level. The purpose of this review is to provide an insightful guideline on how to set up a proteomic investigation for revealing biomolecule mechanisms, protein biomarkers and metabolism networks related to stress response, pollutant recognition, transport and biodegradation, and providing an insightful high-throughput approach for screening functional enzymes and effective microbes based on bioinformatics and functional verification method. Furthermore, the toxicity evaluation of pollutants and byproducts by proteomics approaches provides a scientific insight for early diagnosis of ecological risk and determination of the effectiveness of pollutant treatment techniques.

Characterization of the complete mitochondrial genome of the scale worm, Eunoe nodosa (Phyllodocida; Polynoidae) from the Beaufort Sea

To increase the mitogenome data available for robust phylogeny, we sequenced the complete mitochondrial DNA of the scale worm Eunoe nodosa (Sars, 1861) in the family Polynoidae of the order Phyllodocida. The complete mitogenome has 15,366 bp and has 28.9% A, 13.2% C, 19.0% G, and 38.8% T. Using MITOS and tRNAscan-SE, we identified the 13 typical protein-coding genes (PCGs), 2 ribosomal RNA (rRNA) genes, 22 transfer RNA (tRNA) genes, and a non-coding region. Phylogenomic analysis based on 27 in-group taxa belonging to five families of the subclass Errantia show congruence with the published phylogenetic relationship within the Polynoidae, in which E. nodosa lies in the clade of shallow water species.

Convergent Evolution and Structural Adaptation to the Deep Ocean in the Protein-Folding Chaperonin CCTα

The deep ocean is the largest biome on Earth and yet it is among the least studied environments of our planet. Life at great depths requires several specific adaptations; however, their molecular mechanisms remain understudied. We examined patterns of positive selection in 416 genes from four brittle star (Ophiuroidea) families displaying replicated events of deep-sea colonization (288 individuals from 216 species). We found consistent signatures of molecular convergence in functions related to protein biogenesis, including protein folding and translation. Five genes were recurrently positively selected, including chaperonin-containing TCP-1 subunit α (CCTα), which is essential for protein folding. Molecular convergence was detected at the functional and gene levels but not at the amino-acid level. Pressure-adapted proteins are expected to display higher stability to counteract the effects of denaturation. We thus examined in silico local protein stability of CCTα across the ophiuroid tree of life (967 individuals from 725 species) in a phylogenetically corrected context and found that deep-sea-adapted proteins display higher stability within and next to the substrate-binding region, which was confirmed by in silico global protein stability analyses. This suggests that CCTα displays not only structural but also functional adaptations to deep-water conditions. The CCT complex is involved in the folding of ∼10% of newly synthesized proteins and has previously been categorized as a "cold-shock" protein in numerous eukaryotes. We thus propose that adaptation mechanisms to cold and deep-sea environments may be linked and highlight that efficient protein biogenesis, including protein folding and translation, is a key metabolic deep-sea adaptation.

Comparative transcriptomic analysis of deep- and shallow-water barnacle species (Cirripedia, Poecilasmatidae) provides insights into deep-sea adaptation of sessile crustaceans

BACKGROUND

Barnacles are specialized marine organisms that differ from other crustaceans in possession of a calcareous shell, which is attached to submerged surfaces. Barnacles have a wide distribution, mostly in the intertidal zone and shallow waters, but a few species inhabit the deep-sea floor. It is of interest to investigate how such sessile crustaceans became adapted to extreme deep-sea environments. We sequenced the transcriptomes of a deep-sea barnacle, Glyptelasma gigas collected at a depth of 731 m from the northern area of the Zhongjiannan Basin, and a shallow-water coordinal relative, Octolasmis warwicki. The purpose of this study was to provide genetic resources for investigating adaptation mechanisms of deep-sea barnacles.

RESULTS

Totals of 62,470 and 51,585 unigenes were assembled for G. gigas and O. warwicki, respectively, and functional annotation of these unigenes was made using public databases. Comparison of the protein-coding genes between the deep- and shallow-water barnacles, and with those of four other shallow-water crustaceans, revealed 26 gene families that had experienced significant expansion in G. gigas. Functional annotation showed that these expanded genes were predominately related to DNA repair, signal transduction and carbohydrate metabolism. Base substitution analysis on the 11,611 single-copy orthologs between G. gigas and O. warwicki indicated that 25 of them were distinctly positive selected in the deep-sea barnacle, including genes related to transcription, DNA repair, ligand binding, ion channels and energy metabolism, potentially indicating their importance for survival of G. gigas in the deep-sea environment.

CONCLUSIONS

The barnacle G. gigas has adopted strategies of expansion of specific gene families and of positive selection of key genes to counteract the negative effects of high hydrostatic pressure, hypoxia, low temperature and food limitation on the deep-sea floor. These expanded gene families and genes under positive selection would tend to enhance the capacities of G. gigas for signal transduction, genetic information processing and energy metabolism, and facilitate networks for perceiving and responding physiologically to the environmental conditions in deep-sea habitats. In short, our results provide genomic evidence relating to deep-sea adaptation of G. gigas, which provide a basis for further biological studies of sessile crustaceans in the deep sea.

Insights into high-pressure acclimation: comparative transcriptome analysis of sea cucumber Apostichopus japonicus at different hydrostatic pressure exposures

Background

Global climate change is predicted to force the bathymetric migrations of shallow-water marine invertebrates. Hydrostatic pressure is proposed to be one of the major environmental factors limiting the vertical distribution of extant marine invertebrates. However, the high-pressure acclimation mechanisms are not yet fully understood.

Results

In this study, the shallow-water sea cucumber Apostichopus japonicus was incubated at 15 and 25 MPa at 15 °C for 24 h, and subjected to comparative transcriptome analysis. Nine samples were sequenced and assembled into 553,507 unigenes with a N50 length of 1204 bp. Three groups of differentially expressed genes (DEGs) were identified according to their gene expression patterns, including 38 linearly related DEGs whose expression patterns were linearly correlated with hydrostatic pressure, 244 pressure-sensitive DEGs which were up-regulated at both 15 and 25 MPa, and 257 high-pressure-induced DEGs which were up-regulated at 25 MPa but not up-regulated at 15 MPa.

Conclusions

Our results indicated that the genes and biological processes involving high-pressure acclimation are similar to those related to deep-sea adaptation. In addition to representative biological processes involving deep-sea adaptation (such as antioxidation, immune response, genetic information processing, and DNA repair), two biological processes, namely, ubiquitination and endocytosis, which can collaborate with each other and regulate the elimination of misfolded proteins, also responded to high-pressure exposure in our study. The up-regulation of these two processes suggested that high hydrostatic pressure would lead to the increase of misfolded protein synthesis, and this may result in the death of shallow-water sea cucumber under high-pressure exposure.

Insights into the strategy of micro-environmental adaptation: Transcriptomic analysis of two alvinocaridid shrimps at a hydrothermal vent

Diffusing fluid at a deep-sea hydrothermal vent creates rapid, acute physico-chemical gradients that correlate strongly with the distribution of the vent fauna. Two alvinocaridid shrimps, Alvinocaris longirostris and Shinkaicaris leurokolos occupy distinct microhabitats around these vents and exhibit different thermal preferences. S. leurokolos inhabits the central area closer to the active chimney, while A. longirostris inhabits the peripheral area. In this study, we screened candidate genes that might be involved in niche separation and microhabitat adaptation through comparative transcriptomics. The results showed that among the top 20% of overexpressed genes, gene families related to protein synthesis and structural components were much more abundant in S. leurokolos compared to A. longirostris. Moreover, 15 out of 25 genes involved in cellular carbohydrate metabolism were related to trehalose biosynthesis, versus 1 out of 5 in A. longirostris. Trehalose, a non-reducing disaccharide, is a multifunctional molecule and has been proven to act as a protectant responsible for thermotolerance in Saccharomyces cerevisiae. Putative positively selected genes involved in chitin metabolism and the immune system (lectin, serine protease and antimicrobial peptide) were enriched in S. leurokolos. In particular, one collagen and two serine proteases were found to have experienced strong positive selection. In addition, sulfotransferase-related genes were both overexpressed and positively selected in S. leurokolos. Finally, genes related to structural proteins, immune proteins and protectants were overexpressed or positively selected. These characteristics could represent adaptations of S. leurokolos to its microhabitat, which need to be confirmed by more evidence, such as data from large samples and different development stages of these alvinocaridid shrimps.

The complete mitochondrial genome of the largest amphipod, Alicella gigantea: Insight into its phylogenetic relationships and deep sea adaptive characters

Alicella gigantea (Alicelloidae) is a scavenger with the largest body size among amphipods. It is a participant in the foodweb of deepsea ecosystem and distributed with vast bathymetric and geographic ranges. In this study, the mitochondrial genome of A. gigantea was completely assembled and characterized. The complete sequence has a total length of 16,851 bp, comprising the usual eukaryotic components, with 13 protein-coding genes (PCGs), 2 ribosomal RNA genes (rRNAs), 22 transfer RNA genes (tRNAs), and 2 noncoding control regions (CRs). The gene rearrangement and reverse nucleotide strand bias of its mitochondrial genome are similar to those observed in the deepsea amphipod Eurythenes maldoror (Eurytheneidae), but different from the characters of Halice sp. MT-2017 (Dexaminoidea), an inhabitant of a deeper environment. Phylogenetic analysis indicates that A. gigantea occupies the basal branch of deepsea species-E. maldoror and Hirondellea gigas. This phylogeny supports the hypothesis that the evolution of hadal amphipods has undergone a transition from the abyssal depth. Compared to 41 available shallow water equivalents, the four accessible mitochondrial genomes from the deep sea, including the one produced in this study, show significantly fewer charged amino acids in the 13 PCGs, which suggests an adaption to the deepsea environment.

WEGO 2.0: a web tool for analyzing and plotting GO annotations, 2018 update

WEGO (Web Gene Ontology Annotation Plot), created in 2006, is a simple but useful tool for visualizing, comparing and plotting GO (Gene Ontology) annotation results. Owing largely to the rapid development of high-throughput sequencing and the increasing acceptance of GO, WEGO has benefitted from outstanding performance regarding the number of users and citations in recent years, which motivated us to update to version 2.0. WEGO uses the GO annotation results as input. Based on GO's standardized DAG (Directed Acyclic Graph) structured vocabulary system, the number of genes corresponding to each GO ID is calculated and shown in a graphical format. WEGO 2.0 updates have targeted four aspects, aiming to provide a more efficient and up-to-date approach for comparative genomic analyses. First, the number of input files, previously limited to three, is now unlimited, allowing WEGO to analyze multiple datasets. Also added in this version are the reference datasets of nine model species that can be adopted as baselines in genomic comparative analyses. Furthermore, in the analyzing processes each Chi-square test is carried out for multiple datasets instead of every two samples. At last, WEGO 2.0 provides an additional output graph along with the traditional WEGO histogram, displaying the sorted P-values of GO terms and indicating their significant differences. At the same time, WEGO 2.0 features an entirely new user interface. WEGO is available for free at http://wego.genomics.org.cn.

Transcriptomic analysis reveals insights into deep-sea adaptations of the dominant species, Shinkaia crosnieri (Crustacea: Decapoda: Anomura), inhabiting both hydrothermal vents and cold seeps

Background

Hydrothermal vents and cold seeps are typical deep-sea chemosynthetically-driven ecosystems that allow high abundance of specialized macro-benthos. To gather knowledge about the genetic basis of adaptation to these extreme environments, species shared between different habitats, especially for the dominant species, are of particular interest. The galatheid squat lobster, Shinkaia crosnieri Baba and Williams, 1998, is one of the few dominant species inhabiting both deep-sea hydrothermal vents and cold seeps. In this study, we performed transcriptome analyses of S. crosnieri collected from the Iheya North hydrothermal vent (HV) and a cold seep in the South China Sea (CS) to provide insights into how this species has evolved to thrive in different deep-sea chemosynthetic ecosystems.

Results

We analyzed 5347 orthologs between HV and CS to identify genes under positive selection through the maximum likelihood approach. A total of 82 genes were identified to be positively selected and covered diverse functional categories, potentially indicating their importance for S. crosnieri to cope with environmental heterogeneity between deep-sea vents and seeps. Among 39,806 annotated unigenes, a large number of differentially expressed genes (DEGs) were identified between HV and CS, including 339 and 206 genes significantly up-regulated in HV and CS, respectively. Most of the DEGs associated with stress response and immunity were up-regulated in HV, possibly allowing S. crosnieri to increase its capability to manage more environmental stresses in the hydrothermal vents.

Conclusions

We provide the first comprehensive transcriptomic resource for the deep-sea squat lobster, S. crosnieri, inhabiting both hydrothermal vents and cold seeps. A number of stress response and immune-related genes were positively selected and/or differentially expressed, potentially indicating their important roles for S. crosnieri to thrive in both deep-sea vents and cold seeps. Our results indicated that genetic adaptation of S. crosnieri to different deep-sea chemosynthetic environments might be mediated by adaptive evolution of functional genes related to stress response and immunity, and alterations in their gene expression that lead to different stress resistance. However, further work is required to test these proposed hypotheses. All results can constitute important baseline data for further studies towards elucidating the adaptive mechanisms in deep-sea crustaceans.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-5753-7) contains supplementary material, which is available to authorized users.

Comparative transcriptome analysis of Eogammarus possjeticus at different hydrostatic pressure and temperature exposures

Hydrostatic pressure is an important environmental factor affecting the vertical distribution of marine organisms. Laboratory-based studies have shown that many extant shallow-water marine benthic invertebrates can tolerate hydrostatic pressure outside their known natural distributions. However, only a few studies have focused on the molecular mechanisms of pressure acclimatisation. In the present work, we examined the pressure tolerance of the shallow-water amphipod Eogammarus possjeticus at various temperatures (5, 10, 15, and 20 °C) and hydrostatic pressures (0.1–30 MPa) for 16 h. Six of these experimental groups were used for transcriptome analysis. We found that 100% of E. possjeticus survived under 20 MPa at all temperature conditions for 16 h. Sequence assembly resulted in 138, 304 unigenes. Results of differential expression analysis revealed that 94 well-annotated genes were up-regulated under high pressure. All these findings indicated that the pressure tolerance of E. possjeticus was related to temperature. Several biological processes including energy metabolism, antioxidation, immunity, lipid metabolism, membrane-related process, genetic information processing, and DNA repair are probably involved in the acclimatisation in deep-sea environments.

The transcriptome of the Bermuda fireworm Odontosyllis enopla (Annelida: Syllidae): A unique luciferase gene family and putative epitoky-related genes

The Bermuda fireworm Odontosyllis enopla exhibits an extremely tight circalunar circadian behavior that results in an impressive bioluminescent mating swarm, thought to be due to a conventional luciferase-mediated oxidation of a light-emitting luciferin. In addition, the four eyes become hypertrophied and heavily pigmented, and the nephridial system is modified to store and release gametes and associated secretions. In an effort to elucidate transcripts related to bioluminescence, circadian or circalunar periodicity, as well as epitoky-related changes of the eyes and nephridial system, we examined the transcriptomic profile of three female O. enopla during a bioluminescent swarm in Ferry Reach, Bermuda. Using the well-characterized luciferase gene of the Japanese syllid Odontosyllis undecimdonta as a reference, a complete best-matching luciferase open reading frame (329 amino acids in length) was found in all three individuals analyzed in addition to numerous other paralogous sequences in this new gene family. No photoproteins were detected. We also recovered a predicted homolog of 4-coumarate-CoA ligase (268 amino acids in length) that best matched luciferase of the firefly Luciola with the best predicted template being the crystal structure of luciferase for Photinus pyralis, the common eastern firefly. A wide variety of genes associated with periodicity were recovered including predicted homologs of clock, bmal1, period, and timeless. Several genes corresponding to putative epitoky-related changes of the eyes were recovered including predicted homologs of a phototransduction gene, a retinol dehydrogenase and carotenoid isomerooxygenase as well as a visual perception related retinal rod rhodopsin-sensitive cGMP 3',5'-cyclic phosphodiesterase. Genes correlating to putative epitoky-related changes of the nephridia included predicted homologs of nephrocystin-3 and an egg-release sex peptide receptor.

First comprehensive analysis of lysine acetylation in Alvinocaris longirostris from the deep-sea hydrothermal vents

Background

Deep-sea hydrothermal vents are unique chemoautotrophic ecosystems with harsh conditions. Alvinocaris longirostris is one of the dominant crustacean species inhabiting in these extreme environments. It is significant to clarify mechanisms in their adaptation to the vents. Lysine acetylation has been known to play critical roles in the regulation of many cellular processes. However, its function in A. longirostris and even marine invertebrates remains elusive. Our study is the first, to our knowledge, to comprehensively investigate lysine acetylome in A. longirostris.

Results

In total, 501 unique acetylation sites from 206 proteins were identified by combination of affinity enrichment and high-sensitive-massspectrometer. It was revealed that Arg, His and Lys occurred most frequently at the + 1 position downstream of the acetylation sites, which were all alkaline amino acids and positively charged. Functional analysis revealed that the protein acetylation was involved in diverse cellular processes, such as biosynthesis of amino acids, citrate cycle, fatty acid degradation and oxidative phosphorylation. Acetylated proteins were found enriched in mitochondrion and peroxisome, and many stress response related proteins were also discovered to be acetylated, like arginine kinases, heat shock protein 70, and hemocyanins. In the two hemocyanins, nine acetylation sites were identified, among which one acetylation site was unique in A. longirostris when compared with other shallow water shrimps. Further studies are warranted to verify its function.

Conclusion

The lysine acetylome of A. longirostris is investigated for the first time and brings new insights into the regulation function of the lysine acetylation. The results supply abundant resources for exploring the functions of acetylation in A. longirostris and other shrimps.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-4745-3) contains supplementary material, which is available to authorized users.

Phylogeny, evolution and mitochondrial gene order rearrangement in scale worms (Aphroditiformia, Annelida)

Next-generation sequencing (NGS) has become a powerful tool in phylogenetic and evolutionary studies. Here we applied NGS to recover two ribosomal RNA genes (18S and 28S) from 16 species and 15 mitochondrial genomes from 16 species of scale worms representing six families in the suborder Aphroditiformia (Phyllodocida, Annelida), a complex group of polychaetes characterized by the presence of dorsal elytra or scales. The phylogenetic relationship of the several groups of scale worms remains unresolved due to insufficient taxon sampling and low resolution of individual gene markers. Phylogenetic tree topology based on mitochondrial genomes is comparable with that based on concatenated sequences from two mitochondrial genes (cox1 and 16S) and two ribosomal genes (18S and 28S) genes, but has higher statistical support for several clades. Our analyses show that Aphroditiformia is monophyletic, indicating the presence of elytra is an apomorphic trait. Eulepethidae and Aphroditidae together form the sister group to all other families in this suborder, whereas Acoetidae is sister to Iphionidae. Polynoidae is monophyletic, but within this family the deep-sea subfamilies Branchinotogluminae and Macellicephalinae are paraphyletic. Mitochondrial genomes in most scale-worm families have a conserved gene order, but within Polynoidae there are two novel arrangement patterns in the deep-sea clade. Mitochondrial protein-coding genes in polynoids as a whole have evolved under strong purifying selection, but substitution rates in deep-sea species are much higher than those in shallow-water species, indicating that purifying selection is relaxed in deep-sea polynoids. There are positive selected amino acids for some mitochondrial genes of the deep-sea clade, indicating they may involve in the adaption of deep-sea polynoids. Overall, our study (1) provided more evidence for reconstruction of the phylogeny of Aphroditiformia, (2) provided evidence to refute the assumption that mitochondrial gene order in Errantia is conserved, and (3) indicated that the deep-sea extreme environment may have affected the mitochondrial genome evolution rate and gene order arrangement in Polynoidae.

Evolution of Single-Domain Globins in Hydrothermal Vent Scale-Worms

Hypoxia at deep-sea hydrothermal vents represents one of the most basic challenges for metazoans, which then requires specific adaptations to acquire oxygen to meet their metabolic needs. Hydrothermal vent scale-worms (Polychaeta; Polynoidae) express large amounts of extracellular single- and multi-domain hemoglobins, in contrast with their shallow-water relatives that only possess intracellular globins in their nervous system (neuroglobins). We sequenced the gene encoding the single-domain (SD) globin from nine species of polynoids found in various vent and deep-sea reduced microhabitats (and associated constraints) to determine if the Polynoidae SD globins have been the targets of diversifying selection. Although extracellular, all the SD globins (and multi-domain ones) form a monophyletic clade that clusters within the intracellular globin group of other annelids, indicating that these hemoglobins have evolved from an intracellular myoglobin-like form. Positive selection could not be detected at the major ecological changes that the colonization of the deep-sea and hydrothermal vents represents. This suggests that no major structural modification was necessary to allow the globins to function under these conditions. The mere expression of these globins extracellularly may have been sufficiently advantageous for the polynoids living in hypoxic hydrothermal vents. Among hydrothermal vent species, positively selected amino acids were only detected in the phylogenetic lineage leading to the two mussel-commensal species (Branchipolynoe). In this lineage, the multiplicity of hemoglobins could have lessened the selective pressure on the SD hemoglobin, allowing the acquisition of novel functions by positive Darwinian selection. Conversely, the colonization of hotter environments (species of Branchinotogluma) does not seem to have required additional modifications.