Population subdivision of hydrothermal vent polychaete Alvinella pompejana across equatorial and Easter Microplate boundaries

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.

Balanced Polymorphism at the Pgm-1 Locus of the Pompeii Worm Alvinella pompejana and Its Variant Adaptability Is Only Governed by Two QE Mutations at Linked Sites

The polychaete Alvinella pompejana lives exclusively on the walls of deep-sea hydrothermal chimneys along the East Pacific Rise (EPR), and displays specific adaptations to withstand the high temperatures and hypoxia associated with this highly variable habitat. Previous studies have revealed the existence of a balanced polymorphism on the enzyme phosphoglucomutase associated with thermal variations, where allozymes 90 and 100 exhibit different optimal activities and thermostabilities. Exploration of the mutational landscape of phosphoglucomutase 1 revealed the maintenance of four highly divergent allelic lineages encoding the three most frequent electromorphs over the geographic range of A. pompejana. This polymorphism is only governed by two linked amino acid replacements, located in exon 3 (E155Q and E190Q). A two-niche model of selection, including ‘cold’ and ‘hot’ conditions, represents the most likely scenario for the long-term persistence of these isoforms. Using directed mutagenesis and the expression of the three recombinant variants allowed us to test the additive effect of these two mutations on the biochemical properties of this enzyme. Our results are coherent with those previously obtained from native proteins, and reveal a thermodynamic trade-off between protein thermostability and catalysis, which is likely to have maintained these functional phenotypes prior to the geographic separation of populations across the Equator about 1.2 million years ago.

Phylogeny, biogeography and systematics of Pacific vent, methane seep, and whale-fall Parougia (Dorvilleidae : Annelida), with eight new species

Dorvilleidae is a diverse group of annelids found in many marine environments and also commonly associated with chemosynthetic habitats. One dorvilleid genus, Parougia, currently has 11 described species, of which two are found at vents or seeps: Parougia wolfi and Parougia oregonensis. Eight new Parougia species are recognised and described in this study from collections in the Pacific Ocean, all from whale-falls, hydrothermal vents, or methane seeps at ~600-m depth or greater. The specimens were studied using morphology and phylogenetic analyses of DNA sequences from mitochondrial (cytochrome c oxidase subunit I, 16S rRNA, and cytochrome b) and nuclear (18S rRNA and histone 3) genes. Six sympatric Parougia spp. were found at Hydrate Ridge, Oregon, while three of the Parougia species occurred at different types of chemosynthetic habitats. Two new species were found over wide geographical and bathymetric ranges. Another dorvilleid genus, Ophryotrocha, has previously been highlighted as diversifying in the deep-sea environment. Our results document the hitherto unknown diversity of another dorvilleid genus, Parougia, at various chemosynthetic environments. http://zoobank.org/urn:lsid:zoobank.org:pub:EC7EBBEA-2FB5-43D6-BE53-1A468B541A5C

Amplicon sequencing of 42 nuclear loci supports directional gene flow between South Pacific populations of a hydrothermal vent limpet

In the past few decades, population genetics and phylogeographic studies have improved our knowledge of connectivity and population demography in marine environments. Studies of deep-sea hydrothermal vent populations have identified barriers to gene flow, hybrid zones, and demographic events, such as historical population expansions and contractions. These deep-sea studies, however, used few loci, which limit the amount of information they provided for coalescent analysis and thus our ability to confidently test complex population dynamics scenarios. In this study, we investigated population structure, demographic history, and gene flow directionality among four Western Pacific hydrothermal vent populations of the vent limpet Lepetodrilus aff. schrolli. These vent sites are located in the Manus and Lau back-arc basins, currently of great interest for deep-sea mineral extraction. A total of 42 loci were sequenced from each individual using high-throughput amplicon sequencing. Amplicon sequences were analyzed using both genetic variant clustering methods and evolutionary coalescent approaches. Like most previously investigated vent species in the South Pacific, L. aff. schrolli showed no genetic structure within basins but significant differentiation between basins. We inferred significant directional gene flow from Manus Basin to Lau Basin, with low to no gene flow in the opposite direction. This study is one of the very few marine population studies using >10 loci for coalescent analysis and serves as a guide for future marine population studies.

Population genetic differentiation of the hydrothermal vent crab Austinograea alayseae (Crustacea: Bythograeidae) in the Southwest Pacific Ocean

To understand the origin, migration, and distribution of organisms across disjunct deep-sea vent habitats, previous studies have documented the population genetic structures of widely distributed fauna, such as gastropods, bivalves, barnacles, and squat lobsters. However, a limited number of investigations has been conducted in the Southwest Pacific Ocean, and many questions remain. In this study, we determined the population structure of the bythograeid crab Austinograea alayseae from three adjacent vent systems (Manus Basin, North Fiji Basin, and Tonga Arc) in the Southwest Pacific Ocean using the sequences of two mitochondrial genes (COI and 16S rDNA) and one nuclear gene (28S rDNA). Populations were divided into a Manus clade and a North Fiji-Tonga clade, with sequence divergence values in the middle of the barcoding gap for bythograeids. We inferred that hydrographic and/or physical barriers act on the gene flow of A. alayseae between the Manus and North Fiji basins. Austinograea alayseae individuals interact freely between the North Fiji Basin and the Lau Basin (Tonga Arc). Although further studies of genetic differentiation over a geological time scale, life-history attributes, and genome-based population genetics are needed to improve our understanding of the evolutionary history of A. alayseae, our results contribute to elucidating the phylogeny, evolution, and biogeography of bythograeids.

Deepest and hottest hydrothermal activity in the Okinawa Trough: the Yokosuka site at Yaeyama Knoll

Since the initial discovery of hydrothermal vents in 1977, these 'extreme' chemosynthetic systems have been a focus of interdisciplinary research. The Okinawa Trough (OT), located in the semi-enclosed East China Sea between the Eurasian continent and the Ryukyu arc, hosts more than 20 known vent sites but all within a relatively narrow depth range (600-1880 m). Depth is a significant factor in determining fluid temperature and chemistry, as well as biological composition. However, due to the narrow depth range of known sites, the actual influence of depth here has been poorly resolved. Here, the Yokosuka site (2190 m), the first OT vent exceeding 2000 m depth is reported. A highly active hydrothermal vent site centred around four active vent chimneys reaching 364°C in temperature, it is the hottest in the OT. Notable Cl depletion (130 mM) and both high H₂ and CH₄ concentrations (approx. 10 mM) probably result from subcritical phase separation and thermal decomposition of sedimentary organic matter. Microbiota and fauna were generally similar to other sites in the OT, although with some different characteristics. In terms of microbiota, the H₂-rich vent fluids in Neuschwanstein chimney resulted in the dominance of hydrogenotrophic chemolithoautotrophs such as Thioreductor and Desulfobacterium. For fauna, the dominance of the deep-sea mussel Bathymodiolus aduloides is surprising given other nearby vent sites are usually dominated by B. platifrons and/or B. japonicus, and a sponge field in the periphery dominated by Poecilosclerida is unusual for OT vents. Our insights from the Yokosuka site implies that although the distribution of animal species may be linked to depth, the constraint is perhaps not water pressure and resulting chemical properties of the vent fluid but instead physical properties of the surrounding seawater. The potential significance of these preliminary results and prospect for future research on this unique site are discussed.

Geographical structure of endosymbiotic bacteria hosted by Bathymodiolus mussels at eastern Pacific hydrothermal vents

BACKGROUND

Chemolithoautotrophic primary production sustains dense invertebrate communities at deep-sea hydrothermal vents and hydrocarbon seeps. Symbiotic bacteria that oxidize dissolved sulfur, methane, and hydrogen gases nourish bathymodiolin mussels that thrive in these environments worldwide. The mussel symbionts are newly acquired in each generation via infection by free-living forms. This study examined geographical subdivision of the thiotrophic endosymbionts hosted by Bathymodiolus mussels living along the eastern Pacific hydrothermal vents. High-throughput sequencing data of 16S ribosomal RNA encoding gene and fragments of six protein-coding genes of symbionts were examined in the samples collected from nine vent localities at the East Pacific Rise, Galápagos Rift, and Pacific-Antarctic Ridge.

RESULTS

Both of the parapatric sister-species, B. thermophilus and B. antarcticus, hosted the same numerically dominant phylotype of thiotrophic Gammaproteobacteria. However, sequences from six protein-coding genes revealed highly divergent symbiont lineages living north and south of the Easter Microplate and hosted by these two Bathymodiolus mussel species. High heterogeneity of symbiont haplotypes among host individuals sampled from the same location suggested that stochasticity associated with initial infections was amplified as symbionts proliferated within the host individuals. The mussel species presently contact one another and hybridize along the Easter Microplate, but the northern and southern symbionts appear to be completely isolated. Vicariance associated with orogeny of the Easter Microplate region, 2.5-5.3 million years ago, may have initiated isolation of the symbiont and host populations. Estimates of synonymous substitution rates for the protein-coding bacterial genes examined in this study were 0.77-1.62%/nucleotide/million years.

CONCLUSIONS

Our present study reports the most comprehensive population genetic analyses of the chemosynthetic endosymbiotic bacteria based on high-throughput genetic data and extensive geographical sampling to date, and demonstrates the role of the geographical features, the Easter Microplate and geographical distance, in the intraspecific divergence of this bacterial species along the mid-ocean ridge axes in the eastern Pacific. Altogether, our results provide insights into extrinsic and intrinsic factors affecting the dispersal and evolution of chemosynthetic symbiotic partners in the hydrothermal vents along the eastern Pacific Ocean.

Population subdivision of hydrothermal vent polychaete Alvinella pompejana across equatorial and Easter Microplate boundaries

Background

The Equator and Easter Microplate regions of the eastern Pacific Ocean exhibit geomorphological and hydrological features that create barriers to dispersal for a number of animals associated with deep-sea hydrothermal vent habitats. This study examined effects of these boundaries on geographical subdivision of the vent polychaete Alvinella pompejana. DNA sequences from one mitochondrial and eleven nuclear genes were examined in samples collected from ten vent localities that comprise the species’ known range from 23°N latitude on the East Pacific Rise to 38°S latitude on the Pacific Antarctic Ridge.

Results

Multi-locus genotypes inferred from these sequences clustered the individual worms into three metapopulation segments — the northern East Pacific Rise (NEPR), southern East Pacific Rise (SEPR), and northeastern Pacific Antarctic Ridge (PAR) — separated by the Equator and Easter Microplate boundaries. Genetic diversity estimators were negatively correlated with tectonic spreading rates. Application of the isolation-with-migration (IMa2) model provided information about divergence times and demographic parameters. The PAR and NEPR metapopulation segments were estimated to have split roughly 4.20 million years ago (Mya) (2.42–33.42 Mya, 95 % highest posterior density, (HPD)), followed by splitting of the SEPR and NEPR segments about 0.79 Mya (0.07–6.67 Mya, 95 % HPD). Estimates of gene flow between the neighboring regions were mostly low (2 Nm < 1). Estimates of effective population size decreased with southern latitudes: NEPR > SEPR > PAR.

Conclusions

Highly effective dispersal capabilities allow A. pompejana to overcome the temporal instability and intermittent distribution of active hydrothermal vents in the eastern Pacific Ocean. Consequently, the species exhibits very high levels of genetic diversity compared with many co-distributed vent annelids and mollusks. Nonetheless, its levels of genetic diversity in partially isolated populations are inversely correlated with tectonic spreading rates. As for many other vent taxa, this pioneering colonizer is similarly affected by local rates of habitat turnover and by major dispersal filters associated with the Equator and the Easter Microplate region.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-016-0807-9) contains supplementary material, which is available to authorized users.