# The contrasted evolutionary fates of deep-sea chemosynthetic mussels (Bivalvia, Bathymodiolinae)

Justine Thubaut, Nicolas Puillandre, Jean-Baptiste Faure, Corinne Cruaud, Sarah Samadi

► To cite this version:

Justine Thubaut, Nicolas Puillandre, Jean-Baptiste Faure, Corinne Cruaud, Sarah Samadi. The contrasted evolutionary fates of deep-sea chemosynthetic mussels (Bivalvia, Bathymodiolinae). Ecology and Evolution, 2013, 3 (14), pp.4748-4766. 10.1002/ece3.749. hal-01250918

HAL Id: hal-01250918

https://hal.science/hal-01250918

Submitted on 11 Jan 2016

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished 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_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
# The contrasted evolutionary fates of deep-sea chemosynthetic mussels (Bivalvia, Bathymodiolinae)ª

Justine Thubaut<fnref n="1" />, Nicolas Puillandre<fnref n="1" />, Baptiste Faure<fnref n="2" />,<fnref n="3" />, Corinne Cruaud<fnref n="4" /> & Sarah Samadi<fnref n="1" />

**Keywords**

Bathymodiolinae, chemosynthetic ecosystem, deep-sea, evolution.

**Correspondence**

Département Systématique et Evolution, Muséum National d'Histoire Naturelle, Unité Mixte de Recherche 7138 (UPMC-IRD-MNHN-CNRS), "Systématique, Adaptation et Evolution," 75005 Paris, France.

Tel: +33 (0)1 40 79 37 43; Fax: +33 (0)1 40 79 38 44; E-mail: thubaut@mnhn.fr

**Funding Information**

This study was supported by a grant of the Fondation pour la Recherche sur la Biodiversité (FRB AAP-IN-2009-009), the "Consortium National de Recherche en Génomique", and the "Service de Systématique Moléculaire" of the Muséum National d'Histoire Naturelle (CNRS UMS 2700), part of the agreement n°2005/67 between the Genoscope and the Muséum National d'Histoire Naturelle on the project "Macrophylogeny of life" directed by Guillaume Lecointre and by the network "Bibliothèque du Vivant" funded by the CNRS, the Muséum National d'Histoire Naturelle, the INRA and the CEA (Genoscope)." This work is part of the projects @ SPEED-ID "Accurate SPECiEs Delimitation and IDentification of eukaryotic biodiversity using DNA markers" proposed by F-BoL, the French Barcode of life initiative.

Received: 3 April 2013; Revised: 22 July 2013; Accepted: 24 July 2013

**Ecology and Evolution 2013; 3(14): 4748–4766**

doi: 10.1002/ece3.749

**Abstract**

Bathymodiolinae are giant mussels that were discovered at hydrothermal vents and harboring chemosynthetic symbionts. Due to their close phylogenetic relationship with seep species and tiny mussels from organic substrates, it was hypothesized that they gradually evolved from shallow to deeper environments, and specialized in decaying organic remains, then in seeps, and finally colonized deep-sea vents. Here, we present a multigene phylogeny that reveals that most of the genera are polyphyletic and/or paraphyletic. The robustness of the phylogeny allows us to revise the genus-level classification. Organic remains are robustly supported as the ancestral habitat for Bathymodiolinae. However, rather than a single step toward colonization of vents and seeps, recurrent habitat shifts from organic substrates to vents and seeps occurred during evolution, and never the reverse. This new phylogenetic framework challenges the gradualist scenarios "from shallow to deep." Mussels from organic remains tolerate a large range of ecological conditions and display a spectacular species diversity contrary to vent mussels, although such habitats are yet underexplored compared to vents and seeps. Overall, our data suggest that for deep-sea mussels, the high specialization to vent habitats provides ecological success in this harsh habitat but also brings the lineage to a kind of evolutionary dead end.

<footnote n="1">¹Département Systématique et Evolution, Muséum National d'Histoire Naturelle, Unité Mixte de Recherche 7138 (UPMC-IRD-MNHN-CNRS), "Systématique Adaptation et Evolution", 75005 Paris, France</footnote>

<footnote n="2">²Station Biologique de Roscoff, Unité Mixte de Recherche 7127, Centre National de la Recherche Scientifique, Université Pierre et Marie Curie, 29680 Roscoff, France</footnote>

<footnote n="3">³Biotope, Service Recherche et Développement, BP58 34140 Mèze, France</footnote>

<footnote n="4">⁴Genoscope, CP 5706, 91057 Evry, France</footnote>

<footnote>ªAll the new sequences are submitted to GenBank, and the accession numbers should arrive shortly.</footnote>
**Introduction**

The exploration of the deep sea is relatively recent and many
gaps remain in basic knowledge of marine biodiversity
regarding taxonomy, geographic distribution, or ecology
(Costello et al. 2010). Moreover, marine biologists are
focused on a few emblematic environments such as Antarctic
biota or hydrothermal vents which are characterized by pecu-
liar physical and chemical conditions that represent physio-
logical challenges for organisms. The consequence of such
knowledge gaps is that uncommon and spectacular
phenotypic features or behaviors are easily interpreted, under
the Panglossian paradigm (Gould and Lewontin 1979), as
adaptations to these peculiar environmental conditions. For
example, it was speculated that the gigantism of pycnogonids
from Antarctica is adaptive and stems from high polar oxygen
availability coupled with low metabolic rates (Chapelle and
Peck 1999), but physiological studies have not confirmed
this hypothesis (Woods et al. 2009), which remains contro-
versial (Klok 2009). An alternative explanation might be
that it results from phylogenetic contingency rather than
from adaptive processes. However, to disentangle adaptive
evolution from phylogenetic inertia the diversity of the
related organisms should be covered taxonomically, geo-
graphically, and ecologically (Van Buskirk 2009). For
example, when sampling is biased toward a given habitat,
characters shared by the taxa from this habitat may be mis-
takenly interpreted as an adaptation to it. Thus, because
taxon sampling strongly affects the evolutionary inferences
drawn from phylogenetic trees (Heath et al. 2008), closely
related taxa, from others habitats and/or other regions,
should be included in the phylogenetic analyses.

We here focus on the Bathymodiolinae mussels, often
considered as model organisms in the study of adaptation
to extreme marine environments (Kádár et al. 2005; Kadar
and Powell 2006; Lallier 2006). Bathymodiolinae were first
described from deep-sea hydrothermal vents that, con-
trasting with the more generally oligotrophic deep-sea
habitats, are characterized by extreme physicochemical
conditions associated with an abundant and unique fauna.
Soon after the discovery of this environment, the impres-
sive productivity was explained by the chemosynthetic
process that sustains primary production in absence of
light both through free-living bacteria, grazed by animals,
and through bacteria living in symbiosis with metazoans
(Jannasch 1985). Consequently, among the biological fea-
tures of *Bathymodiolus* (the first genus described in the
Bathymodiolinae to which are attributed most of the giant
mussels sampled at hydrothermal vents), the symbiosis
with chemosynthetic bacteria is often considered as the
key adaptation that explains the ecological success of these
mussels in these harsh habitats. This success was first
explained (Craddock et al. 1995) as resulting from a

gradual evolution toward specialization from shallow
water environments to cold seeps and finally to deep-sea
vents. Then, the evolutionary scenario was refined (Distel
et al. 2000), based on a molecular phylogenetic analysis
that included some small mussels sampled on organic
remains, sunken at the deep-sea floor. This new scenario
introduced a “wooden step” that predates the colonization
of deep-sea seeps and vents. Recently, several studies
(Duperron et al. 2009) revealed the presence of sulfide-
oxidizing bacteria in the gills of all examined deep-sea
mussels sampled on organic remains. Several phylogenetic
studies also confirmed that the evolutionary history of
vent and seep mussels is tightly tied to that of the small
mussels sampled on organic falls of the deep sea and that
the small mussels associated with organic falls should be
included within Bathymodiolinae (Samadi et al. 2007; Lo-
rion et al. 2009, 2010). These results, although still preli-
minary, suggest a more complex evolutionary history than
gradual evolution (i.e., from shallow to deep or from stan-
dard oceanic conditions to “extreme” physical and chemi-
cal conditions). However, many studies (Jones et al. 2006;
Kyuno et al. 2009; Fujiwara et al. 2010; Miyazaki et al.
2010) still suggest that gradual evolutionary scenarios in
which, for example, the presence of intracellular chemo-
synthetic bacteria is interpreted as a “final step” toward
the adaptation to extreme environments (c.f. Miyazaki
et al. 2010).

Presently, most of the studies on Bathymodiolinae,
focusing either on biology, physiology, or behavior, offer
adaptive interpretations of biological features that do not
take the most recent phylogenetic results into account
(Dixon et al. 2004; Hardivillier et al. 2004; Pruski and
Dixon 2007; Serafim et al. 2008; Mestre et al. 2009; Betten-
court et al. 2010). For example, to unravel genes specifically
involved in hydrostatic pressure and chemosynthetic envi-
ronmental adaptations, Bettencourt et al. (2010) compared
the transcriptome profiles of *Bathymodiolus azoricus* and
*Mytilus galloprovincialis*. The underlying assumption justi-
fying the adaptive inferences is that these two species are
"closely related Mytilid family members living in very dis-
tinct marine habitat." However, a recent study (Lorion
et al. 2010) confirmed that many biological features – such
as the presence of chemosynthetic bacteria – are shared by
more closely related mussels that do not inhabit the harsh
vent habitat. Thus, adaptive inferences must be based on
comparative approaches with these closest relatives rather
than with the poorly related shallow water *Mytilus* species.

The main goal of this study was thus to clarify the phylo-
genetic relationships among the Bathymodiolinae *sensus
lato* (i.e., including both Bathymodiolinae *sensus stricto* and
the small mussels from organic falls included in the same
monophyletic clade) from deep-sea environments and to
propose a new classification at the genus level that reflects
the evolutionary history. A barcoding approach was also used, including all species/Evolutionary Significant Units (ESUs) of *Bathymodiolinae sensu lato* available in the literature, in order to attempt to reassess a global classification at the genus level. This new classification aims at clarifying the context in which comparative studies of deep-sea mussels should be conducted to thoroughly test evolutionary hypotheses and thus to interpret biological features as adaptations. The phylogenetic analysis includes, in addition to the mitochondrial marker cytochrome oxydase I (COI) and the nuclear marker 28S rRNA used in the analysis performed by Lorion et al. (2010), one mitochondrial marker (16S rRNA) and four nuclear markers (18S rDNA, H3, HSP70, and Ant), for 51 terminal taxa belonging to Bathymodiolinae *s.l.* We extended the definition of Bathymodiolinae to include all deep-sea mytilids that are associated with chemosynthetic bacteria, but are not restricted to vents and seeps. The more finely resolved and strongly supported phylogeny of the Bathymodiolinae allowed us to explore the paths followed by these mussels during their evolution.

# Materials and Methods

## Taxon sampling

Taxon sampling is mainly based on Lorion et al. (2010) to which five additional terminal taxa were added. First, one specimen of the small mussels collected at the Juan de Fuca hydrothermal vents by McKiness et al. (2005) was included in the data set. These mussels were first attributed to the genus *Bathymodiolus* and labeled *B.* sp JdF (McKiness et al. 2005) and then to *Adipicola* (*Adipicola* MV in [Southward 2008; ] and *Adipicola* sp. JdF [Fontanez and Cavanaugh 2013; ]). Tissue samples of these mussels were provided by E. C. Southward. They are here labeled as *Benthomodiolus* sp. Juan de Fuca following the taxonomic revision proposed in J. Thubaut, R. von Cosel, S. Samadi, and P. Bouchet (unpubl. ms.) based on the phylogenetic results of this study. Second, two specimens collected off Papua New Guinea during the BioPapua cruise (Pante et al. 2012) and corresponding to two distinct morphospecies were added to the data set. These two specimens were collected off the Sepik River in Broken Bay together with many empty shells attributed to the two same morphospecies but also with other organisms typically found at cold seep sites (e.g., tubes of *Lamelibrachia* sp. and *Escarpia* sp., vesycomyid clams, and thyasirids). Both these two morphospecies were attributed to the genus *Gigantidas* (R. von Cosel, pers. comm.) and here referred to as *G.* sp.1 Broken Bay and *G.* sp.2 Broken Bay. We also have obtained four specimens assigned to the species *Idas simpsoni* from the collection of the Museum of London. These were trawled at 162 m in

the northern Atlantic Ocean off Rockall between Iceland and Ireland, associated with whale bones. Last, an additional unstudied morphospecies was found within the collection of the MNHN. These mussels were registered in the MNHN database as collected on vertebrate bones in South Atlantic at 3900 m depth. The shells are similar to those of *Benthomodiolus* mussels (R. von Cosel, pers. comm.). Following Samadi et al. (2007), *Modiolus modiolus* was used as an outgroup. All specimens corresponding to the 51 terminal taxa used in this study are listed in the Table 1.

## Molecular methods

For the five specimens not included in the study of Lorion et al. (2010), DNA was extracted from whole specimens (or gills only for the largest specimens) using the QIAmp® DNA Micro Kit (Qiagen, Valencia, CA). For these specimens, the two gene fragments used in Lorion et al. (2010) were amplified: (1) the COI mitochondrial gene using the universal primer LCO1490 (Folmer et al. 1994) and the reverse primer H691 (5'-GTRTTAAARTGRCGATCAAAAAT-3') designed for deep-sea mussels (Duperron et al. 2008a) and (2) a fragment of the 28S rRNA nuclear gene, covering the D1, D2, and D3 domains (Hassouna et al. 1984) using primers C1' (5'-ACCCGCTGAAT-TTAAGCAT-3') and C4 (5'-TCGGAGGGAACCAGCTACTA-3'). For the whole set of specimens, five additional gene fragments were analyzed: (1) a fragment of the 16S rRNA mitochondrial gene was amplified using the universal 16S primer LRJ-12864 (5'-CTCCGGTTTGAACTCAGATCA-3') and the primer Idas 16S A (5'-GGAR GTASGCCCTGCCCCWATGC-3') designed by Baco-Taylor (2002). Although many data are available for the ND4 mitochondrial gene, we did not use this gene that has been shown to be highly saturated (Samadi et al. 2007). (2) The 18S rDNA nuclear gene was amplified in three overlapping fragments using three pairs of primers: 1F and 5R, 3F and Bi, A2 and 9R (Giribet et al. 1996; Distel 2000; Okusu et al. 2003), (3) a fragment of the histone nuclear gene H3 was amplified using the primers H3F1 (5'-ATGGCTCGTACCAAGCAGACVGC-3') and H3R1 (5'-ATA TCCTTRGGCATRATRGTGAC-3'), (4) a fragment of the 70-kDa heat-shock protein (HSP70) nuclear gene was amplified using primers HSP70F (5'-GGGAAAGTTGA CATTATTGCCAATG-3') and HSP70R (5'-ATTCATAAA TTCTGTCAACATTTTCTGT-3') (developed by BF), and (5) a fragment of the gene encoding the adenosine nucleotide (ADP/ATP) translocase (Ant) was amplified using the pair of primers designed by Audzijonyte and Vrijenhoek (2009). Polymerase chain reaction (PCR) were performed in a 25 µL, final volume, containing approximately 3-ng template DNA, 1.5-mmol/L MgCl₂, 0.26 mmol/L of each dNTPs, 0.3 µmol/L of each primer, 5% dimethyl sulfoxide
<table><caption>**Table 1.** Specimen collection sites, maximum shell size, and GenBank accession numbers (sequences obtained in this study are indicated with stars) for Mytilidae included in phylogenetic analyses.</caption><thead><tr><th>Genus</th><th>Species</th><th>Depth range (m)</th><th>Habitat</th><th>Maximum shell size (mm)</th><th>COI</th><th>28S</th><th>18S</th><th>16S</th><th>H3</th><th>HSP70</th><th>Ant</th></tr></thead><tbody><tr><td><i>Modiolus</i></td><td><i>modiolus</i></td><td>10–300</td><td>Intertidal subtidal</td><td><strong>147</strong></td><td>FJ890501<sup>(1)</sup></td><td>EF526455<sup>(2)</sup></td><td>KF611701*</td><td>KF611732*</td><td>KF720595*</td><td></td><td>KF720516*</td></tr><tr><td rowspan="3"><i>Benthomodiolus</i></td><td rowspan="3"><i>lignocola</i></td><td>1180</td><td>Wood, bone</td><td>18.1</td><td>AY275545<sup>(3)</sup></td><td>AY781131<sup>(4)</sup></td><td>AF221648<sup>(5)</sup></td><td>KF611733*</td><td>KF720596*</td><td></td><td></td></tr><tr><td>2420</td><td>Vent</td><td>16</td><td>KF611694*</td><td>KF611699*</td><td>KF611702*</td><td>KF611734*</td><td>KF720597*</td><td></td></tr><tr><td>3900</td><td>Bone</td><td>8.8</td><td>KF611691*</td><td>KF611698*</td><td>KF611703*</td><td>KF611733*</td><td>KF720593*</td><td>KF720560*</td></tr><tr><td rowspan="5"><i>Vulcanidas</i></td><td>sp. Juan de Fuca</td><td>150–785</td><td>Wood</td><td>36.7</td><td>FJ937079<sup>(1)</sup></td><td>GU065791<sup>(1)</sup></td><td>KF611704*</td><td>KF611736*</td><td>KF720598*</td><td>KF720592*</td><td>KF720517*</td></tr><tr><td>sp. South Atlantic</td><td>150–785</td><td>Wood</td><td>36.7</td><td>FJ937079<sup>(1)</sup></td><td>GU065791<sup>(1)</sup></td><td>KF611704*</td><td>KF611736*</td><td>KF720598*</td><td>KF720592*</td><td>KF720517*</td></tr><tr><td>ESU E</td><td>150–785</td><td>Wood</td><td>36.7</td><td>FJ937079<sup>(1)</sup></td><td>GU065791<sup>(1)</sup></td><td>KF611704*</td><td>KF611736*</td><td>KF720598*</td><td>KF720592*</td><td>KF720517*</td></tr><tr><td>ESU F</td><td>275–560</td><td>Wood</td><td>35</td><td>FJ937127<sup>(1)</sup></td><td>GU065809<sup>(1)</sup></td><td>KF611705*</td><td>KF611737*</td><td>KF720599*</td><td>KF720578*</td><td>KF720518*</td></tr><tr><td>SAL3</td><td>400–1085</td><td>Wood</td><td>38.3</td><td>DQ340772<sup>(6)</sup></td><td>DQ863946<sup>(1)</sup></td><td>DQ340800<sup>(6)</sup></td><td>KF611738*</td><td>KF720600*</td><td>KF720582*</td><td>KF720519*</td></tr><tr><td rowspan="3"><i>Tamu Lignomodiolus</i></td><td>fisheri</td><td>140–504</td><td>Vent</td><td><strong>178</strong></td><td>FJ767936<sup>(1)</sup></td><td>FJ767937<sup>(1)</sup></td><td>KF611706*</td><td>KF611739*</td><td>KF720601*</td><td>KF720558*</td><td>KF720520*</td></tr><tr><td>ESU S'</td><td>546–650</td><td>Seep</td><td><strong>58</strong></td><td>AY649803<sup>(4)</sup></td><td>AY781132<sup>(4)</sup></td><td>AF221642<sup>(5)</sup></td><td></td><td></td><td></td><td></td></tr><tr><td>Western Pacific</td><td>230–380</td><td>Wood</td><td>18.8</td><td>FJ937240<sup>(1)</sup></td><td>GU065816<sup>(1)</sup></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td rowspan="5"><i>Gigantidas</i></td><td>sp.1</td><td>190–400</td><td>Wood</td><td>18.8</td><td>FJ937258<sup>(1)</sup></td><td>GU065829<sup>(1)</sup></td><td>KF611707*</td><td>KF611740*</td><td>KF720602*</td><td>KF720590*</td><td>KF720521*</td></tr><tr><td>Broken Bay</td><td>190–400</td><td>Wood</td><td>18.8</td><td>FJ937258<sup>(1)</sup></td><td>GU065829<sup>(1)</sup></td><td>KF611707*</td><td>KF611740*</td><td>KF720602*</td><td>KF720590*</td><td>KF720521*</td></tr><tr><td>SAL4</td><td>94–886</td><td>Wood</td><td>25</td><td>DQ340776<sup>(6)</sup></td><td>DQ863947<sup>(6)</sup></td><td>DQ340796<sup>(6)</sup></td><td>KF611741*</td><td>KF720603*</td><td>KF720589*</td><td>KF720522*</td></tr><tr><td>ESU R</td><td>490</td><td>Wood</td><td>8.2</td><td>FJ937239<sup>(1)</sup></td><td>GU065877<sup>(1)</sup></td><td>KF611708*</td><td>KF611742*</td><td>KF720604*</td><td>KF720588*</td><td>KF720523*</td></tr><tr><td>ESU Q</td><td>100–130</td><td>Wood</td><td>7.9</td><td>FJ937230<sup>(1)</sup></td><td>GU065875<sup>(1)</sup></td><td>KF611709*</td><td>KF611743*</td><td>KF720605*</td><td>KF720561*</td><td>KF720524*</td></tr><tr><td rowspan="2"><i>Mauritanicus tangaroa</i></td><td>Western Pacific</td><td>150–670</td><td>Wood</td><td>9.6</td><td>FJ937161<sup>(1)</sup></td><td>GU065778<sup>(1)</sup></td><td>KF611710*</td><td>KF611744*</td><td>KF720606*</td><td>KF720591*</td><td>KF720525*</td></tr><tr><td>Western Pacific</td><td>361–750</td><td>Seep</td><td><strong>85</strong></td><td>KF611693*</td><td>KF611696*</td><td></td><td>KF611745*</td><td>KF720607*</td><td>KF720555*</td><td>KF720526*</td></tr><tr><td><i>taiwanensis</i></td><td>Western Pacific</td><td>200–355</td><td>Vent</td><td><strong>56</strong></td><td>GU966638<sup>(1)</sup></td><td>GU966641<sup>(1)</sup></td><td>KF611711*</td><td>KF611746*</td><td>KF720608*</td><td>KF720566*</td><td>KF720527*</td></tr><tr><td rowspan="2"><i>mauritanicus tangaroa</i></td><td>540–2222</td><td>Seep</td><td><strong>127</strong></td><td>FJ890502<sup>(1)</sup></td><td>FJ890504<sup>(1)</sup></td><td>KF611712*</td><td>KF611747*</td><td>KF720609*</td><td>KF720565*</td><td>KF720528*</td><td>KF720528*</td></tr><tr><td>920–1205</td><td>Seep</td><td><strong>200</strong></td><td>AY608439<sup>(3)</sup></td><td>AY781149<sup>(4)</sup></td><td>AY649820<sup>(4)</sup></td><td>KF611748*</td><td>KF720610*</td><td>KF720572*</td><td>KF720529*</td><td>KF720529*</td></tr></tbody><tfoot><tr><td colspan="12" style="text-align: right;">(Continued)</td></tr></tfoot></table>
<table><caption>**Table 1. Continued.**</caption><thead><tr><th>Genus</th><th>Species</th><th>Ocean</th><th>Depth range (m)</th><th>Habitat</th><th>Maximum shell size (mm)</th><th>COI</th><th>28S</th><th>18S</th><th>16S</th><th>H3</th><th>HSP70</th><th>Ant</th></tr></thead><tbody><tr><td rowspan="5">sp.2 Broken Bay<br>gladius</td><td>Western Pacific</td><td>361-750</td><td>Seep</td><td rowspan="5">Wood, bone<br>Wood, bone</td><td rowspan="5">30</td><td rowspan="5">EU702319<sup>(2)</sup><br>EU702315<sup>(2)</sup><br>EU702315<sup>(2)</sup><br>EU702315<sup>(2)</sup></td><td rowspan="5">316</td><td rowspan="5">AY649802<sup>(4)</sup><br>AY649802<sup>(4)</sup><br>AY649802<sup>(4)</sup><br>AY649802<sup>(4)</sup></td><td rowspan="5">KF611713*</td><td rowspan="5">KF611713*</td><td rowspan="5">KF611713*</td><td rowspan="5">KF611713*</td><td rowspan="5">KF611713*</td></tr><tr><td>Western Pacific</td><td>300-460</td><td>Vent</td></tr><tr><td>Western Pacific</td><td>441</td><td></td></tr><tr><td>Western Pacific</td><td>431-493</td><td></td></tr><tr><td>Western Pacific</td><td>400-1767</td><td>Wood</td></tr><tr><td rowspan="5">Nypamodiolus<br>crypta (B')</td><td>Western Pacific</td><td>360-370</td><td>Wood</td><td rowspan="5">Wood</td><td rowspan="5">42</td><td rowspan="5">DQ340773<sup>(6)</sup><br>DQ340773<sup>(6)</sup><br>DQ340773<sup>(6)</sup><br>DQ340773<sup>(6)</sup></td><td rowspan="5">DQ863945<sup>(6)</sup><br>DQ863945<sup>(6)</sup><br>DQ863945<sup>(6)</sup><br>DQ863945<sup>(6)</sup></td><td rowspan="5">DQ340798<sup>(6)</sup><br>DQ340798<sup>(6)</sup><br>DQ340798<sup>(6)</sup><br>DQ340798<sup>(6)</sup></td><td rowspan="5">KF611752*</td><td rowspan="5">KF611752*</td><td rowspan="5">KF611752*</td><td rowspan="5">KF611752*</td><td rowspan="5">KF611752*</td></tr><tr><td>Western Pacific</td><td>190-567</td><td>Wood</td></tr><tr><td>Western Pacific</td><td>220-560</td><td>Wood</td></tr><tr><td>Western Pacific</td><td>162</td><td>Bone</td><td rowspan="2">86.2</td><td rowspan="2">KF611695*</td><td rowspan="2">KF611700*</td><td rowspan="2">KF611731*</td><td rowspan="2">KF611731*</td><td rowspan="2">KF611731*</td><td rowspan="2">KF611731*</td><td rowspan="2">KF611731*</td></tr><tr><td>Western Pacific</td><td>1627</td><td>Vent</td></tr><tr><td rowspan="5">Bathymodiolus<br>simpsoni<br>manusensis</td><td>Western Pacific</td><td>880</td><td>Bone</td><td rowspan="5">Wood</td><td rowspan="5">29</td><td rowspan="5">FJ937033<sup>(1)</sup><br>FJ937033<sup>(1)</sup><br>FJ937033<sup>(1)</sup><br>FJ937033<sup>(1)</sup></td><td rowspan="5">GU065879<sup>(1)</sup><br>GU065879<sup>(1)</sup><br>GU065879<sup>(1)</sup><br>GU065879<sup>(1)</sup></td><td rowspan="5">KF611719*</td><td rowspan="5">KF611719*</td><td rowspan="5">KF611719*</td><td rowspan="5">KF611719*</td><td rowspan="5">KF611719*</td></tr><tr><td>Western Pacific</td><td>800-1060</td><td>Bone</td></tr><tr><td>Western Pacific</td><td>3023-3510</td><td>Vent</td><td rowspan="4">119</td><td rowspan="4">AY649796<sup>(4)</sup><br>AY649795<sup>(4)</sup><br>AY649794<sup>(4)</sup><br>AY649794<sup>(4)</sup></td><td rowspan="4">AF221640<sup>(5)</sup><br>AF221640<sup>(5)</sup><br>AF221639<sup>(5)</sup><br>AF221639<sup>(5)</sup></td><td rowspan="4">KF611758*</td><td rowspan="4">KF611758*</td><td rowspan="4">KF611758*</td><td rowspan="4">KF611758*</td><td rowspan="4">KF611758*</td></tr><tr><td>Atlantic</td><td>866-2330</td><td>Vent</td></tr><tr><td>Atlantic</td><td>3314</td><td>Seep</td><td rowspan="3">230</td><td rowspan="3">AY649794<sup>(4)</sup><br>AY649794<sup>(4)</sup><br>AY649794<sup>(4)</sup><br>AY649794<sup>(4)</sup></td><td rowspan="3">AF221639<sup>(5)</sup><br>AF221639<sup>(5)</sup><br>AF221639<sup>(5)</sup><br>AF221639<sup>(5)</sup></td><td rowspan="3">KF611759*</td><td rowspan="3">KF611759*</td><td rowspan="3">KF611759*</td><td rowspan="3">KF611759*</td><td rowspan="3">KF611759*</td></tr><tr><td>Atlantic</td><td>1000-3170</td><td>Seep</td></tr><tr><td>Atlantic</td><td>3589</td><td>Vent</td><td rowspan="2">140</td><td rowspan="2">AY649799<sup>(4)</sup><br>AY649799<sup>(4)</sup><br>AY649799<sup>(4)</sup><br>AY649799<sup>(4)</sup></td><td rowspan="2">AY649824<sup>(4)</sup><br>AY649824<sup>(4)</sup><br>AY649824<sup>(4)</sup><br>AY649824<sup>(4)</sup></td><td rowspan="2">KF611760*</td><td rowspan="2">KF611760*</td><td rowspan="2">KF611760*</td><td rowspan="2">KF611760*</td><td rowspan="2">KF611760*</td></tr><tr><td>Western Pacific</td><td>2460-2747</td><td>Vent</td><td rowspan="4">200</td><td rowspan="4">GU966639<sup>(1)</sup><br>AF456317<sup>(8)</sup><br>AY781140<sup>(4)</sup></td><td rowspan="4">GU966640<sup>(1)</sup><br>AY781140<sup>(4)</sup><br>AY781140<sup>(4)</sup></td><td rowspan="4">KF611761*</td><td rowspan="4">KF611761*</td><td rowspan="4">KF611761*</td><td rowspan="4">KF611761*</td><td rowspan="4">KF611761*</td></tr><tr><td>thermophilus</td><td>Eastern Pacific</td><td>2331</td><td>Vent</td></tr><tr><td>aff. thermophilus</td><td>Eastern Pacific</td><td>2222-3314</td><td>Seep</td><td rowspan="2">180</td><td rowspan="2">AY649798<sup>(4)</sup><br>DQ340775<sup>(6)</sup></td><td rowspan="2">AY781135<sup>(4)</sup><br>DQ863944<sup>(6)</sup><br>DQ340794<sup>(6)</sup><br>DQ340794<sup>(6)</sup></td><td rowspan="2">KF611761*</td><td rowspan="2">KF611761*</td><td rowspan="2">KF611761*</td><td rowspan="2">KF611761*</td><td rowspan="2">KF611761*</td></tr><tr><td>brooksi</td><td>Atlantic</td><td>408-1356</td><td>Seep</td></tr><tr><td>SAL1</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr></tbody></table>

(Continued)
<table><caption>**Table 1. Continued.**</caption><thead><tr><th>Genus</th><th>Species</th><th>Ocean</th><th>Depth range (m)</th><th>Habitat</th><th>Maximum shell size (mm)</th><th>COI</th><th>28S</th><th>18S</th><th>16S</th><th>H3</th><th>HSP70</th><th>Ant</th></tr></thead><tbody><tr><td></td><td>Western Pacific</td><td>Western Pacific</td><td>180–1390</td><td>Wood, bone</td><td>8.7</td><td>FJ937222<sup>(1)</sup></td><td>GU065846<sup>(1)</sup></td><td>KF611721<sup>*</sup></td><td>KF611762<sup>*</sup></td><td>KF720625<sup>*</sup></td><td>KF720557<sup>*</sup></td><td>KF720543<sup>*</sup></td></tr><tr><td>ESU P</td><td>Western Pacific</td><td>Western Pacific</td><td>473–890</td><td>Wood</td><td>6.7</td><td>FJ937211<sup>(1)</sup></td><td>GU065763<sup>(1)</sup></td><td>KF611722<sup>*</sup></td><td>KF611763<sup>*</sup></td><td>KF720626<sup>*</sup></td><td>KF720581<sup>*</sup></td><td>KF720544<sup>*</sup></td></tr><tr><td>ESU O</td><td>Western Pacific</td><td>Western Pacific</td><td>500–540</td><td>Wood</td><td>6</td><td>FJ937192<sup>(1)</sup></td><td>GU065868<sup>(1)</sup></td><td>KF611723<sup>*</sup></td><td>KF611764<sup>*</sup></td><td>KF720627<sup>*</sup></td><td>KF720568<sup>*</sup></td><td>KF720545<sup>*</sup></td></tr><tr><td>ESU K</td><td>Western Pacific</td><td>Eastern Pacific</td><td>960–1910</td><td>Vent, wood, bone</td><td>8.6</td><td>AY275546<sup>(3)</sup></td><td>AY781146<sup>(4)</sup></td><td>AF221645<sup>(5)</sup></td><td></td><td></td><td></td><td></td></tr><tr><td>Washingtonius</td><td>Eastern Pacific</td><td>Western Pacific</td><td>556–1724</td><td>Wood</td><td>15.6</td><td>EU702357<sup>(2)</sup></td><td>EU683275<sup>(2)</sup></td><td>KF611724<sup>*</sup></td><td>KF611765<sup>*</sup></td><td>KF720628<sup>*</sup></td><td>KF720587<sup>*</sup></td><td>KF720546<sup>*</sup></td></tr><tr><td>sp. D</td><td>Western Pacific</td><td>Western Pacific</td><td>275–1285</td><td>Wood, bone</td><td>16.7</td><td>EU702376<sup>(2)</sup></td><td>EU683260<sup>(6)</sup></td><td>KF611725<sup>*</sup></td><td>KF611766<sup>*</sup></td><td>KF720629<sup>*</sup></td><td>KF720586<sup>*</sup></td><td>KF720547<sup>*</sup></td></tr><tr><td>sp. C</td><td>Western Pacific</td><td>Western Pacific</td><td>590–720</td><td>Wood</td><td>10.8</td><td>FJ937202<sup>(1)</sup></td><td>GU065845<sup>(1)</sup></td><td>KF611726<sup>*</sup></td><td>KF611767<sup>*</sup></td><td>KF720630<sup>*</sup></td><td>KF720583<sup>*</sup></td><td>KF720548<sup>*</sup></td></tr><tr><td>ESU M</td><td>Western Pacific</td><td>Western Pacific</td><td>450–1010</td><td>Wood</td><td>8.6</td><td>FJ937193<sup>(1)</sup></td><td>GU065767<sup>(1)</sup></td><td>KF611727<sup>*</sup></td><td>KF611768<sup>*</sup></td><td>KF720631<sup>*</sup></td><td>KF720574<sup>*</sup></td><td>KF720549<sup>*</sup></td></tr><tr><td>ESU L</td><td>Western Pacific</td><td>Western Pacific</td><td>800–1290</td><td>Wood, bone</td><td>10</td><td>FJ937205<sup>(1)</sup></td><td>GU065843<sup>(1)</sup></td><td>KF611728<sup>*</sup></td><td>KF611769<sup>*</sup></td><td>KF720632<sup>*</sup></td><td>KF720573<sup>*</sup></td><td>KF720552<sup>*</sup></td></tr><tr><td>ESU N</td><td>Western Pacific</td><td>Western Pacific</td><td>441–1866</td><td>Wood, bone</td><td>10</td><td>EU702333<sup>(2)</sup></td><td>EU683288<sup>(2)</sup></td><td>KF611729<sup>*</sup></td><td>KF611770<sup>*</sup></td><td>KF720633<sup>*</sup></td><td>KF720559<sup>*</sup></td><td>KF720550<sup>*</sup></td></tr><tr><td>iwaotakii (A')</td><td>Western Pacific</td><td>Western Pacific</td><td>490–2307</td><td>Wood</td><td>10</td><td>EU702322<sup>(2)</sup></td><td>EU683295<sup>(2)</sup></td><td></td><td></td><td></td><td>KF720634<sup>*</sup></td><td>KF720569<sup>*</sup></td></tr><tr><td>iwaotakii (A')</td><td>Western Pacific</td><td>Western Pacific</td><td>650</td><td>Atlantic</td><td>6.7</td><td>AY649804<sup>(4)</sup></td><td>AY781145<sup>(4)</sup></td><td>AF221647<sup>(5)</sup></td><td></td><td></td><td></td><td></td></tr><tr><td>macdonaldi</td><td>Atlantic</td><td>Mediter- ranean</td><td>2129</td><td>Seep</td><td>13</td><td>FJ158585<sup>(9)</sup></td><td>FJ159555<sup>(1)</sup></td><td>KF611730<sup>*</sup></td><td>KF611772<sup>*</sup></td><td>KF720635<sup>*</sup></td><td>KF720570<sup>*</sup></td><td>KF720553<sup>*</sup></td></tr></tbody></table>

Shell sizes higher than 50 mm are given in bold. References for sequences: (1) Lorion et al. (2010), (2) Lorion et al. (2009), (3) Smith et al. (2004), (4) Jones et al. (2006), (5) Distel et al. (2000), (6) Samadi et al. (2007), (7) Audzijonyte and Vrijenhoek (2009), (8) Won et al. (2003) and (9) Lorion et al. (2012). COI, cytochrome oxydase I.
(DMSO), and 0.75 unit of Taq Polymerase (Qbiogene, Illkirch, France) Amplification products were generated by an initial denaturation step of 4 min at 94°C followed by 30–35 cycles of denaturation at 94°C for 40 sec, annealing at 48°C for ND4 gene, 50°C for COI and first and second fragment of 18S gene, 52°C for 28S gene and third fragment of 18S gene, 55°C for 16S, H3 and HSP70 genes for 50 sec, and extension at 72°C for 1 min. A final elongation was performed at 72°C for 10 min. The gene Ant was amplified according to PCR conditions described in Audzijonyte and Vrijenhoek (2009). All PCR products were purified and sequenced at the Genoscope for both DNA strands. The 203 new sequences were deposited in BOLD and GenBank.

## Phylogenetic analyses

Chromatograms were edited with Sequencher 4.1.4, DNA sequences were aligned with the Clustal W module in Mega 4.0 (Tamura et al. 2007), and ambiguous sites were excluded. For each of the seven genes, one specimen of each species or ESUs (Samadi et al. 2007; Lorion et al. 2009, 2010) was used to evaluate the phylogenetic relationships among them. Best-fitting substitution models were estimated for each gene separately and for a data set that combined the seven single-gene data sets using MrModelTest 3.7 with the Akaike information criterion (AIC) (Posada and Crandall 1998). Each of the seven single-gene data sets and the combined data set were fitted to a general time-reversible model with a proportion of invariant sites and gamma-distributed rates (GTR+I+Γ). Phylogenetic reconstructions were performed, for both single-gene data sets and the combined data set, using maximum likelihood (ML) and Bayesian Inference (BI) approaches. ML analyses were performed using the software Treefinder (Jobb et al. 2004) and robustness of the nodes was assessed using non-parametric bootstrapping (1000 replicates) (Felsenstein 1985). BI analyses were performed using MrBayes (Ronquist and Huelsenbeck 2003) as implemented on the Bioportal cluster (University of Oslo, Bioportal: http://www.bioportal.uio.no). Eight Markov chains and two parallel analyses were run over 50M generations, sampled every 5000th step, and associated with a heating temperature of 0.01. To assess if the two independent runs converged, likelihood curves and standard deviation of split frequencies were analyzed using Tracer v1.4.1 (Rambaud and Drummond 2007) and AWTY website (http://king2.scs.fsu.edu/CEBProjects/awty/awty_start.php) (Nylander et al. 2008).

## Molecular assignation of species to genera

The taxon sampling used for the phylogenetic analysis covers the lineage diversity but not for all lineages the

species diversity. Thus, to compare the species diversity among redefined genera we use a DNA barcoding-like approach (i.e., based on COI data), to assign species to genera. We first gathered from World Register of Marine Species (WORMS; http://www.marinespecies.org) a provisional list of valid names for each genus, and gathered, when available, the associated COI data from GenBank. Second, we gathered the COI data for the species identified in the literature using molecular data, but that have not yet been attributed to species names. This set of species was identified by Miyazaki et al. (2004), Iwasaki et al. (2006), Fujita et al. (2009), and Kyuno et al. (2009) using the COI and ND4 gene fragments. These taxa were labeled according to Kyuno et al. (2009) and were, first, tentatively attributed to either named species or ESUs of our data set. In the literature, the mean intraspecific genetic distance estimated using the Kimura 2-parameters (K2P) distance was about ~1% in Won et al. (2003) and Miyazaki et al. (2004), and of 1.8% in Lorion et al. (2010). We thus used a conservative threshold of 2% K2P genetic distance to attribute a sequence to defined species or ESU. Second, an extended COI data set was constituted to perform a phylogenetic analysis. This data set included a reference data set (i.e., one sequence per named species or ESUs of our data set that are all attributed to a genus name based on the seven-genes phylogeny), and all COI data gathered from GenBank. This extended COI data set was analyzed using BI and ML approaches (see Table 2) to determine among the reference data set the closest relatives of each of the additional species or ESUs gathered from the literature. The attribution to a genus name was based on the sister relationship with a species – or a clade of species – of the reference data set. For one species name (B. elongatus), no COI data were available, although data for another gene fragment (ITS2) were available in GenBank. In this case, the attribution to genus followed the result of the phylogenetic analysis provided by the authors (see Table 2). Last, for named species for which no genetic data were available, species were tentatively placed in genera considering the available morphological data.

## Character evolution

Ancestral character state reconstruction of habitat use and depth was conducted using the ML method implemented in Mesquite 2.75 (Maddison and Maddison 2011). Taxa were assigned to four different habitats (hydrothermal vent, cold seep, organic falls, nonchemosynthetic environment). The character “habitat” was coded as “0” for nonchemosynthetic environment, “1” for hydrothermal vent, “2” for cold seep, and “3” for organic remain. The character
<table><caption>**Table 2.** Specimen collection sites, GenBank accession numbers for Bathymodiolinae included in the Barcode analysis (sequences obtained in this study are indicated with stars), and localization of symbionts.</caption><thead><tr><th>Genus</th><th>Species</th><th>Locality</th><th>Habitat</th><th>Depth (m)</th><th>COI</th><th>ITS2</th><th>Localization of symbionts</th></tr></thead><tbody><tr><td rowspan="3">Benthomodiolus Dell, 1987</td><td><strong>lignocola</strong>* Dell, 1987</td><td>WP</td><td>w, b</td><td><strong>1180</strong></td><td>AY275545<sup>(1)</sup></td><td></td><td></td></tr><tr><td>geikotsucola Okutani & Miyazaki, 2007</td><td>WP</td><td>b</td><td><strong>4051</strong></td><td>AB257513<sup>(2)</sup></td><td></td><td></td></tr><tr><td>abyssicola (Knudsen, 1970)</td><td>EP</td><td>w</td><td><strong>3270–3670</strong></td><td></td><td></td><td></td></tr><tr><td rowspan="4">Vulcanidas Cosel and Marshall, 2010</td><td><strong>sp. Juan de Fuca</strong></td><td>EP</td><td>v</td><td><strong>2420</strong></td><td>KF611694<sup>(3)</sup></td><td></td><td>Extracellular?<sup>(12)</sup></td></tr><tr><td><strong>sp. South Atlantic</strong></td><td>SA</td><td>b</td><td><strong>3900</strong></td><td>KF611691<sup>(3)</sup></td><td></td><td></td></tr><tr><td><strong>insolatus</strong>* Cosel and Marshall, 2010</td><td>WP</td><td>v</td><td>140–504</td><td>FJ767936<sup>(4)</sup></td><td></td><td>Intracellular?<sup>(13)</sup></td></tr><tr><td><strong>SAL3</strong></td><td>WP</td><td>w</td><td>400–1085</td><td>DQ340772<sup>(5)</sup></td><td></td><td></td></tr><tr><td rowspan="3">Tamu Gustafson, Turner, Lutz & Vrijenhoek, 1998</td><td><strong>ESU E</strong></td><td>WP</td><td>w</td><td>150–785</td><td>FJ937079<sup>(4)</sup></td><td></td><td></td></tr><tr><td><strong>ESU F</strong></td><td>WP</td><td>w</td><td>275–560</td><td>FJ937127<sup>(4)</sup></td><td></td><td></td></tr><tr><td><strong>fisheri</strong>* Gustafson, Turner, Lutz & Vrijenhoek, 1998</td><td>GM</td><td>s</td><td>546–650</td><td>AY649803<sup>(6)</sup></td><td></td><td></td></tr><tr><td rowspan="6">Lignomodiolus Cosel & Thubaut in prep.</td><td><strong>miniboomerang</strong>* (ESU G) Cosel & Thubaut in prep.</td><td>WP</td><td>w</td><td>150–670</td><td>FJ937161<sup>(4)</sup></td><td></td><td>Extracellular<sup>(14)</sup></td></tr><tr><td><strong>ESU S'</strong></td><td>WP</td><td>w</td><td>230–380</td><td>FJ937240<sup>(4)</sup></td><td></td><td></td></tr><tr><td><strong>ESU S''</strong></td><td>WP</td><td>w</td><td>190–400</td><td>FJ937245<sup>(4)</sup></td><td></td><td></td></tr><tr><td><strong>SAL4</strong></td><td>WP</td><td>w</td><td>94–886</td><td>DQ340776<sup>(5)</sup></td><td></td><td></td></tr><tr><td><strong>ESU R</strong></td><td>WP</td><td>w</td><td>490</td><td>FJ937238<sup>(4)</sup></td><td></td><td></td></tr><tr><td><strong>ESU Q</strong></td><td>WP</td><td>w</td><td>100–130</td><td>FJ937229<sup>(4)</sup></td><td></td><td></td></tr><tr><td rowspan="11">Gigantidas Cosel & Marshall, 2003</td><td><strong>gladius</strong>* Cosel & Marshall, 2003</td><td>WP</td><td>v</td><td>300–460</td><td>AY649802<sup>(6)</sup></td><td></td><td>Intracellular<sup>(15)</sup></td></tr><tr><td><strong>horikoshii</strong> Hashimoto & Yamane, 2005</td><td>WP</td><td>v</td><td>486</td><td>AB257538<sup>(2)</sup></td><td></td><td></td></tr><tr><td><strong>Nikko G. sp.</strong></td><td>WP</td><td>v</td><td>485</td><td>AB257544<sup>(2)</sup></td><td></td><td></td></tr><tr><td><strong>Sumisu G. sp.</strong></td><td>WP</td><td>v</td><td>676–686</td><td>AB257553<sup>(2)</sup></td><td></td><td></td></tr><tr><td><strong>Ashizuri sp.</strong></td><td>WP</td><td>s</td><td>575</td><td>AB257529<sup>(2)</sup></td><td></td><td></td></tr><tr><td><strong>childressi</strong> Gustafson, Turner, Lutz & Vrijenhoek, 1998</td><td>GM</td><td>s</td><td><strong>1859</strong></td><td>AB257532<sup>(2)</sup></td><td></td><td>Intracellular<sup>(16)</sup></td></tr><tr><td><strong>platifrons</strong> Hashimoto & Okutani, 1994</td><td>WP</td><td>v, s</td><td><strong>1180</strong></td><td>AB101419<sup>(7)</sup></td><td></td><td>Intracellular<sup>(17)</sup></td></tr><tr><td><strong>japonicus</strong> Hashimoto & Okutani, 1994</td><td>WP</td><td>s</td><td><strong>1170</strong></td><td>AB101422<sup>(7)</sup></td><td></td><td>Intracellular<sup>(17)</sup></td></tr><tr><td><strong>mauritanicus</strong> Cosel, 2002</td><td>NA</td><td>s</td><td>540–2222</td><td>FJ890502<sup>(4)</sup></td><td></td><td></td></tr><tr><td><strong>tangaroa</strong> Cosel & Marshall, 2003</td><td>WP</td><td>s</td><td>920–1205</td><td>AY608439<sup>(6)</sup></td><td></td><td></td></tr><tr><td><strong>Sissano B. sp.3</strong></td><td>WP</td><td>s</td><td><strong>1881</strong></td><td>AB257551<sup>(2)</sup></td><td></td><td></td></tr><tr><td rowspan="3">Sissano B. sp.1</td><td><strong>securiformis</strong> Okutani, Fujikura & Sasaki, 2003</td><td>WP</td><td>s</td><td>637–642</td><td>AB170048<sup>(8)</sup></td><td></td><td></td></tr><tr><td><strong>hirtus</strong> Okutani, Fujikura & Sasaki, 2003</td><td>WP</td><td>s</td><td>637–642</td><td>AB170047<sup>(8)</sup></td><td></td><td></td></tr><tr><td><strong>anteumbonatus</strong> Cosel & Janssen, 2008</td><td>WP</td><td>s</td><td><strong>1574–1628</strong></td><td></td><td></td><td></td></tr><tr><td rowspan="4">Sissano B. sp.2</td><td><strong>edisonensis</strong> Cosel & Janssen, 2008</td><td>WP</td><td>s</td><td><strong>1574–1629</strong></td><td></td><td></td><td></td></tr><tr><td><strong>taiwanensis</strong> Cosel, 2008</td><td>WP</td><td>v</td><td>271</td><td>GU966638<sup>(4)</sup></td><td></td><td></td></tr><tr><td><strong>Kikaijima B. sp.</strong></td><td>WP</td><td>s</td><td><strong>1430</strong></td><td>AB257556<sup>(2)</sup></td><td></td><td></td></tr><tr><td><strong>Sissano B. sp.1</strong></td><td>WP</td><td>s</td><td><strong>1646</strong></td><td>AB257548<sup>(2)</sup></td><td></td><td></td></tr><tr><td rowspan="2">Chamorro B. sp.</td><td><strong>Sissano B. sp.2</strong></td><td>WP</td><td>s</td><td><strong>1881</strong></td><td>AB257547<sup>(2)</sup></td><td></td><td></td></tr><tr><td><strong>Chamorro B. sp.</strong></td><td>WP</td><td>s</td><td><strong>2899</strong></td><td>AB257530<sup>(2)</sup></td><td></td><td></td></tr></tbody></table>

(Continued)
<table><caption>**Table 2. Continued.**</caption><thead><tr><th>Genus</th><th>Species</th><th>Locality</th><th>Habitat</th><th>Depth (m)</th><th>COI</th><th>ITS2</th><th>Localization of symbionts</th></tr></thead><tbody><tr><td rowspan="12">Nypamodiolus Cosel & Thubaut in prep</td><td><strong>crypta (B' and B")</strong><br>Dall, Bratsch &amp; Rehder, 1938</td><td>WP</td><td>w, b</td><td>431–493</td><td>EU702315<sup>(9)</sup></td><td></td><td>Intracellular<sup>(9)</sup></td></tr><tr><td><strong>sp.1 Broken Bay</strong></td><td>WP</td><td>s</td><td>361–750</td><td>KF611693<sup>(3)</sup></td><td></td><td></td></tr><tr><td><strong>sp.2 Broken Bay</strong></td><td>WP</td><td>s</td><td>361–750</td><td>KF611692<sup>(3)</sup></td><td></td><td></td></tr><tr><td><strong>Aitape G. sp.</strong></td><td>WP</td><td>s</td><td>470</td><td>AB257524<sup>(2)</sup></td><td></td><td></td></tr><tr><td><strong>longissimus*</strong><br>(Thiele & Jaeckel, 1931)</td><td>WP</td><td>w</td><td>400–1767</td><td>DQ340783<sup>(5)</sup></td><td></td><td>Intracellular<sup>(14)</sup></td></tr><tr><td><strong>ESU J</strong></td><td>WP</td><td>w</td><td>360–370</td><td>FJ937189<sup>(4)</sup></td><td></td><td></td></tr><tr><td><strong>ESU I</strong></td><td>WP</td><td>w</td><td>190–567</td><td>FJ937182<sup>(4)</sup></td><td></td><td></td></tr><tr><td><strong>japonica Habe, 1967 (ESU H)</strong></td><td>WP</td><td>w</td><td>150–110</td><td>FJ937060<sup>(4)</sup></td><td></td><td></td></tr><tr><td><strong>simpsoni (Marshall, 1900)</strong></td><td>NA</td><td>b</td><td>162</td><td>KF611695<sup>(3)</sup></td><td></td><td></td></tr><tr><td><strong>manusensis</strong><br>Hashimoto & Furuta 2007</td><td>WP</td><td>v</td><td><strong>1627</strong></td><td>GU966637<sup>(4)</sup></td><td></td><td></td></tr><tr><td><strong>Lau B. sp.</strong></td><td>WP</td><td>v</td><td><strong>1818</strong></td><td>AB257539<sup>(2)</sup></td><td></td><td></td></tr><tr><td><strong>NZ B. sp.</strong></td><td>WP</td><td>v</td><td>unknown</td><td>AB255739<sup>(2)</sup></td><td></td><td></td></tr><tr><td><strong>aduloides Hashimoto &amp; Okutani 1994</strong></td><td>WP</td><td>v</td><td><strong>1378–1451</strong></td><td>AB170054<sup>(8)</sup></td><td></td><td>Intracellular<sup>(18)</sup></td></tr><tr><td>Terua Dall, Bratsch &amp; Rehder, 1938</td><td><strong>pacifica* Dall, Bratsch &amp; Rehder, 1938</strong></td><td>WP</td><td>v</td><td>229</td><td>AB170040<sup>(8)</sup></td><td></td><td>Extracellular<sup>(19)</sup></td></tr><tr><td rowspan="2">Bathymodiolus Kenk & Wilson, 1985</td><td><strong>arcuatilis Dell, 1995</strong></td><td>WP</td><td>v</td><td>880</td><td>FJ937033<sup>(4)</sup></td><td></td><td></td></tr><tr><td><strong>ESU T</strong></td><td>WP</td><td>v</td><td>800–1060</td><td>FJ937275<sup>(4)</sup></td><td></td><td></td></tr><tr><td rowspan="12">Idas Jeffreys, 1876</td><td><strong>thermophilus*</strong><br>Kenk & Wilson, 1985</td><td>EP</td><td>v</td><td><strong>2460–2747</strong></td><td>GU966639<sup>(4)</sup></td><td></td><td>Intracellular<sup>(20)</sup></td></tr><tr><td><strong>azoricus Cosel & Comtet, 1999</strong></td><td>NA</td><td>v</td><td>866–2330</td><td>AY649795<sup>(6)</sup></td><td></td><td>Intracellular<sup>(21)</sup></td></tr><tr><td><strong>puteoserpentis</strong><br>Cosel, Métivier & Hashimoto, 1994</td><td>NA</td><td>v</td><td><strong>3023–3510</strong></td><td>AY649796<sup>(6)</sup></td><td></td><td>Intracellular<sup>(22)</sup></td></tr><tr><td><strong>elongatus Cosel, Métivier & Hashimoto, 1994</strong></td><td>WP</td><td>v</td><td><strong>2765</strong></td><td></td><td>DQ513471<sup>(11)</sup></td><td></td></tr><tr><td><strong>brevior Cosel, Métivier & Hashimoto, 1994</strong></td><td>WP</td><td>v</td><td><strong>3589</strong></td><td>AY649799<sup>(6)</sup></td><td></td><td>Intracellular<sup>(23)</sup></td></tr><tr><td><strong>marisindicus Hashimoto, 2001</strong></td><td>I</td><td>v</td><td><strong>2454</strong></td><td>AB170042<sup>(8)</sup></td><td></td><td>Intracellular<sup>(24)</sup></td></tr><tr><td><strong>septendierum Hashimoto & Okutani, 1994</strong></td><td>WP</td><td>v</td><td><strong>1288</strong></td><td>AB101424<sup>(2)</sup></td><td></td><td>Intracellular<sup>(17)</sup></td></tr><tr><td><strong>boomerang Cosel & Olu 1998</strong></td><td>NA</td><td>v</td><td><strong>1000–3170</strong></td><td>FJ890503<sup>(4)</sup></td><td></td><td>Intracellular<sup>(25)</sup></td></tr><tr><td><strong>heckerae Gustafson, Turner, Lutz & Vrijenhoek, 1998</strong></td><td>GM</td><td>v</td><td><strong>3314</strong></td><td>AY649794<sup>(6)</sup></td><td></td><td>Intracellular<sup>(26)</sup></td></tr><tr><td><strong>brooksi Gustafson, Turner, Lutz & Vrijenhoek, 1998</strong></td><td>GM</td><td>v</td><td><strong>2222–3314</strong></td><td>AY649798<sup>(6)</sup></td><td></td><td>Intracellular<sup>(26)</sup></td></tr><tr><td><strong>argenteus* Jeffreys, 1876</strong></td><td>NA</td><td>w</td><td><strong>1836</strong></td><td></td><td></td><td></td></tr><tr><td><strong>macdonaldi Gustafson, Turner, Lutz & Vrijenhoek, 1998</strong></td><td>GM</td><td>s</td><td>650</td><td>AY649804<sup>(6)</sup></td><td></td><td>Intracellular<sup>(15)</sup></td></tr><tr><td><strong>modiolaeformis (Sturany, 1896)</strong></td><td>M</td><td>s</td><td><strong>2129</strong></td><td>EF210072<sup>(10)</sup></td><td></td><td>Intracellular<sup>(27)</sup></td></tr><tr><td><strong>iwaotakii (A' and A'') (Habe, 1958)</strong></td><td>WP</td><td>w, b</td><td>441–2307</td><td>EU702322<sup>(9)</sup><br>EU702323<sup>(9)</sup></td><td></td><td></td></tr><tr><td><strong>washingtonius (Bernard, 1978)</strong></td><td>EP</td><td>v, w and b</td><td>960–1910</td><td>AY275546<sup>(1)</sup></td><td></td><td>Intracellular<sup>(28)</sup></td></tr></tbody></table>

(Continued)
<table><caption>**Table 2. Continued.**</caption><thead><tr><th>Genus</th><th>Species</th><th>Locality</th><th>Habitat</th><th>Depth (m)</th><th>COI</th><th>ITS2</th><th>Localization of symbionts</th></tr></thead><tbody><tr><td></td><td>sp. SAL1</td><td>WP</td><td>w, b</td><td>408–1356</td><td>FJ937271<sup>(5)</sup></td><td></td><td></td></tr><tr><td></td><td>ESU P</td><td>WP</td><td>w</td><td>180–1390</td><td>FJ937213<sup>(4)</sup></td><td></td><td></td></tr><tr><td></td><td>ESU O</td><td>WP</td><td>w</td><td>473–890</td><td>FJ937206<sup>(4)</sup></td><td></td><td></td></tr><tr><td></td><td>ESU K</td><td>WP</td><td>w</td><td>500–540</td><td>FJ937191<sup>(4)</sup></td><td></td><td></td></tr><tr><td></td><td>sp. D</td><td>WP</td><td>w</td><td>556–1724</td><td>EU702350<sup>(9)</sup></td><td></td><td>Extracellular<sup>(29)</sup></td></tr><tr><td></td><td>sp. C</td><td>WP</td><td>w, b</td><td>275–1285</td><td>EU702360<sup>(9)</sup></td><td></td><td>Extracellular<sup>(9)</sup></td></tr><tr><td></td><td>ESU M</td><td>WP</td><td>w</td><td>590–720</td><td>FJ937202<sup>(4)</sup></td><td></td><td></td></tr><tr><td></td><td>ESU L</td><td>WP</td><td>w</td><td>450–1010</td><td>FJ937193<sup>(4)</sup></td><td></td><td></td></tr><tr><td></td><td>ESU N</td><td>WP</td><td>w, b</td><td>800–1290</td><td>FJ937203<sup>(4)</sup></td><td></td><td></td></tr></tbody></table>

Species/ESUs used in the multigenic analyses are given in bold. Stars indicated the type species of the various genera. Depth scales higher than 1000 m are given in bold. Species underlined are labeled following (19). For locality: EP, Eastern Pacific; GM, Gulf of Mexico; I, Indian Ocean; M, Mediterranean Sea; NA, Northern Atlantic; SA, Southern Atlantic; and WP, Western Pacific. For habitat: w, wood; b, bone; v, vent; and s, seep. Reference for host sequence and symbiont localization: (1) Smith et al. (2004), (2) Miyazaki et al. (2010), (3) this study, (4) Lorion et al. (2010), (5) Samadi et al. (2007), (6) Jones et al. (2006), (7) Miyazaki et al. (2004), (8) Iwasaki et al. (2006), (9) Lorion et al. (2009), (10) Duperron et al. (2008a), (11) Olu-Le Roy et al. (2007), (12) McKiness et al. (2005), (13) von Cosel and Marshall (2010), (14) Duperron et al. (2008b), (15) Won et al. (2008), (16) Childress et al. (1986), (17) Fujiwara et al. (2000), (18) Yamanaka et al. (2000), (19) Fujiwara et al. (2010), (20) Distel et al. (1988), (21) Fiala-Medioni et al. (2002), (22) Cavanaugh et al. (1992), (23) Dubilier et al. (1998), (24) Yamanaka et al. (2003), (25) Duperron et al. (2005), (26) Duperron et al. (2007), (27) Olu-Le Roy et al. (2004), (28) Deming et al. (1997), and (29) Duperron et al. (2009). COI, cytochrome oxydase I.

"depth" was coded according to Craddock et al. (1995) and Jones et al. (2006) in taking the shallowest recorded depth lesser ("0") and greater than 1000 m ("1"). Ancestral states were reconstructed for all Bayesian trees retained from the analysis of the combined data set and their mean likelihood was then plotted on the maximum clade credibility tree.

## Results

### Phylogenetic relationships among mytilids associated with deep-sea reducing habitats and genus-level classification

Maximum likelihood and BI of each of the individual gene fragments provided similar patterns of relationships, with no incongruence, but poorly resolved phylogenetic trees, whereas the multigene alignment of 51 taxa/ESUs and 5912 bp produced 10 robustly supported clades (Fig. 1).

A first monophyletic lineage (99.7% of bootstraps, posterior probabilities (PP) = 1) was sister group to all the other mussels from chemosynthetic habitats. This lineage included *Be. lignocola* but also two additional undescribed species (sp. Juan de Fuca and sp. South Atlantic). As *Be. lignocola* is the type species for *Benthomodiolus*, both sp. Juan de Fuca and sp. South Atlantic were attributed to this genus.

All the other lineages were included in a well-supported clade (95.2% of bootstraps, PP = 1). The position of

type species was used to name six lineages: *Vulcanidas* Cosel & Marshall 2010, *Tamu* Gustafson, Turner, Lutz & Vrijenhoek 1998, *Lignomodiolus* Cosel & Thubaut in prep., *Nypamodiolus* Cosel & Thubaut in prep., *Gigantidas* Cosel & Marshall 2003, and *Bathymodiolus* Kenk & Wilson 1985 (Fig. 1). The species *Terua pacifica*, type species of the genus *Terua*, clustered with the species *arcuatilis* and the ESU T based on the COI and ND4 genes (as revealed by the DNA barcoding approach – see below), and the name *Terua* was thus attributed to this clade. The sister clade of *Bathymodiolus* mainly included species generally attributed to *Idas*, and we applied this name to the clade. For only one lineage no genus name could be applied. The species "B." *manusensis* is alone and we refrain to describe it as a new genus pending for additional species to be recognized as belonging to this lineage.

We also recovered the lineage, first, highlighted by Jones et al. (2006) and labeled as the "childressi" group: it fell here within the *Gigantidas* lineage. The two morphospecies collected off the Bismarck Sea in the Broken bay (Papua New Guinea) also clustered within the *Gigantidas* lineage: G. sp.1 Broken Bay was robustly included within the "childressi" group (Jones et al. 2006) as the sister species of G. *taiwanensis*, and G. sp.2 Broken Bay was the sister species of G. *gladius*. Finally, the specimens attributed to *Idas simpsoni* from the northern Atlantic clustered in the *Nypamodiolus* lineage and thus this species must be included in this genus.
<figure><img src="image_2.png" /><figcaption>**Figure 1.** Phylogenetic relationships among Bathymodiolinae obtained from Bayesian analysis of the multigene data set (two mitochondrial and five nuclear genes), including species from hydrothermal vents, cold seeps, and organic falls available on GenBank and new samples. Posterior probabilities (PP) and bootstraps values obtained from maximum likelihood (ML) analysis are given above and below nodes, respectively. PP and bootstraps values lower than 0.90 and 50%, respectively, are not shown. Scale bar represents 0.1% estimated base substitution. Environments and depth range inhabited by each species are given. *B.*, *Bathymodiolus*.</figcaption></figure>
<figure><img src="image_3.png" /></figure>

**Figure 2.** Bayesian tree obtained from the analysis of the Cytochrome Oxydase I (COI) mt DNA data set (see Table 2). Only PP (right) higher than 0.50 and bootstraps values (left) obtained from ML analysis higher than 50% are given. Scale bar represents 0.1% estimated base substitution.

## Species diversity within genera

Six species identified in the literature by molecular data were attributed to either a named species or an ESU of our data set using the 2% K2P threshold: NZ *B*. sp. and Lau *B*. sp. both fall with "*B.*" *manusensis*, Nikko *G*. sp. and Sumisu *G*. sp. both fall with *G. horikoshii*, Sissano *B*. sp. 3 falls with *G. tangaroa*, and finally Aitape *G*. sp. falls with *Gigantidas* sp.2 Broken Bay (see Table 2 and Fig. 2). Based on mitochondrial data and using the same threshold, Van Dover et al. (2001) suggested that *B. brevior* and *B. marisindicus* might be conspecific. More recent studies shown that these two species and *B. septemdierum* are closely related and are connected by gene flows (see Kyuno et al. 2009 and Miyazaki et al. 2010). Thus, considering all the data gathered for this three putative species, we agree with Vrijenhoek (2009) that they might be conspecific. The name *B. septemdierum* is the oldest name and has the nomenclatural priority over the two others. Similarly, *B. boomerang* and *B. heckerae* in one hand and *G. childressi* and *G. mauritanicus* in the other hand belong to two complexes of very closely related species that although geographically distant might not be reproductively isolated (Olu-Le Roy et al. 2007) and could thus be conspecific. A similar pattern of weak structure over large geographic has been detected for the small sunken wood mussels (Lorion et al. 2010), but contrary to these studies the pairs of populations geographically very distant but genetically weakly distant have been considered as unique ESUs.

The attribution based on the distance threshold was confirmed for most of the added species by the phylogenetic analysis of the extended COI data set (Fig. 2). The only three exceptions are Chamorro *B*. sp., Kikaijima *B*. sp., and Sissano *B*. sp. 1 for which the sister taxa could not be determined using the COI data set. These three taxa are, however, robustly placed in the "childressi" group by Miyazaki et al. (2010), using a combined phylogenetic analysis with two mitochondrial (COI and ND4 genes), a lineage that branched in *Gigantidas* in our multigenic phylogeny. Last, *B. elongatus* was maintained in *Bathymodiulus* because it falls close to *B. brevior* in the ITS2 data set of Olu-Le Roy et al. (2007). This species might, however, also be conspecific with *B. brevior* considering the very small genetic distance with species for this genetic marker.

Last, for four species no molecular data were available. First, no data are available for *I. argenteus*, but as it is the type species of the genus, it is de facto in this genus. Sec-
<figure><img src="image_4.png" /><figcaption>**Figure 3.** Bayesian cladogram of the combined data set with maximum likelihood estimates of ancestral character state for (A) habitat and (B) depth. Pie charts correspond to average likelihoods of each state. Numeric values of some nodes of interest are given.</figcaption></figure>

ond, *Benthomodiolus abyssicola* was maintained in this genus because contrary to other lineages the shells appeared to be discriminant at the genus level (R. von Cosel, pers. inform.). Indeed, a new taxon added in this study was successfully attributed to this genus using the shell morphology, and the morphological data discussed in Lorion et al. (2010) also suggest that this lineage is morphologically diagnosable. Last, von Cosel and Janssen (2008) described two new species as belonging to the “childressi” group that we also tentatively maintained in *Gigantidas*.

## Character evolution

The analysis confirmed that organic remains were the ancestral habitat for *Bathymodiolinae* (Fig. 3A). Furthermore, the colonization from organic falls to hydrothermal vents and/or cold seeps occurred at least in four lineages: *Benthomodiolus* (*Be.* sp. Juan de Fuca and *Be.* *lignocola* vs. *Be.* sp. South Atlantic), *Vulcanidas* (*V.* *insolatus* vs. the three unidentified species from sunken wood), *Gigantidas* (the “childressi” group, *G.* *gladius*, and *G.* sp. Broken Bay vs. the *G.* *crypta* species complex), and at the divergence

between *Bathymodiolus* and *Idas*. We could not conclude if the ancestral habitat was shallow or deep for bathymodiolin mussels (ancestral state: 0 = 0.53 and 1 = 0.47; Fig. 3B).

## Discussion

Including more taxa and more genes, we significantly improved the resolution of phylogenetic relationships among deep-sea mussels, both for many of the deeper nodes (among lineages) and for some of the internal nodes (within lineages). Most of the genera, as currently defined, are polyphyletic and/or paraphyletic and thus the current classification does not reflect evolutionary relationships. The existence of several lineages within *Bathymodiolinae* from seeps and vents was earlier hypothesized based on anatomical characters. For example, in the original description of *G. childressi* (Gustafson et al. 1998), the unique combination of morphological characters was underlined but a new genus was not erected because the genetic distances were appraised too low to erect a new genus. Similarly, von Cosel and Olu (1998) doubted that
the newly described species *B. boomerang* belong to the same lineage than *B. thermophilus* and suggested that part of the shared features (such as a straight intestine) were independently acquired. The first robust molecular studies (Jones and Vrijenhoek 2006) showed that *G. childressi* indeed belong to a distinct lineage, but that *B. boomerang* belongs to the “thermophilus” lineage. Although the classification used to date, although debated (von Cosel 2002; Génio et al. 2012), does not reflect evolutionary relationships and this has hindered comparative biology studies of deep-sea mussels. For example, Schultz et al. (2011) failed to efficiently cross-amplify microsatellite loci developed for “*B.*” *manusensis* in another species placed in the same genus (*B. heckerae*). We here show that these two species belong to two distant lineages and should be now classified in two distinct genera. The low cross-amplification success is thus explained by the fact that microsatellites only poorly cross-amplify among distant species and should be expected to be successfully transferred only to closely related species (Sharma et al. 2007).

Revising the taxonomy at the genus level, based on robust phylogenetic results, aims at guiding such comparative studies, and to challenge a priori assumptions about the likelihood of “evolutionary steps.” For example, the fact that embryos of the shallow water mussel *Mytilus edulis* can develop within a range of temperature and pressure encountered at vents is used by Mestre et al. (2009) as an argument to reject the “wooden-step” hypothesis. Mestre et al. (2009) assumed that an evolutionary step from eurythermal environments to cold stenothermal environments is *generally assumed highly unlikely*. However, the phylogenetic tree clearly shows that the evolutionary history of mussels from cold stenothermal chemosynthetic environments (i.e., both organic remains and cold seeps) and from the eurythermal environment found around hot deep-sea vents is tightly entangled. Moreover, our phylogenetic reconstruction shows that such evolutionary steps occurred more than once and are thus not so unlikely.

Based on this new phylogeny, organic remains are robustly supported as the ancestral habitat of Bathymodiolinae mussels (Lorion et al. 2010). Indeed, most of internal nodes of each genus are robustly assigned to the habitat “organic falls.” Only the ancestor of the *Bathymodiolus* clade is not statistically resolved. Moreover, whereas most habitat shifts occurred from bone and wood falls to vent and seep sites and a few occurred from seeps to vents, only one transition is inferred from vent to seep (ancestor of the sister-species *G. taiwanensis*/*G. sp.1 Broken Bay*). This result supports the “wooden steps” hypothesis proposed by Distel et al. (2000), in which the ancestors of seep and vent mussels were associated with organic falls. However, our analysis provided evidences

that recurrent colonization events from organic remains to the harsh deep-sea seep and vent habitats occurred during the evolution of the Bathymodiolinae: a “multiple wooden steps” hypothesis would thus be more adequate. If some Bathymodiolinae species are able to inhabit both seeps and vents (e.g., *G. platifrons* and *G. japonicus* [Miyazaki et al. 2010; ]) or both seeps and sunken wood (Thubaut 2012), there is no species known to be able to colonize both organic remains and hydrothermal vents. Cold seeps seem to be an intermediate habitat for the colonization of deep vents, and a switch back from harsh vent habitats to organic falls is never inferred with sampling to date.

This phylogeny also allows a reevaluation of the evolutionary significance of some striking characters of the vent and seep mussels. For example, the repeated evolutionary shifts to vents and seeps are always associated with increase in mussels’ size. The large size seems to be correlated with the habitat and thus could be a convergent evolution in independent lineages. Experimental results also suggest an important role of environmental parameters in the growth of Bathymodiolinae. For example, the size, the body condition, and the growth rate of *G. childressi* is correlated with the concentration of methane (Nix et al. 1995). These results were later confirmed by transplant experiments (Bergquist et al. 2004) in which mussels acquired nearly the characteristics of their host population. However, within the genus *Idas*, all the species either from organic falls or seeps and vents are very small (smaller than 1.5 cm c.f. [Lorion et al. 2010]). The polymorphism of some of the shells described from some species associated with organic falls (i.e., *V. ESU E, L. ESU S*, or *N. longissimus*) suggests that growth rate may also vary with environmental parameters. Additional experimental results are, however, needed to determine to what extent the size and/or growth rate of the small mussels from organic remains, belonging to each evolutionary lineage, may vary with the concentration of reduced compounds (such as hydrogen sulfide). The correlation of gigantism of mussels with seep and vent environments may, at least in part, explain why most of them have been placed in the genus *Bathymodiolus* and have received more attention from marine scientists than the small mussels from organic remains.

The inclusion of more data from various depths and environments shows that the evolution from shallow to deep (Craddock et al. 1995; Jones et al. 2006) is no more supported than the opposite hypothesis and that more samples from the deepest environments are needed to determine whether there is an evolutionary trend linked to the depth ranges. Data for organic falls from abyssal environments remain very scarce due to the difficulty of localizing them at the deep-sea floor. The only mean to obtain
such data is to deploy experimental organic remains in abyssal environments (Voight and Segonzac 2012).

The phylogenetic interminglement of mussels from organic remains with mussels from seeps and vents challenges the scenarios drawn from previous phylogenetic analyses that were mainly based on a small and habitat-biased taxon sampling. For example, Fujiwara et al. (2010) and Miyazaki et al. (2010) suggest an evolutionary scenario in which seep and vent mussels were considered as a monophyletic lineage characterized by chemosynthetic endosymbionts. Based on the basal position of the few taxa associated with organic remains that have ectosymbionts, Fujiwara et al. (2010) suggested that extracellular symbiosis is an “earlier stage” in the evolution of Bathymodiolin than intracellular symbiosis. Although only little data are available on the position of the symbionts in the gill tissues for most of the species (see Table 2) (preventing any robust reconstruction of the ancestral states), our phylogeny shows that (1) several lineages include both types of symbionts and that (2) the symbionts position is not always related to habitats. For example, in the genus *Vulcanidas*, *V. insolatus* is suggested to harbor endosymbiotic bacteria (von Cosel and Marshall 2010), whereas ectosymbiotic bacteria were observed in the gill tissues of *Vulcanidas* sp. E sampled on sunken wood (Duperron et al. 2008b). Conversely, within *Giganitidas*, all species, either from vent and seeps or organic remains, have endosymbionts (Table 2). Finally, the two sister genera *Bathymodiolus* and *Idas*, as redefined here, display contrasted pattern of symbionts positions. Indeed, all *Bathymodiolus* species, from hot vents or cold seeps, have endosymbionts (reviewed in Duperron et al. [2009]), whereas *Idas* species, mainly associated with organic remains – but also some from vents or seeps – have either endosymbionts (Baco-Taylor 2002; Baco and Smith 2003; Won et al. 2008) or ectosymbionts (Duperron et al. 2009). Although data are lacking for many species, our phylogenetic results suggest that endosymbiosis has been acquired more than once. In addition, the endosymbiotic bacteria observed in the gills of *G. crypta*, a species sampled either on sunken wood or bone remains, show that endosymbiosis should not be considered strictly as a specialization to vent or seep habitats.

Hydrothermal vents have been thoroughly explored in the past few decades and only a few new species are still being found, mainly from poorly explored biogeographic areas. The most recently described mussel species from vents that putatively fall into *Bathymodiolus* is *B. marisindicus* from Indian Ocean, an area that remains a yet underexplored biogeographic province even for vent habitats. However, molecular data suggested that *B. brevior* and *B. marisindicus* (Van Dover et al. 2001), but also *B. elongatus* (Olu-Le Roy et al. 2007), may be conspecific.

The other recently described vent species fall putatively in *Gigantidas* (Hashimoto and Yamane 2005), “*Bathymodiolus*” (e.g., “*B.*”. *manusensis*), or in the new genus *Vulcanidas* (von Cosel and Marshall 2010). However, in these genera a lot more new species from seeps have been recently described (von Cosel and Janssen 2008) or have yet to be formally described. Similarly, the recent interest for organic falls added a bulk of species that are mostly not named. The sampling of deep-sea mussels from organic remains is far from being saturated (Lorion and Samadi 2010). Indeed, organic falls are still very poorly explored and we hypothesize that much more new species will be discovered from organic remains, but also from hydrocarbon seeps, than from vents.

The genus *Idas* is not only speciose but also ecologically diverse. It includes *I. modiolaeformis*, the species harboring the highest diversity of symbionts (Duperron et al. 2008a). Within this genus, species display either extra- or intracellular symbionts (Deming et al. 1997; Won et al. 2008; Duperron et al. 2009). The genus is found over wide geographic (Pacific, Atlantic, and Mediterranean) and bathymetric (180–2307 m) ranges and present in all kind of habitats. Generally, the population densities are low, although it may be locally abundant. For the relatively well-sampled species *I. iwaotakii*, the depth range within species covers more than 2000 m and the geographic range spreads over the Western Pacific. In contrast, the genus *Bathymodiolus*, sister lineage of *Idas*, includes species recovered in dense populations mainly from vents and a few from seeps. All species harbor intracellular symbionts and are restricted to depths greater than 800 m, with most of the species apparently restricted to a limited depth range (Fig. 1 and Table 2). Similarly, the “*Bathymodiolus*” lineage includes only two species (“*B.*” *manusensis* and “*B.*” *aduloides*). This lineage appears ecologically restricted: specimens were collected solely at vents or seeps, and only at sites below 1300 m (Iwasaki et al. 2006; Fujita et al. 2009; Kyuno et al. 2009; Miyazaki et al. 2010). In contrast, *Gigantidas* is well diversified, with a large number of species awaiting a formal description, and covers all kind of habitats (wood, bone, seep, and vent) over a quite large depth range (200–2222 m). Finally, *Benthomodiolus*, the sister genus of all other Bathymodiolinae, is still only known by a few species sampled mainly at very deep sites, both in the Atlantic and in the Pacific and both associated with organic falls and vents. Although no molecular data are available, at least another species (*Benthomodiolus abyssicolae*) is mentioned in the Gulf of Panama at 3197 meters attached to wood (Wolff 1979). Wood falls below 1500 m have only been poorly explored (Samadi et al. 2010) and we can guess that exploration of deeper sites will bring more species in this lineage.
Overall, Bathymodiolinae lineages in which mussels are mainly associated with organic remains are more diversified and cover a larger range of ecological conditions, whereas lineages that mainly include vents species are less diverse and ecologically more restricted, but with very dense populations. This trend should be reinforced in the future as organic falls, particularly at depths below 1500 m, are still very poorly explored. Except for the *Bathymodiolus* lineage, the inferred ancestral habitat for the lineages was never the vent habitat. Vent species are ecologically successful, but, because of their low species richness and the apparent difficulty to switch back to other habitats, their evolutionary fate may be compromised. As a consequence, taking into account the species diversity associated with organic falls, we suggest that for mussels, adaptation to vents, although ecologically very efficient, is likely an evolutionary dead end.

## Acknowledgments

New material was collected on board of R/V Alis during the BioPapua 2010 cruise (Principal investigators: S. Samadi and L. Corbari). Sincere thanks to B. Buge and J. Brisset for curation of specimens and the staff of the “Service de Systématique Moléculaire” (UMS2700 CNRS-MNHN). It is part of the agreement no. 2005/67 between the Genoscope and the Museum National d’Histoire Naturelle on the project “Macrophylogeny of life” directed by G. Lecointre. We thank Philippe Bouchet and Patrice David for constructive comments, and Stéphane Hourdez and Eric Pante for discussions and English improvement of this manuscript. This work was supported in part by the grant TF-DeepEvo ANR-12-ISV7-0005, and by the “Action Transversale du MNHN : Taxonomie moléculaire : DNA Barcode et gestion durable des collections”.

## Conflict of Interest

None declared.

## References

Audzijonyte, A., and R. C. Vrijenhoek. 2009. Three nuclear genes for phylogenetic, SNP and population genetic studies of molluscs and other invertebrates. Mol. Ecol. Resour. 10:200–204.

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

Baco-Taylor, A. R. 2002. Food-web structure, succession, and phylogenetics on deep-sea whale skeletons. p. 275 in *Oceanography. The University of Hawaii, Manoa, HI.* [Ph.D. dissertation].

Bergquist, D. C., C. Fleckenstein, E. B. Szalai, J. Knisel, and C. R. Fisher. 2004. Environment drives physiological variability in the cold seep mussel *Bathymodiolus childressi*. Limnol. Oceanogr. 49:706–715.

Bettencourt, R., M. Pinheiro, C. Egas, P. Gomes, M. Afonso, T. Shank, et al. 2010. High-throughput sequencing and analysis of the gill tissue transcriptome from the deep-sea hydrothermal vent mussel *Bathymodiolus azoricus*. BMC Genomics 11:559.

Cavanaugh, C. M., C. O. Wirsen, and H. W. Jannasch. 1992. Evidence for methylotrophic symbionts in a hydrothermal vent mussel (Bivalvia: Mytilidae) from the Mid-Atlantic Ridge. Appl. Environ. Microbiol. 58:3799–3803.

Chapelle, G., and L. S. Peck. 1999. Polar gigantism dictated by oxygen availability. Nature 399:114–115.

Childress, J. J., C. R. Fisher, J. M. Brooks, M. C. Kennicutt, R. Bidigare, and A. E. Anderson. 1986. A methanotrophic marine molluscan (Bivalvia, Mytilidae) symbiosis: mussels fueled by gas. Science 233:1306–1308.

von Cosel, R. 2002. A new species of bathymodioline mussel (Mollusca, Bivalvia, Mytilidae) from Mauritania (West Africa), with comments on the genus *Bathymodiolus* Kenk & Wilson, 1985. Zoosystema 24:259–272.

von Cosel, R., and R. Janssen. 2008. Bathymodioline mussels of the *Bathymodiolus* (s. l.) *childressi* clade from methane seeps near Edison Seamount, New Ireland, Papua New Guinea: (Bivalvia: Mytilidae). Arch. Molluskenkd. 137:195–224.

von Cosel, R., and B. A. Marshall. 2010. A new genus and species of large mussel (Mollusca: Bivalvia: Mytilidae) from the Kermadec Ridge. Nautilus 117:31–46.

von Cosel, R., and K. Olu. 1998. Gigantism in Mytilidae. A new *Bathymodiolus* from cold seep areas on the Barbados accretionary prism. C. R. Acad. Sci. Ser. III Sci. Vie 321:655–663.

Costello, M. J., M. Coll, R. Danovaro, P. Halpin, H. Ojaveer, and P. Miloslavich. 2010. A census of marine biodiversity knowledge, resources, and future challenges. PLoS One 5: e12110.

Craddock, C., W. R. Hoeh, R. G. Gustafson, R. A. Lutz, J. Hashimoto, and R. J. Vrijenhoek. 1995. Evolutionary relationships among deep-sea mytilids (Bivalvia: Mytilidae) from hydrothermal vents and cold-water methane/sulfide seeps. Mar. Biol. 121:477–485.

Deming, J. W., A. L. Reysenbach, S. A. Macko, and C. R. Smith. 1997. Evidence for the microbial basis of a chemoautotrophic invertebrate community at a whale fall on the deep seafloor: bone-colonizing bacteria and invertebrate endosymbionts. Microsc. Res. Tech. 37:162–170.

Distel, D. L. 2000. Phylogenetic relationships among mytilidae (Bivalvia): 18S rRNA data suggest convergence in mytilid body plans. Mol. Phylogenet. Evol. 15:25–33.
Distel, D. L., D. J. Lane, G. J. Olsen, S. J. Giovannoni, B. Pace, N. R. Pace, et al. 1988. Sulfur-oxidizing bacterial endosymbionts: analysis of phylogeny and specificity by 16S rRNA sequences. J. Bacteriol. 170:2506–2510.

Distel, D. L., A. R. Baco, E. Chuang, W. Morrill, C. Cavanaugh, and C. R. Smith. 2000. Do mussels take wooden steps to deep-sea vents? Nature 403: 725–726.

Dixon, D. R., A. M. Pruski, and L. R. J. Dixon. 2004. The effects of hydrostatic pressure change on DNA integrity in the hydrothermal-vent mussel *Bathymodiolus azoricus*: implications for future deep-sea mutagenicity studies. Mutat. Res. 552:235–246.

Dubilier, N., R. Windoffer, and O. Giere. 1998. Ultrastructure and stable carbon isotope composition of the hydrothermal vent mussels *Bathymodiolus brevior* and *B. sp. affinis brevior* from the North Fiji Basin, western Pacific. Mar. Ecol. Prog. Ser. 165:187–193.

Duperron, S., T. Nadalig, J.-C. Caprais, M. Sibuet, A. Fiala-Medioni, R. Amann, et al. 2005. Dual symbiosis in a *Bathymodiolus* sp. mussel from a methane seep on the Gabon continental margin (Southeast Atlantic): 16S rRNA phylogeny and distribution of the symbionts in gills. Appl. Environ. Microbiol. 71:1694–1700.

Duperron, S., M. Sibuet, B. J. MacGregor, M. M. M. Kuypers, C. R. Fisher, and N. Dubilier. 2007. Diversity, relative abundance and metabolic potential of bacterial endosymbionts in three *Bathymodiolus* mussel species from cold seeps in the Gulf of Mexico. Environ. Microbiol. 9:1423–1438.

Duperron, S., S. Halary, J. Lorion, M. Sibuet, and F. Gaill. 2008a. Unexpected co-occurrence of six bacterial symbionts in the gills of the cold seep mussel *Idas* sp. (Bivalvia: Mytilidae). Environ. Microbiol. 10:433–445.

Duperron, S., M. C. Z. Laurent, F. Gaill, and O. Gros. 2008b. Sulphur-oxidizing extracellular bacteria in the gills of Mytilidae associated with wood falls. FEMS Microbiol. Ecol. 63:338–349.

Duperron, S., J. Lorion, S. Samadi, O. Gros, and F. Gaill. 2009. Symbioses between deep-sea mussels (Mytilidae: Bathymodiolinae) and chemosynthetic bacteria: diversity, function and evolution. C. R. Biol. 332:298–310.

Felsenstein, J. 1985. Confidence limits on phylogenies: an approach using the bootstrap. Evolution 39:783–791.

Fiala-Medioni, A., Z. McKiness, P. Dando, J. Boulegue, A. Mariotti, A. Alayse-Danet, et al. 2002. Ultrastructural, biochemical, and immunological characterization of two populations of the mytilid mussel *Bathymodiolus azoricus* from the Mid-Atlantic Ridge: evidence for a dual symbiosis. Mar. Biol. 141:1035–1043.

Folmer, O., M. Black, W. Hoeh, R. Lutz, and R. Vrijenhoek. 1994. DNA primers for amplification of mitochondrial cytochrome c oxidase subunit I from diverse metazoan invertebrates. Mol. Mar. Biol. Biotech. 3:294–299.

Fontanez, K. M., and C. M. Cavanaugh. 2013. Phylogenetic relationships of hydrothermal vent mussels (*Bathymodiolinae*) and their symbionts. Mar. Ecol. Prog. Ser. 474:147–154.

Fujita, Y., H. Matsumoto, Y. Fujiwara, J. Hashimoto, S. V. Galkin, R. Ueshima, et al. 2009. Phylogenetic relationships of deep-sea *Bathymodiolus* mussels to their mytilid relatives from sunken whale carcasses and wood. Venus 67:123–134.

Fujiwara, Y., K. Takai, K. Uematsu, S. Tsuchida, J. C. Hunt, and J. Hashimoto. 2000. Phylogenetic characterization of endosymbionts in three hydrothermal vent mussels: influence on host distributions. Mar. Ecol. Prog. Ser. 208:147–155.

Fujiwara, Y., M. Kawato, C. Noda, G. Kinoshita, T. Yamanaka, Y. Fujita, et al. 2010. Extracellular and mixotrophic symbiosis in the whale-fall mussel *Adipicola pacifica*: a trend in evolution from extra- to intracellular symbiosis. PLoS One 5:e11808.

Génio, L., S. Kiel, M. R. Cunha, J. Grahame, and C. T. S. Little. 2012. Shell microstructures of mussels (Bivalvia: Mytilidae: Bathymodiolinae) from deep-sea chemosynthetic sites: do they have a phylogenetic significance? Deep-Sea Res. I 64:86–103.

Giribet, G., S. Carranza, J. Baguna, M. Riutort, and C. Ribera. 1996. First molecular evidence for the existence of a Tardigrada+ Arthropoda clade. Mol. Biol. Evol. 13:76–84.

Gould, S. J., and R. C. Lewontin. 1979. The spandrels of San Marco and the Panglossian paradigm: a critique of the adaptationist programme. Proc. R. Soc. Lond. B Biol. Sci. 205:581–598.

Gustafson, R. G., R. D. Turner, R. A. Lutz, and R. C. Vrijenhoek. 1998. A new genus and five new species of mussels (Bivalvia, Mytilidae) from deep-sea sulfide/ hydrocarbon seeps in the Gulf of Mexico. Malacologia 40:63–112.

Hardivillier, Y., V. Leignel, F. Denis, G. Uguen, R. Cosson, and M. Laulier. 2004. Do organisms living around hydrothermal vent sites contain specific metallothioneins? The case of the genus *Bathymodiolus* (Bivalvia, Mytilidae). Comp. Biochem. Physiol. C Toxicol. Pharmcol. 139:111–118.

Hashimoto, J., and T. Yamane. 2005. A new species of *Gigantidas* (Bivalvia: Mytilidae) from a vent site on the Kaikata Seamount southwest of the Ogasawara (Bonin) Islands, southern Japan. Venus 64:1–10.

Hassouna, N., B. Mithot, and J. P. Bachellerie. 1984. The complete nucleotide sequence of mouse 28S rRNA gene. Implications for the process of size increase of the large subunit rRNA in higher eukaryotes. Nucleic Acids Res. 12:3563–3583.

Heath, T. A., D. J. Zwickl, J. Kim, and D. M. Hillis. 2008. Taxon sampling affects inferences of macroevolutionary processes from phylogenetic trees. Syst. Biol. 57:160–166.

Iwasaki, H., A. Kyuno, M. Shintaku, Y. Fujita, Y. Fujiwara, K. Fujikura, et al. 2006. Evolutionary relationships of deep-sea mussels inferred by mitochondrial DNA sequences. Mar. Biol. 149:1111–1122.
Jannasch, H. W. 1985. The chemosynthetic support of life and the microbial diversity at deep-sea hydrothermal vents. *Proc. R. Soc. Lond. B Biol. Sci.* 225:277–297.

Jobb, G., A. von Haeseler, and K. Strimmer. 2004. TREEFINDER: a powerful graphical analysis environment for molecular phylogenetics. BMC Evol. Biol. 4:18.

Jones, W. J., and R. C. Vrijenhoek. 2006. Evolutionary relationships within the "Bathymodiolus" childressi group. Cah. Biol. Mar. 47:403–407.

Jones, W. J., Y. J. Won, P. A. Y. Maas, P. J. Smith, R. A. Lutz, and R. C. Vrijenhoek. 2006. Evolution of habitat use by deep-sea mussels. Mar. Biol. 148:841–851.

Kadar, E., and J. Powell. 2006. Post-capture investigations of hydrothermal vent macro-invertebrates to study adaptations to extreme environments. Rev. Environ. Sci. Biotechnol. 5:193–201.

Kádár, E., R. Bettencourt, V. Costa, R. S. Santos, A. Lobo-da-Cunha, and P. Dando. 2005. Experimentally induced endosymbiont loss and re-acquirement in the hydrothermal vent bivalve *Bathymodiolus azoricus*. J. Exp. Mar. Biol. Ecol. 318:99–110.

Klok, C. J. 2009. O₂ levels may not determine gigantism. J. Exp. Biol. 212:iv.

Kyuno, A., M. Shintaku, Y. Fujita, H. Matsumoto, M. Utsumi, H. Watanabe, et al. 2009. Dispersal and differentiation of deep-sea mussels of the genus *Bathymodiolus* (Mytilidae, Bathymodiolinae). J. Mar. Biol. 2009:1–15.

Lallier, F. H. 2006. Thioautotrophic symbiosis: towards a new step in eukaryote evolution. Cah. Biol. Mar. 47:391–396.

Lorion, J., and S. Samadi. 2010. Species richness, sampling bias and phylogenetics in deep-sea mussels. Cah. Biol. Mar. 51:435–539.

Lorion, J., S. Duperron, O. Gros, C. Cruaud, and S. Samadi. 2009. Several deep-sea mussels and their associated symbionts are able to live both on wood and on whale falls. Proc. Biol. Sci. 276:177–185.

Lorion, J., B. Buge, C. Cruaud, and S. Samadi. 2010. New insights into diversity and evolution of deep-sea Mytilidae (Mollusca: Bivalvia). Mol. Phylogenet. Evol. 57:71–83.

Lorion, J., S. Halary, and J. do Nascimento, S. Samadi, A. Couloux, and S. Duperron. 2012. Evolutionary history of Idas sp. Med, (Bivalvia: Mytilidae), a cold seep mussel bearing multiple symbionts. Cah. Biol. Mar. 53:77–87.

Maddison, W. P., and D. R. Maddison. 2011. Mesquite: a modular system for evolutionary analysis. Version 2.75. Available at http://mesquiteproject.org website (accessed October 2012).

McKiness, Z. P., E. R. McMullin, C. R. Fisher, and C. M. Cavanaugh. 2005. A new bathymodioline mussel symbiosis at the Juan de Fuca hydrothermal vents. Mar. Biol. 148:109–116.

Mestre, N. C., S. Thatje, and P. A. Tyler. 2009. The ocean is not deep enough: pressure tolerances during early ontogeny

of the blue mussel *Mytilus edulis*. Proc. Biol. Sci. 276:717–726.

Miyazaki, J. I., M. Shintaku, A. Kyuno, Y. Fujiwara, J. Hashimoto, and H. Iwasaki. 2004. Phylogenetic relationships of deep-sea mussels of the genus *Bathymodiolus* (Bivalvia: Mytilidae). Mar. Biol. 144:527–535.

Miyazaki, J. I., L. de Oliveira Martins, Y. Fujita, H. Matsumoto, and Y. Fujiwara. 2010. Evolutionary process of deep-sea *Bathymodiolus* mussels. PLoS One 5:e10363.

Nix, E. R., C. R. Fisher, J. Vodenichar, and K. M. Scott. 1995. Physiological ecology of a mussel with methanotrophic endosymbionts at three hydrocarbon seep sites in the Gulf of Mexico. Mar. Biol. 122:605–617.

Nylander, J. A. A., J. C. Wilgenbusch, D. L. Warren, and D. L. Swofford. 2008. AWTY (are we there yet?): a system for graphical exploration of MCMC convergence in Bayesian phylogenetics. Bioinformatics 24:581–583.

Okusu, A., E. Schwabe, D. J. Eernisse, and G. Giribet. 2003. Towards a phylogeny of chitons (Mollusca, Polyplacophora) based on combined analysis of five molecular loci. Org. Divers. Evol. 3:281–302.

Olu-Le Roy, K., M. Sibuet, A. Fiala-Médioni, S. Gofas, C. Salas, A. Mariotti, et al. 2004. Cold seep communities in the deep eastern Mediterranean Sea: composition, symbiosis and spatial distribution on mud volcanoes. Deep-Sea Res. I 51:1915–1936.

Olu-Le Roy, K., R. von Cosel, S. Hourdez, S. L. Carney, and D. Jollivet. 2007. Amphi-Atlantic cold-seep *Bathymodiolus* species complexes across the equatorial belt. Deep-Sea Res. I 54:1890–1911.

Pante, E., L. Corbari, J. Thubaut, T. Y. Chan, R. Mana, M. C. Boisselier, et al. 2012. Exploration of the deep-sea fauna of Papua New Guinea. Oceanography 25:214–225.

Posada, D., and K. A. Crandall. 1998. Modeltest: testing the model of DNA substitution. Bioinformatics 14:817–818.

Pruski, A. M., and D. R. Dixon. 2007. Heat shock protein expression pattern (HSP70) in the hydrothermal vent mussel *Bathymodiolus azoricus*. Mar. Environ. Res. 64:209–224.

Rambaud, A., and A. J. Drummond. 2007. Tracer, version 1.4. Available at http://beast.bio.ed.ac.uk/Tracer (accessed December 2012).

Ronquist, F., and J. P. Huelsenbeck. 2003. MrBayes 3: Bayesian phylogenetic inference under mixed models. Bioinformatics 19:1572–1574.

Samadi, S., E. Quemere, J. Lorion, A. Tillier, R. von Cosel, P. Lopez, et al. 2007. Molecular phylogeny in mytilids supports the wooden steps to deep-sea vents hypothesis. C. R. Biol. 330:446–456.

Samadi, S., L. Corbari, J. Lorion, S. Hourdez, T. Haga, J. Dupont, et al. 2010. Biodiversity of deep-sea organisms associated with sunken-wood or other organic remains sampled in the tropical Indo-Pacific. Cah. Biol. Mar. 51:459–466.
Schultz, T. F., P. Y. Hsing, A. Eng, K. A. Zelnio, A. D. Thaler, J. Carlsson, et al. 2011. Characterization of 18 polymorphic microsatellite loci from *Bathymodiolus manusensis* (Bivalvia, Mytilidae) from deep-sea hydrothermal vents. Conserv. Genet. 3:25–27.

Serafim, A., B. Lopes, R. Company, A. M. Ferreira, and M. J. Bebianno. 2008. Comparative petroleum hydrocarbons levels and biochemical responses in mussels from hydrothermal vents (*Bathymodiolus azoricus*) and coastal environments (*Mytilus galloprovincialis*). Mar. Pollut. Bull. 57:529–537.

Sharma, P. C., A. Grover, and G. Kahl. 2007. Mining microsatellites in eukaryotic genomes. Trends Biotechnol. 25:490–498.

Smith, P. J., S. M. McVeagh, Y. Won, and R. C. Vrijenhoek. 2004. Genetic heterogeneity among New Zealand species of hydrothermal vent mussels (*Mytilidae: Bathymodiolus*). Mar. Biol. 144:537–545.

Southward, E. C. 2008. The morphology of bacterial symbioses in the gills of mussels of the genera *Adipicola* and *Idas* (Bivalvia: Mytilidae). J. Shellfish Res. 27:139–146.

Tamura, K., J. Dudley, M. Nei, and S. Kumar. 2007. MEGA4: molecular evolutionary genetics analysis (MEGA) software version 4.0. Mol. Biol. Evol. 24:1596–1599.

Thubaut, J. 2012. Evolutionary history and biology of the populations of mytilids associated with sunken organic substrates. Comparative approach among Bathymodiolinae from deep-sea reducing environments. [Ph.D. thesis], Muséum National d'Histoire Naturelle, Paris.

Van Buskirk, J. 2009. Getting in shape: adaptation and phylogenetic inertia in morphology of Australian anuran larvae. J. Evol. Biol. 22:1326–1337.

Van Dover, C. L., S. E. Humphris, D. Fornari, C. M. Cavanaugh, R. Collier, S. K. Goffredi, et al. 2001.

Biogeography and ecological setting of Indian Ocean hydrothermal vents. Science 294:818.

Voight, J. R., and M. Segonzac. 2012. At the bottom of the deep blue sea: a new wood-boring bivalve (Mollusca, Pholadidae, *Xylophaga*) from the Cape Verde Abyssal Plain (subtropical Atlantic). Zoosystema 34:171–180.

Vrijenhoek, R. C. 2009. Cryptic species, phenotypic plasticity, and complex life histories: assessing deep-sea faunal diversity with molecular markers. Deep-Sea Res. II 56:1713–1723.

Wolff, T. 1979. Macrofaunal utilization of plant remains in the deep sea. Sarsia 64:117–136.

Won, Y., C. R. Young, R. A. Lutz, and R. C. Vrijenhoek. 2003. Dispersal barriers and isolation among deep-sea mussel populations (*Mytilidae: Bathymodiolus*) from eastern Pacific hydrothermal vents. Mol. Ecol. 12:169–184.

Won, Y. J., W. J. Jones, and R. C. Vrijenhoek. 2008. Absence of cospeciation between deep-sea Mytilids and their thiotrophic endosymbionts. J. Shellfish Res. 27:129–138.

Woods, H. A., A. L. Moran, C. P. Arango, L. Mullen, and C. Shields. 2009. Oxygen hypothesis of polar gigantism not supported by performance of Antarctic pycnogonids in hypoxia. Proc. Biol. Sci. 276:1069–1075.

Yamanaka, T., C. Mizota, Y. Maki, K. Fujikura, and H. Chiba. 2000. Sulfur isotope composition of soft tissues of deep-sea mussels, *Bathymodiolus* spp., in Japanese waters. Benthos Res. 55:63–68.

Yamanaka, T., C. Mizota, Y. Fujiwara, H. Chiba, J. Hashimoto, T. Gamo, et al. 2003. Sulphur-isotopic composition of the deep-sea mussel *Bathymodiolus marisindicus* from currently active hydrothermal vents in the Indian Ocean. J. Mar. Biol. Assoc. U. K. 83:841–848.