Dispersal and Differentiation of Deep-Sea Mussels of the Genus Bathymodiolus (Mytilidae, Bathymodiolinae)

1 Institute of Biological Sciences, University of Tsukuba, Tsukuba, Ibaraki 305-8572, Japan 2 Institute of Agricultural and Forest Engineering, University of Tsukuba, Tsukuba, Ibaraki 305-8572, Japan 3 Research Program for Marine Biology and Ecology, Japan Agency for Marine-Earth Science and Technology (JAMSTEC), Natsushima, Yokosuka, Kanagawa 237-0061, Japan 4 Faculty of Education and Human Sciences, University of Yamanashi, Kofu, Yamanashi 400-8510, Japan


Introduction
Deep-sea hydrothermal vents and their attendant dense biological communities were first discovered along the Galapagos Rift [1].Since then, various deep-sea communities surviving under reductive environments rich in sulfide and methane have been discovered in hydrothermal vents on spreading ridges and back-arc basins and in coldwater seeps along subduction zones.These communities contain many endemic species whose primary production is based on bacterial chemosynthesis.Mussels of the genus Bathymodiolus are among the dominant macroorganisms in these communities.They rely primarily on chemoautotrophic endosymbionts for their nutrition similar to the other dominant groups of macroorganisms, such as vesicomyid clams and vestimentiferan tubeworms.The deep-sea mussels belong to one of the subfamilies, Bathymodiolinae, in the family Mytilidae of molluscan Bivalvia.Since the description of B. thermophilus in 1985 [2], 19 species of the genus Bathymodiolus have thus far been described [3][4][5][6][7][8][9][10][11][12].Three bathymodioline species belonging to genera, Tamu and Gigantidas, have been described [6,10,13].
Patchy and ephemeral deep-sea hydrothermal vents and coldwater seeps are separated from each other by various distances, for example, vent sectors are usually separated by a few hundred kilometers and within a vent sector, vent fields including sites which undergo the same temporal variations are separated by hundreds of meters to a few kilometers [14].It is likely for the organisms of chemosynthesis-based communities to be genetically isolated in these discontinuous habitats; however, in Japanese waters, identical Bathymodiolus species are distributed in the Sagami Bay and the Okinawa Trough, which are approximately 1500 km away from each other [4].On the other hand, there is no species shared between the Sagami Bay and the Izu-Ogasawara Islandarc, which are approximately 500 km away from each other.Thus, speciation events do not necessarily depend on the geographical distances between habitats.Genetic differentiation and consequent speciation of deep-sea organisms in the community are caused by a combination of factors shared by diverse taxa (topography, geological histories, and oceanic currents) and those unique to their respective taxa (dispersal ability, physiology, and settlement cues) [15].The dispersal ability of Bathymodiolus mussels is suggested to be high based on the larval shell morphology [16] and small egg size [17], which favors colonization in patchy and ephemeral habitats.Studies on genetic population structures can provide useful information on evolutionary processes such as dispersion, isolation, and speciation of deep-sea macroorganisms.It is tempting to examine the intraspecific relationships of Bathymodiolus mussels to search for factors that lead to speciation and populational differentiation.
Hydrothermal vents and cold-water seeps are driven by different geological processes.Hydrothermal vents are located at spreading centers and back-arc basins and emit water that is heated by the underlying magma chambers.Cold-water seeps are situated in passive margins along subduction zones and supply seawater, which is as cold as the ambient deep-sea water.Seeps are relatively stable, while vents persist for only a few decades [18].Only 3 Bathymodiolus species in Japanese waters are capable of inhabiting both vents and seeps [4], although many species of chemosynthesis-based communities are restricted to either.This study examines whether the seep and vent populations of these Bathymodiolus species are genetically differentiated as a consequence of adaptation to highly different environments.
Dispersal ability and adaptability to the deep-sea environments have been found to be associated with speciation and thus the evolutionary process of deep-sea organisms [15,19,20].Few studies of genetic population structures aimed at gaining an insight into dispersal ability have been done so far.Exceptions are for the northern and southern Bathymodiolus species of the East Pacific Rise [21] and the Mid-Atlantic Ridge [22].There was no evidence of dispersal of northern species to the territory of southern species and vice versa [21,22], with hybrid zones on the boundary of the territories in the case of the Atlantic mussels [22].Genetic population structures of neoverrucid barnacles showed that they are unable to migrate between the Izu-Ogasawara Island-arc and the Okinawa Trough [23].A study using East Pacific annelids showed that genetic population structures differed among species and suggested that those annelids had their own distinct abilities to disperse [24].
Only some Bathymodiolus species from restricted areas have been subjects in earlier molecular phylogenetic studies [25][26][27][28].Then, given increasing sequence data, molecular phylogenetics searched for the phylogeny of about ten species [29,30], and mytilid relatives from sunken whale carcasses and wood were included to trace the origin of Bathymodiolus mussels [31,32].In our previous studies [33], we showed, by sequencing of the mitochondrial COI (cytochrome c oxidase subunit I) gene of more than 15 nominal and cryptic bathymodioline species, that mussels in the subfamily Bathymodiolinae comprised 3 groups with one exception of Tamu fisheri.The first group (Group 1) consisted of the West Pacific and Atlantic Bathymodiolus (Group 1-1) and Gigantidas mussels (Group 1-2).The second group (Group 2) consisted of Bathymodiolus mussels, which were subdivided into 3 subclusters (the Indo-West Pacific, Atlantic, and East Pacific species).The third group (Group 3) consisted of the West Pacific Bathymodiolus mussels.In the present study, we employed the faster-evolving mitochondrial ND4 (NADH dehydrogenase subunit 4) gene to investigate the genetic population structure and assess the dispersal ability and adaptability to deep-sea environments of Bathymodiolus mussels.We also investigated the phylogenetic relationships of deep-sea Bathymodiolus mussels and their mytilid relatives to understand the evolutionary processes of deep-sea animals.

Materials.
The specimens used in this study are listed in Table 1, and the collection sites are mapped in Figure 1 sp.) was collected as described previously [27].Adipicola pacifica, A. crypta, and Benthomodiolus geikotsucola (the subfamily Modiolinae) were collected from sunken whale carcasses during dives by submersibles from JAMSTEC.The mussels attached to sunken wood (modioline A. iwaotakii and Idasola japonica) were obtained by trawling.All the mussels collected for this study were frozen and preserved at −80 • C or preserved in 100% ethanol, and deposited in JAMSTEC.

Analysis.
The DNA sequences were edited and aligned using DNASIS (Hitachi Software Engineering Co., Ltd., Tokyo) and MEGA 3.1 [37] and were corrected by visual inspection.We used 423-bp ND4 sequences and constructed dendrograms by the neighbor-joining (NJ) and maximum parsimony (MP) methods using PAUP * 4.0 beta10 [38].Genetic distances were computed according to Kimura's two-parameter method [39].The reliability of the trees was evaluated by producing 1,000 bootstrap replicates.The majorityrule consensus MP tree was constructed by conducting a heuristic search based on the 1,000 bootstrap replicates with an unweighted ts/tv ratio.The Bayesian tree was constructed using MrBayes version 3.1 [40] based on the model evaluated by the Mrmodel test 2.2 [41].The Monte Carlo Markov chain (MCMC) length was 1 × 10 6 generations, and we sampled the chain after every 100 generations.MCMC convergence was assessed by calculating the potential scale reduction factor, and the first 2 × 10 3 generations were discarded.We used Modiolus nipponicus (the subfamily, Modiolinae) as an outgroup species.We estimated the genetic divergences (F st ) and the bidirectional mean rates of gene flow (Nm; the virtual average number of migrants exchanged per generation) between the populations using Arlequin 3.1 [42].We evaluated the significance of F st by calculating 1 × 10 6 values.We also calculated the mismatch distribution [43] and constructed

Phylogenetic Relationships among the Bathymodiolus
Mussels.The partial DNA fragments of the mitochondrial ND4 gene (423 bp) were sequenced from 1 to 5 specimens (if available, more than 5) of Bathymodiolus species and their mytilid relatives.The sequence data were deposited in the DDBJ, EMBL, and GenBank databases under accession numbers AB478422-AB478475.A part of the sequence data was previously reported (AB175280-AB175326, [29]).No deletions or insertions were found in the sequences after excluding the sequences of Modiolus and Mytilus.Sequences encoding 2 amino acids were deleted or inserted when Modiolus and Mytilus species were included.The deep-sea mussels of the subfamily Bathymodiolinae formed a poorly supported cluster consisting of 3 major groups (Figure 2).The ND4 trees constructed by the NJ, MP, and Bayesian methods yielded fundamentally the same topology.Previously published COI trees [33] also presented essentially the same topology (Figure 3).The first group (Group 1) was marginally supported (73, 57, and 0.97 for NJ, MP, and Bayes, resp.) and was subdivided into 2 clades.Group  unidentified mussels from Aitape (Aitape G. sp.) and off Ashizuri Cape (Ashizuri G. sp.) in the West Pacific.We regarded the mussels from the Sumisu Caldera (Sumisu G. sp.) and the Nikko Seamount (Nikko G. sp.) as conspecific with G. horikoshii, because they were very closely related to each other.Adipicola crypta belonging to the subfamily Modiolinae formed a marginally supported cluster together with Group 1-2.On the other hand, A. crypta was a sister taxon to the cluster including Groups 1 and 2 in the COI tree (Figure 3).The species of the subfamily Modiolinae obtained from sunken wood and whale carcasses and seeps, namely, Benthomodiolus lignicola, Benthomodiolus geikotsucola, Idas macdonaldi, I. washingtonia, Adipicola iwaotakii, A. pacifica, Idasola japonica, and 1 unidentified mussels from Macauley Cone (B.sp.NZ3), were outside the cluster including the bathymodioline species and modioline A. crypta.The relationships were also shown in the COI tree (Figure 3).B. sp.NZ3 was reported to belong to the genus Bathymodiolus [30].However, it was distantly related to Bathymodiolus as in the COI tree, and its phylogenetic position remains to be studied.

Genetic Population Structure.
A minimum spanning tree was constructed based on the ND4 sequences of a total of 60 specimens of B. platifrons and B. japonicus (Figure 4(a)).The new sequence data were deposited in the DDBJ, EMBL, and GenBank databases under accession numbers AB480561-AB480578.As expected, no haplotype was shared by B. platifrons and B. japonicus.Significantly high F st (0.968 = (0.965 + 0.972)/2) and small Nm (0.016 = (0.018 + 0.014)/2) were estimated between the 2 species (Table 2).On the other hand, the major haplotype was shared by 24 (60%) specimens of B. platifrons from the seeps of the Sagami Bay and the vents of the Okinawa Trough.Negative F st (−0.007) and Nm of infinity were     estimated between the 2 populations.These results showed high gene flow between the seep and vent populations.
Neither the seep nor the vent population was monophyletic, and the members of both populations were intermingled (Figure 2), indicating that the type of environment is not the primary factor responsible for habitat segregation.
A minimum spanning tree was also constructed based on the ND4 sequences of a total of 70 specimens of B. septemdierum, B. brevior, and B. marisindicus (Figure 4(b)).The new sequence data were deposited in the DDBJ, EMBL, and GenBank databases under accession numbers AB485606-AB485629.The haplotype of the greatest majority was shared by 21 (30%) specimens of B. septemdierum from the Myojin Knoll and the Suiyo Seamount of the Izu-Ogasawara Island-arc and B. brevior from the North Fiji Basin of the Southwest Pacific.One of the 2 major haplotypes was possessed exclusively by 6 specimens (8.6%) of B. marisindicus from the Kairei field of the Southern Central Indian Ridge, and the other was shared by 7 specimens (10%) of the 3 species, although they are distinct species.Low F st (0.017) and large Nm (29.619) were estimated between B. septemdierum and B. brevior, while significantly high F st (0.254 = (0.252 + 0.257)/2) and small Nm (1.466 = (1.486+ 1.445)/2) were estimated between B. marisindicus and the 2 West Pacific species (Table 2).
Mismatch distribution showed that the τ values were 2.197 for B. japonicus and 2.111 and 1.967 for B. platifrons from the Sagami Bay and the Okinawa Trough, respectively (Figures 5(a

Phylogenetic Relationships of the Bathymodiolus Mussels and their
Relatives.The present study based on the ND4 sequences presented fundamentally the same phylogenetic relationships of the Bathymodiolus mussels and their relatives as those reported by our previous study based on the COI sequences [33].Although the discrepancy was found in the positions of Tamu fisheri and Adipicola crypta, it did not affect the major conclusions presented in our previous study.The bathymodioline Tamu fisheri was closely related to the Bathymodiolus mussels of Group 2 in the ND4 study, while it was shown to be distantly related to the other bathymodioline species in the COI study.Although the modioline Adipicola crypta was closely related to Group1-2 in the ND4 study, it was more closely related to the cluster consisting of Groups 1-1, 1-2, and 3 in the COI study.Both studies suggested that the subfamily Bathymodiolinae and the genus Bathymodiolus were not monophyletic because the monophyly of the former and that of the latter were refuted by the existence of A. crypta and two Gigantidas species, respectively.More extensive morphological investigations are needed to reevaluate the classification.The branching orders in Groups 1 to 3 and in the 3 subclusters of Group 2 differed between the ND4 and COI studies.However, their divergences appeared trichotomous because of the short branch lengths between the nodes leading to the groups and clusters.

Adaptation of the Mussels to the Abyssal Environment.
The present study supports the "Evolutionary stepping stone hypothesis" [25,45].According to this hypothesis, the ancestors of Bathymodiolus mussels exploited sunken wood and whale carcasses in their progressive adaptation to the deepsea environment with regard to nutrition and tolerance to high pressure, cold seawater, and toxicity of hydrogen sulfide.Both ND4 and COI [33] trees showed that species from sunken wood and whale carcasses assumed the outgroup position to the Bathymodiolus and Gigantidas mussels from the vents and cold seeps, with only the exception of A. crypta from the whale carcasses.Shallow water mytilid mussels such as Modiolus nipponicus, Mytilus edulis, M. gallloprovincialis, and M. trossulus were positioned more distantly to the vent/seep mussels.The findings indicate an evolutionary transition from the shallow water to vent/seep sites via the wood/whale carcass sites, and a reversion to the whale carcass sites from the vents or seeps in the case of A. crypta.The studies also suggest independent invasion into the seeps in case of Idas macdonaldi and into the vents in case of I. washingtonia and B. sp.NZ3.
Most species of the deep-sea chemosynthesis-based communities are restricted either to seeps or vents.Only three known Bathymodiolus species endemic to Japanese waters can inhabit both seeps and vents.Bathymodiolus japonicus and B. platifrons live in the seeps in the Sagami Bay and the vents in the Okinawa Trough.We analyzed the genetic population structure to examine whether the seep and vent populations of B. platifrons were genetically differentiated owing to their adaptation to the highly different habitats.
The present results showed no significant genetic differentiation between the seep and vent populations, indicating a high adaptability of the species to the abyssal environment.Genetic similarity between populations from the Sagami Bay and the Okinawa Trough was also shown in deep-sea bresiliid shrimp Alvinocaris longirostris [46].
Habitat segregation is not caused by habit types (seep versus vent), but depth in some species of the chemosynthesisbased community.The vestimentiferan tubeworm Escarpia sp.inhabits the seeps in Japanese waters and the vents in the Manus Basin, and no genetic differentiation was detected between their seep and vent populations, although their habitats are approximately 4400 km away from each other [47].The populations from the seeps in the Nankai Trough and the vents in the Lau Basin were not distinguished genetically in the vestimentiferan Lamellibrachia columna, although their habitats are 7500 km away from each other [48].Habitation of Calyptogena clams in Japanese waters is constrained by depth, but not restricted by the type of environment, and their colonization is not limited by geographical distances between habitats [49].Habitat segregation by depth in the Calyptogena clams is probably ascribed to the differences in their physiological tolerance to pressure [49,50].
It is unlikely that habitation is constrained by depth [26] or colonization is limited by geographical distances (as described below) in some species of the genus Bathymodiolus.Instead habitat segregation and colonization of deep-sea mussels can be ascribed to their preference to one (or some) specific ambient condition(s).Bathymodiolus species that harbor only methanotrophic endosymbionts occur in coldwater seeps and hydrothermal vents with higher methane concentrations, and the deep-sea mussels that depend on thioautotrophic endosymbionts for their nutrition occur in vents with lower methane concentrations [51].This suggests that the chemical environment in their habitats is one of the factors restricting the distribution of the Bathymodiolus species.

Dispersal of the Deep-Sea Mussels.
High gene flow (Nm = infinity) was detected between the populations from the Sagami Bay and the Okinawa Trough in B. platifrons, although the sites are more than 1500 km away from each other.An Nm value of more than 1 is indicated to be sufficient to maintain genetic continuity among populations [42].
The species in Cluster A were very closely related to one another, although they are distributed over vast distances, including B. marisindicus in the Kairei field, B. brevior in the North Fiji Basin, and B. septemdierum in the Myojin Knoll and the Suiyo Seamount.Our previous study [29] showed that their interspecific genetic distances were considerably smaller than those of species except the Cluster A species and approximated to intraspecific genetic distances of the latter.High gene flow (Nm = ca.30) existed between B. septemdierum from the Myojin Knoll and the Suiyo Seamount and B. brevior from the North Fiji Basin.The localities of the two species are approximately 5,000 km away.Furthermore, the gene flow between B. septemdierum and B. marisindicus from the Kairei field (Nm = ca.1.5) and that between B. brevior and B. marisindicus (Nm = ca.1.4) were not negligible despite significantly high genetic divergences (F st = ca.0.25 and 0.26, resp.).The locality of B. marisindicus is approximately 10,000 km away from those of B. septemdierum and B. brevior.Therefore, the present results showed that (1)  Bathymodiolus mussels are suggested to have the high dispersal ability, based on their larval shell morphology [16] and small egg size [17] that are indicative of planktotrophic (actively feeding planktonic larval) development.Developmental arrest at cold temperatures also appears to play an important role in extension of the planktonic stage and increasing the dispersal distance [52].
Some deep-sea animals are known to have an ability to disperse their larvae over very long distances.The East Pacific deep-sea mussel B. thermophilus and clam Calyptogena magnifica are estimated to have dispersal capabilities of at least 2,370 km and 3,340 km, respectively [53].The vestimentiferans Escarpia sp. and Lamellibrachia sp. have larvae with long-distance dispersal capability [47,48].Kojima et al. [54] showed the existence of active gene flow between the populations in the Manus Basin and the North Fiji Basin (3500 km away from each other) in the vent-endemic gastropod Alviniconcha sp.Riftia pachyptila had a larval stage of approximately 38 d under conditions similar to the in situ environment, suggesting that larvae can disperse over 100 km albeit influenced by deep-water circulation regimes [55].The larvae of the verrucomorph barnacle Neoverruca sp. had a planktonic period of over 70 d at 4 • C under 1 atm, suggesting its high dispersal ability [56].Since Riftia pachyptila and Neoverruca sp. have nonplanktotrophic larvae, it is conceivable that Bathymodiolus mussels can disperse their planktotrophic larvae over longer distances.
Mismatch distribution showed that the τ values decreased in the order of Sagami Bay B. japonicus, Sagami Bay B. platifrons, and Okinawa Trough B. platifrons (Figures 5(a)-5(c)).The results suggest the immigration of ancestral B. platifrons into the Okinawa Trough from the Sagami Bay, which is consistent with the history of the Japanese archipelago.The Okinawa Trough has been habitable for animals in the chemosynthesis-based community since ca.200 MYA, while the Sagami Bay since more than 500 MYA [57].However, this immigration event appears unlikely because the difference in the τ values was small (2.111 versus 1.967), and the intense stream, due to the existence of the Kuroshio Current, runs from the Okinawa Trough to the Sagami Bay down to a depth of 1,000 m.Therefore, we suppose that immigration might have occurred to the Okinawa Trough and the Sagami Bay from an unknown home location somewhere in the West Pacific.The τ values decreased in the order of B. marisindicus, B. septemdierum, and B. brevior (Figures 5(d)-5(f)), proposing that the ancestor of the 3 species might have migrated from the Southern Central Indian Ridge to the Izu-Ogasawara Islandarc via the Southwest Pacific.However, the present definite gene flow between B. septemdierum and B. brevior revealed that they were conspecific or sibling species that had recently differentiated, thus suggesting a greater probability that the Southern Central Indian Ridge might be the more ancient residence rather than the West Pacific.The possible routes for dispersal appear to be along the seeps and vents localized to north and south of Australia or through sunken wood and whale carcasses.The subduction zones reel along the South Asian Islands, and the Southeast Indian Ridge runs from the east to the west.
We suggest collecting novel or more specimens to answer the questions regarding the dispersal routes and barriers.Such future work will provide more information on the sea floor topography, oceanic circulation, distribution of unknown seeps and vents, and larval development.
. All deep-sea mussels of the genus Bathymodiolus and Gigantidas (the subfamily Bathymodiolinae), except the East Pacific species and 3 Atlantic species, were collected during dives by submersibles from the Japan Agency for Marine-Earth Science and Technology (JAMSTEC).The East Pacific species B. thermophilus and the 2 Atlantic species B. puteoserpentis and B. azoricus were collected during the cruise of the scientific research vessel Akademik Mistislav Keldysh belonging to the Institute of Oceanology of the Russian Academy of Sciences.The Atlantic species B. childressi was collected in an oil-seep in the Gulf of Mexico during R/V Seward Johnson cruise (dive number 4568).The undescribed West Pacific species from off New Zealand (herein referred to as NZ B.

Figure 1 :
Figure 1: The sampling sites of the deep-sea Bathymodiolus mussels and their relatives used in this study.Refer to Table 1 for the details of the sampling sites (1 to 42).(a) worldwide map; (b) the magnified map showing the areas around Japanese waters., hydrothermal vent;, cold-water seep; , wood/whale bone; , shallow.

Group 2
was well supported (94 and 88 for NJ and MP, respectively), but the topology of the Bayesian tree was different from those of the NJ and MP trees.This group consisted of the 8 nominal species, namely B. septemdierum and B. brevior from the West Pacific, B. marisindicus from the Indian Ocean, B. azoricus, B. puteoserpentis, B. heckerae, and B. brooksi from the Atlantic, and B. thermophilus and 1 undescribed (morphologically examined but not described yet) Bathymodiolus species from the East Pacific (East Pacific B. sp.).Bathymodiolus septemdierum, B. brevior, and B. marisindicus comprised the closely related species group (Cluster A).Mussels from the Eifuku Seamount were included in Cluster A. Group 2 was subdivided into 3 well-supported clades comprising the Indo-West Pacific, Atlantic, and East Pacific species.The only exception was the Atlantic species B. brooksi, which diverged basally to the Indo-West Pacific clade and the Atlantic clade including the other Atlantic species.In the COI tree (Figure3), B. brooksi diverged basally to the Indo-West Pacific, Atlantic, and East Pacific clades.Tamu fisheri was found to be distantly related to the other bathymodioline species in our previous study based on COI sequences (Figure3); however, the species was closely related to the Bathymodiolus species of Group 2 although the alliance was poorly supported.Group 3 was well supported (100, 99, 1.00) and consisted of the 2 nominal species B. aduloides and B. manusensis from the West Pacific.Mussels from the Lau Basin (Lau B. sp.), North Fiji Basin (NF B. sp.), and off New Zealand (NZ B. sp.) should be conspecific with B. manusensis, because they were very closely related to each other.The distribution of Group 3 was restricted to the West Pacific.

Figure 3 :
Figure3: The phylogenetic relationships of the deep-sea Bathymodiolus mussels and their relatives based on the 401-bp mitochondrial COI sequences.The tree was modified from our published data[33] for reference., hydrothermal vent; , cold-water seep; , wood/whale bone; , shallow.(c) Malacological Society of Japan.

Figure 4 :
Figure 4: The minimum spanning trees revealing the genetic population structure based on the 423-bp mitochondrial ND4 sequences.The tree (a) was constructed using total 60 specimens of B. platifrons from the seeps of the Sagami Bay (PH1∼20), B. platifrons from the vents of the Okinawa Trough (PT1∼10, 12∼15, PI1∼4, PY1, 2), and B. japonicus from the seeps of the Sagami Bay (JH1, 2, 4∼21).Black, haplotypes possessed by B. platifrons specimens from the vents; gray, haplotypes shared by B. platifrons specimens from the seeps and vents; white, haplotypes possessed by B. platifrons or B. japonicus specimens from the seeps.The tree (b) was constructed using a total of 70 specimens of B. septemdierum from the Myojin Knoll (SM1∼10) and the Suiyo Seamount (SS1∼11), B. brevior from the North Fiji Basin (BN1∼17, 19∼30), and B. marisindicus from the Kairei field (MK1∼19, MK from a database).Black, haplotypes possessed by B. marisindicus specimens; gray, haplotypes shared by the three species; white with the bold outline, haplotypes possessed by B. brevior; white with the thin outline, haplotypes possessed by B. septemdierum or those shared by B. brevior and B. septemdierum.

Figure 5 :
Figure 5: Mismatch distribution based on the 423-bp mitochondrial ND4 sequences.(a) B. japonicus from the Sagami Bay; B. platifrons from the Sagami Bay; B. platifrons from the Okinawa Trough; (d) B. marisindicus from the Kairei field; (e) B. septemdierum from the Myojin Knoll and the Suiyo Seamount; (f) B. brevior from the North Fuji Basin.
gene flow was present between B. septemdierum and B. brevior and (2) although B. marisindicus was not isolated from B. septemdierum and B. brevior, gene flow was relatively limited.These results indicate the high dispersal ability of deep-sea mussels, although various factors such as oceanic circulation patterns, water temperature, and sea-floor topography may change the actual dispersal distances.

Table 1
[44]mum spanning tree[44]using Arlequin 3.1.We performed goodness-of-fit tests to evaluate discrepancy between the observed and model values of the mismatch distribution.The age of demographic expansion (τ = 2ut) is proportional to the number of generations (t) since a population at equilibrium of size entered a demographic expansion phase, although the mutation rate (u) of mytilid mussels is unknown.For analysis of the genetic population structures, we determined the ND4 sequences of 20 specimens each of B. platifrons from the Sagami Bay and the Okinawa Trough and B. japonicus from the Sagami Bay.We also determined the ND4 sequences of 20 specimens of B. marisindicus from the Kairei field, 29 specimens of B. brevior from the North Fiji Basin, and 21 specimens of B. septemdierum from the Myojin Knoll and the Suiyo Seamount.

Table 2 :
F st (bold values) and Nm (light values).JH, B. japonicus from the Sagami Bay; PH, B. platifrons from the Sagami Bay; PIPTPY, B. platifrons from the Okinawa Trough; MK, B. marisindicus from the Kairei Field; SMSS, B. septemdierum from the Myojin Knoll and the Suiyo Seamount; BN, B. brevior from the North Fiji Basin.