Skip to main content

ORIGINAL RESEARCH article

Front. Microbiol., 28 June 2023
Sec. Microbiological Chemistry and Geomicrobiology
This article is part of the Research Topic Insights in Microbiological Chemistry and Geomicrobiology : 2022/2023 View all 10 articles

Variation in epibiotic bacteria on two squat lobster species of Munidopsidae

  • 1CAS Key Laboratory of Marine Ecology and Environmental Sciences, Institute of Oceanology, Chinese Academy of Sciences, Qingdao, China
  • 2Center of Deep Sea Research, Institute of Oceanology, Chinese Academy of Sciences, Qingdao, China
  • 3University of Chinese Academy of Sciences, Beijing, China
  • 4South China Sea Institute of Oceanology, Chinese Academy of Sciences, Guangzhou, China
  • 5Department of Marine Organism Taxonomy and Phylogeny, Institute of Oceanology, Chinese Academy of Sciences, Qingdao, China
  • 6College of Life Sciences, Zaozhuang University, Zaozhuang, China

The relationships between epibiotic bacteria on deep-sea hosts and host lifestyle factors are of particular interest in the field of deep-sea chemoautotrophic environmental adaptations. The squat lobsters Shinkaia crosnieri and Munidopsis verrilli are both dominant species in cold-seep ecosystems, and they have different distributions and feeding behaviors. These species may have evolved to have distinct epibiotic microbiota. Here, we compared the epibiotic bacterial communities on the M. verrilli carapace (MVcarapace), S. crosnieri carapace (SCcarapace), and S. crosnieri ventral plumose setae (SCsetae). The epibiotic bacteria on SCsetae were dense and diverse and had a multi-layer configuration, while those on MVcarapace and SCcarapace were sparse and had a monolayer configuration. Chemoautotrophic bacteria had the highest relative abundance in all epibiotic bacterial communities. The relative abundance of amplicon sequence variant 3 (ASV3; unknown species in order Thiotrichales), which is associated with sulfide oxidation, was significantly higher in SCsetae than MVcarapace and SCcarapace. Thiotrichales species seemed to be specifically enriched on SCsetae, potentially due to the synthetic substrate supply, adhesion preference, and host behaviors. We hypothesize that the S. crosnieri episymbionts use chemical fluxes near cold seeps more efficiently, thereby supporting the host’s nutrient strategies, resulting in a different distribution of the two species of squat lobster.

1. Introduction

Deep-Sea chemoautotrophic ecosystems, including cold seeps, are often characterized by extreme environmental conditions and endemic fauna. Decapod crustaceans, such as alvinocaridid shrimps (Komai et al., 2005), lithodid crabs (Niemann et al., 2013), and squat lobsters (Chevaldonne and Olu, 1996), are common macrofauna in cold seeps. Squat lobsters, as one of the dominant microbenthic crustaceans, are key components of deep-sea chemosynthetic ecosystems, considering their wide distribution and high abundance (Martin and Haney, 2005). Squat lobsters belong to the superfamilies Galatheoidea and Chirostyloidea, with over 1,000 species (Baba and Fujita, 2017). They are found at all depths, and the Munidopsidae (Galatheoidea) species almost completely dominate at extreme depths, with typical adaptations such as eye degeneration (Warrant and Locket, 2004). The Munidopsidae species occupy various ecological niches in the deep sea, including in cold seeps (Martin and Haney, 2005). Most of those observed at cold seeps [such as Munidopsis andamanica (Hoyoux et al., 2012) and Munidopsis sp. (Macavoy et al., 2008)] have a wide distribution around the vents, with an omnivorous diet. Others [such as Kiwa puravida (Thurber et al., 2011) and Kiwa hirsute (Macpherson et al., 2005)] are concentrated around vents, with distinctive dense setae on their body surfaces. Numerous invertebrate species have formed symbiotic associations with their microbiomes, which can exert a range of effects on their hosts, predominantly beneficial ones (Mioduchowska et al., 2018). Goffredi et al. proposed that the bacteria that clustered on the plumose setae of Kiwa hirsute could be used as a nutrient resource to meet the energetic demands of the host (Goffredi et al., 2008). The provision of nutrients is a major feature of host–microbe interactions (Schwartzman et al., 2015). Dattagupta et al. found that Niphargus ictus, a freshwater cave amphipod, developed a nutritional strategy involving scraping epibiotic bacteria (by grooming itself with its gnathopods) to use as food, and it was adapted to environments with high concentrations of sulfides (Dattagupta et al., 2009). Specialized nutritional strategies might have driven the evolution of some Munidopsidae species (Hoyoux et al., 2009). Munidopsidae species’ epibiotic bacterial communities, as a source of nutrients, might be drivers of speciation in squat lobsters.

Shinkaia crosnieri (Galatheoidea: Munidopsidae) was first recorded in hydrothermal vents in the Okinawa Trough, mainly dwelling in chemosynthetic ecosystems in the west Pacific Ocean (Baba and Williams, 1998). Munidopsis verrilli (Galatheoidea: Munidopsidae) has been widely reported in various deep-sea ecosystems around the world, such as the San Clemente Basin (Williams et al., 2000), Tasmania Island (Baba and Poore, 2002), eastern Hokkaido (Komai and Matsuzaki, 2016), and northeastern South China Sea (Dong and Li, 2015). These two species are closely evolutionarily related and belong to different genera of the same family (Cheng et al., 2020). The plumose setae on S. crosnieri (SCsetae), which harbors dense epibiotic bacterial communities, are the main basis for distinguishing between the two squat lobsters (Wehrtmann and Acuna, 2011). According to a previous study conducted by Watsuji, the host can use the abundant episymbionts that attach to the setae as dense biofilm to obtain nutrients through in vitro digestion experiments (Watsuji et al., 2015). Some of these epibiotic bacteria were found to be thioautotrophic based on FISH and Nano-SIMS analysis, as suggested by Watsuji et al. (2012). On the other hand, M. verrilli’s underside was smooth without dense plumose setae. Leinberger et al. found evidence for a possible beneficial relationship between Munidopsis alvisca and its epibiotic bacteria (Leinberger et al., 2022). However, Munidopsis are not generally considered to have a close symbiotic association with their microbiomes but an opportunist. The two species also differ in their distribution in the chemosynthetic ecosystems. Shinkaia crosnieri forms impressive aggregations around cold seeps (refs) and hydrothermal vents in the Okinawa Trough (Tsuchida et al., 2003). In contrast, M. verrilli, associated with whale carcasses, appears to be an opportunist with a much broader geological range (Williams et al., 2000). It is evident that the concentration of reductive materials is higher in the core regions of vents and seeps, which can lead to higher primary productivity. The distance from the seepages may have influenced the microbiological communities in different squat lobsters. The microbiome’s composition and abundance may, in turn, influence the behavior and evolution of squat lobsters. Thus, we hypothesized that the establishment of symbiosis in SC may contribute to the formation of the current distribution pattern. A systematic comparison of epibiotic bacterial communities associated with different Munidopsidae species would help to explore this hypothesis.

Site F is one of the active cold seeps in the South China Sea and supports several macrobenthic communities, including squat lobsters (Xiao et al., 2020). Both S. crosnieri and M. verrilli are dominant species in Site F. They have distinct differences in their distribution in relation to the cold seep. Zhao et al. found that S. crosnieri only congregated in the immediate vicinity of the cold seep at Site F, with very high abundances, reaching up to 300 individuals per square meter (Zhao et al., 2020). Munidopsis verrilli, on the other hand, is distributed over many regions of Site F, but not at high densities (Wang X. et al., 2022). These findings reflect the differential distributions of S. crosnieri and M. verrilli within a relatively small area. The co-occurrences of the two species provide an excellent opportunity for a comparative study of the link between the epibiotic bacterial communities and the distribution of squat lobsters.

The aim of this study was to investigate the epibiotic bacterial communities on squat lobsters with different ecological niches. We used scanning electron microscopy (SEM) to visualize the morphology of the epibiotic bacteria, quantitative real-time PCR to determine the difference in the absolute density of the bacteria, and high-throughput sequencing of 16S rRNA gene to explore the differences in the composition of the epibiotic bacterial communities.

2. Materials and methods

2.1. Sample collection

Squat lobsters were collected from Site F in the South China Sea (N119°17′8.22″ E, 22°06′55.26” N) during a cruise in May, 2021. Shinkaia crosnieri was caught around the cold seep, where conspicuous methane bubbles can be observed. Munidopsis verrilli was mainly caught at peripheral mussel beds around 10 m away from the cold seep. Specimens were collected using a carousel suction sampler on the remotely operated vehicle (ROV) Faxian on the research vessel “Kexue.” They were transported in a homothermic macrofauna carrier which maintained a temperature similar to that at the cold seep, at around 4°C.

All animals were identified as S. crosnieri or M. verrilli via morphological classification, according to the descriptions given by Baba and Poore (2002) and Baba and Williams (1998). The carapaces of the two squat lobster species (MVcarapace and SCcarapace) and the plumose setae of S. crosnieri (SCsetae) were rinsed with sterilized seawater and then immediately dissected with scissors on board. Samples for molecular analysis were snap-frozen in liquid nitrogen and stored at −80°C until DNA extraction. Cold seep interface water associated with the community of S. crosnieri (designated environment [ENV]), which was used to investigate the bacterioplankton, was filtered using an in-situ microbe sampler with a 0.22-μm filter membrane and stored at −80°C.

2.2. SEM

The animal samples for SEM were preserved in glutaraldehyde, underwent replacement of glutaraldehyde after 24 h, and were then stored at 4°C. In the laboratory, the samples were dehydrated in a graded ethanol series, underwent critical-point drying, were coated with gold (sputter/carbon thread; EM ACE200, Leica, Germany), and were then observed using SEM (S-3400 N, Hitachi, Japan).

2.3. Nucleic acid extraction and 16S rRNA gene amplicon sequencing

For 16S amplicon sequencing, we obtained ten different samples. These included three carapace samples of M. verrilli, three carapace samples of S. crosnieri, three setae samples of S. crosnieri, and one environmental water sample for bacterioplankton. We used a DNeasy PowerWater Kit (Qiagen, United States) to extract total DNA from these samples, following the manufacturer’s instructions. The DNA was then diluted to 1 ng/μl for PCR amplification. The PCR amplification was performed using the primer pair 341F (5′-CCTAYGGGRBGCASCAG-3′) and 806R (5′-GGACTACNNGGGTATCTAAT-3′), which target the hypervariable regions between V3 and V4 of the prokaryotic 16S rRNA gene. All PCRs were performed with 10 μl Phusion® HighFidelity PCR Master Mix (NEB, United States), 1 μM forward primer, 1 μM reverse primer, and 1 μL template DNA. The thermal cycling program involved initial denaturation at 98°C for 1 min, 30 cycles of denaturation at 98°C for 10 s, annealing at 50°C for 30 s, and extension at 72°C for 30 s, and a final elongation at 72°C for 5 min. Electrophoresis was then performed using 2% agarose gel. Sequencing libraries were generated using an NEBNext® Ultra™ II DNA Library Prep Kit and were quantified by Qubit and qPCR. After that, NovaSeq 6000 (Illumina) was used for on-machine sequencing.

2.4. Quantification of the epibiotic bacteria

Real-time PCR was performed on QuantStudio 1 (Thermo Fisher Scientific, USA) in triplicate. 331F (5′-TCCTACGGGAGGCAGCAGT-3′) and 797R (5′-GGACTACCAGGGTATCTAATCCTGTT-3′) primers were used and each reaction (20 μL) contained 1× SYBR® Premix Ex Taq™ (TaKaRa), 1 μl template DNA, 1 μM forward primer, and 1 μM reverse primer (Nadkarni et al., 2002). The standard curve was constructed based on known amounts of purified PCR product (102 to 108 copies/μl) obtained from E. coli genomic DNA by using the bacterial 16S rRNA gene-specific primers 27F (5′-AGAGTT TGATCMTGGCTCAG-3′) and 1492R (5′-TACGGYTACCTTGTTACGACTT-3′) (Weisburg et al., 1991). The thermal cycling program involved 95°C for 15 min, followed by 35 cycles of 95°C for 15 s, 58°C for 30 s, and 72°C for 30 s. The wet weight of each sample was recorded using an electronic balance and the number of copies per unit mass (wet weight) was calculated using the standard curve.

2.5. Bioinformatic analyses

Bioinformatic analyses were performed using QIIME 2 (v2022.2) (Bolyen et al., 2019). Raw sequence data were demultiplexed and quality filtered using the q2-demux plugin followed by denoising with DADA2 (Callahan et al., 2016). All amplicon sequence variants (ASVs) were aligned using MAFFT (Katoh and Standley, 2013) and were used to construct a phylogeny using FastTree2 (Price et al., 2010). To minimize differences in sequencing depth among the samples for the alpha and beta diversity analyses, all samples were normalized by subsample at a depth of 32,549. Alpha diversity analysis, beta diversity analysis, and principal coordinate analysis (PCoA) were conducted using q2-diversity. Taxonomy was assigned to ASVs using the q2-feature-classifier (classify-sklearn naïve Bayes taxonomy classifier; Bokulich et al., 2018) and the Silva138 database (Quast et al., 2012). Analysis of variance (ANOVA) and least significant difference (LSD) post-hoc comparisons were applied in SPSS v26 to show the significance of the difference in the alpha diversity index (Shannon index; Fahy et al., 2015). STAMP v2.1.3 was used with Welch’s t-test to compare the dominant bacterial orders and ASVs in different groups (Parks et al., 2014). Analysis of similarities (ANOSIM) based on the Bray–Curtis distance was conducted using the ‘vegan’ package in R v3.5.1 (Liu and Tong, 2017). A bacterial co-occurrence network consisting of the top 200 ASVs was constructed based on a Spearman correlation matrix using the ‘ggClusterNet’ package (Wen et al., 2022).

3. Results

3.1. Observation of the microbiota on squat lobsters

According to Figures 1C,I, the ventral surface of S. crosnieri was covered with numerous plumose setae, whereas the ventral parts of M. verrilli was much smoother with only nolittle setae. The SEM images showed that the percentage of surfaces covered by biofilm was significantly higher for SCsetae (Figure 1G) than SCcarapace and MVcarapace (Figures 1F,H). In terms of the morphology, the biofilm on the carapace of both species of squat lobster tended to be dominated by groups <5 μm in length and had a monolayer configuration. In contrast, the SCsetae epibiotic bacterial community consisted of diverse morphological groups, ranging from 1 to 200 μm in length, with rod or filament shapes, and the community had a multi-layer configuration. In summary, in the setae, the biofilm was much denser and groups with slender linear shapes were more common (Figure 1).

FIGURE 1
www.frontiersin.org

Figure 1. Epibiotic microbiota. (A) Photographs of squat lobsters at Site F. The dominant squat lobsters were Shinkaia crosnieri (B,C) and Munidopsis verrilli (D,E). Scanning electron microscopy (SEM) images of the morphology of bacteria attached to the S. crosnieri carapace (SCcarapace) (F), S. crosnieri ventral setae (SCsetae) (G), and M. verrilli carapace (MVcarapace) (H). Dissecting microscope images showing the ventral view of M. verrilli, without the setae (I). White scale bar = 100 μm (F–H).

3.2. Quantification of the microbiota on squat lobsters

We established the standard curve with the R2 of 0.996 to estimate the abundance of the bacteria for each sample (Supplementary Figure S1). Real-time PCR confirmed the significant difference in the 16S rRNA gene copy number of epibiotic bacteria on the two species of squat lobster. The 16S rRNA gene copy number of SCsetae exceeded that of MVcarapace and SCcarapace by around three orders of magnitude, while the values for MVcarapace and SCcarapace were close. The results suggest that SCsetae might provide epibiotic bacteria with an excellent attachment substrate, leading to the enrichment of various groups, with a multi-layer configuration (Figure 2).

FIGURE 2
www.frontiersin.org

Figure 2. (A) Box plots of the 16S rRNA gene copy number of epibiotic bacteria from the M. verrilli carapace (MVcarapace), S. crosnieri carapace (SCcarapace), and S. crosnieri setae (SCsetae). (***: p < 0.01). (B) Alpha diversity, based on the Shannon index and Chao1 index, of bacterioplankton in cold seep interface water (environment [ENV]) and epibiotic bacteria on M. verrilli carapace (MVcarapace), S. crosnieri carapace (SCcarapace), and S. crosnieri setae (SCsetae). (*: p < 0.05).

3.3. Bacterial diversity of different groups

A total of 673,093 reads were generated from the bacterioplankton and epibiotic bacteria associated with MVcarapace, SCcarapace, and SCsetae, and 2,100 ASVs were retained. All sequences found in this study were bacteria due to the use of bacteria-specific primers. The rarefaction curves based on the sequencing data were stable, indicating that sufficient sequencing depth was achieved (Supplementary Figure S2). Based on comparisons with the Silva 138 bacterial database, 52 microbial phyla were detected in the 10 samples, comprising 143 classes, 354 orders, 591 families, and 1,123 genera.

The alpha diversity of the epibiotic bacterial and bacterioplankton communities were compared. The Shannon index varied from 5.13 to 8.68, and differed significantly among the three groups of epibiotic bacteria. The Shannon index was highest for MVcarapace and lowest for SCsetae. And samples at MVcarapace showed highest Chao1 indexes. The difference in the epibiotic bacterial communities between carapace and setae was significant (p < 0.05). However, the differences between the bacterioplankton (in the single ENV sample) and the epibiotic bacteria were insignificant.

3.4. Comparison of bacterial orders and ASVs among groups

Regarding the three epibiotic bacterial communities, Methylococcales (10.25 to 36.75%), Campylobacterales (4.32 to 34.32%), and Thiotrichales (6.28 to 51.29%) were all in the top three orders in each community, although their ranks varied. In the SCsetae epibiotic bacterial communities, Thiotrichales accounted for the highest relative abundance (42.25%). In contrast, in the MVcarapace and SCcarapace epibiotic bacterial communities, the relative abundance of Thiotrichales (14.60 and 10.74%, respectively) was lower than that of Methylococcales (24.16% and 29.54, respectively) or Campylobacterales (17.37 and 29.26%, respectively). Furthermore, the ranking of the top three orders was similar between the epibiotic bacterial communities of MVcarapace and SCcarapace. The bacterioplankton community in the ENV sample wasdominated by Alteromonadales (15.99%), Campylobacterales (15.03%), and Methylococcales (14.77%). Notably, Alteromonadales was only exclusively dominated the bacterioplankton community, and the proportion of Thiotrichales was much lower (3.87%; Figure 3).

FIGURE 3
www.frontiersin.org

Figure 3. Stacked bar chart of relative abundances of bacterial orders. The top 20 orders in the M. verrilli carapace (MVcarapace), S. crosnieri carapace (SCcarapace), and S. crosnieri setae (SCsetae) samples and the bacterioplankton community in the cold seep interface water (environment [ENV]) sample are shown.

To further explore the compositional differences at the ASV level, PCoA, heatmap analysis, and bacterial co-occurrence network analysis were used (Figure 4). The dominant order Campylobacterales mainly involved ASV1 (genus Sulfurovum), Thiotrichales mainly involved ASV3 (unknown species in order Thiotrichales), and Methylococcales mainly involved ASV2 (genus Marine_Methylotrophic_Group_2). ANOSIM revealed significant differences in the epibiotic bacterial communities within each of the three groups (p < 0.01). The SCsetae epibiotic bacterial community was highly separated from the others based on PCoA and heatmap clustering (Figures 4A,B). According to the STAMP analysis, the high abundance of ASV3 (order Thiotrichales), associated with sulfide oxidation, in SCsetae seemed to be a major source of variation among the groups (Supplementary Figure S3). There was a high degree of similarity between the MVcarapace and SCcarapace epibiotic bacterial communities, reflected by partial overlaps based on PCoA and heatmap clustering. There were slight differences between individuals in the MVcarapace group, as evidenced by the dispersion shown by our results (Figure 4A). The bacterial co-occurrence network, containing 144 nodes and 533 edges, was parsed into various modules, which are groups of tightly connected ASVs. The relative abundance of the module graph_3 (Figure 4C; Supplementary Figure S4), which mainly contained ASV3 (order Thiotrichales) associated with sulfur oxidation, was much higher in SCsetae than MVcarapace and SCcarapace, which might indicate an enrichment of specific sulfur-oxidizing bacteria in SCsetae.

FIGURE 4
www.frontiersin.org

Figure 4. (A) PCoA plot of M. verrilli carapace (MVcarapace), S. crosnieri carapace (SCcarapace), and S. crosnieri setae (SCsetae) epibiotic bacterial communities and the bacterioplankton community in cold seep interface water (environment [ENV]), showing the differences in the epibiotic bacterial and bacterioplankton communities at the ASV level. The SCsetae samples were clustered far apart from the others. The rank of the dissimilarities within groups of epibiotic bacteria was estimated using analysis of similarity (ANOSIM), revealing significant differences in the epibiotic bacterial communities within each of the three groups. (B) Heatmap of the top 10 ASVs associated with sulfur and methane metabolism. (C) Bacterial co-occurrence network based on the top 200 ASVs. Each node represents a bacterial ASV. The size of each node is proportional to the number of correlations and the links represent statistically significant correlations (|ρ| > 0.7, p < 0.05). The background color around the nodes indicates the modules (groups of tightly connected ASVs), i.e., graph_2, graph_3, graph_14, and graph_17, that the ASVs are affiliated with.

4. Discussion

4.1. Comparison of bacterial orders and ASVs among groups

Cold seeps harbor thriving communities of unique organisms, supported mainly by chemosynthetic microbes (such as Thiotrichales, Methylococcales, and Campylobacterales; Sogin et al., 2021). These microbes produce biomass for the higher trophic levels using energy from reducing gases such as hydrogen sulfide and methane. In this study, we present the first report of the epibiotic bacterial communities associated with M. verrilli and the first comparative study of the epibiotic bacterial communities on Munidopsidae species in cold seeps.

Despite the wide distribution of M. verrilli in the global deep-sea environment, previous studies on it were limited to morphological identification (Baba and Poore, 2002). Our results showed that chemosynthetic bacteria, including Thiotrichales, Methylococcales, and Campylobacterales, dominated the MVcarapace epibiotic bacterial community. Campylobacterales species often play a crucial role in the sulfur, nitrogen, and hydrogen cycles, as oxidizers or reducers (Wang Y. et al., 2022). Methylococcales belongs to facultative anaerobic bacteria and exhibits prominent capacity for the anaerobic oxidation of methane (Bhattarai et al., 2019). Research on serpulids at a seep known as Jaco Scar showed that Methylococcales dominated the epibiotic bacterial communities (Goffredi et al., 2020). Consistent with our results, epibiotic growth of chemoautotrophic bacteria on invertebrates is quite common in cold seeps (Martin and Haney, 2005). Similar to our results on the MVcarapace epibiotic bacterial communities, research on Munidopsis alvisca showed that chemosynthetic bacterial families were also the dominant families (Leinberger et al., 2022). Research on Alvinocarididae crustaceans at Site F showed that the epibiotic microbes had a similar composition, including Methylococcales, Campylobacterales, and Thiotrichales (Hui et al., 2022). The high relative abundance of chemosynthetic bacteria in MVcarapace in our study was in line with the extreme environment involving high hydrogen sulfide and methane concentrations (Cao et al., 2021). Most of the reported mutualistic relationships in deep-sea chemosynthetic ecosystems involve chemosynthetic bacteria with sulfide-oxidizing or methane-oxidizing activity (Levin et al., 2016). Epibiotic microbes have been reported to contribute to the Munidopsis alvisca diet, with bacterial biofilms being a relevant food source for Munidopsis sp. (Hoyoux et al., 2012). Based on the co-occurrence on M. verrilli of methanotrophic and thiotrophic bacteria with high relative abundances, we hypothesize that the epibiotic bacteria might provide nutrition to their hosts, although only a small amount, given their low absolute density.

As for S. crosnieri, its epibiotic bacteria have been intensively studied. Watsuji et al. concluded that there is a close relationship between the epibiotic bacteria of S. crosnieri and the host in terms of nutrient supply, as the epibiotic bacteria adhering to the setae (including those with sulfide-oxidizing or methane-oxidizing activity) could be scraped by the host to use as food (Watsuji et al., 2017). Our result on the SCsetae epibiotic bacterial community composition concurred with previous studies, in which Thiotrichales, Methylococcales, and Campylobacterales made up a large proportion (Xu et al., 2022).

4.2. Absence or presence of setae is the main cause of variation in epibiotic bacterial communities

The significant differences in epibiotic bacterial community morphology, density, and composition of SCsetae compared to MVcarapace and SCcarapace hinted at the function of setae in shaping the epibiotic bacterial communities.

The morphology of the dense epibiotic bacteria on S. crosnieri that we observed is consistent with previous studies (Zhang et al., 2018). Microbes often live in consortia bound to substrates, forming biofilms, which are quite common in chemosynthetic systems, e.g., euglenozoans in cold seeps are entirely covered by biofilms of closely packed epibiotic bacteria (Buck et al., 2000). Bacterial communities within a biofilm often have a multi-layered structure (Hansen et al., 2007). In our study, the SCsetae bacteria, forming a biofilm, were quite different from other epibiotic bacterial communities, with diverse taxa and a multi-layer configuration. A previous study suggested that a multi-layer biofilm configuration involving metabolic cooperation exhibited increased stability and productivity, increasing biomass accumulation due to optimal nutrient utilization (Elias and Banin, 2012). In addition, pereopods densely covered with specialized setae maximized the adhesion space for epibiotic bacteria and the dense structure intensified the adhesion to the surfaces, thereby reducing loss of the bacteria (Christensen et al., 2002; Federle, 2006). Moreover, chemical intermediates produced by reactions within the chemosynthetic microbiota maintained the stability of the biofilm in the face of temporal heterogeneity (Dolinšek et al., 2016). Research has shown that the stable intermediates of the SCsetae epibiotic bacterial community play essential roles in metabolic interactions between the hosts and bacteria in the face of temporal oscillations in resource availability (Xu et al., 2022).

The PCoA analysis revealed slight differences in the epibiotic bacteria on the MVcarapace and SCcarapace, which is possibly caused by their different distribution patterns. However, such difference was not supported by the co-occurrence network results and was not visible in the heatmaps, indicating that the microbiome difference is not as remarkable as the difference between the setae and the carapaces. More research is necessary to fully comprehend whether the setae selectively attach to the biofilms of the SC-symbionts. In addition, the overlapping distribution of some individuals of the two squat lobster species around the vents could reduce the significance of the difference between two carapaces in the statistical analysis. The similarity between the MVcarapace and SCcarapace groups suggested a weak selection of the carapaces by bacteria.

Epibiotic bacterial communities often reflect the environment and can exchange members with the bacterioplankton community in the surrounding water via direct contact (He et al., 2020). A study of Chinese mitten crabs indicated that the bacterial diversity in the water somewhat shaped the epibiotic bacterial communities (Zhang et al., 2016). In our study, almost all the orders present in the bacterioplankton community were present in the epibiotic bacterial community, The presence of plume setae with a soft texture may increase the three-dimensional space for chemosynthetic bacteria to interact with each other and may provide a more suitable adhesion substrate for certain groups. In addition, the extremely low alpha diversity in the SCsetae epibiotic bacterial community in contrast to its high density indicates an enrichment of specific bacterial groups.

To identify the enriched taxa, we performed various analyses at the ASV level. Studies have shown that mutualistic host–microbe relationships likely exist for S. crosnieri and M. verrilli (Watsuji et al., 2012; Leinberger et al., 2022). Additionally, chemosynthetic bacteria played a crucial role in the adaptation of crustaceans to cold seeps by forming close symbiotic relationships with their hosts (Cowart et al., 2017). Therefore, we focused on these chemosynthetic groups in epibiotic bacterial communities to investigate whether these groups were specifically enriched and thus may be involved in differential adaptive strategies. We found that the chemosynthetic orders (Thiotrichales, Methylococcales, and Campylobacterales) accounted for a high proportion of the epibiotic bacterial communities (>60%), while the proportions of these orders were less in the bacterioplankton community (<40%). These differences indicated that the epibiotic bacteria were not random recruited from the bacterioplankton.

The difference in the relative abundance of ASV3 (order Thiotrichales) among the epibiotic bacterial communities is striking. Thiotrichales, which are often involved in primary production in chemosynthetic communities (Cambon-Bonavita et al., 2020), can exhibit sulfur-oxidizing activity in the presence of hydrogen sulfide and its spontaneously oxidized compounds and has been widely reported in the deep-sea research field, as sulfur-oxidizing symbionts of invertebrate hosts (Hourdez et al., 2021). Dissolved sulfate and elemental sulfur are regarded as indicators of biochemical processes (Feng et al., 2016). Du et al. detected strong Raman peak of the sulfate and elemental sulfur in the fluids in the chemosynthetic communities at Site F (Du et al., 2018). The presence of these chemicals in the vicinity of cold seep could sustain the high biomass production in epibiotic bacterial communities. The relative abundance of ASV3 (order Thiotrichales) was strikingly higher for SCsetae than MVcarapace and SCcarapace, which might indicate that sulfide oxidation was quite important in the formation of epibiotic bacterial communities on S. crosnieri setae.Notably, a recent study found that chemoautotrophic bacteria belonging to Gammaproteobacteria (including the order Thiotrichales) could coevolve with their hosts, whereas Campylobacteria (such as the order Campylobacterales) had a weak preference for hosts (Lee et al., 2021). This indicates the possibility of enrichment of ASV3 (order Thiotrichales) in SCsetae and explains the high relative abundance, across all samples, of ASV1 (order Campylobacterales), which is also a sulfur-oxidizing bacteria like ASV3 (order Thiotrichales).

S. crosnieri cluster near the most active cold seeps, where fluids are expelled as a result of geological processes such as sediment compaction (Feng et al., 2018). The chemical fluxes are higher near the seeps with the most active bubble outflows (Sibuet and Olu, 1998). The multi-layer biofilm configuration and high density of the SCsetae epibiotic bacterial community helps the host to efficiently utilize the high chemical fluxes around it. Additionally, certain physical behaviors (such as dancing behavior) of the host can indirectly deepen the symbiotic relationship to increase the productivity of its episymbionts and the adaptation of the host, such as the yeti crab Kiwa puravida (Thurber et al., 2011). Research suggested that S. crosnieri exhibited similar behaviors, allowing it to drive water through the endogenous water flow around its setae (Watsuji et al., 2018). At Site F, there is a vertical gradient regarding the concentration of hydrogen sulfide, with the highest concentration at the bottom interface of the chemosynthetic communities (Cao et al., 2021). Higher concentrations of reducing substances could be produced from the bottom interface due to the endogenous water flow of the hosts. This improves the energy metabolism and the assimilation of chemosynthetic inorganic carbon under the water flow, contributing to an enrichment of sulfur-oxidizing bacteria (such as Thiotrichales species) on the setae to support host survival.

As for M. verrilli, although the relative abundance of chemoautotrophic bacteria was high, their body surfaces are lacking dense setae, and endogenous water flow and dancing behavior have not been reported or observed. These chemoautotrophic bacteria could not form as dense community as S. crosnieri. Previous studies indicated that M. verrilli tended to be opportunists and had broader feeding targets (Martin and Haney, 2005). This might indicate that, compared to S. crosnieri, a relatively small proportion of the nutrient supply of M. verrilli is derived from chemosynthesis by epibiotic microbes utilizing reducing substances (such as sulfide). We hypothesize that the presence of setae shaped the structure of the epibiotic bacterial communities and that the numerous chemosynthetic epibiotic bacteria tightly adhering to SCsetae form an ideal biofilm, in contrast to the relatively sparse and simple epibiotic bacterial communities on carapaces, improving metabolic efficiency and the overall primary productivity.

5. Conclusion

In light of the unique composition of the SCsetae group and the higher relative abundance of ASV3 (order Thiotrichales), we suggest that the morphology of S. crosnieri might allow the epibiotic bacterial community to form a bushy jungle in the plumose setae, which improves primary production and provides S. crosnieri with sufficient nutrition by selectively obtaining sulfur-oxidizing episymbionts from the environment. Consistently, S. crosnieri clusters around the vents in order to use the high chemical fluxes for chemosynthetic reactions. In contrast, M. verrilli has a relatively low dependence on epibiotic bacterial communities for its nutrient supply, so they are more dispersed. Our research showed that the different distributions of S. crosnieri and M. verrilli relate to the differences in their epibiotic bacterial communities, involving differential adaptation mechanisms around cold seeps.

Data availability statement

The data presented in the study are deposited in the NCBI repository, accession number PRJNA947936 and PRJNA954131.

Ethics statement

The studies involving animals were reviewed and approved by the Science and Technology Ethics Committee of the Institute of Oceanology Chinese Academy of Sciences.

Author contributions

WF performed the statistical analysis and wrote the first draft of the manuscript. MW and CL contributed to conception and design of the study. ZX, DD, MH, ZZ, HZ, and LF contributed to sampling. All authors contributed to the article and approved the submitted version.

Funding

This work was supported by the National Natural Science Foundation of China, grant/award number: 42076091 and 41906124, the NSFC Innovative Group Grant (No. 42221005), the National Key R&D Program of China (2022YFC2804003), the Strategic Priority Research Program of the Chinese Academy of Sciences (XDA22050303 and XDB42020401).

Acknowledgments

We thank Weiyue Liu, Longzhao Li, Yujie Yan and Yan Sun for the helping of project implementation.

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Supplementary material

The Supplementary material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmicb.2023.1197476/full#supplementary-material

References

Baba, K., and Fujita, Y. (2017). Squat lobsters of the genus Galathea associated with comatulid crinoids (Decapoda: Anomura: Galatheidae) from the Ryukyu Islands, Japan. Crustacean Res 37, 43–62.

Google Scholar

Baba, K., and Poore, G. C. B. (2002). Munidopsis (Decapoda, Anomura) from south-eastern Australia. Crustaceana 75, 231–252. doi: 10.1163/156854002760095363

CrossRef Full Text | Google Scholar

Baba, K., and Williams, A. B. (1998). New Galatheoidea (Crustacea, Decapoda, Anomura) from hydrothermal systems in the West Pacific Ocean: Bismarck archipelago and Okinawa trough. Zoosystema 20, 143–156.

Google Scholar

Bhattarai, S., Cassarini, C., and Lens, P. N. L. (2019). Physiology and distribution of Archaeal Methanotrophs that couple anaerobic oxidation of methane with sulfate reduction. Microbiol. Mol. Biol. Rev. 83:e00074-18. doi: 10.1128/MMBR.00074-18

PubMed Abstract | CrossRef Full Text | Google Scholar

Bokulich, N. A., Kaehler, B. D., Rideout, J. R., Dillon, M., Bolyen, E., Knight, R., et al. (2018). Optimizing taxonomic classification of marker-gene amplicon sequences with QIIME 2’s q2-feature-classifier plugin. Microbiome 6:90. doi: 10.1186/s40168-018-0470-z

CrossRef Full Text | Google Scholar

Bolyen, E., Rideout, J., Dillon, M., Bokulich, N., Abnet, C., Al-Ghalith, G., et al. (2019). Reproducible, interactive, scalable and extensible microbiome data science using QIIME 2. Nat. Biotechnol. 37:1091. doi: 10.1038/s41587-019-0252-6

CrossRef Full Text | Google Scholar

Buck, K. R., Barry, J. P., and Simpson, A. G. B. (2000). Monterey Bay cold seep biota: Euglenozoa with chemoautotrophic bacterial epibionts. Eur. J. Protistol. 36, 117–126. doi: 10.1016/S0932-4739(00)80029-2

CrossRef Full Text | Google Scholar

Callahan, B. J., Mcmurdie, P. J., Rosen, M. J., Han, A. W., Johnson, A. J., and Holmes, S. P. (2016). DADA2: high resolution sample inference from Illumina amplicon data. Nat. Methods 13, 581–583. doi: 10.1038/nmeth.3869

PubMed Abstract | CrossRef Full Text | Google Scholar

Cambon-Bonavita, M.-A., Aubé, J., Cueff-Gauchard, V., and Reveillaud, J. (2020). Niche partitioning in the Rimicaris exoculata holobiont: the case of the first symbiotic Zetaproteobacteria. Microbiome 9:87. doi: 10.1186/s40168-022-01257-4

CrossRef Full Text | Google Scholar

Cao, L., Lian, C., Zhang, X., Zhang, H., Wang, H., Zhou, L., et al. (2021). In situ detection of the fine scale heterogeneity of active cold seep environment of the Formosa ridge, the South China Sea. J. Mar. Syst. 218:103530. doi: 10.1016/j.jmarsys.2021.103530

CrossRef Full Text | Google Scholar

Cheng, J., Hui, M., Li, Y. L., and Sha, Z. L. (2020). Genomic evidence of population genetic differentiation in deep-sea squat lobster Shinkaia crosnieri (crustacea: Decapoda: Anomura) from northwestern Pacific hydrothermal vent and cold seep. Deep Sea Res I Oceanogr Res Pap. 156:103188. doi: 10.1016/j.dsr.2019.103188

CrossRef Full Text | Google Scholar

Chevaldonne, P., and Olu, K. (1996). Occurrence of anomuran crabs (Crustacea: Decapoda) in hydrothermal vent and cold-seep communities: a review. Proc. Biol. Soc. Wash. 109, 286–298.

Google Scholar

Christensen, B. B., Haagensen, J. A. J., Heydorn, A., and Molin, S. (2002). Metabolic commensalism and competition in a two-species microbial consortium. Appl. Environ. Microbiol. 68, 2495–2502. doi: 10.1128/AEM.68.5.2495-2502.2002

PubMed Abstract | CrossRef Full Text | Google Scholar

Cowart, D. A., Durand, L., Cambon-Bonavita, M.-A., and Arnaud-Haond, S. (2017). Investigation of bacterial communities within the digestive organs of the hydrothermal vent shrimp Rimicaris exoculata provide insights into holobiont geographic clustering. PLoS One 12:e0172543. doi: 10.1371/journal.pone.0172543

PubMed Abstract | CrossRef Full Text | Google Scholar

Dattagupta, S., Schaperdoth, I., Montanari, A., Mariani, S., Kita, N., Valley, J. W., et al. (2009). A novel symbiosis between chemoautotrophic bacteria and a freshwater cave amphipod. ISME J. 3, 935–943. doi: 10.1038/ismej.2009.34

PubMed Abstract | CrossRef Full Text | Google Scholar

Dolinšek, J., Goldschmidt, F., and Johnson, D. R. (2016). Synthetic microbial ecology and the dynamic interplay between microbial genotypes. FEMS Microbiol. Rev. 40, 961–979. doi: 10.1093/femsre/fuw024

CrossRef Full Text | Google Scholar

Dong, D., and Li, X. (2015). Galatheid and chirostylid crustaceans (Decapoda: Anomura) from a cold seep environment in the northeastern South China Sea. Zootaxa 4057, 91–105. doi: 10.11646/zootaxa.4057.1.5

PubMed Abstract | CrossRef Full Text | Google Scholar

Du, Z., Zhang, X., Luan, Z.-D., Wang, M., Xi, S., Li, L., et al. (2018). In situ Raman quantitative detection of the cold seep vents and fluids in the chemosynthetic communities in the South China Sea. Geochemistry 19, 2049–2061. doi: 10.1029/2018GC007496

CrossRef Full Text | Google Scholar

Elias, S., and Banin, E. (2012). Multi-species biofilms: living with friendly neighbors. FEMS Microbiol. Rev. 36, 990–1004. doi: 10.1111/j.1574-6976.2012.00325.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Fahy, B. G., Vasilopoulos, T., Ford, S., Gravenstein, D., and Enneking, F. K. (2015). A single consent for serial anesthetics in burn surgery. Anesth. Analg. 121, 219–222. doi: 10.1213/ANE.0000000000000780

PubMed Abstract | CrossRef Full Text | Google Scholar

Federle, W. (2006). Why are so many adhesive pads hairy? J. Exp. Biol. 209, 2611–2621. doi: 10.1242/jeb.02323

PubMed Abstract | CrossRef Full Text | Google Scholar

Feng, D., Peng, Y., Bao, H., Peckmann, J., Roberts, H. H., and Chen, D. (2016). A carbonate-based proxy for sulfate-driven anaerobic oxidation of methane. Geology 44, 999–1002. doi: 10.1130/G38233.1

CrossRef Full Text | Google Scholar

Feng, D., Qiu, J.-W., Hu, Y., Peckmann, J., Guan, H., Tong, H., et al. (2018). Cold seep systems in the South China Sea: an overview. J. Asian Earth Sci. 168, 3–16. doi: 10.1016/j.jseaes.2018.09.021

CrossRef Full Text | Google Scholar

Goffredi, S. K., Jones, W. J., Erhlich, H., Springer, A., and Vrijenhoek, R. C. (2008). Epibiotic bacteria associated with the recently discovered yeti crab, Kiwa hirsuta. Environ. Microbiol. 10, 2623–2634. doi: 10.1111/j.1462-2920.2008.01684.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Goffredi, S. K., Tilic, E., Mullin, S. W., Dawson, K. S., Keller, A., Lee, R. W., et al. (2020). Methanotrophic bacterial symbionts fuel dense populations of deep-sea feather duster worms (Sabellida, Annelida) and extend the spatial influence of methane seepage. Sci. Adv. 6:6. doi: 10.1126/sciadv.aay8562

CrossRef Full Text | Google Scholar

Hansen, S. K., Haagensen, J., Gjermansen, M., Jørgensen, T. M., Tolker-Nielsen, T., and Molin, S. (2007). Characterization of a Pseudomonas putida rough variant evolved in a mixed-species biofilm with Acinetobacter sp. strain C6. J. Bacteriol. 189, 4932–4943. doi: 10.1128/JB.00041-07

PubMed Abstract | CrossRef Full Text | Google Scholar

He, D., Ren, L., and Wu, Q. L. (2020). Growing season drives the compositional changes and assembly processes of epiphytic bacterial communities of two submerged macrophytes in Taihu Lake. FEMS Microbiol. Ecol. 96:fiaa025. doi: 10.1093/femsec/fiaa025

PubMed Abstract | CrossRef Full Text | Google Scholar

Hourdez, S., Boidin-Wichlacz, C., Jollivet, D., Massol, F., Rayol, M. C., Bruno, R., et al. (2021). 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. Sci. Total Environ. 798:149149. doi: 10.1016/j.scitotenv.2021.149149

PubMed Abstract | CrossRef Full Text | Google Scholar

Hoyoux, C., Zbinden, M., Samadi, S., Gaill, F., and Compère, P. (2009). Wood-based diet and gut microXora of a galatheid crab associated with PaciWc deep-sea wood falls. Mar. Biol. 156, 2421–2439. doi: 10.1007/s00227-009-1266-2

CrossRef Full Text | Google Scholar

Hoyoux, C., Zbinden, M., Samadi, S., Gaill, F., and Compère, P. (2012). Diet and gut microorganisms of Munidopsis squat lobsters associated with natural woods and mesh-enclosed substrates in the deep South Pacific. Mar. Biol. Res. 8, 28–47. doi: 10.1080/17451000.2011.605144

CrossRef Full Text | Google Scholar

Hui, M., Wang, A., Cheng, J., and Sha, Z. (2022). Full-length 16S rRNA amplicon sequencing reveals the variation of epibiotic microbiota associated with two shrimp species of alvinocarididae: possibly co-determined by environmental heterogeneity and specific recognition of hosts. PeerJ 10:e13758. doi: 10.7717/peerj.13758

PubMed Abstract | CrossRef Full Text | Google Scholar

Katoh, K., and Standley, D. M. (2013). MAFFT multiple sequence alignment software version 7: improvements in performance and usability. Mol. Biol. Evol. 30, 772–780. doi: 10.1093/molbev/mst010

PubMed Abstract | CrossRef Full Text | Google Scholar

Komai, T., and Matsuzaki, K. (2016). Two deep-sea decapod crustaceans collected off eastern Hokkaido, Japan: Sclerocrangon rex n. sp. (Caridea: Crangonidae) and Munidopsis verrilli Benedict, 1902 (Anomura: Munidopsidae). Zootaxa 4162, 92–106. doi: 10.11646/zootaxa.4162.1.4

PubMed Abstract | CrossRef Full Text | Google Scholar

Komai, T., Shank, T. M., and Dover, C. L. V. (2005). A new species of Alvinocaris (Crustacea: Decapoda: Caridea: Alvinocarididae) and a new record of A. muricola from methane seeps on the Blake ridge Diapir, Northwestern Atlantic. Zootaxa 1019, 27–42. doi: 10.11646/zootaxa.1019.1.2

CrossRef Full Text | Google Scholar

Lee, W.-K., Juniper, S. K., Perez, M., Ju, S.-J., and Kim, S.-J. (2021). Diversity and characterization of bacterial communities of five co-occurring species at a hydrothermal vent on the Tonga arc. Ecol. Evol. 11, 4481–4493. doi: 10.1002/ece3.7343

PubMed Abstract | CrossRef Full Text | Google Scholar

Leinberger, J., Milke, F., Christodoulou, M., Poehlein, A., Caraveo-Patino, J., Teske, A., et al. (2022). Microbial epibiotic community of the deep-sea galatheid squat lobster Munidopsis alvisca. Sci. Rep. 12:2675. doi: 10.1038/s41598-022-06666-x

PubMed Abstract | CrossRef Full Text | Google Scholar

Levin, L. A., Baco, A. R., Bowden, D. A., Colaço, A., Cordes, E., Cunha, M. R., et al. (2016). Hydrothermal vents and methane seeps: rethinking the sphere of influence. Front. Mar. Sci. 3:72. doi: 10.3389/fmars.2016.00072

CrossRef Full Text | Google Scholar

Liu, D., and Tong, C. (2017). Bacterial community diversity of traditional fermented vegetables in China. LWT Food Sci. Technol. 86, 40–48. doi: 10.1016/j.lwt.2017.07.040

CrossRef Full Text | Google Scholar

Macavoy, S. E., Carney, R. S., Morgan, E., and Macko, S. A. (2008). Stable isotope variation among the mussel Bathymodiolus childressi and associated heterotrophic Fauna at four cold-seep communities in the Gulf of Mexico. J. Shellfish Res. 27, 147–151. doi: 10.2983/0730-8000(2008)27[147:SIVATM]2.0.CO;2

CrossRef Full Text | Google Scholar

Macpherson, E., Jones, W., and Segonzac, M. (2005). A new squat lobster family of Galatheoidea (Crustacea, Decapoda: Anomura) from the hydrothermal vents of the Pacific-Antarctic ridge. Zoosystema 27, 709–723.

Google Scholar

Martin, J. W., and Haney, T. A. (2005). Decapod crustaceans from hydrothermal vents and cold seeps: a review through 2005. Zool. J. Linnean Soc. 145, 445–522. doi: 10.1111/j.1096-3642.2005.00178.x

CrossRef Full Text | Google Scholar

Mioduchowska, M., Czyż, M. J., Gołdyn, B., Kilikowska, A., Namiotko, T., Pinceel, T., et al. (2018). Detection of bacterial endosymbionts in freshwater crustaceans: the applicability of non-degenerate primers to amplify the bacterial 16S rRNA gene. Peer J. 6:e6039. doi: 10.7717/peerj.6039

CrossRef Full Text | Google Scholar

Nadkarni, M. A., Martin, F. E., Jacques, N. A., and Hunter, N. (2002). Determination of bacterial load by real-time PCR using a broad-range (universal) probe and primers set. Microbiology 148, 257–266. doi: 10.1099/00221287-148-1-257

PubMed Abstract | CrossRef Full Text | Google Scholar

Niemann, H., Linke, P., Knittel, K., Macpherson, E., Boetius, A., Ckmann, W. B., et al. (2013). Methane-carbon flow into the benthic food web at cold seeps – a case study from the Costa Rica Subduction zone. PLoS One 8:10. doi: 10.1371/journal.pone.0074894

CrossRef Full Text | Google Scholar

Parks, D. H., Tyson, G. W., Hugenholtz, P., and Beiko, R. G. (2014). STAMP: statistical analysis of taxonomic and functional profiles. Bioinformatics 30, 3123–3124. doi: 10.1093/bioinformatics/btu494

CrossRef Full Text | Google Scholar

Price, M. N., Dehal, P. S., and Arkin, A. P. (2010). FastTree 2 – approximately maximum-likelihood trees for large alignments. PLoS One 5:e9490. doi: 10.1371/journal.pone.0009490

PubMed Abstract | CrossRef Full Text | Google Scholar

Quast, C., Pruesse, E., Yilmaz, P., Gerken, J., Schweer, T., Yarza, P., et al. (2012). The SILVA ribosomal RNA gene database project: improved data processing and web-based tools. Nucleic Acids Res. 41, D590–D596. doi: 10.1093/nar/gks1219

PubMed Abstract | CrossRef Full Text | Google Scholar

Schwartzman, J. A., Koch, E., Heath-Heckman, E. A. C., Zhou, L., Kremer, N., Mcfall-Ngai, M. J., et al. (2015). The chemistry of negotiation: rhythmic, glycan-driven acidification in a symbiotic conversation. Proc. Natl. Acad. Sci. U. S. A. 112, 566–571. doi: 10.1073/pnas.1418580112

PubMed Abstract | CrossRef Full Text | Google Scholar

Sibuet, M., and Olu, K. (1998). Biogeography, biodiversity and fluid dependence of deep-sea cold-seep communities at active and passive margins. Deep Sea Res II 45, 517–567. doi: 10.1016/S0967-0645(97)00074-X

CrossRef Full Text | Google Scholar

Sogin, M. E., Kleiner, M., Borowski, C., Gruber-Vodicka, H. R., and Dubilier, N. (2021). Life in the dark: phylogenetic and physiological diversity of chemosynthetic symbioses. Ann. Rev. Microbiol. 75, 695–718. doi: 10.1146/annurev-micro-051021-123130

CrossRef Full Text | Google Scholar

Thurber, A. R., Jones, W. J., and Schnabel, K. (2011). Dancing for food in the Deep Sea: bacterial farming by a new species of yeti crab. PLoS One 6:e26243. doi: 10.1371/journal.pone.0026243

PubMed Abstract | CrossRef Full Text | Google Scholar

Tsuchida, S., Fujiwara, Y., and Fujikura, K. (2003). Distribution and population structure of the Galatheid crab Shinkaia crosnieri (Decapoda: Anomura: Galatheidae) in the southern Okinawa trough. Plankon Benthos Res. 58, 84–88. doi: 10.5179/benthos.58.84

CrossRef Full Text | Google Scholar

Wang, Y., Bi, H.-Y., Chen, H., Zheng, P., Zhou, Y.-L., and Li, J. (2022). Metagenomics reveals dominant unusual sulfur oxidizers inhabiting active hydrothermal chimneys from the southwest Indian ridge. Front. Microbiol. 13:861795. doi: 10.3389/fmicb.2022.861795

CrossRef Full Text | Google Scholar

Wang, X., Guan, H., Qiu, J.-W., Xu, T., Peckmann, J., Chen, D., et al. (2022). Macro-ecology of cold seeps in the South China Sea. Geosyst. Geoenviron. 1:100081. doi: 10.1016/j.geogeo.2022.100081

CrossRef Full Text | Google Scholar

Warrant, E. J., and Locket, N. A. (2004). Vision in the deep sea. Biol. Rev. 79, 671–712. doi: 10.1017/S1464793103006420

CrossRef Full Text | Google Scholar

Watsuji, T.-O., Motoki, K., Hada, E., Nagai, Y., Takaki, Y., Yamamoto, A., et al. (2018). Compositional and functional shifts in the Epibiotic bacterial Community of Shinkaia crosnieri Baba & Williams (a squat lobster from hydrothermal vents) during methane-fed rearing. Microbes Environ. 33, 348–356. doi: 10.1264/jsme2.ME18072

CrossRef Full Text | Google Scholar

Watsuji, T.-O., Nishizawa, M., Morono, Y., Hirayama, H., Kawagucci, S., Takahata, N., et al. (2012). Cell-specific thioautotrophic productivity of epsilon-proteobacterial epibionts associated with Shinkaia crosnieri. PLoS One :7:e46282. doi: 10.1371/journal.pone.0046282

CrossRef Full Text | Google Scholar

Watsuji, T.-O., Tsubaki, R., Chen, C., Nagai, Y., Nakagawa, S., Yamamoto, M., et al. (2017). Cultivation mutualism between a deep-sea vent galatheid crab and its chemosynthetic epibionts. Deep Sea Res. I Oceanogr. Res. Pap. 127, 13–20. doi: 10.1016/j.dsr.2017.04.012

CrossRef Full Text | Google Scholar

Watsuji, T.-O., Yamamoto, A., Motoki, K., Ueda, K., Hada, E., Takaki, Y., et al. (2015). Molecular evidence of digestion and absorption of epibiotic bacterial community by deep-sea crab Shinkaia crosnieri. ISME J. 9, 821–831. doi: 10.1038/ismej.2014.178

CrossRef Full Text | Google Scholar

Wehrtmann, I. S., and Acuna, E. (2011). “Squat Lobsters Fisheries,” in The Biology of Squat Lobsters. eds. G. C. B. Poore, S. Ahyong and J. Taylor (Museum Victoria: CSIRO Publishing), 297–322.

Google Scholar

Weisburg, W. G., Barns, S. M., Pelletier, D. A., and Lane, D. J. W. (1991). 16S ribosomal DNA amplification for phylogenetic study. J. Bacteriol. 173, 697–703. doi: 10.1128/jb.173.2.697-703.1991

PubMed Abstract | CrossRef Full Text | Google Scholar

Wen, T., Xie, P., Yang, S., Niu, G., Liu, X., Ding, Z., et al. (2022). ggClusterNet: an R package for microbiome network analysis and modularity-based multiple network layouts. iMeta. 1:e32. doi: 10.1002/imt2.32

CrossRef Full Text | Google Scholar

Williams, A. B., Smith, C. R., and Baco, A. R. (2000). New species of Paralomis (Decapoda, Anomura, Lithodidae) from a sunken whale carcass in the San Clemente Basin off southern California. J. Crustac. Biol. 20, 281–285. doi: 10.1163/1937240X-90000030

CrossRef Full Text | Google Scholar

Xiao, Y., Xu, T., Sun, J., Wang, Y., Wong, W. C., Kwan, Y. H., et al. (2020). Population genetic structure and gene expression plasticity of the Deep-Sea vent and seep squat lobster Shinkaia crosnieri. Front. Mar. Sci. 7:587686. doi: 10.3389/fmars.2020.587686

CrossRef Full Text | Google Scholar

Xu, Z., Wang, M., Zhang, H., He, W., Cao, L., Lian, C., et al. (2022). Metabolism interactions promote the overall functioning of the episymbiotic chemosynthetic community of Shinkaia crosnieri of cold seeps. mSystems 7:4. doi: 10.1128/msystems.00320-22

CrossRef Full Text | Google Scholar

Zhang, M. L., Sun, Y., Chen, L., Cai, C., Qiao, F., Du, Z. Y., et al. (2016). Symbiotic Bacteria in gills and guts of Chinese mitten crab (Eriocheir sinensis) differ from the free-living Bacteria in water. PLoS One 11:e0148135. doi: 10.1371/journal.pone.0148135

CrossRef Full Text | Google Scholar

Zhang, J., Zeng, Z.-G., Chen, S., and Sun, L. (2018). Bacterial communities associated with Shinkaia crosnieri from the Iheya north, Okinawa trough: microbial diversity and metabolic potentials. J. Mar. Syst. 180, 228–236. doi: 10.1016/j.jmarsys.2017.02.011

CrossRef Full Text | Google Scholar

Zhao, Y., Xu, T., Law, Y. S., Feng, D., Li, N., Xin, R. J., et al. (2020). Ecological characterization of cold-seep epifauna in the South China Sea. Deep Sea Res. I Oceanogr. Res. Pap. 163:103361. doi: 10.1016/j.dsr.2020.103361

CrossRef Full Text | Google Scholar

Keywords: cold seep, squat lobster, epibiotic microbiome, sulfide, adaptation

Citation: Feng W, Wang M, Dong D, Hui M, Zhang H, Fu L, Zhong Z, Xu Z and Li C (2023) Variation in epibiotic bacteria on two squat lobster species of Munidopsidae. Front. Microbiol. 14:1197476. doi: 10.3389/fmicb.2023.1197476

Received: 31 March 2023; Accepted: 25 May 2023;
Published: 28 June 2023.

Edited by:

David Emerson, Bigelow Laboratory for Ocean Sciences, United States

Reviewed by:

Jiwen Liu, Ocean University of China, China
Hongmei Jing, Institute of Deep-Sea Science and Engineering, China

Copyright © 2023 Feng, Wang, Dong, Hui, Zhang, Fu, Zhong, Xu and Li. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Chaolun Li, lcl@qdio.ac.cn

Disclaimer: All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.