Genetic sex determination in three closely related hydrothermal vent gastropods, including one species with intersex individuals

Molluscs have undergone many transitions between separate sexes and hermaphroditism, which is of interest for studying the evolution of sex determination and differentiation. Here we combined multi-locus genotypes obtained from RAD sequencing with anatomical observations of the gonads for three deep-sea hydrothermal vent gastropods of the genus Alviniconcha living in the southwest Pacific. We found that all three species (A. boucheti, A. strummeri, and A. kojimai) share the same male-heterogametic XY sex determination system, but that the gonads of XX A. kojimai individuals are invaded by a variable proportion of male reproductive tissue. The identification of Y-specific RAD loci (found only in A. boucheti) and the phylogenetic analysis of three sex-linked loci shared by all species suggested that X-Y recombination has evolved differently within each species. This situation of three species showing variation in gonadal development around a common sex determination system provides new insights into the reproductive mode of poorly known deep-sea species and opens up an opportunity to study the evolution of recombination suppression on sex chromosomes and its association with mixed or transitory sexual systems.

Highly structured populations of copepods at risk to deep-sea mining: Integration of genomic data with demogenetic and biophysical modelling

Copepoda is the most abundant taxon in deep-sea hydrothermal vents, where hard substrate is available. Despite the increasing interest in seafloor massive sulphides exploitation, there have been no population genomic studies conducted on vent meiofauna, which are known to contribute over 50% to metazoan biodiversity at vents. To bridge this knowledge gap, restriction-site-associated DNA sequencing, specifically 2b-RADseq, was used to retrieve thousands of genome-wide single-nucleotide polymorphisms (SNPs) from abundant populations of the vent-obligate copepod Stygiopontius lauensis from the Lau Basin. SNPs were used to investigate population structure, demographic histories and genotype-environment associations at a basin scale. Genetic analyses also helped to evaluate the suitability of tailored larval dispersal models and the parameterization of life-history traits that better fit the population patterns observed in the genomic dataset for the target organism. Highly structured populations were observed on both spatial and temporal scales, with divergence of populations between the north, mid, and south of the basin estimated to have occurred after the creation of the major transform fault dividing the Australian and the Niuafo'ou tectonic plate (350 kya), with relatively recent secondary contact events (<20 kya). Larval dispersal models were able to predict the high levels of structure and the highly asymmetric northward low-level gene flow observed in the genomic data. These results differ from most studies conducted on megafauna in the region, elucidating the need to incorporate smaller size when considering site prospecting for deep-sea exploitation of seafloor massive sulphides, and the creation of area-based management tools to protect areas at risk of local extinction, should mining occur.

Highly structured populations of deep-sea copepods associated with hydrothermal vents across the Southwest Pacific, despite contrasting life history traits

Hydrothermal vents are extreme environments, where abundant communities of copepods with contrasting life history traits co-exist along hydrothermal gradients. Here, we discuss how these traits may contribute to the observed differences in molecular diversity and population genetic structure. Samples were collected from vent locations across the globe including active ridges and back-arc basins and compared to existing deep-sea hydrothermal vent and shallow water data, covering a total of 22 vents and 3 non-vent sites. A total of 806 sequences of mtDNA from the Cox1 gene were used to reconstruct the phylogeny, haplotypic relationship and demography within vent endemic copepods (Dirivultidae, Stygiopontius spp.) and non-vent-endemic copepods (Ameiridae, Miraciidae and Laophontidae). A species complex within Stygiopontius lauensis was studied across five pacific back-arc basins at eight hydrothermal vent fields, with cryptic species being restricted to the basins they were sampled from. Copepod populations from the Lau, North Fiji and Woodlark basins are undergoing demographic expansion, possibly linked to an increase in hydrothermal activity in the last 10 kya. Highly structured populations of Amphiascus aff. varians 2 were also observed from the Lau to the Woodlark basins with populations also undergoing expansion. Less abundant harpacticoids exhibit little to no population structure and stable populations. This study suggests that similarities in genetic structure and demography may arise in vent-associated copepods despite having different life history traits. As structured meta-populations may be at risk of local extinction should major anthropogenic impacts, such as deep-sea mining, occur, we highlight the importance of incorporating a trait-based approach to investigate patterns of genetic connectivity and demography, particularly regarding area-based management tools and environmental management plans.

The diversification of the antimicrobial peptides from marine worms is driven by environmental conditions

Antimicrobial peptides (AMPs) play a key role in the external immunity of animals, offering an interesting model for studying the influence of the environment on the diversification and evolution of immune effectors. Alvinellacin (ALV), arenicin (ARE) and polaricin (POL, a novel AMP identified here), characterized from three marine worms inhabiting contrasted habitats ('hot' vents, temperate and polar respectively), possess a well conserved BRICHOS domain in their precursor molecule despite a profound amino acid and structural diversification of the C-terminal part containing the core peptide. Data not only showed that ARE, ALV and POL display an optimal bactericidal activity against the bacteria typical of the habitat where each worm species lives but also that this killing efficacy is optimal under the thermochemical conditions encountered by their producers in their environment. Moreover, the correlation between species habitat and the cysteine contents of POL, ARE and ALV led us to investigate the importance of disulfide bridges in their biological efficacy as a function of abiotic pressures (pH and temperature). The construction of variants using non-proteinogenic residues instead of cysteines (α-aminobutyric acid variants) leading to AMPs devoid of disulfide bridges, provided evidence that the disulfide pattern of the three AMPs allows for a better bactericidal activity and suggests an adaptive way to sustain the fluctuations of the worm's environment. This work shows that the external immune effectors exemplified here by BRICHOS AMPs are evolving under strong diversifying environmental pressures to be structurally shaped and more efficient/specific under the ecological niche of their producer.

Contrasted phylogeographic patterns of hydrothermal vent gastropods along South West Pacific: Woodlark Basin, a possible contact zone and/or stepping-stone

Understanding drivers of biodiversity patterns is essential to evaluate the potential impact of deep-sea mining on ecosystems resilience. While the South West Pacific forms an independent biogeographic province for hydrothermal vent fauna, different degrees of connectivity among basins were previously reported for a variety of species depending on their ability to disperse. In this study, we compared phylogeographic patterns of several vent gastropods across South West Pacific back-arc basins and the newly-discovered La Scala site on the Woodlark Ridge by analysing their genetic divergence using a barcoding approach. We focused on six genera of vent gastropods widely distributed in the region: Lepetodrilus, Symmetromphalus, Lamellomphalus, Shinkailepas, Desbruyeresia and Provanna. A wide-range sampling was conducted at different vent fields across the Futuna Volcanic Arc, the Manus, Woodlark, North Fiji, and Lau Basins, during the CHUBACARC cruise in 2019. The Cox1-based genetic structure of geographic populations was examined for each taxon to delineate putative cryptic species and assess potential barriers or contact zones between basins. Results showed contrasted phylogeographic patterns among species, even between closely related species. While some species are widely distributed across basins (i.e. Shinkailepas tollmanni, Desbruyeresia melanioides and Lamellomphalus) without evidence of strong barriers to gene flow, others are restricted to one (i.e. Shinkailepas tufari complex of cryptic species, Desbruyeresia cancellata and D. costata). Other species showed intermediate patterns of isolation with different lineages separating the Manus Basin from the Lau/North Fiji Basins (i.e. Lepetodrilus schrolli, Provanna and Symmetromphalus spp.). Individuals from the Woodlark Basin were either endemic to this area (though possibly representing intermediate OTUs between the Manus Basin and the other eastern basins populations) or, coming into contact from these basins, highlighting the stepping-stone role of the Woodlark Basin in the dispersal of the South West Pacific vent fauna. Results are discussed according to the dispersal ability of species and the geological history of the South West Pacific.

Genomic patterns of divergence in the early and late steps of speciation of the deep-sea vent thermophilic worms of the genus Alvinella

Background

The transient and fragmented nature of the deep-sea hydrothermal environment made of ridge subduction, plate collision and the emergence of new rifts is currently acting to separate of vent populations, promoting local adaptation and contributing to bursts of speciation and species specialization. The tube-dwelling worms Alvinella pompejana called the Pompeii worm and its sister species A. caudata live syntopically on the hottest part of deep-sea hydrothermal chimneys along the East Pacific Rise. They are exposed to extreme thermal and chemical gradients, which vary greatly in space and time, and thus represent ideal candidates for understanding the evolutionary mechanisms at play in the vent fauna evolution.

Results

We explored genomic patterns of divergence in the early and late stages of speciation of these emblematic worms using transcriptome assemblies and the first draft genome to better understand the relative role of geographic isolation and habitat preference in their genome evolution. Analyses were conducted on allopatric populations of Alvinella pompejana (early stage of separation) and between A. pompejana and its syntopic species Alvinella caudata (late stage of speciation). We first identified divergent genomic regions and targets of selection as well as their position in the genome over collections of orthologous genes and, then, described the speciation dynamics by documenting the annotation of the most divergent and/or positively selected genes involved in the isolation process. Gene mapping clearly indicated that divergent genes associated with the early stage of speciation, although accounting for nearly 30% of genes, are highly scattered in the genome without any island of divergence and not involved in gamete recognition or mito-nuclear incompatibilities. By contrast, genomes of A. pompejana and A. caudata are clearly separated with nearly all genes (96%) exhibiting high divergence. This congealing effect however seems to be linked to habitat specialization and still allows positive selection on genes involved in gamete recognition, as a possible long-duration process of species reinforcement.

Conclusion

Our analyses highlight the non-negligible role of natural selection on both the early and late stages of speciation in the iconic thermophilic worms living on the walls of deep-sea hydrothermal chimneys. They shed light on the evolution of gene divergence during the process of speciation and species specialization over a very long period of time.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12862-022-02057-y.

Origin, diversity, and biogeography of Antarctic scale worms (Polychaeta: Polynoidae): a wide-scale barcoding approach

The Antarctic marine environment hosts diversified and highly endemic benthos owing to its unique geologic and climatic history. Current warming trends have increased the urgency of understanding Antarctic species history to predict how environmental changes will impact ecosystem functioning. Antarctic benthic lineages have traditionally been examined under three hypotheses: (1) high endemism and local radiation, (2) emergence of deep‐sea taxa through thermohaline circulation, and (3) species migrations across the Polar Front. In this study, we investigated which hypotheses best describe benthic invertebrate origins by examining Antarctic scale worms (Polynoidae). We amassed 691 polynoid sequences from the Southern Ocean and neighboring areas: the Kerguelen and Tierra del Fuego (South America) archipelagos, the Indian Ocean, and waters around New Zealand. We performed phylogenetic reconstructions to identify lineages across geographic regions, aided by mitochondrial markers cytochrome c oxidase subunit I (Cox1) and 16S ribosomal RNA (16S). Additionally, we produced haplotype networks at the species scale to examine genetic diversity, biogeographic separations, and past demography. The Cox1 dataset provided the most illuminating insights into the evolution of polynoids, with a total of 36 lineages identified. Eunoe sp. was present at Tierra del Fuego and Kerguelen, in favor of the latter acting as a migration crossroads. Harmothoe fuligineum, widespread around the Antarctic continent, was also present but isolated at Kerguelen, possibly resulting from historical freeze–thaw cycles. The genus Polyeunoa appears to have diversified prior to colonizing the continent, leading to the co‐occurrence of at least three cryptic species around the Southern and Indian Oceans. Analyses identified that nearly all populations are presently expanding following a bottleneck event, possibly caused by habitat reduction from the last glacial episodes. Findings support multiple origins for contemporary Antarctic polynoids, and some species investigated here provide information on ancestral scenarios of (re)colonization. First, it is apparent that species collected from the Antarctic continent are endemic, as the absence of closely related species in the Kerguelen and Tierra del Fuego datasets for most lineages argues in favor of Hypothesis 1 of local origin. Next, Eunoe sp. and H. fuligineum, however, support the possibility of Kerguelen and other sub‐Antarctic islands acting as a crossroads for larvae of some species, in support of Hypothesis 3. Finally, the genus Polyeunoa, conversely, is found at depths greater than 150 m and may have a deep origin, in line with Hypothesis 2. These “non endemic” groups, nevertheless, have a distribution that is either north or south of the Antarctic Polar Front, indicating that there is still a barrier to dispersal, even in the deep sea.

Global 16S rRNA diversity of provannid snail endosymbionts from Indo-Pacific deep-sea hydrothermal vents

Symbioses between invertebrate animals and chemosynthetic bacteria build the foundation of deep-sea hydrothermal ecosystems worldwide. Despite the importance of these symbioses for ecosystem functioning, the diversity of symbionts within and between host organisms and geographic regions is still poorly understood. In this study we used 16S rRNA amplicon sequencing to determine the diversity of gill endosymbionts in provannid snails of the genera Alviniconcha and Ifremeria, which are key species at deep-sea hydrothermal vents in the Indo-Pacific Ocean. Our analysis of 761 snail samples across the distributional range of these species confirms previous findings that symbiont lineages are strongly partitioned by host species and broad-scale geography. Less structuring was observed within geographic regions, probably due to insufficient strain resolution of the 16S rRNA gene. Symbiont richness in individual hosts appeared to be unrelated to host size, suggesting that provannid snails might acquire their symbionts only during a permissive time window in early developmental stages in contrast to other vent molluscs that obtain their symbionts throughout their lifetime. Despite the extent of our dataset, symbiont accumulation curves did not reach saturation, highlighting the need for increased sampling efforts to uncover the full diversity of symbionts within these and other hydrothermal vent species.

Subtle limits to connectivity revealed by outlier loci within two divergent metapopulations of the deep-sea hydrothermal gastropod Ifremeria nautilei

Hydrothermal vents form archipelagos of ephemeral deep-sea habitats that raise interesting questions about the evolution and dynamics of the associated endemic fauna, constantly subject to extinction-recolonization processes. These metal-rich environments are coveted for the mineral resources they harbor, thus raising recent conservation concerns. The evolutionary fate and demographic resilience of hydrothermal species strongly depend on the degree of connectivity among and within their fragmented metapopulations. In the deep sea, however, assessing connectivity is difficult and usually requires indirect genetic approaches. Improved detection of fine-scale genetic connectivity is now possible based on genome-wide screening for genetic differentiation. Here, we explored population connectivity in the hydrothermal vent snail Ifremeria nautilei across its species range encompassing five distinct back-arc basins in the Southwest Pacific. The global analysis, based on 10 570 single nucleotide polymorphism (SNP) markers derived from double digest restriction-site associated DNA sequencing (ddRAD-seq), depicted two semi-isolated and homogeneous genetic clusters. Demo-genetic modeling suggests that these two groups began to diverge about 70 000 generations ago, but continue to exhibit weak and slightly asymmetrical gene flow. Furthermore, a careful analysis of outlier loci showed subtle limitations to connectivity between neighboring basins within both groups. This finding indicates that migration is not strong enough to totally counterbalance drift or local selection, hence questioning the potential for demographic resilience at this latter geographical scale. These results illustrate the potential of large genomic datasets to understand fine-scale connectivity patterns in hydrothermal vents and the deep sea.

Investigation of Capitella spp. symbionts in the context of varying anthropic pressures: First occurrence of a transient advantageous epibiosis with the giant bacteria Thiomargarita sp. to survive seasonal increases of sulfides in sediments

Capitella spp. is considered as an important ecological indicator of eutrophication due to its high densities in organic-rich, reduced, and sometimes polluted coastal ecosystems. We investigated whether such ability to cope with adverse ecological contexts might be a response to the microorganisms these worms are associated with. In populations from the French Atlantic coast (Roscoff, Brittany), we observed an epibiotic association covering the tegument of 20-30% specimens from an anthropized site while individuals from a reference, non-anthropized site were devoid of any visible epibionts. Using RNAseq, molecular and microscopic analyses, we described and compared the microbial communities associated with the epibiotic versus the non-epibiotic specimens at both locations. Interestingly, data showed that the epibiosis is characterized by sulfur-oxidizing bacteria among which the giant bacterium Thiomargarita sp., to date only described in deep sea habitats. Survey of Capitella combined with the geochemical analysis of their sediment revealed that epibiotic specimens are always found in muds with the highest concentration of sulfides, mostly during the summer. Concomitantly, tolerance tests demonstrated that the acquisition of epibionts increased survival against toxic level of sulfides. Overall, the present data highlight for the first time a peculiar plastic adaptation to seasonal variations of the habitat based on a transcient epibiosis allowing a coastal species to survive temporary harsher conditions.

Comparative genomic and transcriptomic analyses of transposable elements in polychaetous annelids highlight LTR retrotransposon diversity and evolution

Background

With the expansion of high throughput sequencing, we now have access to a larger number of genome-wide studies analyzing the Transposable elements (TEs) composition in a wide variety of organisms. However, genomic analyses often remain too limited in number and diversity of species investigated to study in depth the dynamics and evolutionary success of the different types of TEs among metazoans. Therefore, we chose to investigate the use of transcriptomes to describe the diversity of TEs in phylogenetically related species by conducting the first comparative analysis of TEs in two groups of polychaetes and evaluate the diversity of TEs that might impact genomic evolution as a result of their mobility.

Results

We present a detailed analysis of TEs distribution in transcriptomes extracted from 15 polychaetes depending on the number of reads used during assembly, and also compare these results with additional TE scans on associated low-coverage genomes. We then characterized the clades defined by 1021 LTR-retrotransposon families identified in 26 species. Clade richness was highly dependent on the considered superfamily. Copia elements appear rare and are equally distributed in only three clades, GalEa, Hydra and CoMol. Among the eight BEL/Pao clades identified in annelids, two small clades within the Sailor lineage are new for science. We characterized 17 Gypsy clades of which only 4 are new; the C-clade largely dominates with a quarter of the families. Finally, all species also expressed for the majority two distinct transcripts encoding PIWI proteins, known to be involved in control of TEs mobilities.

Conclusions

This study shows that the use of transcriptomes assembled from 40 million reads was sufficient to access to the diversity and proportion of the transposable elements compared to those obtained by low coverage sequencing. Among LTR-retrotransposons Gypsy elements were unequivocally dominant but results suggest that the number of Gypsy clades, although high, may be more limited than previously thought in metazoans. For BEL/Pao elements, the organization of clades within the Sailor lineage appears more difficult to establish clearly. The Copia elements remain rare and result from the evolutionary consistent success of the same three clades.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13100-021-00252-0.

The genome of Ectocarpus subulatus - A highly stress-tolerant brown alga

Brown algae are multicellular photosynthetic stramenopiles that colonize marine rocky shores worldwide. Ectocarpus sp. Ec32 has been established as a genomic model for brown algae. Here we present the genome and metabolic network of the closely related species, Ectocarpus subulatus Kützing, which is characterized by high abiotic stress tolerance. Since their separation, both strains show new traces of viral sequences and the activity of large retrotransposons, which may also be related to the expansion of a family of chlorophyll-binding proteins. Further features suspected to contribute to stress tolerance include an expanded family of heat shock proteins, the reduction of genes involved in the production of halogenated defence compounds, and the presence of fewer cell wall polysaccharide-modifying enzymes. Overall, E. subulatus has mainly lost members of gene families down-regulated in low salinities, and conserved those that were up-regulated in the same condition. However, 96% of genes that differed between the two examined Ectocarpus species, as well as all genes under positive selection, were found to encode proteins of unknown function. This underlines the uniqueness of brown algal stress tolerance mechanisms as well as the significance of establishing E. subulatus as a comparative model for future functional studies.

Worms' Antimicrobial Peptides

Antimicrobial peptides (AMPs) are natural antibiotics produced by all living organisms. In metazoans, they act as host defense factors by eliminating microbial pathogens. But they also help to select the colonizing bacterial symbionts while coping with specific environmental challenges. Although many AMPs share common structural characteristics, for example having an overall size between 10-100 amino acids, a net positive charge, a γ-core motif, or a high content of cysteines, they greatly differ in coding sequences as a consequence of multiple parallel evolution in the face of pathogens. The majority of AMPs is specific of certain taxa or even typifying species. This is especially the case of annelids (ringed worms). Even in regions with extreme environmental conditions (polar, hydrothermal, abyssal, polluted, etc.), worms have colonized all habitats on Earth and dominated in biomass most of them while co-occurring with a large number and variety of bacteria. This review surveys the different structures and functions of AMPs that have been so far encountered in annelids and nematodes. It highlights the wide diversity of AMP primary structures and their originality that presumably mimics the highly diverse life styles and ecology of worms. From the unique system that represents marine annelids, we have studied the effect of abiotic pressures on the selection of AMPs and demonstrated the promising sources of antibiotics that they could constitute.

Mollusc genomes reveal variability in patterns of LTR-retrotransposons dynamics

Background

The three superfamilies of Long Terminal Repeat (LTR) retrotransposons are a widespread kind of transposable element and a major factor in eukaryotic genome evolution. In metazoans, recent studies suggested that Copia LTR-retrotransposons display specific dynamic compared to the more abundant and diverse Gypsy elements. Indeed, Copia elements show a relative scarcity and the prevalence of only a few clades in specific hosts. Thus, BEL/Pao seems to be the second most abundant superfamily. However, the generality of these assumptions remains to be assessed. Therefore, we carried out the first large-scale comparative genomic analysis of LTR-retrotransposons in molluscs. The aim of this study was to analyse the diversity, copy numbers, genomic proportions and distribution of LTR-retrotransposons in a large host phylum.

Results

We compare nine genomes of molluscs and further added LTR-retrotransposons sequences detected in databases for 47 additional species. We identified 1709 families, which enabled us to define 31 clades. We show that clade richness was highly dependent on the considered superfamily. We found only three Copia clades, including GalEa and Hydra which appear to be widely distributed and highly dominant as they account for 96% of the characterised Copia elements. Among the seven BEL/Pao clades identified, Sparrow and Surcouf are characterised for the first time. We find no BEL or Pao elements, but the rare clades Dan and Flow are present in molluscs. Finally, we characterised 21 Gypsy clades, only five of which had been previously described, the C-clade being the most abundant one. Even if they are found in the same number of host species, Copia elements are clearly less abundant than BEL/Pao elements in copy number or genomic proportions, while Gypsy elements are always the most abundant ones whatever the parameter considered.

Conclusions

Our analysis confirms the contrasting dynamics of Copia and Gypsy elements in metazoans and indicates that BEL/Pao represents the second most abundant superfamily, probably reflecting an intermediate dynamic. Altogether, the data obtained in several taxa highly suggest that these patterns can be generalised for most metazoans. Finally, we highlight the importance of using database information in complement of genome analyses when analyzing transposable element diversity.

Electronic supplementary material

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

Proteome Evolution of Deep-Sea Hydrothermal Vent Alvinellid Polychaetes Supports the Ancestry of Thermophily and Subsequent Adaptation to Cold in Some Lineages

Temperature, perhaps more than any other environmental factor, is likely to influence the evolution of all organisms. It is also a very interesting factor to understand how genomes are shaped by selection over evolutionary timescales, as it potentially affects the whole genome. Among thermophilic prokaryotes, temperature affects both codon usage and protein composition to increase the stability of the transcriptional/translational machinery, and the resulting proteins need to be functional at high temperatures. Among eukaryotes less is known about genome evolution, and the tube-dwelling worms of the family Alvinellidae represent an excellent opportunity to test hypotheses about the emergence of thermophily in ectothermic metazoans. The Alvinellidae are a group of worms that experience varying thermal regimes, presumably having evolved into these niches over evolutionary times. Here we analyzed 423 putative orthologous loci derived from 6 alvinellid species including the thermophilic Alvinella pompejana and Paralvinella sulfincola. This comparative approach allowed us to assess amino acid composition, codon usage, divergence, direction of residue changes and the strength of selection along the alvinellid phylogeny, and to design a new eukaryotic thermophilic criterion based on significant differences in the residue composition of proteins. Contrary to expectations, the alvinellid ancestor of all present-day species seems to have been thermophilic, a trait subsequently maintained by purifying selection in lineages that still inhabit higher temperature environments. In contrast, lineages currently living in colder habitats likely evolved under selective relaxation, with some degree of positive selection for low-temperature adaptation at the protein level.

HYDROTHERMAL-VENT ALVINELLID POLYCHAETE DISPERSAL IN THE EASTERN PACIFIC. 2. A METAPOPULATION MODEL BASED ON HABITAT SHIFTS

Marine organisms typically fall into two main categories: those with a high level of population structuring and those with a low one. The first are often found to be poor dispersers, following isolation by distance or stepping-stone theoretical predictions. The second are commonly associated with high-dispersal taxa and are best described by the island model. Deep-sea hydrothermal vent systems represent a good model for studying one-dimensional metapopulations. Whereas isolation by distance might be expected to be the rule in such a system for species with limited dispersal capabilities, a biological paradox can be observed: an apparent genetic homogeneity in some vent species with short-scale dispersal potential, in a one-dimensional fragmented habitat. This can be explained if one key assumption of the existing models is not met: gene flow between populations and genetic drift may not have the time to equilibrate. Geophysical models revealed that hydrothermal convection is intrinsically unstable, inducing processes of coalescence or splitting of venting areas in a chaotic manner. This is likely to generate frequent extinctions and recolonizations. Theoretical genetic predictions derived from extinctions/recolonizations cannot satisfactorily model a situation where habitat shifts are frequent and constantly affect the metapopulation equilibrium. Because neither the island and the stepping-stone models nor the classical metapopulation models resemble the hydrothermal vent reality, we present here a realistic model developed to provide a compromise between existing conceptual models and what is currently known of the biology and ecology of one of the most peculiar and best-studied vent species, the polychaete Alvinella pompejana. This model allows us to define the boundaries between which the metapopulation is evolutionary stable in an unstable context. Simulations show different patterns in which metapopulation size and recolonization vary but reach an equilibrium despite chaotic vent extinctions. In contrast, the model also shows that displacing habitat continuously affects the equilibrium between gene flow and drift. As a consequence, the time required to balance these evolutionary forces can never be attained, leading to chaotic fluctuations in F-statistics. Those fluctuations are mainly due to stochastic changes of the interpatch distance which affect migration rates. The shifting of active zones of venting can episodically counterbalance differentiation and allow a long-term genetic homogenization at the ridge scale. These findings lead to a new concept in which the exchanges between populations would mainly depend on the habitat's movements along the ridge axis rather than the organim's dispersal. We therefore propose a new model based on patch-network displacements in which transient contact zones allow low levels of gene flow throughout the metapopulation.

Antagonistic evolution of an antibiotic and its molecular chaperone: how to maintain a vital ectosymbiosis in a highly fluctuating habitat

Evolution of antimicrobial peptides (AMPs) has been shown to be driven by recurrent duplications and balancing/positive selection in response to new or altered bacterial pathogens. We use Alvinella pompejana, the most eurythermal animal known on Earth, to decipher the selection patterns acting on AMP in an ecological rather than controlled infection approach. The preproalvinellacin multigenic family presents the uniqueness to encode a molecular chaperone (BRICHOS) together with an AMP (alvinellacin) that controls the vital ectosymbiosis of Alvinella. In stark contrast to what is observed in the context of the Red queen paradigm, we demonstrate that exhibiting a vital and highly conserved ecto-symbiosis in the face of thermal fluctuations has led to a peculiar selective trend promoting the adaptive diversification of the molecular chaperone of the AMP, but not of the AMP itself. Because BRICHOS stabilizes beta-stranded peptides, this polymorphism likely represents an eurythermal adaptation to stabilize the structure of alvinellacin, thus hinting at its efficiency to select and control the epibiosis across the range of temperatures experienced by the worm; Our results fill some knowledge gaps concerning the function of BRICHOS in invertebrates and offer perspectives for studying immune genes in an evolutionary ecological framework.

Current hypotheses to explain genetic chaos under the sea

Chaotic genetic patchiness (CGP) refers to surprising patterns of spatial and temporal genetic structure observed in some marine species at a scale where genetic variation should be efficiently homogenized by gene flow via larval dispersal. Here we review and discuss 4 mechanisms that could generate such unexpected patterns: selection, sweepstakes reproductive success, collective dispersal, and temporal shifts in local population dynamics. First, we review examples where genetic differentiation at specific loci was driven by diversifying selection, which was historically the first process invoked to explain CGP. Second, we turn to neutral demographic processes that may drive genome-wide effects, and whose effects on CGP may be enhanced when they act together. We discuss how sweepstakes reproductive success accelerates genetic drift and can thus generate genetic structure, provided that gene flow is not too strong. Collective dispersal is another mechanism whereby genetic structure can be maintained regardless of dispersal intensity, because it may prevent larval cohorts from becoming entirely mixed. Theoretical analyses of both the sweepstakes and the collective dispersal ideas are presented. Finally, we discuss an idea that has received less attention than the other ones just mentioned, namely temporal shifts in local population dynamics.

Biophysical and Population Genetic Models Predict the Presence of "Phantom" Stepping Stones Connecting Mid-Atlantic Ridge Vent Ecosystems

Deep-sea hydrothermal vents are patchily distributed ecosystems inhabited by specialized animal populations that are textbook meta-populations. Many vent-associated species have free-swimming, dispersive larvae that can establish connections between remote populations. However, connectivity patterns among hydrothermal vents are still poorly understood because the deep sea is undersampled, the molecular tools used to date are of limited resolution, and larval dispersal is difficult to measure directly. A better knowledge of connectivity is urgently needed to develop sound environmental management plans for deep-sea mining. Here, we investigated larval dispersal and contemporary connectivity of ecologically important vent mussels (Bathymodiolus spp.) from the Mid-Atlantic Ridge by using high-resolution ocean modeling and population genetic methods. Even when assuming a long pelagic larval duration, our physical model of larval drift suggested that arrival at localities more than 150 km from the source site is unlikely and that dispersal between populations requires intermediate habitats ("phantom" stepping stones). Dispersal patterns showed strong spatiotemporal variability, making predictions of population connectivity challenging. The assumption that mussel populations are only connected via additional stepping stones was supported by contemporary migration rates based on neutral genetic markers. Analyses of population structure confirmed the presence of two southern and two hybridizing northern mussel lineages that exhibited a substantial, though incomplete, genetic differentiation. Our study provides insights into how vent animals can disperse between widely separated vent habitats and shows that recolonization of perturbed vent sites will be subject to chance events, unless connectivity is explicitly considered in the selection of conservation areas.

Divergent ecological histories of two sister Antarctic krill species led to contrasted patterns of genetic diversity in their heat-shock protein (hsp70) arsenal

The Arctic and the Antarctic Peninsula are currently experiencing some of the most rapid rates of ocean warming on the planet. This raises the question of how the initial adaptation to extreme cold temperatures was put in place and whether or not directional selection has led to the loss of genetic variation at key adaptive systems, and thus polar species’ (re)adaptability to higher temperatures. In the Southern Ocean, krill represents the most abundant fauna and is a critical member at the base of the Antarctic food web. To better understand the role of selection in shaping current patterns of polymorphisms, we examined genetic diversity of the cox‐1 and hsp70 genes by comparing two closely related species of Euphausiid that differ in ecology. Results on mtcox‐1 agreed with previous studies, indicating high and similar effective population sizes. However, a coalescent‐based approach on hsp70 genes highlighted the role of positive selection and past demographic changes in their recent evolution. Firstly, some form of balancing selection was acting on the inducible isoform C, which reflected the maintenance of an ancestral adaptive polymorphism in both species. Secondly, E. crystallorophias seems to have lost most of its hsp70 diversity because of a population crash and/or directional selection to cold. Nonsynonymous diversities were always greater in E. superba, suggesting that it might have evolved under more heterogeneous conditions. This can be linked to species’ ecology with E. superba living in more variable pelagic conditions, while E. crystallorophias is strictly associated with continental shelves and sea ice.

Selective forces acting during multi-domain protein evolution: the case of multi-domain globins

Multi-domain proteins form the majority of proteins in eukaryotes. During their formation by tandem duplication or gene fusion, new interactions between domains may arise as a result of the structurally-forced proximity of domains. The proper function of the formed proteins likely required the molecular adjustment of these stress zones by specific amino acid replacements, which should be detectable by the molecular signature of selection that governed their changes. We used multi-domain globins from three different invertebrate lineages to investigate the selective forces that acted throughout the evolution of these molecules. In the youngest of these molecules [Branchipolynoe scaleworm; original duplication ca. 60 million years (Ma)], we were able to detect some amino acids under positive selection corresponding to the initial duplication event. In older lineages (didomain globin from bivalve mollusks and nematodes), there was no evidence of amino acid positions under positive selection, possibly the result of accumulated non-adaptative mutations since the original duplication event (165 and 245 Ma, respectively). Some amino acids under positive selection were sometimes detected in later branches, either after speciation events, or after the initial duplication event. In Branchipolynoe, the position of the amino acids under positive selection on a 3D model suggests some of them are located at the interface between two domains; while others are locate in the heme pocket.

Diversification, evolution and sub-functionalization of 70kDa heat-shock proteins in two sister species of antarctic krill: differences in thermal habitats, responses and implications under climate change

BACKGROUND

A comparative thermal tolerance study was undertaken on two sister species of Euphausiids (Antarctic krills) Euphausia superba and Euphausia crystallorophias. Both are essential components of the Southern Ocean ecosystem, but occupy distinct environmental geographical locations with slightly different temperature regimes. They therefore provide a useful model system for the investigation of adaptations to thermal tolerance.

METHODOLOGY/PRINCIPAL FINDING

Initial CTmax studies showed that E. superba was slightly more thermotolerant than E. crystallorophias. Five Hsp70 mRNAs were characterized from the RNAseq data of both species and subsequent expression kinetics studies revealed notable differences in induction of each of the 5 orthologues between the two species, with E. crystallorophias reacting more rapidly than E. superba. Furthermore, analyses conducted to estimate the evolutionary rates and selection strengths acting on each gene tended to support the hypothesis that diversifying selection has contributed to the diversification of this gene family, and led to the selective relaxation on the inducible C form with its possible loss of function in the two krill species.

CONCLUSIONS

The sensitivity of the epipelagic species E. crystallorophias to temperature variations and/or its adaptation to cold is enhanced when compared with its sister species, E. superba. These results indicate that ice krill could be the first of the two species to be impacted by the warming of coastal waters of the Austral ocean in the coming years due to climate change.

Transcriptomic response of the hydrothermal mussel Bathymodiolus azoricus in experimental exposure to heavy metals is modulated by the Pgm genotype and symbiont content

Hydrothermal vent mussels belonging to the genus Bathymodiolus dominate communities at hydrothermal sites of the Mid-Atlantic Ridge. The mussel Bathymodiolus azoricus harbors thiotrophic and methanotrophic symbiotic bacteria in its gills and evolves in naturally highly metal contaminated environments. In the context of investigations on metal tolerance/effect in B. azoricus, we focused our work on the short-term adaptive response (15days) of mussels to different metals exposure at a molecular level using metal concentrations chosen to mimic natural situations at three vents sites. The expression of a set of 38 genes involved in different steps of the metal uptake, detoxication and various metabolisms was analysed by qPCR. Mussels were also genotyped at 10 enzyme loci to explore the relationships among natural genetic variation and gene expression. Relation between symbiont content (both sulfur-oxidizing and methanogen bacteria) and gene expression was also analysed. Our study demonstrated the influence of metal cocktail composition and time exposure on the transcriptome regulation with a specific pattern of regulation observed for the three metal cocktail tested. We also evidenced the significant influence of some specific Pgm genotype on the global gene expression in our experimental populations and a general trend of a higher gene expression in individuals carrying a high symbiont content.

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.

A new barrier to dispersal trapped old genetic clines that escaped the Easter Microplate tension zone of the Pacific vent mussels

Comparative phylogeography of deep-sea hydrothermal vent species has uncovered several genetic breaks between populations inhabiting northern and southern latitudes of the East Pacific Rise. However, the geographic width and position of genetic clines are variable among species. In this report, we further characterize the position and strength of barriers to gene flow between populations of the deep-sea vent mussel Bathymodiolus thermophilus. Eight allozyme loci and DNA sequences of four nuclear genes were added to previously published sequences of the cytochrome c oxidase subunit I gene. Our data confirm the presence of two barriers to gene flow, one located at the Easter Microplate (between 21°33′S and 31°S) recently described as a hybrid zone, and the second positioned between 7°25′S and 14°S with each affecting different loci. Coalescence analysis indicates a single vicariant event at the origin of divergence between clades for all nuclear loci, although the clines are now spatially discordant. We thus hypothesize that the Easter Microplate barrier has recently been relaxed after a long period of isolation and that some genetic clines have escaped the barrier and moved northward where they have subsequently been trapped by a reinforcing barrier to gene flow between 7°25′S and 14°S.

Thermal limit for metazoan life in question: in vivo heat tolerance of the Pompeii worm

The thermal limit for metazoan life, expected to be around 50°C, has been debated since the discovery of the Pompeii worm Alvinella pompejana, which colonizes black smoker chimney walls at deep-sea vents. While indirect evidence predicts body temperatures lower than 50°C, repeated in situ temperature measurements depict an animal thriving at temperatures of 60°C and more. This controversy was to remain as long as this species escaped in vivo investigations, due to irremediable mortalities upon non-isobaric sampling. Here we report from the first heat-exposure experiments with live A. pompejana, following isobaric sampling and subsequent transfer in a laboratory pressurized aquarium. A prolonged (2 hours) exposure in the 50-55°C range was lethal, inducing severe tissue damages, cell mortalities and triggering a heat stress response, therefore showing that Alvinella's upper thermal limit clearly is below 55°C. A comparison with hsp70 stress gene expressions of individuals analysed directly after sampling in situ confirms that Alvinella pompejana does not experience long-term exposures to temperature above 50°C in its natural environment. The thermal optimum is nevertheless beyond 42°C, which confirms that the Pompeii worm ranks among the most thermotolerant metazoans.

Deep transcriptome-sequencing and proteome analysis of the hydrothermal vent annelid Alvinella pompejana identifies the CvP-bias as a robust measure of eukaryotic thermostability

BACKGROUND

Alvinella pompejana is an annelid worm that inhabits deep-sea hydrothermal vent sites in the Pacific Ocean. Living at a depth of approximately 2500 meters, these worms experience extreme environmental conditions, including high temperature and pressure as well as high levels of sulfide and heavy metals. A. pompejana is one of the most thermotolerant metazoans, making this animal a subject of great interest for studies of eukaryotic thermoadaptation.

RESULTS

In order to complement existing EST resources we performed deep sequencing of the A. pompejana transcriptome. We identified several thousand novel protein-coding transcripts, nearly doubling the sequence data for this annelid. We then performed an extensive survey of previously established prokaryotic thermoadaptation measures to search for global signals of thermoadaptation in A. pompejana in comparison with mesophilic eukaryotes. In an orthologous set of 457 proteins, we found that the best indicator of thermoadaptation was the difference in frequency of charged versus polar residues (CvP-bias), which was highest in A. pompejana. CvP-bias robustly distinguished prokaryotic thermophiles from prokaryotic mesophiles, as well as the thermophilic fungus Chaetomium thermophilum from mesophilic eukaryotes. Experimental values for thermophilic proteins supported higher CvP-bias as a measure of thermal stability when compared to their mesophilic orthologs. Proteome-wide mean CvP-bias also correlated with the body temperatures of homeothermic birds and mammals.

CONCLUSIONS

Our work extends the transcriptome resources for A. pompejana and identifies the CvP-bias as a robust and widely applicable measure of eukaryotic thermoadaptation.

Proteome adaptation to high temperatures in the ectothermic hydrothermal vent Pompeii worm

Taking advantage of the massive genome sequencing effort made on thermophilic prokaryotes, thermal adaptation has been extensively studied by analysing amino acid replacements and codon usage in these unicellular organisms. In most cases, adaptation to thermophily is associated with greater residue hydrophobicity and more charged residues. Both of these characteristics are positively correlated with the optimal growth temperature of prokaryotes. In contrast, little information has been collected on the molecular 'adaptive' strategy of thermophilic eukaryotes. The Pompeii worm A. pompejana, whose transcriptome has recently been sequenced, is currently considered as the most thermotolerant eukaryote on Earth, withstanding the greatest thermal and chemical ranges known. We investigated the amino-acid composition bias of ribosomal proteins in the Pompeii worm when compared to other lophotrochozoans and checked for putative adaptive changes during the course of evolution using codon-based Maximum likelihood analyses. We then provided a comparative analysis of codon usage and amino-acid replacements from a greater set of orthologous genes between the Pompeii worm and Paralvinella grasslei, one of its closest relatives living in a much cooler habitat. Analyses reveal that both species display the same high GC-biased codon usage and amino-acid patterns favoring both positively-charged residues and protein hydrophobicity. These patterns may be indicative of an ancestral adaptation to the deep sea and/or thermophily. In addition, the Pompeii worm displays a set of amino-acid change patterns that may explain its greater thermotolerance, with a significant increase in Tyr, Lys and Ala against Val, Met and Gly. Present results indicate that, together with a high content in charged residues, greater proportion of smaller aliphatic residues, and especially alanine, may be a different path for metazoans to face relatively 'high' temperatures and thus a novelty in thermophilic metazoans.

Insights into metazoan evolution from Alvinella pompejana cDNAs

Background

Alvinella pompejana is a representative of Annelids, a key phylum for evo-devo studies that is still poorly studied at the sequence level. A. pompejana inhabits deep-sea hydrothermal vents and is currently known as one of the most thermotolerant Eukaryotes in marine environments, withstanding the largest known chemical and thermal ranges (from 5 to 105°C). This tube-dwelling worm forms dense colonies on the surface of hydrothermal chimneys and can withstand long periods of hypo/anoxia and long phases of exposure to hydrogen sulphides. A. pompejana specifically inhabits chimney walls of hydrothermal vents on the East Pacific Rise. To survive, Alvinella has developed numerous adaptations at the physiological and molecular levels, such as an increase in the thermostability of proteins and protein complexes. It represents an outstanding model organism for studying adaptation to harsh physicochemical conditions and for isolating stable macromolecules resistant to high temperatures.

Results

We have constructed four full length enriched cDNA libraries to investigate the biology and evolution of this intriguing animal. Analysis of more than 75,000 high quality reads led to the identification of 15,858 transcripts and 9,221 putative protein sequences. Our annotation reveals a good coverage of most animal pathways and networks with a prevalence of transcripts involved in oxidative stress resistance, detoxification, anti-bacterial defence, and heat shock protection. Alvinella proteins seem to show a slow evolutionary rate and a higher similarity with proteins from Vertebrates compared to proteins from Arthropods or Nematodes. Their composition shows enrichment in positively charged amino acids that might contribute to their thermostability. The gene content of Alvinella reveals that an important pool of genes previously considered to be specific to Deuterostomes were in fact already present in the last common ancestor of the Bilaterian animals, but have been secondarily lost in model invertebrates. This pool is enriched in glycoproteins that play a key role in intercellular communication, hormonal regulation and immunity.

Conclusions

Our study starts to unravel the gene content and sequence evolution of a deep-sea annelid, revealing key features in eukaryote adaptation to extreme environmental conditions and highlighting the proximity of Annelids and Vertebrates.

Determining gene flow and the influence of selection across the equatorial barrier of the East Pacific Rise in the tube-dwelling polychaete Alvinella pompejana

Background

Comparative phylogeography recently performed on the mitochondrial cytochrome oxidase I (mtCOI) gene from seven deep-sea vent species suggested that the East Pacific Rise fauna has undergone a vicariant event with the emergence of a north/south physical barrier at the Equator 1-2 Mya. Within this specialised fauna, the tube-dwelling polychaete Alvinella pompejana showed reciprocal monophyly at mtCOI on each side of the Equator (9°50'N/7°25'S), suggesting potential, ongoing allopatric speciation. However, the development of a barrier to gene flow is a long and complex process. Secondary contact between previously isolated populations can occur when physical isolation has not persisted long enough to result in reproductive isolation between genetically divergent lineages, potentially leading to hybridisation and subsequent allelic introgression. The present study evaluates the strength of the equatorial barrier to gene flow and tests for potential secondary contact zones between A. pompejana populations by comparing the mtCOI gene with nuclear genes.

Results

Allozyme frequencies and the analysis of nucleotide polymorphisms at three nuclear loci confirmed the north/south genetic differentiation of Alvinella pompejana populations along the East Pacific Rise. Migration was oriented north-to-south with a moderate allelic introgression between the two geographic groups over a narrow geographic range just south of the barrier. Multilocus analysis also indicated that southern populations have undergone demographic expansion as previously suggested by a multispecies approach. A strong shift in allozyme frequencies together with a high level of divergence between alleles and a low number of 'hybrid' individuals were observed between the northern and southern groups using the phosphoglucomutase gene. In contrast, the S-adenosylhomocysteine hydrolase gene exhibited reduced diversity and a lack of population differentiation possibly due to a selective sweep or hitch-hiking.

Conclusions

The equatorial barrier leading to the separation of East Pacific Rise vent fauna into two distinct geographic groups is still permeable to migration, with a probable north-to-south migration route for A. pompejana. This separation also coincides with demographic expansion in the southern East Pacific Rise. Our results suggest that allopatry resulting from ridge offsetting is a common mechanism of speciation for deep-sea hydrothermal vent organisms.

Origin and evolution of the unique tetra-domain hemoglobin from the hydrothermal vent scale worm Branchipolynoe

Hemoglobin is the most common respiratory pigment in annelids. It can be intra or extracellular, and this latter type can form large multimeric complexes. The hydrothermal vent scale worms Branchipolynoe symmytilida and Branchipolynoe seepensis express an extracellular tetra-domain hemoglobin (Hb) that is unique in annelids. We sequenced the gene for the single-domain and tetra-domain globins in these two species. The single-domain gene codes for a mature protein of 137 amino acids, and the tetra-domain gene codes for a mature protein of 552 amino acids. The single-domain gene has a typical three exon/two intron structure, with introns located at their typical positions (B12.2 and G7.0). This structure is repeated four times in the tetra-domain gene, with no bridge introns or linker sequences between domains. The phylogenetic position of Branchipolynoe globins among known annelid globins revealed that, although extracellular, they cluster within the annelid intracellular globins clade, suggesting that the extracellular state of these Hbs is the result of convergent evolution. The tetra-domain structure likely resulted from two tandem duplications, domain 1 giving rise to domain 2 and after this the two-domain gene duplicated to produce domains 3 and 4. The high O(2) affinity of Branchipolynoe extracellular globins may be explained by the two key residues (B10Y and E7Q) in the heme pocket in each of the domains of the single and tetra-domain globins, which have been shown to be essential in the oxygen-avid Hb from the nematode Ascaris suum. This peculiar globin evolutionary path seems to be very different from other annelid extracellular globins and is most likely the product of evolutionary tinkering associated with the strong selective pressure to adapt to chronic hypoxia that characterizes hydrothermal vents.

Speciation in the deep sea: multi-locus analysis of divergence and gene flow between two hybridizing species of hydrothermal vent mussels

BACKGROUND

Reconstructing the history of divergence and gene flow between closely-related organisms has long been a difficult task of evolutionary genetics. Recently, new approaches based on the coalescence theory have been developed to test the existence of gene flow during the process of divergence. The deep sea is a motivating place to apply these new approaches. Differentiation by adaptation can be driven by the heterogeneity of the hydrothermal environment while populations should not have been strongly perturbed by climatic oscillations, the main cause of geographic isolation at the surface.

METHODOLOGY/PRINCIPAL FINDING

Samples of DNA sequences were obtained for seven nuclear loci and a mitochondrial locus in order to conduct a multi-locus analysis of divergence and gene flow between two closely related and hybridizing species of hydrothermal vent mussels, Bathymodiolus azoricus and B. puteoserpentis. The analysis revealed that (i) the two species have started to diverge approximately 0.760 million years ago, (ii) the B. azoricus population size was 2 to 5 time greater than the B. puteoserpentis and the ancestral population and (iii) gene flow between the two species occurred over the complete species range and was mainly asymmetric, at least for the chromosomal regions studied.

CONCLUSIONS/SIGNIFICANCE

A long history of gene flow has been detected between the two Bathymodiolus species. However, it proved very difficult to conclusively distinguish secondary introgression from ongoing parapatric differentiation. As powerful as coalescence approaches could be, we are left by the fact that natural populations often deviates from standard assumptions of the underlying model. A more direct observation of the history of recombination at one of the seven loci studied suggests an initial period of allopatric differentiation during which recombination was blocked between lineages. Even in the deep sea, geographic isolation may well be a crucial promoter of speciation.

Global depression in gene expression as a response to rapid thermal changes in vent mussels

Hydrothermal vent mussels belonging to the genus Bathymodiolus are distributed worldwide and dominate communities at shallow Atlantic hydrothermal sites. While organisms inhabiting coastal ecosystems are subjected to predictable oscillations of physical and chemical variables owing to tidal cycles, the vent mussels sustain pronounced temperature changes over short periods of time, correlated to the alternation of oxic/anoxic phases. In this context, we focused on the short-term adaptive response of mussels to temperature change at a molecular level. The mRNA expression of 23 genes involved in various cell functions of the vent mussel Bathymodiolus azoricus was followed after heat shocks for either 30 or 120 min, at 25 and 30 degrees C over a 48 h recovery period at 5 degrees C. Mussels were genotyped at 10 enzyme loci to explore a relationship between natural genetic variation, gene expression and temperature adaptation. Results indicate that the mussel response to increasing temperature is a depression in gene expression, such a response being genotypically correlated at least for the Pgm-1 locus. This suggests that an increase in temperature could be a signal triggering anaerobiosis for B. azoricus or this latter alternatively behaves more like a 'cold' stenotherm species, an attribute more related to its phylogenetic history, a cold seeps/wood fall origin.

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.

Evidence for a slightly deleterious effect of intron polymorphisms at the EF1alpha gene in the deep-sea hydrothermal vent bivalve Bathymodiolus

A multilocus analysis was initiated in order to infer the general effect of demography and the indirect effect of positive selection on some chromosome segments in Bathymodiolus. Mussels of the genus Bathymodiolus inhabit the very hostile, fragmented and variable environment of deep-sea hydrothermal vents which is thought to cause recurrent population bottlenecks via extinction/colonisation processes and adaptation to new environmental conditions. In the course of this work we discovered that the assumption of neutrality of non-coding polymorphisms usually made in genome scan experiments was likely to be violated at one of the loci we analysed. The direct effect of slight purifying selection on non-coding polymorphisms shares many resemblances with the indirect effect of positive selection through genetic hitchhiking. Combining polymorphism with divergence data for several closely related species allowed us to obtain different expectations for the direct effect of negative selection and the indirect effect of positive selection. We observed a strong excess of rare non-coding polymorphisms at the second intron of the EF1alpha gene in the two species Bathymodiolus azoricus and Bathymodiolus thermophilus, while two other loci, the mitochondrial COI gene and an intron of the Lysozyme gene, did not exhibit such a deviation. In addition, the divergence rate of the EF1alpha intron was estimated to be unexpectedly low when calibrated using the closure of the Panama Isthmus that interrupted gene flow between the two species. The polymorphism to divergence ratio was similar to the one observed for the other two loci, in accordance to the hypothesis of purifying selection. We conclude that slight purifying selection is likely to act on polymorphic intronic mutations of the EF1alpha second intron and discuss the possible relationship with the specific biology of Bathymodiolus mussels.

Incomplete cryptic speciation between intertidal and subtidal morphs of Acrocnida brachiata (Echinodermata: Ophiuroidea) in the Northeast Atlantic

The brittle-star Acrocnida brachiata (Montagu) lives in sandy-bottom habitat of both intertidal and subtidal zones along the coasts of the northwestern Europe. An allozyme frequency-based survey (five enzyme loci) was combined with a mitochondrial (mt) COI haplotype analysis (598-bp sequences) on 17 populations to trace back past colonization pathways from the actual population structure of the species. Both genetic markers display a sharp genetic break between intertidal (clade I) and subtidal populations (clade S). This break corresponds to an allele frequency inversion at three enzyme loci (Hk, Pgm and Pgi) and a deep divergence of about 20% in mtCOI sequences between most of the intertidal populations and other samples. The geographic distribution of clade I seems to be more restricted than clade S as it is absent from the intertidal of the eastern English Channel and North Sea and may be replaced by clade S in south Brittany. Applying previously published rates of mutation, divergence between the two clades is estimated to pre-date 5 million years ago and may be due to allopatric speciation processes at the Mio-Pliocene transition. The occurrence of putative hybrids in a few localities, however, suggests incomplete cryptic speciation with secondary contact zones. The relative importance of colonization history vs. habitat specialization are discussed in the light of neutral evolution as tested from mtCOI gene sequences. While differential selection seems to have contributed little to the separation of the lineages, it may have played a role in the emergence of adaptive polymorphisms in the hybrid zone. Furthermore, congruent spatial patterns of differentiation were observed in both clades suggesting a recent increase in population size. These findings are in agreement with a recent expansion of the populations during or after the formation of the English Channel, from a southern refuge for the subtidal clade whereas the intertidal clade may have persisted further north. As previously suspected for a species with a very short pelagic larval phase, contemporary gene flow between distant or adjacent populations appears to be extremely reduced or even absent.

Comparative phylogeography of two coastal polychaete tubeworms in the Northeast Atlantic supports shared history and vicariant events

The historic processes which have led to the present-day patterns of genetic structure in the marine coastal fauna of the Northeast Atlantic are still poorly understood. While tectonic uplifts and changes in sea level may have caused large-scale vicariance, warmer conditions during glacial maxima may have allowed pockets of diversity to persist to a much wider extent than in the Northwestern Atlantic. The large-scale geographic distribution of deeply divergent lineages of the coastal polychaete tubeworms Pectinaria koreni (two clades) and Owenia fusiformis (three clades) were compared using a fragment of the mitochondrial cytochrome oxidase I gene (mtCOI). All lineages were present along the biogeographic transition zone on the north coast of Brittany (France) and we found evidence pointing towards congruence in the timing of cladogenic events between Pectinaria sp. (P. auricoma/P. belgica and P. koreni) and Owenia sp., suggesting a shared history of vicariant events. More conserved 16SrRNA sequences obtained from four species of Pectinariidae together with mtCOI sequences of P. koreni seem consistent with an initial establishment of pectinariids in the north, and a southward colonization of the Northeast Atlantic. Phylogeographic patterns in O. fusiformis were also consistent with a north/south pattern of lineage splitting and congruent levels of divergence were detected between lineages of both species. We observed signatures of both persistence in small northern glacial refugia, and of northwards range expansion from regions situated closer to the Mediterranean. However, whether the recolonization of the Northeast Atlantic by both species actually reflects separate interglacial periods is unclear with regards to the lack of molecular clock calibration in coastal polychaete species.

Temperature resistance studies on the deep-sea vent shrimpMirocaris fortunata

The shrimp Mirocaris fortunata is a hydrothermal vent species that is found at most vent-sites along the Mid-Atlantic Ridge. This endemic species is found across a hydrothermal gradient, with thermal conditions ranging from 2–9°C in ambient seawater to fairly warm values of about 25°C. We performed in vivo experiments on M. fortunata specimens originating from different sites and depths (850 m to 2300 m), both at atmospheric pressure and in pressurized aquaria, to characterise the upper thermal limits of this species. Atmospheric pressure results show that thermal physiology should be studied at each population's native pressure. At in situ pressure, shrimps from Menez Gwen (850 m depth) and Lucky Strike(1700 m depth) do not survive temperatures of 39°C, and the `loss of equilibrium' response suggests that their critical thermal maximum(Ctmax), is about 36±1°C for both sites. This value is similar to those found for another vent shrimp, Rimicaris exoculata, which is thought to be a more temperature-resistant organism,so temperature resistance does not appear to be a crucial factor for explaining differences in distribution of shrimp species in a given vent site. Finally, the data for both vent shrimps are also comparable to those of other non-vent tropical caridean species.

Sharp genetic break between Atlantic and English Channel populations of the polychaete Pectinaria koreni, along the North coast of France

This study uses enzymatic and mitochondrial genes to infer the relative importance of historical processes and contemporary hydrodynamic features on the observed patterns of genetic structure in subdivided populations of Pectinaria koreni (Polychaeta: Pectinariidae) along the coasts of Brittany and the English Channel. Nucleotide sequence variation of a 603-bp fragment of the mtDNA cytochrome oxidase subunit I gene revealed a surprisingly deep phylogeographic break of about 16% divergence separating the Brittany and Channel populations, which coincides with a biogeographic boundary along the western coast of Brittany. Deep sequence divergence with fixed haplotype differences and the inversion of allele frequencies at two enzyme loci suggests the occurrence of potential cryptic or sibling species of P. koreni. The two clades showed opposite features. Channel populations exhibited bimodal match-mismatch curves due to two highly divergent haplotypes occurring at high frequencies and no overall heterozygote deficiencies at enzyme loci, suggesting respectively, a historic secondary contact between two differentiated populations followed by contemporary panmixia. On the contrary, Brittany populations displayed unimodal curves with low nucleotide diversity and highly significant heterozygote deficiencies, probably reminiscent of a recent population expansion and recolonisation of Brittany with contemporary admixture of divergent populations.

Thermal selection of PGM allozymes in newly founded populations of the thermotolerant vent polychaete Alvinella pompejana

Alvinella pompejana lives on the top of chimneys at deep-sea hydrothermal vents of the East Pacific Rise. It is thought to be one of the most thermotolerant and eurythermal metazoans. Our experimental approach combines methods of population genetics and biochemistry, considering temperature as a potential selective factor. Phosphoglucomutase (Pgm-1 locus) is one of the most polymorphic loci of A. pompejana and exhibits four alleles, from which alleles 90 and 100 dominate with frequencies of approximately 0.5 in populations. Results from previous studies suggested that allele 90 might be more thermostable than allele 100. Significant genetic differentiation was found by comparing contrasted microhabitats, especially the young, still hot, versus older and colder chimneys, with allele 90 being at highest frequency on young chimneys. Moreover the frequency of allele 90 was positively correlated with mean temperature at the opening of Alvinella tubes. In parallel, thermostability and thermal optimum experiments demonstrated that allele 90 is more thermostable and more active at higher temperatures than allele 100. This dataset supports an additive model of diversifying selection in which allele 90 is favoured on young hot chimneys but counterbalanced over the whole metapopulation by the dynamics of the vent ecosystem.

Escarpia southwardae sp. nov., a new species of vestimentiferan tubeworm (Annelida, Siboglinidae) from West African cold seeps

A new species of vestimentiferan tubeworm belonging to the genus Escarpia is described from cold seeps off the western coast of Africa. The description is based on two collections (one of 180 animals, the other of 30 animals) using both morphological and molecular techniques. Morphologically, the African tubeworms are very similar to Escarpia laminata Jones, 1985 but differ from all other escarpids by the lack of branchial pinnules, a unique feature among vestimentiferans. Molecular evidence from sequences of the cytochrome-c oxidase subunit I gene places the species in the escarpid clade, closely related to E. laminata and Escarpia spicata Jones, 1985, but fails to discriminate among the three species. Four morphotypes are identified in the African species, corresponding to the four permutations of the following characters: presence or absence of an axial rod on the obturaculum and presence or absence of a split on the posterior ventral margin of the vestimentum. However, molecular data could not distinguish them as separate species. We suggest that the lack of an axial rod reflects predation. Biometrical data indicate a discontinuous recruitment period, as is known for other vestimentiferan species. Sex ratios are balanced, but females tend to be larger than males. We hypothesize that the males grow more slowly or die younger than the females.

Heat-shock response and temperature resistance in the deep-sea vent shrimp Rimicaris exoculata

The shrimp Rimicaris exoculata swarms around hydrothermal black smoker chimneys at most vent sites along the Mid-Atlantic Ridge. This species maintains close proximity to the hydrothermal fluid, where temperatures can reach 350 degrees C and steep thermal and chemical gradients are expected. We performed in vivo experiments in pressurized aquaria to determine the upper thermal limit [critical thermal maximum (CT(max))] of R. exoculata and to investigate some characteristics of the shrimp stress response to heat exposure. These experiments showed that the shrimp does not tolerate sustained exposure to temperatures in the 33-37 degrees C range (CT(max)). A heat-inducible stress protein belonging to the hsp70 family was identified in R. exoculata, and its synthesis threshold induction temperature is below 25 degrees C. The R. exoculata optimal thermal habitat may thus be restricted to values lower than previously expected (<25 degrees C).

The loss of the hemoglobin H2S-binding function in annelids from sulfide-free habitats reveals molecular adaptation driven by Darwinian positive selection

The hemoglobin of the deep-sea hydrothermal vent vestimentiferan Riftia pachyptila (annelid) is able to bind toxic hydrogen sulfide (H(2)S) to free cysteine residues and to transport it to fuel endosymbiotic sulfide-oxidising bacteria. The cysteine residues are conserved key amino acids in annelid globins living in sulfide-rich environments, but are absent in annelid globins from sulfide-free environments. Synonymous and nonsynonymous substitution analysis from two different sets of orthologous annelid globin genes from sulfide rich and sulfide free environments have been performed to understand how the sulfide-binding function of hemoglobin appeared and has been maintained during the course of evolution. This study reveals that the sites occupied by free-cysteine residues in annelids living in sulfide-rich environments and occupied by other amino acids in annelids from sulfide-free environments, have undergone positive selection in annelids from sulfide-free environments. We assumed that the high reactivity of cysteine residues became a disadvantage when H(2)S disappeared because free cysteines without their natural ligand had the capacity to interact with other blood components, disturb homeostasis, reduce fitness and thus could have been counterselected. To our knowledge, we pointed out for the first time a case of function loss driven by molecular adaptation rather than genetic drift. If constraint relaxation (H(2)S disappearance) led to the loss of the sulfide-binding function in modern annelids from sulfide-free environments, our work suggests that adaptation to sulfide-rich environments is a plesiomorphic feature, and thus that the annelid ancestor could have emerged in a sulfide-rich environment.

Introgression patterns in the mosaic hybrid zone between Mytilus edulis and M. galloprovincialis

Hybrid zones are fascinating systems to investigate the structure of genetic barriers. Marine hybrid zones deserve more investigation because of the generally high dispersion potential of planktonic larvae which allows migration on scales unrivalled by terrestrial species. Here we analyse the genetic structure of the mosaic hybrid zone between the marine mussels Mytilus edulis and M. galloprovincialis, using three length-polymorphic PCR loci as neutral and diagnostic markers on 32 samples along the Atlantic coast of Europe. Instead of a single genetic gradient from M. galloprovincialis on the Iberian Peninsula to M. edulis populations in the North Sea, three successive transitions were observed in France. From South to North, the frequency of alleles typical of M. galloprovincialis first decreases in the southern Bay of Biscay, remains low in Charente, then increases in South Brittany, remains high in most of Brittany, and finally decreases again in South Normandy. The two enclosed patches observed in the midst of the mosaic hybrid zone in Charente and Brittany, although predominantly M. edulis-like and M. galloprovincialis-like, respectively, are genetically original in two respects. First, considering only the various alleles typical of one species, the patches show differentiated frequencies compared to the reference external populations. Second, each patch is partly introgressed by alleles of the other species. When introgression is taken into account, linkage disequilibria appear close to their maximum possible values, indicating a strong genetic barrier within all transition zones. Some pre- or postzygotic isolation mechanisms (habitat specialization, spawning asynchrony, assortative fertilization and hybrid depression) have been documented in previous studies, although their relative importance remains to be evaluated. We also provided evidence for a recent migratory 'short-cut' connecting M. edulis-like populations of the Charente patch to an external M. edulis population in Normandy and thought to reflect artificial transfer of spat for aquaculture.

Evolution of the sulfide-binding function within the globin multigenic family of the deep-sea hydrothermal vent tubeworm Riftia pachyptila

The giant extracellular hexagonal bilayer hemoglobin (HBL-Hb) of the deep-sea hydrothermal vent tube worm Riftia pachyptila is able to transport simultaneously O(2) and H(2)S in the blood from the gills to a specific organ: the trophosome that harbors sulfide-oxidizing endosymbionts. This vascular HBL-Hb is made of 144 globins from which four globin types (A1, A2, B1, and B2) coevolve. The H(2)S is bound at a specific location (not on the heme site) onto two of these globin types. In order to understand how such a function emerged and evolved in vestimentiferans and other related annelids, six partial cDNAs corresponding to the six globins known to compose the multigenic family of R. pachyptila have been identified and sequenced. These partial sequences (ca. 120 amino acids, i.e., 80% of the entire protein) were used to reconstruct molecular phylogenies in order to trace duplication events that have led to the family organization of these globins and to locate the position of the free cysteine residues known to bind H(2)S. From these sequences, only two free cysteine residues have been found to occur, at positions Cys + 1 (i.e., 1 a.a. from the well-conserved distal histidine) and Cys + 11 (i.e., 11 a.a. from the same histidine) in globins B2 and A2, respectively. These two positions are well conserved in annelids, vestimentiferans, and pogonophorans, which live in sulfidic environments. The structural comparison of the hydrophobic environment that surrounds these cysteine residues (the sulfide-binding domain) using hydrophobic cluster analysis plots, together with the cysteine positions in paralogous strains, suggests that the sulfide-binding function might have emerged before the annelid radiation in order to detoxify this toxic compound. Moreover, globin evolutionary rates are highly different between paralogous strains. This suggests that either the two globin subfamilies involved in the sulfide-binding function (A2 and B2) have evolved under strong directional selective constraints (negative selection) and that the two other globins (A1 and B1) have accumulated more substitutions through positive selection or have evolved neutrally after a relaxation of selection pressures. A likely scenario on the evolution of this multigenic family is proposed and discussed from this data set.