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

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

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

Shining light on a deep-sea bacterial symbiont population structure with CRISPR

Many foundation species in chemosynthesis-based ecosystems rely on environmentally acquired symbiotic bacteria for their survival. Hence, understanding the biogeographic distributions of these symbionts at regional scales is key to understanding patterns of connectivity and predicting resilience of their host populations (and thus whole communities). However, such assessments are challenging because they necessitate measuring bacterial genetic diversity at fine resolutions. For this purpose, the recently discovered clustered regularly interspaced short palindromic repeats (CRISPR) constitutes a promising new genetic marker. These DNA sequences harboured by about half of bacteria hold their viral immune memory, and as such, might allow discrimination of different lineages or strains of otherwise indistinguishable bacteria. In this study, we assessed the potential of CRISPR as a hypervariable phylogenetic marker in the context of a population genetic study of an uncultured bacterial species. We used high-throughput CRISPR-based typing along with multi-locus sequence analysis (MLSA) to characterize the regional population structure of the obligate but environmentally acquired symbiont species Candidatus Endoriftia persephone on the Juan de Fuca Ridge. Mixed symbiont populations of Ca. Endoriftia persephone were sampled across individual Ridgeia piscesae hosts from contrasting habitats in order to determine if environmental conditions rather than barriers to connectivity are more important drivers of symbiont diversity. We showed that CRISPR revealed a much higher symbiont genetic diversity than the other housekeeping genes. Several lines of evidence imply this diversity is indicative of environmental strains. Finally, we found with both CRISPR and gene markers that local symbiont populations are strongly differentiated across sites known to be isolated by deep-sea circulation patterns. This research showed the high power of CRISPR to resolve the genetic structure of uncultured bacterial populations and represents a step towards making keystone microbial species an integral part of conservation policies for upcoming mining operations on the seafloor.

Contrasting influences on bacterial symbiont specificity by co-occurring deep-sea mussels and tubeworms

Relationships fueled by sulfide between deep-sea invertebrates and bacterial symbionts are well known, yet the diverse overlapping factors influencing symbiont specificity are complex. For animals that obtain their symbionts from the environment, both host identity and geographic location can impact the ultimate symbiont partner. Bacterial symbionts were analysed for three co-occurring species each of Bathymodiolus mussels and vestimentiferan tubeworms, from three deep methane seeps off the west coast of Costa Rica. The bacterial internal transcribed spacer gene was analysed via direct and barcoded amplicon sequencing to reveal fine-scale symbiont diversity. Each of the three mussel species (B. earlougheri, B. billschneideri and B. nancyschneideri) hosted genetically distinct thiotrophic endosymbionts, despite living nearly side-by-side in their habitat, suggesting that host identity is crucial in driving symbiont specificity. The dominant thiotrophic symbiont of co-occurring tubeworms Escarpia spicata and Lamellibrachia (L. barhami and L. donwalshi), on the other hand, was identical regardless of host species or sample location, suggesting lack of influence by either factor on symbiont selectivity in this group of animals. These findings highlight the specific, yet distinct, influences on the environmental acquisition of symbionts in two foundational invertebrates with similar lifestyles, and provide a rapid, precise method of examining symbiont identities.

High‐contiguity genome assembly of the chemosynthetic gammaproteobacterial endosymbiont of the cold seep tubeworm Lamellibrachia barhami

Symbiotic relationships between vestimentiferan tubeworms and chemosynthetic Gammaproteobacteria build the foundations of many hydrothermal vent and hydrocarbon seep ecosystems in the deep sea. The association between the vent tubeworm Riftia pachyptila and its endosymbiont Candidatus Endoriftia persephone has become a model system for symbiosis research in deep‐sea vestimentiferans, while markedly fewer studies have investigated symbiotic relationships in other tubeworm species, especially at cold seeps. Here we sequenced the endosymbiont genome of the tubeworm Lamellibrachia barhami from a cold seep in the Gulf of California, using short‐ and long‐read sequencing technologies in combination with Hi‐C and Dovetail Chicago libraries. Our final assembly had a size of ~4.17 MB, a GC content of 54.54%, 137X coverage, 4153 coding sequences, and a checkm completeness score of 97.19%. A single scaffold contained 99.51% of the genome. Comparative genomic analyses indicated that the L. barhami symbiont shares a set of core genes and many metabolic pathways with other vestimentiferan symbionts, while containing 433 unique gene clusters that comprised a variety of transposases, defence‐related genes and a lineage‐specific CRISPR/Cas3 system. This assembly represents the most contiguous tubeworm symbiont genome resource to date and will be particularly valuable for future comparative genomic studies investigating structural genome evolution, physiological adaptations and host‐symbiont communication in chemosynthetic animal‐microbe symbioses.