Range extensions of Pacific bone-eating worms (Annelida, Siboglinidae, Osedax)

Abstract First described in 2004 off California, Osedax worms are now known from many of the world's oceans, ranging from 10 to over 4000 m in depth. Currently, little is known about species ranges, since most descriptions are from single localities. In this study, we used new sampling in the north-eastern Pacific and available GenBank data from off Japan and Brazil to report expanded ranges for five species: Osedaxfrankpressi, O.knutei, O.packardorum, O.roseus and O.talkovici. We also provided additional DNA sequences from previously reported localities for two species: Osedaxpriapus and O.randyi. To assess the distribution of each species, we used cytochrome c oxidase subunit I (COI) sequences to generate haplotype networks and assess connectivity amongst localities where sampling permitted. Osedaxfrankpressi, O.packardorum, O.priapus, O.roseus and O.talkovici all had one or more dominant COI haplotypes shared by individuals at multiple localities, suggesting high connectivity throughout some or all of their ranges. Low ΦST values amongst populations for O.packardorum, O.roseus and O.talkovici confirmed high levels of gene flow throughout their known ranges. High ΦST values for O.frankpressi between the eastern Pacific and the Brazilian Atlantic showed little gene flow, reflected by the haplotype network, which had distinct Pacific and Atlantic haplotype clusters. This study greatly expands the ranges and provides insights into the phylogeography for these nine species.


Introduction
Osedax (Rouse et al. 2004), part of Siboglinidae, secrete acid to dissolve sunken bone and teeth as a habitat and, aided by symbiotic bacteria, feed on the organic matrix (Rouse et al. 2004, Goffredi et al. 2005, Tresguerres et al. 2013, Rouse and Goffredi 2023).Osedax can exploit the remains of diverse vertebrates, from sharks to teleost fishes to mammals, which, together with high fecundity and lecithotrophic larvae may enable them to span extensive ranges (Rouse et al. 2009, Rouse et al. 2018, Zhou et al. 2020, Rouse and Goffredi 2023).To date, 29 Osedax species have been formally named, with several others yet to be described (Rouse et al. 2004, Rouse et al. 2018, Fujiwara et al. 2019, McClain et al. 2019, Shimabukuro and Sumida 2019, Eilertsen et al. 2020, Georgieva et al. 2023a); see Suppl.material 1.Most Osedax species have only been collected from their type localities (Rouse et al. 2018), though there are a few exceptions (Figs 1, 2); see Suppl.material 1.For example, Osedax rubiplumus (Rouse et al. 2004), originally described from Monterey Bay in central California at 2891 m depth, has subsequently been found in the eastern Pacific, Indian Ocean and Southern Ocean (Zhou et al. 2020).Osedax deceptionensis was originally described from Deception Island near the Antarctic Peninsula (Glover et al. 2013) and was subsequently recorded from near South Georgia Island in the Subantarctic (Taboada et al. 2015).Osedax docricketts ( Rouse et al. 2018), O. randyi (Rouse et al. 2018), O. roseus (Rouse et al. 2008) and O. westernflyer (Rouse et al. 2018) were all originally described from the eastern Pacific, but are also found in the western Pacific, in Japanese waters (Rouse et al. 2018).Osedax priapus (Rouse et al. 2015) was originally described from Monterey Bay and Oregon (Rouse et al. 2015).Finally, O. frankpressi (Rouse et al. 2004) is known from the eastern Pacific and the western Atlantic (Rouse et al. 2018, Shimabukuro andSumida 2019).Much is still unknown about Osedax species distributions and the genetic structure across their ranges.In this study, we noted expanded ranges for five Osedax species, most of which were previously only known from single localities.We used haplotype networks, based on mitochondrial cytochrome oxidase subunit I (COI), to document range extensions and compare phylogeography amongst Osedax species.

Material and methods
We aligned all available mitochondrial cytochrome oxidase subunit I (COI) sequence data for Osedax from GenBank with new sequences generated from specimens collected from naturally occurring animal falls and experimentally sunken bones off California and Oregon (USA) and off the Pacific coast of Costa Rica (Tables 1, 2).DNA extractions and PCR products were amplified, purified and sequenced following previous protocols (Vrijenhoek et al. 2008, Vrijenhoek et al. 2009).Depth ranges and regions of occurrence for all Osedax species reported to date, including undescribed species referenced under informal names.Details and sources are in Suppl.material 1.

Data resources
All COI sequences in this paper are available on NCBI GenBank, see (Table 2).

Results
We  2).Osedax roseus's range, previously known from Sagami Bay (Japan) and Monterey Bay, was extended southwards to off San Diego (Fig. 1).
Osedax frankpressi, previously recorded from Monterey Bay and the Brazilian Atlantic, was found off Oregon, establishing a new northern record and also south to Costa Rica's Pacific coast (Fig. 1).The Oregon record of O. frankpressi was found at 642 m, expanding the species' minimum known depth by 787 m for a total depth range of 2249 m (642 -2891 m), representing the widest known range for any Osedax species (Fig. 2, Suppl.material 1).An additional sequence was provided from Monterey Bay (type locality) for O. randyi.New sequences were provided for O. priapus from the two previously reported localities of Oregon and Monterey Bay.
Uncorrected maximum intraspecific pairwise distances ranged from 4.5% for O. knutei and 3.9% for O. frankpressi to as low as 0.9% for O. randyi (Table 3).Osedax talkovici, O. roseus and O. packardorum had the largest sample sizes, but not the largest intraspecific pairwise distances (Table 3).Maximum pairwise distances for O. frankpressi were 1% amongst samples from the Pacific and 1.7% for the Brazilian Atlantic ( We used TCS haplotype networks of COI to visualise the diversity and biogeography of the nine species of Osedax.The geographical distribution of O. frankpressi was the largest examined, spanning from the Pacific to Atlantic Oceans (Fig. 1).The network for O. frankpressi revealed two divergent haplotype clusters, one from Brazil and the other from Oregon, California and Costa Rica (Fig. 3).Osedax frankpressi differed across its range by nearly 3.9% (uncorrected pairwise distance) and by a minimum of 3% between the Pacific and Brazilian sequences (Fig. 3, Table 3).In the eastern Pacific, one haplotype of O. frankpressi was shared from Oregon to Costa Rica and the maximum intraspecific distance was less than 1% (Fig. 3).
Uncorrected maximum intraspecific COI pairwise distance matrices for Osedax in this study.Four species had trans-Pacific distributions.Osedax roseus was found off Japan (Sagami Bay) and California (Fig. 1).Intraspecific diversity was high with three distinct subnetworks, but limited geographic divergence was observed (Fig. 4).Several haplotypes were shared between Japan and California, although a distinct subnetwork was found in Sagami Bay (Fig. 4).Though O. docricketts, O. randyi and O. westernflyer had trans-Pacific distributions (Fig. 1), the limited samples available revealed no shared haplotypes (Figs 5, 6, 7).Haplotype diversity in western Pacific samples of O. docrickets was high compared to samples from Monterey and haplotypes were divergent (Fig. 5).9).Both species had had similar network topologies with one or two predominant haplotypes and many singleton haplotypes which were somewhat divergent (Figs 8,9).Osedax packardorum and O. talkovici were distributed from Oregon to San Diego, California (Fig. 1).Both species had many individual haplotypes as well as several haplotypes shared amongst several localities (Figs 10,11).Each showed some predominant haplotypes shared across most localities (Figs 10,11).Osedax talkovici had  the largest sample size with 116 sequences and the highest levels of haplotype variability along the eastern Pacific (Fig. 11).

Discussion
The data added in this study revealed that many Osedax species tend to exhibit higher intraspecific divergence than other siboglinid taxa with comparable ranges (Table 3).For example, the iconic vent vestimentiferan tubeworm Riftia pachyptila has a range spanning > 7000 km along the East Pacific Rise, Galapagos  3).The relatively high value obtained for O. knutei suggests that the taxon might contain cryptic species and needs further investigation with data from additional genes and samples.
Amongst annelids, the siboglinid clade Vestimentifera appears to be an extreme case of low interspecific distances, as evidenced by the nominal species Escarpia laminata, E. southwardae, E. spicata and E. tritentaculata, which actually share a COI haplotype, though data from morphology and other genes suggest that they are vaild species on present evidence (Cowart et al. 2013, Georgieva et al. 2023b).Other vestimentiferan  ' and 7.4% between O. lehmani and O. packardorum (Rouse et al. 2018).Other annelid genera and species with comparable interspecific distances include the dorvilleid Parougia, which has minimum interspecific distances of 7% or more (Yen and Rouse 2020), the phyllodocid Eumida sanguinea with minimum interspecific distance of 5.5% (Teixeira et al. 2022) and the amphinomid Eurythoe complanata cryptic species complex, with an interspecific distance of 10% in the Atlantic (Barroso et al. 2009).However, there is no clear standard when it comes to species delimitations in annelids.For example, a 5% intraspecific distance was sufficient to split the dorvilleids Ophryotroca japonica and O. glandulata (Paxton and Akesson 2010).Nygren ( 2013) found that minimal interspecific distances of 2 -23% have been used to delineate cryptic annelid species and distances of ~ 7% are often typical for named congeneric species.This places Osedax within the normal minimum interspecific ranges for annelids and makes Vestimentifera somewhat exceptional.
Large geographic ranges in Osedax did not always correspond with large intraspecific distances (Table 3).While eastern Pacific samples of O. frankpressi differed by up to 3.9% from Brazil Atlantic samples, O. knutei had greater intraspecific distances (up to 4.5%) across a range spanning only the eastern Pacific from Monterey to Costa Rica.Similarly, O. packardorum, O. priapus and O. talkovici had relatively high intraspecific distances (2% to 3%) amongst samples from the western margin of the United States.Osedax docricketts (up to 3.5%) and O. roseus (up to 2.4%) both had high intraspecific distances though they have trans-Pacific ranges.Osedax randyi and O. westernflyer also had trans-Pacific ranges, but intraspecific distances were low (≤ 1%).Osedax rubiplumus had the largest known range of any Osedax, spanning from Antarctica, across the eastern and western Pacific and the Indian Ocean; yet, its maximum COI distance has been recorded at 1.39% (GTR corrected) between California and the Indian Ocean (Zhou et al. 2020).
Despite exhibiting some relatively large geographical distances, O. packardorum, O. priapus, O. roseus and O. talkovici exhibited evidence for connectivity across their known ranges.For example, O. roseus spans > 8000 km from Sagami Bay and Monterey Bay, as demonstrated by Φ values ≤ 0.191.Φ for O. roseus was 0.00 between Monterey Bay and San Diego, suggesting that the populations might be effectively panmictic.The moniliferan siboglinid Sclerolinum contortum also has a large range, but relatively large sampling has revealed no shared haplotypes between geographical populations (Eilertsen et al. 2018).On the other hand, widely distributed Osedax species (O.packardorum, O. priapus, O. roseus and O. talkovici) had haplotypes shared across multiple localities further indicating either good dispersal potential across their respective ranges or considerable intermediate habitat (i.e.bones).
Eight Osedax species had no haplotypes shared across multiple localities.For O. randyi and O. westernflyer, the lack of shared haplotypes was likely due to very small sample sizes.Conversely O. docricketts and O. knutei might encompass cryptic species ST ST complexes.For example, nine divergent O. docricketts COI sequences occurred in the Sagami Bay, suggesting that cryptic species may occur in Japanese waters, while the real O. docricketts may occur in both Sagami Bay and Monterey (Fig. 5).The most divergent O. docricketts sequence exhibited 55 nucleotide substitutions from the holotype sequence (asterisk in Fig. 5) (Rouse et al. 2018).In contrast, the O. talkovici sample included 116 sequences and had a maximum intraspecific distance of 2.3% (Fig. 11).Osedax knutei had the largest intraspecific distance of any Osedax species at 4.5%.The haplotype network for O. knutei showed (Fig. 8) that many individuals share a haplotype in Monterey Bay, but there were also divergent haplotypes in Monterey, San Diego and Costa Rica.The large intraspecific distance and the absence of shared haplotypes amongst the three localities suggested that O. knutei could be a cryptic species complex, though in sympatry in Monterey Bay.
Osedax frankpressi and O. rubiplumus have the broadest known geographic and depth ranges in this genus (Fig. 2, Suppl.material 1).Osedax frankpressi also had the largest Φ values and one of the greatest intraspecific distances reported in this study (Tables 3, 4 ).No COI haplotypes were shared between the Brazilian Atlantic and eastern Pacific samples; however, one common haplotype was shared amongst Oregon, Monterey Bay and Costa Rica samples.A prior study found ~ 3% divergence between Atlantic (Brazil) and the Pacific (California to Costa Rica) samples, with maximum distances of 0.7% within the Brazil population and 0.3% in the Pacific (Shimabukuro and Sumida 2019).Adding in the new sequences from Oregon, California and Costa Rica samples raised the intraspecific pairwise distances to nearly 3.9%, though the minimum distance between the Brazilian Atlantic and the Pacific remained ~ 3%.The Φ value of 0.86 for O. frankpressi clearly demonstrated population subdivision between Pacific and Atlantic populations.Although one haplotype was shared amongst samples from Oregon to Costa Rica, a distance of over 6,000 km, further sampling of bones along the east and west coasts of South America might reveal evidence of historical connectivity between Atlantic and Pacific populations, as previously suggested (Shimabukuro and Sumida 2019).

Miyamoto et al. 2013).
While it is clear that many Osedax species are known to be widely dispersed, the large number of species found in Monterey Bay is interesting (Rouse et al. 2004, Vrijenhoek et al. 2009, Rouse et al. 2015, Rouse et al. 2018).Perhaps other deepocean canyons will reveal comparable species diversity as exploration and sampling increase worldwide.Osedax's life history traits make them well suited to wide oceanic dispersal and ecological success.As this study demonstrates, a number of Osedax species are as widely distributed as other deep-sea invertebrates that experience little population subdivision across their ranges.

Figure 1 .
Figure 1.Map of geographic distributions of Osedax species analysed in this work.This map was generated using the R package marmap (Pante and Simon-Bouhet 2013).

Figure 4 .
Figure 4. Osedax roseus COI haplotype network coloured by sampling locality.Cross-hatches and black circles represent missing mutations.Network made with alignment of 730 bp.

Figure 9 .
Figure 9. Osedax priapus COI haplotype network coloured by sampling locality.Cross-hatches and black circles represent missing mutations.Network made with alignment of 891 bp.

Table 1 .
Number of COI sequences of Osedax used in this study and number of samples from each locality.Range extension = *.

Table 2 .
GenBank accession numbers used for the Osedax species in this study.Alternative names listed on GenBank are also listed.New sequences are in bold.A total of 258 new sequences were included in this study and released on GenBank.