# Colonization of plant substrates at hydrothermal vents and cold seeps in the northeast Atlantic and Mediterranean and occurrence of symbiont-related bacteria

Kamil M. Szafranski, Philippe Deschamps, Marina R. Cunha, Sylvie Marylène Gaudron, Sébastien Duperron

<figure><img src="image_1.png" /></figure>

## To cite this version:

Kamil M. Szafranski, Philippe Deschamps, Marina R. Cunha, Sylvie Marylène Gaudron, Sébastien Duperron. Colonization of plant substrates at hydrothermal vents and cold seeps in the northeast Atlantic and Mediterranean and occurrence of symbiont-related bacteria. Frontiers in Microbiology, 2015, 6, pp.162. 10.3389/fmicb.2015.00162. hal-01221485

HAL Id: hal-01221485

https://hal.sorbonne-universite.fr/hal-01221485

Submitted on 28 Oct 2015

HAL is a multi-disciplinary open access archive for the deposit and dissemination of scientific research documents, whether they are published or not. The documents may come from teaching and research institutions in France or abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est destinée au dépôt et à la diffusion de documents scientifiques de niveau recherche, publiés ou non, émanant des établissements d'enseignement et de recherche français ou étrangers, des laboratoires publics ou privés.

<figure><img src="image_2.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_3.png" /></figure>

# Colonization of plant substrates at hydrothermal vents and cold seeps in the northeast Atlantic and Mediterranean and occurrence of symbiont-related bacteria

Kamil M. Szafranski<sup>1,2</sup>, Philippe Deschamps<sup>3</sup>, Marina R. Cunha<sup>4</sup>, Sylvie M. Gaudron<sup>1,2†</sup> and Sébastien Duperron<sup>1,2,5*</sup>

**Edited by:**

Peter R. Girguis, Harvard University,
USA

**Reviewed by:**

Sang-Jin Kim, Korea Ocean
Research and Development
Institute, South Korea
Julie A. Huber, Marine Biological
Laboratory, USA

**\*Correspondence:**

Sébastien Duperron, UMR 7208
BOREA, Adaptation aux Milieux
Extrêmes, Université Pierre et Marie
Curie, Bât A, 4ème étage, 7 quai
St. Bernard, 75005 Paris, France
e-mail: sebastien.duperron@
snv.jussieu.fr

**†Present address:**

Sylvie M. Gaudron, UMR8187
Laboratoire d'Océanologie et de
Géosciences (UL1 Centre National
de la Recherche Scientifique,
University of the Littoral Opal
Coast), Station Marine de
Wimereux, Wimereux, France

Reducing conditions with elevated sulfide and methane concentrations in ecosystems such as hydrothermal vents, cold seeps or organic falls, are suitable for chemosynthetic primary production. Understanding processes driving bacterial diversity, colonization and dispersal is of prime importance for deep-sea microbial ecology. This study provides a detailed characterization of bacterial assemblages colonizing plant-derived substrates using a standardized approach over a geographic area spanning the North-East Atlantic and Mediterranean. Wood and alfalfa substrates in colonization devices were deployed for different periods at 8 deep-sea chemosynthesis-based sites in four distinct geographic areas. Pyrosequencing of a fragment of the 16S rRNA-encoding gene was used to describe bacterial communities. Colonization occurred within the first 14 days. The diversity was higher in samples deployed for more than 289 days. After 289 days, no relation was observed between community richness and deployment duration, suggesting that diversity may have reached saturation sometime in between. Communities in long-term deployments were different, and their composition was mainly influenced by the geographical location where devices were deployed. Numerous sequences related to horizontally-transmitted chemosynthetic symbionts of metazoans were identified. Their potential status as free-living forms of these symbionts was evaluated based on sequence similarity with demonstrated symbionts. Results suggest that some free-living forms of metazoan symbionts or their close relatives, such as Epsilon proteobacteria associated with the shrimp *Rimicaris exoculata*, are efficient colonizers of plant substrates at vents and seeps.

**Keywords:** cold seeps, colonization, deep-sea, symbiont, hydrothermal vents, wood falls

## INTRODUCTION

Hydrothermal vents and other chemosynthesis-based ecosystems like cold seeps are deep-sea hotspots of primary production, which is ensured by chemoautotrophic prokaryotes, many of which live in symbiosis with various metazoans (Van Dover, 2000; Dubilier et al., 2008). The deep-sea ecosystem can also benefit from significant input of exogenous sources of carbon such as plant remains of different origins, like sunken wood (Wolff, 1979). Wood falls contain cellulose, hemicellulose and lignin, i.e., energy-rich polysaccharides that can be degraded by microbial assemblages. This results in the production of reduced compounds such as hydrogen sulfide, dihydrogen and methane, providing a niche for thiotrophic and methanotrophic bacteria (Leschine, 1995; Palacios et al., 2006). Wood falls, but also whale skeletons and others large inputs of organic material, are quickly localized and efficiently colonized

by opportunistic fauna (Turner, 1973; Bennett et al., 1994; Baco and Smith, 2003; Smith and Baco, 2003; Laurent et al., 2013). Reducing conditions with elevated sulfide concentrations around degrading organic falls also constitute a prime environment for metazoans with sulfide-oxidizing symbioses (Duperron et al., 2008; Lorion et al., 2009; reviewed in: Dubilier et al., 2008; Laurent et al., 2009; Treude et al., 2009). Seep and vent ecosystems are distributed worldwide but are often separated by large distances, and various authors postulated that wood and whale falls could act as stepping stones between habitats (Vrijenhoek, 2010). Nevertheless, the most successful animal colonizers remain specialists for one substrate type, remarkable examples including *Xylophaga* spp. or *Xyloredo* spp. bivalves on wood (Turner, 1973; Distel and Roberts, 1997) and *Osedax* spp. polychaetes on bones (Rouse et al., 2004; Glover et al., 2005).

<footnote n="1">¹ Sorbonne Universités, UPMC Univ. Paris 06, UMR 7208, Adaptation aux Milieux Extrêmes, Paris, France</footnote>

<footnote n="2">² UMR MNHN UPMC CNRS IRD UCBN 7208, Biologie des Organismes Aquatiques et Ecosystèmes, Paris, France</footnote>

<footnote n="3">³ UMR8079 Unité d'Ecologie, Systématique et Evolution, CNRS Université Paris-Sud 11, Orsay, France</footnote>

<footnote n="4">⁴ Departamento de Biologia and CESAM, Universidade de Aveiro, Aveiro, Portugal</footnote>

<footnote n="5">⁵ Institut Universitaire de France, Paris, France</footnote>
Despite the significance of sunken plant substrates as locale organic enrichments, studies on the associated microbial assemblages remain scarce. The first trial to investigate and compare free-living and attached prokaryotic assemblages from sunken wood at several sites and water depths was carried out by Palacios et al. (2009). Authors combined electron microscopy and molecular fingerprinting (capillary electrophoresis single-stranded conformation polymorphisms—CE-SSCP) to test the influence of depth of immersion, geographic location, deployment time, and wood type on the structure of microbial assemblages (Palacios et al., 2009). The phylogeny and diversity of Bacteria and Archaea associated with wood falls was recently investigated using clone libraries (Fagervold et al., 2012), and authors demonstrated the occurrence of various free-living sulfate-reducing and sulfur-oxidizing bacteria (SRB and SOX respectively) and methanogenic archaea. This microbial community was modified with time of immersion because of changes in wood chemistry. The combination of ARISA fingerprinting and 454 pyrosequencing recently allowed Bienhold et al. (2013) to gain a more detailed overview of bacterial assemblages in a pine wood deployment lasting 1 year and surrounding sediments at a cold seep site in the Eastern Mediterranean Sea. They established that wood-boring bivalves, cellulolytic and sulfate-reducing bacteria colonized the organic substrate first, and then attracted chemosymbiotic fauna (Bienhold et al., 2013). Finally, using pyrosequencing, Fagervold et al. (2013) also demonstrated that the microbial diversity in oak cubes deployed at the Blanes Canyon (Western Mediterranean Sea) was significantly higher compared to cubes deployed in the adjacent open slope (Fagervold et al., 2013). To date, these high-throughput diversity studies include a single region, use only one substrate type, and are all based on different experimental designs, and as a consequence knowledge of deep-sea wood-associated microbial assemblages is still scarce.

Organic fall-associated and vent/seep faunas include many metazoans harboring chemosynthetic symbionts, and many of these symbionts are environmentally acquired anew at each generation (Gros et al., 1996; McFall-Ngai, 2000; Won et al., 2003; Nussbaumer et al., 2006; Duperron et al., 2008; Lorion et al., 2009; reviewed in: Dubilier et al., 2008; Bright and Bulgheresi, 2010; Vrijenhoek, 2010). This mode of symbiont transmission implies the existence of free-living forms of symbionts which may infect the appropriate host, and for which symbiotic lifestyle may be facultative (Bright and Bulgheresi, 2010). Only few studies successfully documented the occurrence of such free-living symbionts in marine shallow waters (Lee and Ruby, 1992; Gros et al., 2003; Aida et al., 2008), even less in the deep-sea (Harmer et al., 2008). In which habitat and numbers do these free-living forms exist thus remains an open question.

Since 2006, colonization devices (CHEMECOLIs) filled with wood cubes and alfalfa grass have been deployed for various periods of time at several chemosynthesis-based sites located in the Eastern Mediterranean, the Gulf of Cadiz, the Mid-Atlantic Ridge and the Haakon Mosby Mud Volcano (Gaudron et al., 2010; Cunha et al., 2013). In this study, we characterize the bacterial communities colonizing these organic substrates. The aims are to identify bacteria colonizing the substrates, to evaluate the effect of region, type of substrate, depth, duration of deployment,

and temperature on the bacterial assemblages, and to screen for potential free-living forms or relatives of metazoan-associated chemosynthetic symbionts. For this, a 454 pyrosequencing-based approach on a 16S rRNA-encoding gene fragment is employed.

# MATERIALS AND METHODS

## SAMPLING

Thirteen sets of standardized colonization devices (CHEMECOLI: CHEMosynthetic Ecosystem Colonization by Larval Invertebrates, described by Gaudron et al., 2010) were deployed at 8 hydrothermal vents and cold seeps sites located in 4 distinct geographic areas (Figure 1): the Mid-Atlantic Ridge (MAR, 3 sites), the Gulf of Cadiz (GoC, 3 sites), the Eastern Mediterranean (EM, 1 site), and the Norwegian Sea (Haakon Mosby Mud Volcano, HMMV, 1 site). Each set consisted of two CHEMECOLIs, one filled with dried alfalfa grass (A) and the second with 2 cm pine wood cubes (W). Deployment sites were located at depths ranging from 354 to 2300 m. CHEMECOLIs were deployed a few meters away from any visible fluid escape and recovered in individual sterile-water filled hermetic boxes after periods of 10 days to 3 years. Operations were performed by various ROVs (Remotely Operated Vehicles) over the course of 13 cruises spanning over years 2006–2013 (Table 1). Immediately after recovery, pieces of alfalfa and wood cubes were randomly selected in a cold room, fixed in 96% ethanol and stored at 4°C until DNA extraction. Substrates from a short-term deployment at Menez Gwen were also frozen in liquid nitrogen and stored at –80°C in order to compare different sample fixation procedures (MAR-A-V-C, MAR-W-V-C). Non-deployed substrates were used as controls (Contr-A and Contr-W).

## DNA EXTRACTION AND 454 PYROSEQUENCING

Samples were screened under a dissecting microscope to check for the absence of animal material (tissue, eggs, and larvae). Small pieces (2–3 g) of wood cubes were cut in thin slices and homogenized (3 × 30 s) in a metal jar with beads (12 mm diameter) using

<figure><img src="image_4.png" /><figcaption>FIGURE 1 | Map of sampling sites and their depths, located in 4 geographical regions.</figcaption></figure>
<table><caption>Table 1 | Main characteristics of samples used for sequence analyses.</caption><thead><tr><th>Alfalfa</th><th>Wood</th><th>Sample name</th><th>Region</th><th>Site</th><th>Depth</th><th>Days</th><th>Temperature</th><th>Site localization</th><th>Negative control</th><th>Deployment date</th><th>Cruise</th><th>Recovery date</th><th>Cruise</th></tr></thead><tbody><tr><td colspan="14">Contr-A -</td></tr><tr><td>GoC-A-I</td><td>GoC-W-I</td><td>Gulf of Cadiz</td><td>Mercator</td><td>354 m</td><td>289</td><td>12</td><td>35°-17.916°N</td><td>6°-38.709°W</td><td>19/05/2007</td><td>JC10</td><td>03/03/2008</td><td>64PE284</td></tr><tr><td>GoC-A-II</td><td>GoC-W-II</td><td>Gulf of Cadiz</td><td>Mercator</td><td>354 m</td><td>731</td><td>12</td><td>35°-17.916°N</td><td>6°-38.709°W</td><td>19/05/2007</td><td>JC10</td><td>19/05/2009</td><td>Belgica cruise ST0914</td></tr><tr><td>GoC-A-III</td><td>GoC-W-III</td><td>Gulf of Cadiz</td><td>Meknes</td><td>698 m</td><td>445</td><td>11</td><td>34°-59.09°N</td><td>7°-04.42°W</td><td>01/03/2008</td><td>64PE284</td><td>20/05/2009</td><td>Belgica cruise ST0914</td></tr><tr><td>GoC-A-IV</td><td>GoC-W-IV</td><td>Gulf of Cadiz</td><td>Darwin</td><td>1110 m</td><td>729</td><td>10</td><td>35°-23.523°N</td><td>7°-11.513°W</td><td>21/05/2007</td><td>JC10</td><td>19/05/2009</td><td>Belgica cruise ST0914</td></tr><tr><td>EM-A-I</td><td>EM-W-I</td><td>Eastern Mediterranean</td><td>NDSF</td><td>1693 m</td><td>357</td><td>13</td><td>32°-31.9772°N</td><td>30°-21.1779°E</td><td>18/11/2006</td><td>Bionil</td><td>10/11/2007</td><td>Medeco</td></tr><tr><td>MAR-A-I</td><td>MAR-W-I</td><td>Mid-Atlantic Ridge</td><td>Rainbow</td><td>2300 m</td><td>328</td><td>3.7</td><td>36°-13.766°N</td><td>33°-54.192°W</td><td>21/08/2006</td><td>Momareto</td><td>15/07/2007</td><td>Momar dream</td></tr><tr><td>MAR-A-II</td><td>MAR-W-II</td><td>Mid-Atlantic Ridge</td><td>Rainbow</td><td>2279 m</td><td>414</td><td>3.7</td><td>36°-13.7454°N</td><td>33°-54.0513°W</td><td>14/07/2007</td><td>Momardream</td><td>31/08/2008</td><td>MOMAR08</td></tr><tr><td>MAR-A-III</td><td>MAR-W-III</td><td>Mid-Atlantic Ridge</td><td>Lucky Strike</td><td>1696 m</td><td>1112</td><td>4.6</td><td>37°-17.339°N</td><td>32°-16.533°W</td><td>04/09/2006</td><td>Momareto</td><td>20/09/2009</td><td>Bathyluck</td></tr><tr><td>HMMV-A-I</td><td>HMMV-W-I</td><td>Haakon-Mosby Mud Volcano</td><td>HMMV</td><td>1257 m</td><td>388</td><td>2</td><td>72°-00.1443°N</td><td>14°-43.2282°E</td><td>06/06/2006</td><td>Vicking</td><td>29/06/2007</td><td>ARK XXII/1b</td></tr><tr><td>HMMV-A-II</td><td>HMMV-W-II</td><td>Haakon-Mosby Mud Volcano</td><td>HMMV</td><td>1257 m</td><td>752</td><td>2</td><td>72°-00.3390°N</td><td>14°-43.2209°E</td><td>03/07/2007</td><td>ARK XXII/1b</td><td>24/07/2009</td><td>ARKXXIV/2</td></tr><tr><td>MAR-A-IV</td><td>MAR-W-IV</td><td>Mid-Atlantic Ridge</td><td>Rainbow</td><td>2275 m</td><td>10</td><td>3.7</td><td>36°-13.7553°N</td><td>33°-54.11°W</td><td>31/08/2008</td><td>MOMAR08</td><td>10/09/2008</td><td>MOMAR08</td></tr><tr><td>-</td><td>EM-W-II</td><td>Eastern Mediterranean</td><td>NDSF</td><td>1697 m</td><td>14</td><td>13</td><td>32°-31.993°N</td><td>30°-21.193°E</td><td>04/11/2006</td><td>Bionil</td><td>18/11/2006</td><td>Bionil</td></tr><tr><td>MAR-A-V</td><td>MAR-W-V</td><td>Mid-Atlantic Ridge</td><td>Menez Gwen</td><td>860 m</td><td>13</td><td>4</td><td>37°-50.6786°N</td><td>31°-31.1493°W</td><td>05/08/2013</td><td>BioBaz</td><td>18/08/2013</td><td>BioBaz</td></tr></tbody></table>

Environmental data for each site were used in redundancy analysis (RDA). Sample nomenclature is as follows: region-substrate-deployment reference. Sites are: Contr, control; GoC, Gulf of Cadiz; EM, eastern Mediterranean; MAR, mid-Atlantic Ridge; HMMV, Haakon Mosby mud volcano. Substrate: alfalfa (A) or wood (W). Deployment reference from I to V.
a Retch MM301 BeadBeater (Retch, Germany). Approximately equal volume of alfalfa samples were ground to powder in liquid nitrogen. A modified Doyle and Doyle (1987) DNA extraction protocol was applied. About 200–300 mg of powder was mixed with 1 mL of preheated CTAB isolation buffer containing 2% hexadecyltrimethylammonium bromide, 1.4 M NaCl, 0.2% β-mercaptoethanol, 20 mM EDTA, 100 mM Tris-HCl, 1% PVP, and 0.2 mg × mL-1 proteinase K. Samples were agitated (1 h, 60°C) then mixed with one volume of chloroform-isoamyl alcohol (24:1). Nucleic acids from the upper aqueous phase were precipitated by isopropanol and centrifuged. Dry DNA pellets were re-suspended in sterile water. Concentration and quality were checked. DNA extraction was performed in duplicate, using another randomly selected wood cube or alfalfa grass sample from the same CHEMECOLI, resulting in two replicates (region-substrate-nb-1 and region-substrate-nb-2).

Primers V5V6_F (CAAACAGGATTAGATACCCTG) and V5V6_R (CGTTRCGGGACTTAACCCAACA) (designed by GENOSCREEN) were used to amplify the V5-V6 hypervariable region corresponding to regions 770–1094 in the *E.coli* 16S rRNA-encoding gene. Short (~300 bp) pyrotags were sequenced by 454 pyrosequencing (GsFLX-Roche Diagnostics, GENOSCREEN, France). A 10-bp molecular identifier (MID) tag was inserted between the GS-FLX adapter and the specific primer to facilitate further sequence binning. PCR products were purified and quantified by Picogreen (Invitrogen, USA) and the same amount of amplicons was mixed prior to 454 pyrosequencing. 5/8 of a plate was used to sequence 52 samples (including replicates), expecting a minimum of 7700 reads per sample. Sequences have been submitted to MG-RAST (http://metagenomics.anl.gov/linkin.cgi?project=9917) with MG-RAST ID's ranging from 4571074.3 to 4571125.3 (Table S1).

## DATA ANALYSES

Resulting binary “.sff” files were extracted using Mothur (Schloss et al., 2009). Sequences shorter than 250 bp, longer than 350 bp and containing Ns were eliminated from further analysis. Filtered sequences were sorted by their MID sequences into separate Fasta files. Typical 454 sequencing errors and PCR single base errors were screened using the PyroNoise and SeqNoise modules of the AmpliconNoise software (Quince et al., 2011) with default parameters. Sequences with MIDs and primers removed were used to generate a Needleman-Wunsch distance matrix (function NDist) and clustered into operational taxonomic units (OTUs, function Fcluster) using the same software. The matrix of sequence abundances per OTU was generated using Python scripts using 97% identity threshold for OTU definition (Dataset_1 in Table S1). The 10 best hits (maximal e-value = 1 e⁻¹⁰) for each OTU were found using BLAST (Altschul et al., 1997) using the SSURef_NR99_115 Silva database (Quast et al., 2013), and added to the OTU abundance matrix.

## ANALYSES OF OTU ABUNDANCES AND SAMPLES COMPARISON

The microbial α-diversity indices (Shannon and Invsimpson), community microbial richness indices (Chao1 and ACE) and rarefaction data were computed with Mothur on Dataset_1 (Schloss et al., 2009). The second dataset was prepared by eliminating all

sequences present in less than 5 copies in the total number of reads per sample (rare OTUs eliminated) and was screened for symbiont-related sequences (Dataset_2 in Table S1). For further community abundance and statistical analyses, OTUs containing in total less than 40 sequences, which represents ~1% of the total number of sequences in the smallest sample (4543) were eliminated, thus generating a reduced data matrix with abundant sequences (Dataset_3 in Table S1). All statistical analyses were performed using R (R Development Core Team, 2013). Non-Metric Multidimensional Scaling (NMDS) was performed using “metaMDS” function in the MASS package. Average agglomerative clustering dendrograms were generated using Unweighted Pair-Group Method with arithmetic Average (UPGMA) on a Bray-Curtis distance matrix. Abundance data were transformed using the Hellinger method (Legendre and Gallagher, 2001) and then used for transformation-based redundancy analyses (tb-RDA). Substrate type and region were used as factors and water temperature, depth, and deployment time were used as variables into the constrained RDA, in order to estimate their contribution to the global variance. Significance was assessed using ANOVA permutation tests.

## IDENTIFICATION OF OTUs RELATED TO METAZOAN-ASSOCIATED BACTERIAL SYMBIONTS

The Dataset_2 was used for screening for symbiont-related sequences. Sequences from this dataset were blasted (with 97% identity threshold) against a local database (Supplementary File), containing sequences of bacteria living in symbiosis with deep-sea fauna (Potential_symbionts in Table S1). These sequences were then blasted against public database and their blast hits definitions were then filtered by scanning for the “symbiont” key word (Confirmed_symbionts in Table S1). In parallel, a phylogenetic analysis was computed in order to evaluate the phylogenetic relationship of sequences to known symbionts. For this, 16S rRNA sequences of documented symbiotic bacteria from vent and seep metazoans and several best blast hits from our OTUs were aligned with SINA Web Aligner (Pruesse et al., 2012) and truncated to the V5-V6 region. Symbiont-related OTUs from this study were added to this dataset, aligned with ClustalX (Larkin et al., 2007), and alignments were manually checked. Phylogenetic relationships among sequences were estimated from a 250-bp alignment with MEGA6 (Tamura et al., 2013) using distance methods and neighbor-joining. Bootstrap values were computed on 1000 replicates.

## RESULTS

### PREPARATION OF DATASETS

Bacterial diversity was investigated on 25 individual CHEMECOLIs filled with either pine wood cubes or alfalfa and deployed for periods of 10–1112 days in four areas in the North East Atlantic and Mediterranean. A total of 564925 V5-V6 reads was obtained. After length filtering, elimination of PCR and random sequencing errors, a raw dataset containing 364633 sequences distributed in 22721 OTUs (97% cut-off) was obtained. Screening for symbiont-related sequences was made on a reduced dataset including 332621 sequences representing 2641 OTUs. Community composition and comparisons were
performed on a dataset containing only abundant sequences, overall 306716 sequences in 658 OTUs. On average, dataset sizes were reduced by 35, 41, and 46% of total read numbers and contained 7012, 6397, and 5898 sequences per sample (2 replicate samples per CHEMECOLI). Details on sequence abundances in each sample and dataset, and the mean number of reads per sample are summarized in Table S1.

## TAXONOMIC RICHNESS AND DIVERSITY

Rarefaction curves based on a similarity threshold of 97% were generated for each substrate type separately (Figure S1). Generally, alfalfa samples showed more diverse (up to 1.6-fold) bacterial communities than their wood equivalents. This was confirmed by bacterial community richness and diversity indices (Table 2). Among alfalfa samples, all long-term ones (>289 days of deployment) from Gulf of Cadiz cold seeps (GoC), Eastern Mediterranean cold seep (EM) and Lucky Strike (L-S) hydrothermal vent on Mid Atlantic Ridge (MAR) (GoC-A-I to GoC-A-IV, EM-A-I and MAR-A-III) presented more diverse bacterial communities than other samples (Table 2). In samples mentioned above the Chao1 index was higher than 1500 and the Shannon index had a value higher than 4.5 compared to wood samples. On any given set of CHEMECOLIs, the bacterial community was less diverse in wood than alfalfa samples (Table 2). Within wood samples GoC-W-II, GoC-W-III and MAR-W-III had the highest values of all considered diversity indices (Table 2). Replicate samples for both alfalfa and wood displayed comparable values for all analyzed diversity indices (Table 2), with two exceptions, namely the replicates in a GoC alfalfa sample (GoC-A-II-1 and GoC-A-II-2), and a short-term wood replicates from MAR (MAR-W-IV-1 and MAR-W-IV-2) (Table 3). Although Shannon and Invsimpson diversity indices showed similar values for the two GoC-A-II replicates, replicates from MAR-W-IV displayed a 6-fold difference in their Invsimpson index (Table 2).

## TAXONOMIC COMPOSITION OF BACTERIAL COMMUNITIES

At the division or sub-division level, 10 taxa included altogether more than 96% of all sequences present in Dataset_3 (Figure 2; Table S2). Proteobacteria was the most abundant division in this dataset and alone represented above 85% of the total number of reads (Figure 2; Table S2). Within them three taxa (Gamma-, Delta-, and Alphaproteobacteria) corresponded to 77% of the total number of amplicons (Table S2). Gammaproteobacteria were remarkably over-represented in short-term samples (MAR-A/W-IV; EM-W-II and MAR-A-V, Figure 2; Table S2). Epsilonproteobacteria were also highly abundant in short-term samples (MAR-A/W-IV and MAR-A/W-V), while not exceeding 2% of total sequence numbers in long term samples. They were particularly dominant (98%) in the short-term wood sample from M-G on MAR. While almost absent from short-term samples, Delta- and Alphaproteobacteria represented a significant fraction of reads in all long-term samples. Other taxa represented below 30% of reads per sample.

For further community analyses at the OTU level, we focused on the most abundant OTUs, which represented above 1% of the total number of sequences in the dataset (>3000

amplicons; Table S3). 15 OTUs within the Proteobacteria and 1 within the Bacteroidetes matched this criterion, representing 131712 sequences (almost 43% of the total sequences) (Figure 3; Table S3). All were related to marine bacteria. Four major OTUs (OTU_00201, OTU_11277, OTU_00411 and OTU_00023) were abundant in almost all long-term samples and constituted altogether 38 to 86% of the total number of sequences in each sample (Table S3). Communities from short-term experiments were different from the long-term ones, comprising 6 highly abundant OTUs that were rare to absent in long-term deployments (Figure 3; Table S3).

## COMPARISON AMONG BACTERIAL COMMUNITIES

Out of the 24 paired replicates of samples, 21 were grouped together in the UPGMA-based cluster dendrogram based on the Bray-Curtis dissimilarity matrix (Figure 4). The exceptions were two GoC samples (GoC-W-I and GoC-W-IV) and one short-term MAR sample (MAR-W-IV), for which replicates were not the closest relatives. This indicates congruency within most replicates, except for three out of 24 samples. Samples that had been fixed following two different procedures (ethanol fixation for MAR-A/W-V-1/2 and deep freezing for MAR-A/W-V-C) were also on the same branches, indicating that the type of fixation had limited influence on the results. Some groups could be easily distinguished. Both substrates for long-term EM samples (EM-A/W-I) clustered together. Alfalfa substrates from GoC (GoC-A-I to GoC-A-IV) and L-S (MAR-A-III) grouped together; as did wood substrates from GoC (GoC-W-I to GoC-W-IV). A larger group included samples corresponding to both substrates from MAR and HMMV samples (MAR-A/W-I to MAR-A/W-III, HMMV-A/W-I, and HMMV-A/W-II). Finally, all short-term samples formed a clearly distant group (MAR-A/W-IV, MAR-A/W-V, and EM-A/W-II). Interestingly, CHEMECOLIS deployed for different periods of time in the long-term experiment clustered within the same larger groups but do not with each other. These include the GoC Mercator deployments for 289 (GoC-A/W-I) and 731 days (GoC-A/W-II), the MAR Rainbow deployments for 328 and 414 days (MAR-A/W-I and MAR-A/W-II), and the HMMV deployments for 388 and 752 days (HMMV-A/W-I and HMMV-A/W-II). Similar groupings were found on the NMDs graph (Figure S2).

The influence of factors (substrate type, deployment area) and variables (water temperature, duration of deployment and depth) on the community composition was evaluated by constrained RDA (Figure 5 and Figure S4). On the two dimensional graph, clear segregation of samples by region was visible, samples clustering around the centroids representing their region of origin. The influence of the type of substrate was less marked, as centroids representing pine wood and alfalfa were not very distant from one another (Figure 5 and Figure S4). Variables were represented as three vectors with their lengths proportional to the variance explained. Samples from MAR were for example aligned parallel to the vector "DAYS," according to the duration of deployment: short-term deployments were in the upper left corner of the graph (MAR-A/W-IV and MAR-A/W-V), 1-year in the center left (MAR-A/W-I and MAR-A/W-II), and longest deployments (MAR-A/W-III, 1112 days) in the center bottom
<table><caption>Table 2 | Species richness and diversity indices calculated for each sample in Mothur (Schloss et al., 2009) for sequences with OTU identity threshold at 97%.</caption><thead><tr><th colspan="2">Threshold</th><th colspan="4">Index</th></tr><tr><th colspan="2">0.03</th><th colspan="2">Richness</th><th colspan="2">Diversity</th></tr><tr><th>Sample</th><th>Sobs</th><th>Chao1</th><th>ACE</th><th>Shannon</th><th>Invsimpson</th></tr><tr><th>nb_seq</th><th></th><th colspan="4">4543</th></tr></thead><tbody><tr><td>Contr-A</td><td>279.11</td><td>526.96</td><td>733.44</td><td>3.39</td><td>11.18</td></tr><tr><td>GoC-A-I-1</td><td>1196.32</td><td>2784.75</td><td>4440.27</td><td>5.64</td><td>52.92</td></tr><tr><td>GoC-A-I-2</td><td>899.13</td><td>2002.70</td><td>3446.10</td><td>4.99</td><td>28.71</td></tr><tr><td>GoC-A-II-1</td><td>589.27</td><td>851.89</td><td>912.04</td><td>5.22</td><td>85.76</td></tr><tr><td>GoC-A-II-2</td><td>1678.54</td><td>4108.65</td><td>7117.76</td><td>6.39</td><td>158.41</td></tr><tr><td>GoC-A-III-1</td><td>740.86</td><td>1558.90</td><td>2554.81</td><td>4.90</td><td>40.70</td></tr><tr><td>GoC-A-III-2</td><td>1108.94</td><td>2297.72</td><td>3608.17</td><td>5.30</td><td>32.70</td></tr><tr><td>GoC-A-IV-1</td><td>1265.55</td><td>2957.10</td><td>5319.83</td><td>5.71</td><td>65.85</td></tr><tr><td>GoC-A-IV-2</td><td>1396.16</td><td>3102.94</td><td>4999.33</td><td>5.90</td><td>73.61</td></tr><tr><td>EM-A-I-1</td><td>916.56</td><td>1984.98</td><td>3157.16</td><td>5.23</td><td>45.95</td></tr><tr><td>EM-A-I-2</td><td>1058.51</td><td>2370.38</td><td>3736.46</td><td>5.22</td><td>35.52</td></tr><tr><td>MAR-A-I-1</td><td>154.85</td><td>207.59</td><td>234.98</td><td>3.02</td><td>9.74</td></tr><tr><td>MAR-A-I-2</td><td>216.49</td><td>421.47</td><td>617.40</td><td>2.69</td><td>6.24</td></tr><tr><td>MAR-A-II-1</td><td>644.48</td><td>1238.09</td><td>1883.15</td><td>4.42</td><td>21.70</td></tr><tr><td>MAR-A-II-2</td><td>460.93</td><td>969.89</td><td>1496.23</td><td>4.02</td><td>18.40</td></tr><tr><td>MAR-A-III-1</td><td>1041.82</td><td>2283.52</td><td>3766.65</td><td>5.20</td><td>29.86</td></tr><tr><td>MAR-A-III-2</td><td>858.28</td><td>2079.23</td><td>3740.74</td><td>4.41</td><td>11.43</td></tr><tr><td>HMMV-A-I-1</td><td>598.65</td><td>1149.13</td><td>1601.53</td><td>4.37</td><td>21.37</td></tr><tr><td>HMMV-A-I-2</td><td>499.62</td><td>862.40</td><td>1135.98</td><td>4.15</td><td>18.58</td></tr><tr><td>HMMV-A-II-1</td><td>625.12</td><td>1315.23</td><td>2082.32</td><td>4.49</td><td>25.48</td></tr><tr><td>HMMV-A-II-2</td><td>673.01</td><td>1283.04</td><td>1924.22</td><td>4.62</td><td>31.34</td></tr><tr><td>MAR-A-IV-1</td><td>46.91</td><td>126.08</td><td>192.93</td><td>0.81</td><td>1.42</td></tr><tr><td>MAR-A-IV-2</td><td>46.19</td><td>102.66</td><td>174.43</td><td>0.98</td><td>1.60</td></tr><tr><td>MAR-A-V-1</td><td>276.45</td><td>604.08</td><td>830.39</td><td>3.43</td><td>11.72</td></tr><tr><td>MAR-A-V-2</td><td>279.00</td><td>508.53</td><td>852.31</td><td>3.43</td><td>11.92</td></tr><tr><td>MAR-A-V-C</td><td>237.16</td><td>503.26</td><td>763.51</td><td>3.15</td><td>9.49</td></tr><tr><td>GoC-W-I-1</td><td>282.77</td><td>433.37</td><td>535.60</td><td>4.11</td><td>27.14</td></tr><tr><td>GoC-W-I-2</td><td>379.09</td><td>656.86</td><td>900.68</td><td>3.84</td><td>12.01</td></tr><tr><td>GoC-W-II-1</td><td>782.17</td><td>1722.28</td><td>2878.57</td><td>4.61</td><td>19.39</td></tr><tr><td>GoC-W-II-2</td><td>709.60</td><td>1440.33</td><td>2247.49</td><td>4.39</td><td>18.23</td></tr><tr><td>GoC-W-III-1</td><td>529.85</td><td>1069.75</td><td>1658.48</td><td>4.53</td><td>30.35</td></tr><tr><td>GoC-W-III-2</td><td>536.81</td><td>1123.29</td><td>1637.86</td><td>4.71</td><td>39.90</td></tr><tr><td>GoC-W-IV-1</td><td>391.71</td><td>714.70</td><td>1096.37</td><td>4.48</td><td>42.43</td></tr><tr><td>GoC-W-IV-2</td><td>282.05</td><td>460.62</td><td>510.23</td><td>3.67</td><td>12.68</td></tr><tr><td>EM-W-I-1</td><td>367.52</td><td>636.47</td><td>824.42</td><td>4.40</td><td>39.49</td></tr><tr><td>EM-W-I-2</td><td>408.57</td><td>741.91</td><td>1013.12</td><td>4.48</td><td>36.14</td></tr><tr><td>MAR-W-I-1</td><td>178.22</td><td>308.73</td><td>435.47</td><td>2.75</td><td>5.51</td></tr><tr><td>MAR-W-I-2</td><td>97.21</td><td>152.68</td><td>168.77</td><td>2.28</td><td>4.06</td></tr><tr><td>MAR-W-II-1</td><td>178.81</td><td>264.22</td><td>308.09</td><td>3.62</td><td>20.69</td></tr><tr><td>MAR-W-II-2</td><td>257.03</td><td>439.33</td><td>542.02</td><td>3.26</td><td>8.00</td></tr><tr><td>MAR-W-III-1</td><td>691.66</td><td>1391.93</td><td>2370.47</td><td>4.60</td><td>28.79</td></tr><tr><td>MAR-W-III-2</td><td>580.78</td><td>1105.29</td><td>1658.46</td><td>4.43</td><td>25.42</td></tr><tr><td>HMMV-W-I-1</td><td>291.77</td><td>568.01</td><td>833.71</td><td>3.54</td><td>13.81</td></tr><tr><td>HMMV-W-I-2</td><td>407.70</td><td>794.03</td><td>1054.50</td><td>4.12</td><td>20.23</td></tr><tr><td>HMMV-W-II-1</td><td>347.54</td><td>625.72</td><td>947.47</td><td>3.59</td><td>10.12</td></tr><tr><td>HMMV-W-II-2</td><td>434.71</td><td>795.56</td><td>1031.26</td><td>3.99</td><td>13.96</td></tr><tr><td>MAR-W-IV-1</td><td>77.80</td><td>106.13</td><td>109.30</td><td>1.48</td><td>2.54</td></tr><tr><td>MAR-W-IV-2</td><td>354.33</td><td>687.95</td><td>995.05</td><td>3.65</td><td>14.37</td></tr><tr><td>EM-W-II-1</td><td>154.90</td><td>266.93</td><td>321.80</td><td>2.97</td><td>10.54</td></tr><tr><td>MAR-W-V-1</td><td>139.41</td><td>291.52</td><td>569.19</td><td>2.60</td><td>7.75</td></tr><tr><td>MAR-W-V-2</td><td>116.56</td><td>238.70</td><td>383.18</td><td>2.43</td><td>6.92</td></tr><tr><td>MAR-W-V-C</td><td>201.96</td><td>442.01</td><td>695.70</td><td>2.67</td><td>6.64</td></tr></tbody></table>

Sample IDs are described in Table 1.
(Figure S4). The RDA overall explained 38% of the total variance. The deployment region accounted for 21% of the variance; days, substrate, and depth between 5 to 4%; and the water temperature was responsible for less than 4% of the total variance (Table S4). The ANOVA test confirmed that all factors and variables were statistically significant, based on permutation tests (p < 0.001, 99 permutations; Table S4).

IDENTIFICATION AND DISTRIBUTION OF BACTERIA RELATED TO CHEMOSYNTHETIC SYMBIONTS

A total of 30607 reads within 34 OTUs (~9.2% of the Dataset_2) have been identified as potential symbionts (>97% of identity with sequences in local database of symbionts) while 21612 reads, distributed in 16 OTUs, had the “symbiont” keyword in at least one of their ten best blast hit definitions (~6.5

<table><caption>Table 3 | Symbionts and animal hosts occurrence in geographical regions analyzed in this study.</caption><thead><tr><th rowspan="2">Host family</th><th colspan="5">Region</th><th rowspan="2">Total nb seqs</th><th colspan="4">Symbiont occurrence</th><th colspan="4">Host occurrence</th><th rowspan="2">References</th></tr><tr><th>C-</th><th>GoC</th><th>EM</th><th>MAR</th><th>HMMV</th><th>GoC</th><th>EM</th><th>MAR</th><th>HMMV</th><th>GoC</th><th>EM</th><th>MAR</th><th>HMMV</th></tr></thead><tbody><tr><td>Alvinocarididae</td><td>0</td><td>11</td><td>13</td><td>9669</td><td>0</td><td>9693</td><td>±</td><td>±</td><td>+</td><td>–</td><td>–</td><td>–</td><td>+</td><td>–</td><td>WoRMS Website</td></tr><tr><td>Mytilidae</td><td>0</td><td>2778</td><td>174</td><td>4766</td><td>1728</td><td>9446</td><td>+</td><td>+</td><td>+</td><td>+</td><td>+</td><td>+</td><td>+</td><td>+</td><td>Duperron et al., 2013</td></tr><tr><td>Teredinidae</td><td>0</td><td>1378</td><td>14</td><td>0</td><td>0</td><td>1392</td><td>+</td><td>±</td><td>–</td><td>–</td><td>–</td><td>–</td><td>–</td><td>–</td><td>Borges et al., 2014</td></tr><tr><td>Tubificidae</td><td>0</td><td>542</td><td>34</td><td>10</td><td>3</td><td>589</td><td>+</td><td>+</td><td>±</td><td>±</td><td>–</td><td>–</td><td>–</td><td>–</td><td>Ruehland and Dubilier, 2010</td></tr><tr><td>Siboglinidae</td><td>0</td><td>393</td><td>2</td><td>55</td><td>37</td><td>487</td><td>+</td><td>±</td><td>+</td><td>+</td><td>+</td><td>+</td><td>–</td><td>+</td><td>The World Polychaeta Database Website</td></tr><tr><td>Thyasiridae</td><td>0</td><td>5</td><td>0</td><td>0</td><td>0</td><td>5</td><td>±</td><td>–</td><td>–</td><td>–</td><td>+</td><td>+</td><td>–</td><td>–</td><td>Rodrigues and Duperron, 2011</td></tr><tr><td>Total nb seqs</td><td>0</td><td>5107</td><td>237</td><td>14500</td><td>1768</td><td>21612</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr></tbody></table>

Symbiont occurrence is a synthesis of data from Table S5, host occurrence is based on available literature study. Sample IDs are described in Table 1.

<figure><img src="image_5.png" /></figure>
<figure><img src="image_6.png" /><figcaption>**FIGURE 3** | Percentage of pyrosequencing reads representing the most abundant OTUs (more than 1% of the total number of sequences, i.e., >3000 reads) and assigned to different bacterial phylotypes. Taxonomic</figcaption></figure>

classification was based on comparison of sequences within the SILVA database with OTU identity threshold at 97%. More details on taxonomic assignment of OTUs analyzed is provided in **Table S3**. Refer to **Table 1** for sample ID.

% of the Dataset_2). Putatively related symbionts were distributed within 9 potential metazoan host families. Above 22000 belonged to relatives of symbionts associated with *Rimicaris* shrimps, *Idas* mussels and *Lyrodus* shipworms. A phylogenetic tree was computed using representative sequences of potential symbiont OTU, one best blast hit from the local database of symbionts and four best blast hits from public databases (**Figure S3**). Only the OTUs displaying the keyword “symbiont” within at least one out of ten best blast hits (public database) and above 97% sequence identity with a symbiont sequence were retained as potential candidates for free-living forms of symbionts (**Table S5**). We obtained 21612 sequences representing 16 OTUs related to symbionts (**Table S5**). OTUs related to *Rimicaris*-associated Epsilonproteobacteria (97.7–99% identity) represented 9693 reads, followed by *Idas*-associated Symbiont G and Bacteroidetes (6888 and 2558 reads with above 97.4 and 99% identities, respectively), and *Lyrodus* symbionts-related OTUs (1392, above 97.4% identity). Symbionts of *Tubificoides benedii* (589 sequences with 99.7% of identity), relatives of bacteria associated with Siboglinidae polychaetes (487, 97.6–99.7% identity) and Thyasiridae bivalves (5 sequences only with 97.9% identity) were far less abundant (**Table S5**).

The *Rimicaris* symbiont-related OTUs were highly abundant in short-term deployments in MAR, very rare in EM and GoC, and absent everywhere else. Relatives of *Idas* Bacteroidetes and Symbiont G, Siboglinidae symbionts and those of *Tubificidae*

were ubiquitous (GoC, EM, MAR, HMMV). Those of *Lyrodus* were present only in GoC and EM, and those of Thyasiridae bivalves only in GoC. The occurrence of symbiont-related OTUs was related to the host presence in the corresponding regions (**Table 3**). Shrimps and their symbionts co-occured on the MAR (WoRMS Website<fnref n="1" />,) with only few symbiont related sequences in GoC and EM. Thyasirids (*Rodrigues and Duperron, 2011*) and their symbionts co-occurred in GoC, while symbionts were absent in MAR, HMMV and EM. Symbionts of *Idas* mussels and siboglinids were present in all regions, while the hosts were absent at HMMV (*Duperron et al., 2013*) and MAR (World Polychaeta Database Website<fnref n="2" />) respectively. Symbionts of *Tubificoides benedii* were also ubiquitous while those of teredinids could be detected mainly in GoC and EM. Metazoans of these two families occur in shallow water.

## DISCUSSION

### METHODOLOGY

All CHEMECOLI samples were analyzed in duplicates. Estimated diversity and richness indices of most duplicates were comparable, save for two exceptions in which Chao1 and ACE indices varied 5–9-fold between replicates for a long-term sample from

<footnote n="1">¹WORMS, World Register of Marine Species. http://www.marinespecies.org/index.php</footnote>

<footnote n="2">²The World Polychaeta Database http://www.marinespecies.org/polychaeta/</footnote>
<figure><img src="image_7.png" /><figcaption>**FIGURE 4** | Dendrogram for average agglomerative clustering generated in R (R Development Core Team, 2013) using Unweighted Pair-Group Method with arithmetic Average (UPGMA) on a Bray-Curtis distance</figcaption></figure>

matrix ("hclust()" function in stats package). The dendrogram is rooted on the negative alfalfa control (Contr-A-I). Sample IDs are described in **Table 1**. Color code is the same as that used in **Figure 1**.

GoC and for short-term one from MAR. These two exceptions emphasize the need for sample replication. The UPGMA dendrograms evidenced the very short distances between most replicated samples. Similar groupings were observed on the RDA (Figure 5) graph, in which replicates were most often in very close vicinity to each other. Finally, high similarity was observed between samples stored in 96% ethanol and deep frozen at -80°C. Ethanol treatment is thus an appropriate method for sample storage and conservation in view of the difficult conditions experienced sometimes on board oceanographic cruises and during shipment.

## BACTERIAL COLONIZATION OF PINE WOOD AND ALFALFA SUBSTRATES IN DEEP-SEA CHEMOSYNTHESIS-BASED ECOSYSTEMS

A variety of devices have been designed to study bacterial colonization in deep-sea ecosystems such as hydrothermal vents, in most cases using relatively inert surfaces (examples in Reysenbach et al., 2000; López-García et al., 2003; Rassa et al., 2009). Wood substrates have been introduced more recently (Fagervold

et al., 2012, 2013; Bienhold et al., 2013). In our study, bacterial communities colonizing alfalfa were usually more diverse than those colonizing wood samples (Figure S1; Table 2). This is not unexpected because wood displays high content of lignin, which consists in refractory hydrophobic polymers harder to degrade than cellulose, and because alfalfa offers greater surface per volume for bacterial colonization. The overall OTU diversity in short-term deployments was lower than in samples deployed for more than 289 days. After 289 days, no correlation was observed between community richness and deployment duration (data not shown), suggesting that diversity may have reached saturation sometime in between 14 and 289 days. Short-term samples were usually recovered during distinct legs of the same cruise, while long-term samples were recovered during follow-up cruises, which occurred on a yearly basis, explaining the lack of intermediate deployment times. These would however be crucial to understand how fast saturation is reached.

The bacterial community in long-term samples included only 10 bacterial subdivisions (Figure 2; Table S2). Three
<figure><img src="image_8.png" /></figure>

**FIGURE 5** | RDA biplot of the Hellinger-transformed OTU abundance data constrained by all environmental variables, scaling 3. Substrate type [factor(SUBS)] and region [factor(REG)] were used as factors and water temperature (TEMP), depth (DEPT), and deployment time (DAYS) were used as variables into the constrained RDA, generated in R using "rda()" and "plot()" functions in "vegan" package. Geographic areas are

indicated by ellipses with color corresponding to region of origin. Color code follows that of **Figures 1, 4**. Sample IDs of individual points can be found in **Figure S4** and environmental variables and factors are described in **Table 1**. The influence of each variable and factor and their significance were assessed using ANOVA permutation tests—details in **Table S4**.

of them (Alpha-, Delta-, and Gammaproteobacteria) were over-represented if compared to others, as shown before (Fagervold et al., 2012, 2013; Bienhold et al., 2013). Among the most dominant OTUs (Figure 3; Table S3) were members of the Oceanospirillales and Alteromonadales, facultative anaerobe growing on cellobiose and previously found on sunken wood (Fagervold et al., 2012). Within the Deltaproteobacteria, OTUs related to Mycococcales and Desulfobacteriales, and two groups of SRBs detected by Bienhold et al. (2013) (OTU refs in Bienhold et al., 2013: Deltaproteobacteria_03, _24, _27, _50, _55, _283, _421, and _737), were identified. Finally, Bacteroidetes were found, which are often chemoorganotrophs specialized in degrading various biopolymers such as cellulose (Kirchman, 2002). The composition of short-term communities was different, and except for the wood sample from M-G, was dominated by Gammaproteobacteria. This is comparable with the 1-day deployment control sample in Bienhold's study (2013). Epsilonproteobacteria were over-represented in the short-term M-G wood sample thanks to the high abundance of *Rimicaris exoculata* ectosymbiont (Table S3), which is discussed below. The abundance of Epsilonproteobacteria is consistent with their dominance at seep and vent ecosystems and their documented quick colonization of newly formed habitats including wood

(Reysenbach et al., 2000; López-García et al., 2003; Huber et al., 2010; Fagervold et al., 2012, 2013; Bienhold et al., 2013).

## MAJOR FACTORS INFLUENCING THE MICROBIAL COMMUNITY COMPOSITION

The RDA explains less than 40% of the total variance observed in the community composition. It means that other factors and variables, likely including environmental parameters not measured during sampling, play a major role. These are hard to obtain because only time-series acquired in the immediate vicinity of CHEMECOLIs and over the course of the experiment would make sense. Nevertheless, physico-chemical data could explain another significant fraction of the total variance (Lee et al., 2014). Among factors considered here, the deployment area contributed 21% of the variance. This is evident from the UPGMA, NMDS and RDA graphs which cluster samples according to their area of origin, with few exceptions (i.e., L-S sample MAR-A/W-III). EM samples are grouped and GoC alfalfa and wood samples build two rather distinct groups. Interestingly, MAR and HMMV samples appear more mixed. This suggests that communities from vent MAR and cold seep HMMV deployments are less distinct, despite the great geographical distance and very different latitudes. This observation reminds of a previous study in which deep
ocean microbial communities from poles and low latitudes dif-
fered less than corresponding samples from the euphotic surface
waters (Ghiglione et al., 2012). Other factors and variables each
contribute 5% or less to the variance. The fact that substrate does
not explain much variance is probably due to the similar com-
position of alfalfa and wood, both being plant-derived material.
The consequence is that although alfalfa communities are more
diverse, the most abundant OTUs tend to be shared. Alfalfa may
host a greater number of rare OTUs. More caution is needed when
looking at other variables. Duration of deployment for example
only explains 5% of the total variance, yet short-term samples
are remarkably less diverse and cluster very far from long-term
samples in the various analyses. This apparent lack of congruency
between the results of RDA and comparisons is clearly due to a
sampling bias. Short-term samples (10–13 days) represented only
5 of the 25 deployed CHEMECOLIs, the other being deployed
for 289–1112 days. Their quantitative contribution to the vari-
ance in the overall datset is thus low, and only if we had as
many short terms as long terms samples would we have a correct
estimate of the percentage of variance explained by deployment
duration. Long-term CHEMECOLIS deployed at the same site
for different durations do not display the most similar commu-
nities, but their communities do not cluster very far from one
another, suggesting a certain level of stability at this stage. If we
had better coverage of intermediate times and a more balanced
design, duration of deployment could become a major explana-
tory factor (Palacios et al., 2009; Fagervold et al., 2012). Although
the temperature variable is well equilibrated with 7 sites display-
ing typical deep-sea temperatures (2–4.6°C) against six others
ranging from 10 to 13°C, it explains less than 4% of the total vari-
ance. Depths cover a broad bathymetric range (354 to 2300 m)
and samples are distributed rather homogeneously, yet it explains
only 4% of the variance. Depth is usually not considered a major
factor in the aphotic zone of the ocean (Fagervold et al., 2013;
Zhu et al., 2013). Neither latitude nor water temperature showed
significant correlations with estimated microbial diversities in the
previous large scale study by Ghiglione et al. (2012). Overall, the
geographical area of deployment seems to be the most impor-
tant factor influencing the microbial colonization in this study.
Short term samples are different from long terms in both their
overall diversity as well as community composition. Both sub-
strates share similar dominant bacteria and alfalfa displays greater
diversity thanks to rare OTUs. Communities become saturated
in term of both diversity and composition between 14 and 289
days after deployment. For longer periods, bacterial composition
appears to be relatively stable, as samples from the same sites
group together. Only marine microbes colonize both substrates
in CHEMECOLIs, no typical contaminants could be detected.
The bacterial taxonomic composition identified in colonization
devices compared to other studies (Fagervold et al., 2012, 2013;
Bienhold et al., 2013) confirms that the degradation of organic
matter starts to take place within the first 2 weeks after immersion.

COLONIZATION OF SUBSTRATES BY SYMBIONT-RELATED BACTERIA

The sequence dataset was screened for sequences related to
chemosynthetic symbionts of metazoans. Based on percent-
age similarity (97%) criterion with documented symbionts, we

retained sequences closely related to documented symbiotic bac-
teria as reasonably supported candidates. Despite that the short
length of reads hampers the quality of phylogenetic recon-
structions, the combination of sequence similarity and phyloge-
netic relatedness appears a rather conservative approach. Almost
70% of symbiont-related sequences were found on pine wood
(Table S1). The higher diversity within alfalfa samples may indeed
"dilute" symbiont-related sequences in higher numbers of other
bacterial OTUs, making their discovery less likely. Putative free-
living forms of symbionts or symbiont-related bacteria from
five host metazoan families were identified. The most abundants
were related and 99% identical to Epsilonproteobacteria ectosym-
bionts of *Rimicaris exculata*. In this species, bacteria quickly
re-colonize the shrimp's gill chamber after each molt, estimated
to occur every 10 days (Zbinden et al., 2004; Corbari et al., 2008).
Our results are in good agreement with this constraint. Shrimp
symbiont-related OTUs were indeed the most abundant in the
short-term MAR samples (MAR-A/W-IV and MAR-A/W-V, 10-
14 days), and absent from all long-term deployments. They thus
seem to be pioneer and efficient colonizers of available surfaces at
vents. This questions the specificity of the symbiotic association.
It may well be that bacteria behave as opportunists colonizing the
new cuticle when it appears, as they would colonize any surface.
The molt cycle would favor iterative colonization by the same
bacteria simply by renewing the surface available every few weeks.

Relatives of gill endosymbionts of the mytilid *Idas* sp. were also
abundant, mostly the Symbiont-G group and the Bacteroidetes,
two symbiont groups which were recently discovered in species
from the eastern Mediterranean and the Gulf of Cadiz (Duperron
et al., 2008; Rodrigues et al., 2013). The usual sulfur-oxidizers are
much rarer, though a few were found. Finding free-living relatives
of mussel symbionts was not unexpected, given that environmen-
tal acquisition of symbionts is supported in *Idas modiolaeformis*
based on the absence of bacteria within gonads and gametes
(Gaudron et al., 2012). In a recent study, Laming et al. (2014)
showed that environmental symbiont acquisition took place at the
early juvenile stage, after settlement.

The occurrence of *Lyrodus* and *Osedax* symbionts-related OTUs with identities between 97.4 and almost 100% was more unexpected in our samples, as shipworms live in shallow coastal waters (Borges et al., 2014) and *Osedax* is a bone-eating specialist (Rouse et al., 2004; Glover et al., 2005). *Lyrodus* gill endosymbionts have been proposed to produce cellulolytic enzymes that contribute to the host's ability to digest wood (Waterbury et al., 1983), and multiple bacterial phylotypes can co-occur in its gill bacteriocytes (Distel et al., 2002). Although the *Lyrodus pedicellatus* symbiont transmission mode remains still unknown, vertical transmission has been evidenced for its relative Bankia setacea (Sipe et al., 2000). But others wood-eating bivalves (*Xylophaga* spp. and *Xyloredo* spp.) were abundant in most samples and although we carefully avoided processing any animal tissue in our analyses, we cannot rule out that some bacterial symbionts were released during processing. Knowing that Bienhold et al. (2013) suggested a close phylogenetic relationship between symbionts of the shallow and deep-sea wood-boring bivalves we cannot objectively conclude whether the identified OTUs are actually free-living forms of *Lyrodus pedicellatus* or *Xylophaga* spp. symbionts.
Presence of an OTU related to *Tubificoides benedii* symbionts is also surprising as these oligochaete worms are found in eutrophic coastal sediments, but at the same time their ectosymbionts have been shown to belong to clades that consist almost exclusively of bacteria associated with invertebrates from deep-sea hydrothermal vents (Ruehland and Dubilier, 2010). The co-occurrence of thyasirid bivalves and their symbionts, represented in our study by only 5 sequences in GoC, may be explained by the mode of life of these animals burrowing deeply in sediments (Dufour and Felbeck, 2003).

## SYMBIONT TRANSMISSION AND COLONIZATION

All symbiont-related OTUs corresponded to environmentally-transmitted bacteria. The existence of free-living forms is thus expected, although their natural habitat is not documented (Gros et al., 1996; Harmer et al., 2008). Symbiont-related sequences were quite frequent in a few samples, and some were present on several sites. The abundance of symbionts on short-term wood samples indicates that they are probably efficient colonizers, although little is known about the mechanism of symbiont acquisition from the environment. The *Idas* sp. symbiont-G related-OTUs were on the other hand found in moderate numbers, but in almost all samples. Most of symbionts detected shared the same distribution pattern as their hosts (Table 3). An exception was the shipworms' symbiont in GoC and in EM, while their host occurs only in coastal shallow waters (Borges et al., 2014), but this may be because of close relatedness with symbionts of *Xylophaga* spp. which are present. We also did not expect the occurrence of Siboglinidae polychaetes symbionts on MAR, nor those of *Idas* sp. in HMMV. The detection of unexpected symbiont-related phylotypes (Table 3; Table S5) fuels the debate on the connectivity between different deep-sea sites (Vrijenhoek, 2010). Moreover, it confirms that short- and long-term colonization devices and high-throughput sequencing are good tools to further investigate the habitat and densities of free-living relatives of the symbionts.

The next step will be to distinguish between free-living close relatives of symbionts and true free-living forms of symbionts, but this will be tricky. Symbiont-specific fluorescence in situ hybridization (FISH) probes may not be able to properly discriminate, and low abundances of free living relatives may render difficult their detection and quantification on substrata, which show a high level of background autofluorescence (data not shown). Finally, the analyzed hypervariable region of 16S rRNA seems to be too short to be a good target for highly specific probes. Standardized devices are currently being deployed in other geographical regions and in other habitats. There is a high potential to improve the protocol we applied in our study with the increasing length of pyrosequencing reads, thus making the identification of microbes more reliable. For better understanding of the influence of environment on the microbial communities, the chemistry of water should be monitored precisely in the proximity of our colonization devices.

## ACKNOWLEDGMENTS

We thank captains, crew and chief scientists of RV's, Pourquoi pas?, Meteor, Atalante, Polarstern, James Cook, Pelagia and Belgica and teams operating ROVs Victor 6000 (Ifremer, France),

Quest 4000 (MARUM, Bremen, Germany) and submersible Nautile (Ifremer, France) and also people on board from all cruises participating in the deployment and recovery of CHEMECOLIs. We thank Françoise Gaill who initiated the CHEMECO EuroDEEP ESF program. This research was supported by UPMC, ANR DeepOases, CHEMECO EuroDEEP ESF, GDR DIWOOD, ITN Symbiomics and in Portugal by European funds (COMPETE) and by national funds through the Portuguese Science Foundation (FCT-EURODEEP/0001/2007 and FCT-PEst-C/MAR/LA0017/2013). Kamil Szafranski was funded through a Ph.D. grant from the Marie Curie Actions Initial Training Network (ITN) SYMBIOMICS (contract number 264774).

## SUPPLEMENTARY MATERIAL

The Supplementary Material for this article can be found online at: http://www.frontiersin.org/journal/10.3389/fmicb.2015.00162/abstract

**Figure S1 | Rarefaction analysis for the alfalfa grass (A) and pine wood (W) samples.** The curves were generated for 97% levels of OTU using Mothur (Schloss et al., 2009). Sample IDs are described in **Table 1**.

**Figure S2 | NMDS biplot of a Bray-Curtis dissimilarity matrix of the OTU abundance data generated in R (R Development Core Team, 2013) using "metaMDS()" and "plot()" functions in MASS package.** Sample data has been added using weighted averages. The short-term samples and the negative control are far from the long-term samples (green rectangle), the latter forming two cloud-like groups. Along the horizontal axis (NMDS1), long-term samples are roughly separated by substrate type, with alfalfa on the left and wood on the right (blue ellipse). The vertical axis displays a clear separation between samples from GoC and EM (red ellipse) and samples from MAR and HMMV (upper part), with the alfalfa sample from L-S (MAR-A-III) in between. Sample IDs are described in **Table 1**.

**Figure S3 | Phylogenetic tree based on sequences of V5-V6 hypervariable region of 16S rDNA gene obtained from the "potential symbiont" dataset and on respective sequences of known symbionts in public database.** Phylogenetic relationships among sequences were estimated with MEGA6 (Tamura et al., 2013) from a 250-bp ClustalX alignment (Larkin et al., 2007) using distance methods and neighbor-joining. Bootstrap values were computed on 1000 replicates. Each OTU is represented by sample ID followed by representative OTU number and total number of sequences in the respective OTU. Reference sequences from public database are represented by the accession number followed by the strain name. Based on the tree, OTUs related to environmental sequences were considered as doubtful candidates were not used for further analysis.

**Figure S4 | RDA biplot of the Hellinger-transformed OTU abundance data constrained by all environmental variables, scaling 3.** The figure is the same as **Figure 5** but samples IDs are detailed.

**Table S1 | Detailed distribution of sequencing data in datasets used in this study.** Refer to **Table 1** for sample IDs.

**Table S2 | The total number of V5-V6 sequences at OTU identity threshold of 97% in each bacterial division or sub-division and relative percentage in total effective bacterial sequences per sample.** Refer to **Table 1** for sample ID.
**Table S3 |** The total number of V5-V6 sequences at OTU identity threshold of 97% in the most abundant OTUs, represented by more than 1% of the total number of sequences in the dataset (>3000 amplicons). OTUs were assigned to different bacterial phylotypes by Blast on SILVA 115 database with OTU identity threshold at 97%. The relative percentage was normalized by the total number of the most abundant bacterial sequences per site. Refer to **Table 1** for sample ID.

**Table S4 |** Estimation of the influence of each variable and factor and their respective significance in redundancy analysis (RDA), assessed using ANOVA permutation tests in R (R Development Core Team, 2013). Sample IDs are described in **Table 1**.

**Table S5 |** The total number of sequences of symbiont-related OTUs identified by the "symbio" key word in at least one of their four best BLAST hit definitions or high similarity with a sequence from the symbiont database and confirmed by the phylogenetic proximity. Only OTUs displaying a symbiont sequence as their closest relative on the phylogenetic tree were retained. Sample IDs are described in **Table 1**.

## REFERENCES

Aida, M., Kanemori, M., Kubota, N., Matada, M., Sasayama, Y., and Fukumori, Y. (2008). Distribution and population of free-living cells related to endosymbiont A Harbored in *Oligobrachia mashikoi* (a Siboglinid Polychaete) Inhabiting Tsukumo Bay. *Microbes. Environ.* 23, 81–88. doi: 10.1264/jsme2.23.81

Altschul, S. F., Madden, T. L., Schäffer, A. A., Zhang, J., Zhang, Z., Miller, W., et al. (1997). Gapped BLAST and PSI-BLAST: a new generation of protein database search programs. *Nucleic Acids Res.* 25, 3389–3402. doi: 10.1093/nar/25.17.3389

Baco, A. R., and Smith, C. R. (2003). High species richness in deep-sea chemoautotrophic whale skeleton communities. *Mar. Ecol. Prog. Ser.* 260, 109–114. doi: 10.3354/meps26019

Bennett, B., Smith, C., Glaser, B., and Maybaum, H. (1994). Faunal community structure of a chemoautotrophic assemblage on whale bones in the deep Northeast Pacific-Ocean. *Mar. Ecol. Prog. Ser.* 108, 205–223. doi: 10.3354/meps108205

Bienhold, C., Pop Ristova, P., Wenzhöfer, F., Dittmar, T., and Boetius, A. (2013). How deep-sea wood falls sustain chemosynthetic life Kirchman, DL (ed). *PLoS ONE* 8:e53590. doi: 10.1371/journal.pone.0053590

Borges, L. M., Merckelbach, L. M., Sampaio, Í., and Cragg, S. M. (2014). Diversity, environmental requirements, and biogeography of bivalve wood-borers (Teredinidae) in European coastal waters. *Front. Zool.* 11:13. doi: 10.1186/1742-9994-11-13

Bright, M., and Bulgheresi, S. (2010). A complex journey: transmission of microbial symbionts. *Nat. Rev. Microbiol.* 8, 218–230. doi: 10.1038/nrmicro2262

Corbari, L., Zbinden, M., CambonBonavita, M., Gaill, F., and Compre, P. (2008). Bacterial symbionts and mineral deposits in the branchial chamber of the hydrothermal vent shrimp Rimicaris exoculata: relationship to moulting cycle. *Aquat. Biol.* 1, 225–238. doi: 10.3354/ab00024

Cunha, M. R., Rodrigues, C. F., Génio, L., Hilário, A., Ravara, A., and Pfannkuche, O. (2013). Macrofaunal assemblages from mud volcanoes in the Gulf of Cadiz: abundance, biodiversity and diversity partitioning across spatial scales. *Biogeosciences* 10, 2553–2568. doi: 10.5194/bg-10-2553-2013

Distel, D. L., Beaudoin, D. J., and Morrill, W. (2002). Coexistence of multiple proteobacterial endosymbionts in the gills of the wood-boring Bivalve Lyrodus pedicellatus (Bivalvia: Teredinidae). *Appl. Environ. Microbiol.* 68, 6292–6299. doi: 10.1128/AEM.68.12.6292-6299.2002

Distel, D. L., and Roberts, S. J. (1997). Bacterial endosymbionts in the gills of the deep-sea wood-boring bivalves Xylophaga atlantica and Xylophaga washingtona. *Biol. Bull.* 192, 253–261. doi: 10.2307/1542719

Doyle, J., and Doyle, J. (1987). Isolation of DNA from fresh plant tissue. *Focus* 12, 13–15.

Dubilier, N., Bergin, C., and Lott, C. (2008). Symbiotic diversity in marine animals: the art of harnessing chemosynthesis. *Nat. Rev. Microbiol.* 6, 725–740. doi: 10.1038/nrmicro1992

Dufour, S. C., and Felbeck, H. (2003). Sulphide mining by the superextensile foot of symbiotic thyasirid bivalves. *Nature* 426, 65–67. doi: 10.1038/nature02095

Duperron, S., Gaudron, S. M., Rodrigues, C. F., Cunha, M. R., Decker, C., and Olu, K. (2013). An overview of chemosynthetic symbioses in bivalves from the North Atlantic and Mediterranean Sea. *Biogeosciences* 10, 3241–3267. doi: 10.5194/bg-10-3241-2013

Duperron, S., Laurent, M. C. Z., Gaill, F., and Gros, O. (2008). Sulphur-oxidizing extracellular bacteria in the gills of Mytilidae associated with wood falls. *FEMS Microbiol. Ecol.* 63, 338–349. doi: 10.1111/j.1574-6941.2008.00438.x

Fagervold, S. K., Bessette, S., Romano, C., Martin, D., Plyuscheva, M., Le Bris, N., et al. (2013). Microbial communities associated with the degradation of oak wood in the Blanes submarine canyon and its adjacent open slope (NW Mediterranean). *Prog. Oceanogr.* 118, 137–143. doi: 10.1016/j.pocean.2013.07.012

Fagervold, S. K., Galand, P. E., Zbinden, M., Gaill, F., Lebaron, P., and Palacios, C. (2012). Sunken woods on the ocean floor provide diverse specialized habitats for microorganisms. *FEMS Microbiol. Ecol.* 82, 616–628. doi: 10.1111/j.1574-6941.2012.01432.x

Gaudron, S. M., Demoyencourt, E., and Duperron, S. (2012). Reproductive traits of the cold-seep symbiotic mussel idas modiolaeformis: gametogenesis and larval biology. *Biol. Bull.* 222, 6–16. doi: 10.2307/23188032

Gaudron, S. M., Pradillon, F., Pailleret, M., Duperron, S., Le Bris, N., and Gaill, F. (2010). Colonization of organic substrates deployed in deep-sea reducing habitats by symbiotic species and associated fauna. *Mar. Environ. Res.* 70, 1–12. doi: 10.1016/j.marenvres.2010.02.002

Ghiglione, J.-F., Galand, P. E., Pommier, T., Pedrós-Alió, C., Maas, E. W., Bakker, K., et al. (2012). Pole-to-pole biogeography of surface and deep marine bacterial communities. *Proc. Natl. Acad. Sci. U.S.A.* 109, 17633–17638. doi: 10.1073/pnas.1208160109

Glover, A. G., Kallstrom, B., Smith, C. R., and Dahlgren, T. G. (2005). World-wide whale worms? A new species of Osedax from the shallow north Atlantic. *Proc. R. Soc. B Biol. Sci.* 272, 2587–2592. doi: 10.1098/rspb.2005.3275

Gros, O., Darrasse, A., Durand, P., Frenkiel, L., and Moueza, M. (1996). Environmental transmission of a sulfur-oxidizing bacterial gill endosymbiont in the tropical lucinid bivalve codakia orbicularis. *Appl. Environ. Microbiol.* 62, 2324–2330.

Gros, O., Liberge, M., Heddi, A., Khatchadourian, C., and Felbeck, H. (2003). Detection of the free-living forms of sulfide-oxidizing gill endosymbionts in the lucinid habitat (Thalassia testudinum environment). *Appl. Environ. Microbiol.* 69, 6264–6267. doi: 10.1128/AEM.69.10.6264-6267.2003

Harmer, T. L., Rotjan, R. D., Nussbaumer, A. D., Bright, M., Ng, A. W., DeChaine, E. G., et al. (2008). Free-living tube worm endosymbionts found at deep-sea vents. *Appl. Environ. Microbiol.* 74, 3895–3898. doi: 10.1128/AEM.02470-07

Huber, J. A., Cantin, H. V., Huse, S. M., Welch, D. B., Sogin, M. L., and Butterfield, D. A. (2010). Isolated communities of Epsilonproteobacteria in hydrothermal vent fluids of the Mariana Arc seamounts. *FEMS Microbiol. Ecol.* 73, 538–549. doi: 10.1111/j.1574-6941.2010.00910.x

Kirchman, D. L. (2002). The ecology of *Cytophaga-Flavobacteria* in aquatic environments. *FEMS Microbiol. Ecol.* 39, 91–100. doi: 10.1111/j.1574-6941.2002.tb00910.x

Laming, S. R., Duperron, S., Cunha, M. R., and Gaudron, S. M. (2014). Settled, symbiotic, then sexually mature: adaptive developmental anatomy in the deep-sea, chemosymbiotic mussel idas modiolaeformis. *Mar. Biol.* 161, 1319–1333. doi: 10.1007/s00227-014-2421-y

Larkin, M. A., Blackshields, G., Brown, N. P., Chenna, R., McGettigan, P. A., McWilliam, H., et al. (2007). Clustal W and Clustal X version 2.0. *Bioinforma. Oxf. Engl.* 23, 2947–2948. doi: 10.1093/bioinformatics/btm404

Laurent, M. C. Z., Gros, O., Brulport, J.-P., Gaill, F., and Le Bris, N. (2009). Sunken wood habitat for thiotrophic symbiosis in mangrove swamps. *Mar. Environ. Res.* 67, 83–88. doi: 10.1016/j.marenvres.2008.11.006

Laurent, M. C. Z., Le Bris, N., Gaill, F., and Gros, O. (2013). Dynamics of wood fall colonization in relation to sulfide concentration in a mangrove swamp. *Mar. Environ. Res.* 87–88, 85–95. doi: 10.1016/j.marenvres.2013.03.007

Lee, K.-H., and Ruby, E. G. (1992). Detection of the light organ symbiont, Vibrio fischeri, in Hawaiian seawater by Using lux Gene Probes. *Appl. Environ. Microbiol.* 58, 942–947.

Lee, O. O., Wang, Y., Tian, R., Zhang, W., Shek, C. S., Bougouffa, S., et al. (2014). *In situ* environment rather than substrate type dictates microbial community structure of biofilms in a cold seep system. *Sci. Rep.* 4:3587. doi: 10.1038/srep03587
Legendre, P., and Gallagher, E. D. (2001). Ecologically meaningful transformations for ordination of species data. *Oecologia* 129, 271–280. doi: 10.1007/s004420100716

Leschine, S. B. (1995). Cellulose degradation in anaerobic environments. *Annu. Rev. Microbiol.* 49, 399–426. doi: 10.1146/annurev.mi.49.100195.002151

López-García, P., Duperron, S., Philippot, P., Foriel, J., Susini, J., and Moreira, D. (2003). Bacterial diversity in hydrothermal sediment and epsilonproteobacterial dominance in experimental microcolonizers at the Mid-Atlantic Ridge. *Environ. Microbiol.* 5, 961–976. doi: 10.1046/j.1462-2920.2003.00495.x

Lorion, J., Duperron, S., Gros, O., Cruaud, C., and Samadi, S. (2009). Several deep-sea mussels and their associated symbionts are able to live both on wood and on whale falls. *Proc. R. Soc. B Biol. Sci.* 276, 177–185. doi: 10.1098/rspb.2008.1101

McFall-Ngai, M. J. (2000). Negotiations between animals and bacteria: the 'diplo-macy' of the squid-vibrio symbiosis. *Comp. Biochem. Physiol. Part A Mol. Integr. Physiol.* 126, 471–480. doi: 10.1016/S1095-6433(00)00233-6

Nussbaumer, A. D., Fisher, C. R., and Bright, M. (2006). Horizontal endosymbiont transmission in hydrothermal vent tubeworms. *Nature* 441, 345–348. doi: 10.1038/nature04793

Palacios, C., Zbinden, M., Baco, A. R., Treude, T., Smith, C. R., Gaill, F., et al. (2006). Microbial ecology of deep-sea sunken wood: quantitative measurements of bacterial biomass and cellulolytic activities. *Cah. Biol. Mar.* 47, 415–420.

Palacios, C., Zbinden, M., Pailleret, M., Gaill, F., and Lebaron, P. (2009). Highly similar prokaryotic communities of sunken wood at shallow and deep-sea sites across the oceans. *Microb. Ecol.* 58, 737–752. doi: 10.1007/s00248-009-9538-4

Pruesse, E., Peplies, J., and Glöckner, F. O. (2012). SINA: accurate high-throughput multiple sequence alignment of ribosomal RNA genes. *Bioinformatics* 28, 1823–1829. doi: 10.1093/bioinformatics/bts252

Quast, C., Pruesse, E., Yilmaz, P., Gerken, J., Schweer, T., Yarza, P., et al. (2013). 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

Quince, C., Lanzen, A., Davenport, R. J., and Turnbaugh, P. J. (2011). Removing noise from pyrosequenced amplicons. *BMC Bioinform.* 12:38. doi: 10.1186/1471-2105-12-38

R Development Core Team (2013). *R: A Language and Environment for Statistical Computing*. Vienna: R Foundation for Statistical Computing. ISBN 3-900051-07-0. Available online at: http://www.R-project.org.

Rassa, A. C., McAllister, S. M., Safran, S. A., and Moyer, C. L. (2009). Zeta-Proteobacteria dominate the colonization and formation of microbial mats in low temperature hydrothermal vents at Loihi Seamount, Hawaii. *Geomicrobiol. J.* 8, 623–638. doi: 10.1080/01490450903263350

Reysenbach, A. L., Longnecker, K., and Kirshtein, J. (2000). Novel bacterial and archaeal lineages from an in situ growth chamber deployed at Mid-Atlantic Ridge hydrothermal vent. *Appl. Environ. Microbiol.* 66, 3798–3806. doi: 10.1128/AEM.66.9.3798-3806.2000

Rodrigues, C. F., Cunha, M. R., Génio, L., and Duperron, S. (2013). A complex picture of associations between two host mussels and symbiotic bacteria in the Northeast Atlantic. *Naturwissenschaften* 100, 21–31. doi: 10.1007/s00114-012-0985-2

Rodrigues, C. F., and Duperron, S. (2011). Distinct symbiont lineages in three thyasirid species (Bivalvia: Thyasiridae) from the eastern Atlantic and Mediterranean Sea. *Naturwissenschaften* 98, 281–287. doi: 10.1007/s00114-011-0766-3

Rouse, G. W., Goffredi, S. K., and Vrijenhoek, R. C. (2004). Osedax: bone-eating Marine Worms with Dwarf males. *Science* 305, 668–671. doi: 10.1126/science.1098650

Ruehland, C., and Dubilier, N. (2010). Gamma- and epsilonproteobacterial ectosymbionts of a shallow-water marine worm are related to deep-sea hydrothermal vent ectosymbionts. *Environ. Microbiol.* 12, 2312–2326. doi: 10.1111/j.1462-2920.2010.02256.x

Schloss, P. D., Westcott, S. L., Ryabin, T., Hall, J. R., Hartmann, M., Hollister, E. B., et al. (2009). Introducing mothra: open-source, platform-independent, community-supported software for describing and comparing microbial communities. *Appl. Environ. Microbiol.* 75, 7537–7541. doi: 10.1128/AEM.01541-09

Sipe, A. R., Wilbur, A. E., and Cary, S. C. (2000). Bacterial symbiont transmission in the wood-boring Shipworm Bankia setacea (Bivalvia: Teredinidae). *Appl. Environ. Microbiol.* 66, 1685–1691. doi: 10.1128/AEM.66.4.1685-1691.2000

Smith, C. R., and Baco, A. R. (2003). Ecology of whale falls at the deep-sea floor. *Oceanogr. Mar. Biol.* 41, 311–354.

Tamura, K., Stecher, G., Peterson, D., Filipski, A., and Kumar, S. (2013). MEGA6: molecular evolutionary genetics analysis version 6.0. *Mol. Biol. Evol.* 30, 2725–2729. doi: 10.1093/molbev/mst197

Treude, T., Smith, C. R., Wenzöhfer, F., Carney, E., Bernardino, A. F., Hannides, A. K., et al. (2009). FEATURE ARTICLE Biogeochemistry of a deep-sea whale fall: sulfate reduction, sulfide efflux and methanogenesis. *Mar. Ecol. Prog. Ser.* 382, 1–21. doi: 10.3354/meps07972

Turner, R. D. (1973). Wood-boring bivalves, opportunistic species in the deep sea. *Science* 180, 1377–1379. doi: 10.1126/science.180.4093.1377

Van Dover, C. (2000). *The Ecology of Deep-Sea Hydrothermal Vents*. Princeton, NJ: Princeton University Press.

Vrijenhoek, R. C. (2010). Genetic diversity and connectivity of deep-sea hydrothermal vent metapopulations. *Mol. Ecol.* 19, 4391–4411. doi: 10.1111/j.1365-294X.2010.04789.x

Waterbury, J. B., Calloway, C. B., and Turner, R. D. (1983). A cellulolytic nitrogen-fixing bacterium cultured from the gland of deshayes in shipworms (bivalvia: teredinidae). *Science* 221, 1401–1403. doi: 10.1126/science.221.4618.1401

Wolff, T. (1979). Macrofaunal utilization of plant remains in the deep-sea. *Sarsia* 64, 117–136.

Won, Y. J., Hallam, S. J., O'Mullan, G. D., Pan, I. L., Buck, K. R., and Vrijenhoek, R. C. (2003). Environmental acquisition of thiotrophic endosymbionts by deep-sea mussels of the genus *Bathymodiolus*. *Appl. Environ. Microbiol.* 69, 6785–6792. doi: 10.1128/AEM.69.11.6785-6792.2003

Zbinden, M., Bris, N. L., Gaill, F., and Compère, P. (2004). Distribution of bacteria and associated minerals in the gill chamber of the vent shrimp Rimicaris exculata and related biogeochemical processes. *Mar. Ecol. Prog. Ser.* 284, 237–251. doi: 10.3354/meps284237

Zhu, D., Tanabe, S.-H., Yang, C., Zhang, W., and Sun, J. (2013). Bacterial community composition of South China Sea sediments through pyrosequencing-based analysis of 16S rRNA genes. *PLoS ONE* 8:e78501. doi: 10.1371/journal.pone.0078501

**Conflict of Interest Statement:** 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.

*Received: 21 November 2014; accepted: 12 February 2015; published online: 27 February 2015.*

*Citation: Szafranski KM, Deschamps P, Cunha MR, Gaudron SM and Duperron S (2015) Colonization of plant substrates at hydrothermal vents and cold seeps in the northeast Atlantic and Mediterranean and occurrence of symbiont-related bacteria. *Front. Microbiol.* 6:162. doi: 10.3389/fmicb.2015.00162*

*This article was submitted to Aquatic Microbiology, a section of the journal Frontiers in Microbiology.*

Copyright © 2015 Szafranski, Deschamps, Cunha, Gaudron and Duperron. 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) or licensor 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.