A step in the deep evolution of Alvinellidae (Annelida: Polychaeta): A phylogenomic comparative approach based on transcriptomes

The Alvinellidae are a family of worms that are endemic to deep-sea hydrothermal vents in the Pacific and Indian Oceans. These annelid worms, a sister group to the Ampharetidae, occupy a wide range of thermal habitats. The family includes the most thermotolerant marine animals described to date such as the Pompeii worm Alvinella pompejana, and other species living at much lower temperatures such as Paralvinella grasslei or Paralvinella pandorae. The phylogeny of this family has not been studied extensively. It is, however, a complex case where molecular phylogenies have given conflicting results, especially concerning the monophyletic or polyphyletic nature of the genus Paralvinella. We carried out a comprehensive study of the phylogeny of this family using the best molecular data currently available from RNAseq datasets, leading to several hundred orthologous transcripts for 11 of the 14 species currently described or in description. The results obtained by the two most popular phylogenetic inference models (using either gene concatenation with maximum likelihood, or a coalescent-based model from gene trees) were compared using a series of ampharetid and terebellid outgroups. Our study shows that the global phylogenetic signal favors the hypothesis of paraphyly for the Paralvinella genus, with P. pandorae being sister species of the other Alvinellidae. However, a high number of gene trees also supports the hypothesis of alternative trees in which the monophyly of the Paralvinella genus, as initially proposed by Desbruyères and Laubier, is valid with the species P. pandorae and Paralvinella unidentata being sister species. According to molecular dating, the radiation of the Alvinellidae was rapid and took place in a short period of time between 80 and 110 million years ago. This is reflected at the genomic scale by gene trees equally separated between different phylogenetic hypothesis, showing high rates of incomplete lineage sorting between the first lineages of the Alvinellidae and probable gene transfers. Although different genomic regions seem to have different phylogenetic stories in the early step of the alvinellid radiation, our study supports the view that the two P. pandorae species can be grouped into a separate genus (possibly Nautalvinella) and that the Miralvinella subgenus, defined by Desbruyères and Laubier, is not valid anymore.

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.

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 complex evolution of the metazoan HSP70 gene family

The metazoan 70-kDa heat shock protein (HSP70) family contains several members localized in different subcellular compartments. The cytosolic members have been classified into inducible HSP70s and constitutive heat shock cognates (HSC70s), but their distinction and evolutionary relationship remain unclear because of occasional reports of “constitutive HSP70s” and the lack of cross-phylum comparisons. Here we provide novel insights into the evolution of these important molecular chaperones. Phylogenetic analyses of 125 full-length HSP70s from a broad range of phyla revealed an ancient duplication that gave rise to two lineages from which all metazoan cytosolic HSP70s descend. One lineage (A) contains a relatively small number of genes from many invertebrate phyla, none of which have been shown to be constitutively expressed (i.e., either inducible or unknown). The other lineage (B) included both inducible and constitutive genes from diverse phyla. Species-specific duplications are present in both lineages, and Lineage B contains well-supported phylum-specific clades for Platyhelminthes, Rotifera, Nematoda, Porifera/Cnidaria, and Chordata. Some genes in Lineage B have likely independently acquired inducibility, which may explain the sporadic distribution of “HSP70” or “HSC70” in previous phylogenetic analyses. Consistent with the diversification history within each group, inducible members show lower purifying selection pressure compared to constitutive members. These results illustrate the evolutionary history of the HSP70 family, encouraging us to propose a new nomenclature: “HSP70 + subcellular localization + linage + copy number in the organism + inducible or constitutive, if known.” e.g., HSP70cA1i for cytosolic Lineage A, copy 1, inducible.

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.

Characterization and function of the first antibiotic isolated from a vent organism: the extremophile metazoan Alvinella pompejana

The emblematic hydrothermal worm Alvinella pompejana is one of the most thermo tolerant animal known on Earth. It relies on a symbiotic association offering a unique opportunity to discover biochemical adaptations that allow animals to thrive in such a hostile habitat. Here, by studying the Pompeii worm, we report on the discovery of the first antibiotic peptide from a deep-sea organism, namely alvinellacin. After purification and peptide sequencing, both the gene and the peptide tertiary structures were elucidated. As epibionts are not cultivated so far and because of lethal decompression effects upon Alvinella sampling, we developed shipboard biological assays to demonstrate that in addition to act in the first line of defense against microbial invasion, alvinellacin shapes and controls the worm's epibiotic microflora. Our results provide insights into the nature of an abyssal antimicrobial peptide (AMP) and into the manner in which an extremophile eukaryote uses it to interact with the particular microbial community of the hydrothermal vent ecosystem. Unlike earlier studies done on hydrothermal vents that all focused on the microbial side of the symbiosis, our work gives a view of this interaction from the host side.

Exploring the limit of metazoan thermal tolerance via comparative proteomics: thermally induced changes in protein abundance by two hydrothermal vent polychaetes

Temperatures around hydrothermal vents are highly variable, ranging from near freezing up to 300°C. Nevertheless, animals thrive around vents, some of which live near the known limits of animal thermotolerance. Paralvinella sulfincola, an extremely thermotolerant vent polychaete, and Paralvinella palmiformis, a cooler-adapted congener, are found along the Juan de Fuca Ridge in the northwestern Pacific. We conducted shipboard high-pressure thermotolerance experiments on both species to characterize the physiological adaptations underlying P. sulfincola's pronounced thermotolerance. Quantitative proteomics, expressed sequence tag (EST) libraries and glutathione assays revealed that P. sulfincola (i) exhibited an upregulation in the synthesis and recycling of glutathione with increasing temperature, (ii) downregulated nicotinamide adenine dinucleotide (NADH) and succinate dehydrogenases (key enzymes in oxidative phosphorylation) with increasing temperature, and (iii) maintained elevated levels of heat shock proteins (HSPs) across all treatments. In contrast, P. palmiformis exhibited more typical responses to increasing temperatures (e.g. increasing HSPs at higher temperatures). These data reveal differences in how a mesotolerant and extremely thermotolerant eukaryote respond to thermal stress, and suggest that P. sulfincola's capacity to mitigate oxidative stress via increased synthesis of antioxidants and decreased flux through the mitochondrial electron transport chain enable pronounced thermotolerance. Ultimately, oxidative stress may be the key factor in limiting all metazoan thermotolerance.

Deep-sea echinoderm oxygen consumption rates and an interclass comparison of metabolic rates in Asteroidea, Crinoidea, Echinoidea, Holothuroidea and Ophiuroidea

Echinoderms are important components of deep-sea communities because of their abundance and the fact that their activities contribute to carbon cycling. Estimating the echinoderm contribution to food webs and carbon cycling is important to our understanding of the functioning of the deep-sea environment and how this may alter in the future as climatic changes take place. Metabolic rate data from deep-sea echinoderm species are, however, scarce. To obtain such data from abyssal echinoderms, a novel in situ respirometer system, the benthic incubation chamber system (BICS), was deployed by remotely operated vehicle (ROV) at depths ranging from 2200 to 3600 m. Oxygen consumption rates were obtained in situ from four species of abyssal echinoderm (Ophiuroidea and Holothuroidea). The design and operation of two versions of BICS are presented here, together with the in situ respirometry measurements. These results were then incorporated into a larger echinoderm metabolic rate data set, which included the metabolic rates of 84 echinoderm species from all five classes (Asteroidea, Crinoidea, Echinoidea, Holothuroidea and Ophiuroidea). The allometric scaling relationships between metabolic rate and body mass derived in this study for each echinoderm class were found to vary. Analysis of the data set indicated no change in echinoderm metabolic rate with depth (by class or phylum). The allometric scaling relationships presented here provide updated information for mass-dependent deep-sea echinoderm metabolic rate for use in ecosystem models, which will contribute to the study of both shallow water and deep-sea ecosystem functioning and biogeochemistry.

Identification of differentially expressed genes in the hydrothermal vent shrimp Rimicaris exoculata exposed to heat stress

The deep-sea vent shrimp Rimicaris exoculata dominates the vagile megafauna at most vent sites along the Mid-Atlantic Ridge. This shrimp swarms around the hot end of the hydrothermal biotope where temperature can exceed its critical maximal temperature (33-38.5 ± 2°C). It may therefore be subjected to a thermal regime that is assumed to be stressful for animals. In this study, we used a global transcriptomic approach by constructing suppression subtractive hybridization cDNA libraries in order to identify specific up- and down-regulated genes in R. exoculata exposed to a severe heat stress (1h at 30°C). A total of 218 sequences representing potentially highly expressed genes in thermally stressed shrimp were obtained. Expression of 11 genes involved in various cell functions was quantified in control and heat shocked specimens using real-time PCR. Differential expression was observed for some specific genes such as mannose receptor C1, metalloprotease, histone H1, and hemocyanin with a strong up-regulation of several genes encoding heat shock proteins. These results suggest that R. exoculata is affected at both cellular and molecular levels by sustained exposure at 30°C. The sequenced ESTs presented here will provide an excellent basis for future thermal stress studies on deep-sea vent fauna.

Molecular identification of differentially regulated genes in the hydrothermal-vent species Bathymodiolus thermophilus and Paralvinella pandorae in response to temperature

BACKGROUND

Hydrothermal vents and cold seeps represent oases of life in the deep-sea environment, but are also characterized by challenging physical and chemical conditions. The effect of temperature fluctuations on vent organisms in their habitat has not been well explored, in particular at a molecular level, most gene expression studies being conducted on coastal marine species. In order to better understand the response of hydrothermal organisms to different temperature regimes, differentially expressed genes (obtained by a subtractive suppression hybridization approach) were identified in the mussel Bathymodiolus thermophilus and the annelid Paralvinella pandorae irlandei to characterize the physiological processes involved when animals are subjected to long term exposure (2 days) at two contrasting temperatures (10 degrees versus 20 degrees C), while maintained at in situ pressures. To avoid a potential effect of pressure, the experimental animals were initially thermally acclimated for 24 hours in a pressurized vessel.

RESULTS

For each species, we produced two subtractive cDNA libraries (forward and reverse) from sets of deep-sea mussels and annelids exposed together to a thermal challenge under pressure. RNA extracted from the gills, adductor muscle, mantle and foot tissue were used for B. thermophilus. For the annelid model, whole animals (small individuals) were used. For each of the four libraries, we sequenced 200 clones, resulting in 78 and 83 unique sequences in mussels and annelids (about 20% of the sequencing effort), respectively, with only half of them corresponding to known genes. Real-time PCR was used to validate differentially expressed genes identified in the corresponding libraries. Strong expression variations have been observed for some specific genes such as the intracellular hemoglobin, the nidogen protein, and Rab7 in P. pandorae, and the SPARC protein, cyclophilin, foot protein and adhesive plaque protein in B. thermophilus.

CONCLUSION

Our results indicate that mussels and worms are not responding in the same way to temperature variations. While the results obtained for the mussel B. thermophilus seem to indicate a metabolic depression (strong decrease in the level of mRNA expression of numerous genes) when temperature increased, the annelid P. pandorae mainly displayed a strong regulation of the mRNA encoding subunits and linkers of respiratory pigments and some proteins involved in membrane structure. In both cases, these regulations seem to be partly due to a possible cellular oxidative stress induced by the simulated thermal environment (10 degrees C to 20 degrees C). This work will serve as a starting point for studying the transcriptomic response of hydrothermal mussels and annelids in future experiments in response to thermal stress at various conditions of duration and temperature challenge.