Genetic connectivity from the Arctic to the Antarctic: Sclerolinum contortum and Nicomache lokii (Annelida) are both widespread in reducing environments

The paradigm of large geographic ranges in the deep sea has been challenged by genetic studies, which often reveal putatively widespread species to be several taxa with more restricted ranges. Recently, a phylogeographic study revealed that the tubeworm Sclerolinum contortum (Siboglinidae) inhabits vents and seeps from the Arctic to the Antarctic. Here, we further test the conspecificity of the same populations of S. contortum with additional mitochondrial and nuclear markers. We also investigate the genetic connectivity of another species with putatively the same wide geographic range - Nicomache lokii (Maldanidae). Our results support the present range of S. contortum, and the range of N. lokii is extended from vents and seeps in the Nordic Seas to mud volcanoes in the Barbados Trench and Antarctic vents. Sclerolinum contortum shows more pronounced geographic structure than N. lokii, but whether this is due to different dispersal capacities or reflects the geographic isolation of the sampled localities is unclear. Two distinct mitochondrial lineages of N. lokii are present in the Antarctic, which may result from two independent colonization events. The environmental conditions inhabited by the two species and implications for their distinct habitat preference is discussed.


Specimens of
. A map of the sampling localities was generated using QGIS 34 .
For Sclerolinum contortum, the same DNA extractions were used as in Georgieva et al. 25 , and sampling and DNA extraction of S. contortum is described in that paper. DNA of Nicomache lokii was extracted using the QIAGEN DNeasy Blood and Tissue Kit, following the manufacturer's protocol (spin-column protocol). PCR and sequencing. Four markers were selected for population genetic analyses: the mitochondrial Cytochrome C oxidase subunit I (COI), 16S rRNA and Cytochrome B (CytB), and the nuclear 28S rRNA (for primers see Supplementary Table S1). COI had already been sequenced for the three populations of S. contortum for a previous paper 25 , but the other three markers were sequenced for this study. The same four markers were used for species delimitation analyses of both N. lokii and S. contortum. For phylogenetic analysis of Nicomache, we used three of the population level markers, COI, 16S and 28S, and also the more conservative nuclear marker 18S rRNA. Three additional species of Nicomache were sequenced (N. lumbricalis, N. minor and N. quadrispinata), and also one species of the closely related genus Petaloproctus (Petaloproctus tenuis) as outgroups. All sequences used for the phylogenetic analyses, including representative sequences of all haplotypes were submitted to GenBank, see Supplementary Tables S2-S4 for accession numbers. PCR reactions and cycling profiles are listed in Supplementary methods.
Quality and quantity of amplicons were assessed by gel electrophoresis imaging using a FastRuler DNA Ladder (Life Technologies) and GeneSnap and GeneTools (SynGene) for image capture and band quantification. When multiple bands were present, the total PCR product was run on a new gel and the target band was extracted from the gel using MinElute Gel Extraction Kit (QIAGEN) following the manufacturer's protocol. In cases where there were multiple weak bands, a small piece of the desired band was cut out from the gel and re-amplified using the original PCR protocol. Successful PCRs were purified using Exonuclease 1 (EXO, 10 U mL −1 ) and Shrimp 90 Alkaline Phosphatase (SAP, 10 U mL −1 , USB Europe, Germany) in 10 μL reactions (0.1 mL EXO, 1 μL SAP, 0.9 μL ddH 2 O, and 8 μL PCR product). Samples were incubated at 37 °C for 15 min followed by an inactivation step at 80 °C for 15 min. The purified PCR products were sequenced using BigDye v3.1 (Life Technologies) and run on an Automatic Sequencer 3730XL at the sequencing facility of the Institute of Molecular Biology, University of Bergen. COI, 16S and 18S were sequenced with forward and reverse primers, while CytB and 28S were mainly sequenced with forward primers only.
Sequence assembly and alignments. Forward and reverse sequences were assembled in Geneious v.6.1.8 (Biomatters Ltd.) and checked for contamination using BLAST 35 . For the protein coding genes (COI and CytB), sequences were translated to amino acids in Geneious to check for stop codons. Sequences were aligned in Geneious using MUSCLE 36 , except the 16S dataset for phylogenetic analysis of N. lokii which was aligned with MAFFT v7 using the online server 37 . Alignments for population genetic analyses were trimmed to the shortest sequence, while for the phylogenetic analysis missing data at the ends were coded with question marks.

Population genetic analyses.
To obtain a good number of samples per population, samples from multiple sampling stations in the Barbados Trench and eastern Scotia Sea were pooled for population genetic analyses (see Fig. 1 and Table 1). The two sampling stations in the Barbados Trench (Atalante East and Manon) are closely located, while the Scotia Sea stations are further apart (approximately 400 km between E2 and E9, see Table 1 for coordinates). However, the same haplotypes were mainly found in all the Scotia Sea stations (see Results). Since HMMV is a cold seep, samples from this locality were indicated with a darker shade of green in the haplotype networks to separate it from the samples from the Loki's Castle vent field. However, the HMMV and Loki's Castle samples did not show much genetic variation, if any, in the haplotype networks, and these samples were pooled for the subsequent analyses. Between 10-35 specimens were sequenced from each population for each of the four markers (see Table 2). Six heterozygous individuals were identified for 28S by double peaks in the chromatograms, and out of these two were heterozygous for more than one position (Antarctic haplotypes A and B). Although the inference of haplotypes in individuals with multiple heterozygous positions is not straightforward, in this case the heterozygous positions matched the polymorphic positions between the only two haplotypes present in the population (as seen in homozygotes), and it was assumed that the individuals carried one copy of each of these haplotypes. For homozygotes, the sequences were duplicated to generate a dataset with two sequences per individual.
Haplotype networks were generated in TCS v1.21 38 with gaps as a fifth state, drawn in PopART (http://popart. otago.ac.nz), and final graphical adjustments were made in Adobe Illustrator CS6 (Adobe Systems, San Jose, CA, USA). Genetic diversity within each population (Arctic, Gulf of Mexico/Barbados and Antarctic) was calculated in DnaSP 39 . Pairwise p-distance (uncorrected p) matrixes for COI haplotypes and within group mean p-distances were calculated in MEGA7 40 . GenAlEx was used to perform AMOVA (analysis of molecular variance) and to calculate pairwise PhiPT values 41 . A permutation test with 999 permutations was performed for both the AMOVA and pairwise PhiPT values to test if the observed values were significantly different from the null distribution.
Phylogenetic analyses and species delimitation. Species delimitation for S. contortum and N. lokii was performed under the multi-species coalescent model (MSC) using STACEY 42 , a package for BEAST2 43 . Firstly, species delimitation for both species were performed without any outgroups, using the same four markers as in the population genetic analyses. Secondly, a species tree was reconstructed for Nicomache, also using STACEY, based on COI, 16S, 28S and 18S. This analysis included three additional species of Nicomache (N. lumbricalis, N. minor and N. quadrispinata), and also one species of the closely related genus Petaloproctus (Petaloproctus tenuis), which was used to root the tree. Species tree reconstruction under the multispecies coalescent model requires sequences from multiple specimens per species, and this was not available for other species of Sclerolinum apart from S. contortum, so therefore no phylogeny was reconstructed for this genus.
Datasets for species delimitation and phylogenetic reconstruction were assembled so that all haplotypes and combinations of haplotypes present in each population were represented (see Supplementary Tables 2-4 Table S4). A concatenated matrix of all genes was generated using Sequence Matrix 44 , which automatically generates empty sequences of question marks for individuals with missing data for one or more genes. Substitution saturation was tested for the third position of COI and CytB using the Xia method implemented in DAMBE6 45 , but no saturation was detected. The best partition scheme and the best fitting model of evolution for each partition was found using Partition Finder v2.1.1 with the greedy algorithm and PhyML 46,47 . The input file was partitioned with each gene and each codon position of COI and CytB as separate partitions. The best-fit partition schemes and site models for each dataset can be found in Supplementary Table S5. Due to statistical concerns regarding the co-estimation of the gamma and invariant-site parameters (discussed in the RAxML manual 48 ) we chose to use only the gamma model for rate heterogeneity for the partitions where Partition Finder suggested to use both. STACEY incorporates species delimitation and phylogenetic reconstruction in the same MCMC run, and therefore it does not require a priori species designations. All specimens were defined as different species in the initial settings, leaving delimitation to the analysis. For the phylogenetic analysis of Nicomache the tree was rooted by defining the ingroup (all specimens of Nicomache spp.) as monophyletic. For all analyses, the site models were linked and substitution models defined as suggested by Partition Finder, but the phylogenetic analysis of Nicomache would not reach convergence with this model. Therefore, the model for this analysis was simplified by setting all site models to HKY, but with site heterogeneity models as before. For all analyses, the clock models were unlinked, while the tree model was linked for all the mitochondrial partitions as these are inherited as a unit,  34 . The two stations in the Barbados Trench were located so close together that these are indicated with only one circle. There were also multiple sampling stations on the E2 ridge in the eastern Scotia Sea (see Table 1). and thus expected to share the same phylogenetic history. It can be argued that 28S and 18S are not completely independent either, as these genes are located close together in the nuclear genome, and therefore we ran two separate analyses for the phylogeny of Nicomache, with the tree models for these genes linked and unlinked. For comparison, species delimitation was also performed for the phylogenetic analyses of Nicomache. Priors and settings for the STACEY analyses can be found in Supplementary Methods.

Results
Population genetic analyses. Between 55-116 sequences per gene for N. lokii and 18-50 for S. contortum were included in the population genetic analyses (counting two sequences per individual for 28S, see Table 2). Alignment lengths varied between 217 bp (CytB for S. contortum) to 656 bp (COI for Nicomache lokii, see Table 2). COI was the most variable marker in both datasets, with the highest number of haplotypes and polymorphic sites (see Table 2 and Fig. 2). The highest pairwise p-distance for COI in S. contortum was 1.7%, between haplotype E from Loki's Castle and haplotype H from GoM (see Supplementary Table S6). For N. lokii the highest pairwise p-distance was 4.1% between haplotype A in the Antarctic and haplotypes F and I from Loki's Castle and haplotype J in the Barbados Trench (see Supplementary Table S7). Not considering haplotype A, the highest COI p-distance for N. lokii was 1.9% between haplotype B from the Barbados Trench and haplotypes I from Loki's Castle and J from the Barbados Trench.
The network topology for 16S of S. contortum was very similar to that found previously by Georgieva et al. 25 for COI, but with lower variation (see Fig. 2). For COI, CytB and 28S of S. contortum, the genetic diversity was very similar for the Arctic and the GoM populations, while the lowest diversity was found in the Antarctic population, which had only a single haplotype for all markers (see Table 2). For 16S the highest genetic diversity was found in the GoM population (Table 2). Both network topology and pairwise PhiPT values for COI and 16S showed a pattern of decreasing genetic similarity with increasing geographic distance, as both the Arctic and Antarctic populations were more genetically similar to the GoM population than to each other (Table 3). However, for COI and CytB pairwise PhiPTs were lower between the GoM and Arctic populations than between the GoM and Antarctic populations (Table 3). For 28S the variation in S. contortum was low in all populations with only 4 haplotypes in total, and one of these were shared between all populations (Fig. 2). Pairwise PhiPTs for 28S were similar for all populations (Table 3). Results from the AMOVA shows that most of the genetic variation in S. contortum was between populations for all markers (between 82-93%, see Table 4).
The patterns of genetic diversity for N. lokii were more complex. The number of haplotypes was similar between the Arctic and Barbados Trench populations (same number of haplotypes for COI and 28S, higher number in Barbados Trench for 16S and CytB), while the Antarctic population had fewer haplotypes for most of  Table 2. Genetic diversity within populations. Arc = Arctic, Bar = Barbados Trench, Ant = Antarctic, GoM = Gulf of Mexico, n = sample size (for 28S this is given as total number of sequences; two per specimen), k = number of different haplotypes, H = haplotype diversity, π = nucleotide diversity, SD = standard deviation, P = polymorphic sites, S/N = synonymous/non-synonymous substitutions, In-Dels = insertions/deletions, M = mean p-distance. *One substitution was outside the reading frame, and it could thus not be determined if it was synonymous or not. **The partition contains an in-del, which was excluded by DnaSP. This position is therefore not included in the calculations of haplotype diversity and nucleotide diversity, but all other statistics include the gap as a 5th state, and the sequences including a gap as a distinct haplotype. the markers (see Table 2). However, for the Barbados Trench and Antarctic populations, the haplotypes did not cluster together in the haplotype networks, and this was especially apparent in the mitochondrial markers (Fig. 2). The high level of divergence between mitochondrial haplotypes in the Barbados Trench and Antarctic populations was also reflected in the measures of nucleotide diversity (π) and mean p-distance within populations (Table 2).
In the Antarctic population, there were only two haplotypes for COI, 16S and 28S, while for CytB there was also a third haplotype which was one basepair different from haplotype F (Fig. 2). For all the mitochondrial markers one of the Antarctic haplotypes was closely related to haplotypes from the Arctic or the Barbados Trench (one or two mutation steps different in the networks, see Fig. 2), while the other haplotype (haplotype A for COI, 16S and CytB) was very different. For COI haplotype A was 22 mutation steps away from the nearest haplotypes (either haplotype G from the Antarctic or haplotype B from the Barbados Trench), which corresponds to 3.5% genetic distance (see Supplementary material, Table S6 for pairwise p-distances between COI haplotypes). However, for 28S the Antarctic group of specimens with the most divergent mitochondrial haplotypes (henceforth referred to as haplogroup A) had two haplotypes, and one of these (28S haplotype A) is shared with the remaining Antarctic specimens, and the other populations (Fig. 2). The two haplotypes for 28S found in the Antarctic population were found both as homozygotes (14 individuals homozygotic for A and one individual homozygotic for B) and heterozygotes (two individuals). Haplogroup A co-occurs with the second Antarctic lineage of N. lokii at both E2 and   (Table 3). However, much like the pattern found in S. contortum, the highest PhiPT values were found between the Arctic and Antarctic populations for all markers, and for COI and 28S the lowest values were found between the Arctic and Barbados Trench populations. The high genetic variation within populations of N. lokii was reflected in the AMOVA results (Table 4). For the mitochondrial markers, between 41-47% of the variation was found within populations, while for 28S the number was even higher, with 75% of the variation found within populations. This was much higher than the values found in S. contortum, where within population variation accounted for between 7-18% of the variation (Table 4).

Phylogenetic analyses and species delimitation.
Based on the criteria outlined in Methods, 21 specimens were included in the species delimitation analysis of S. contortum, of which one specimen was a "mock specimen" (see Supplementary Table S4). 27 specimens were included in the species delimitation analysis of N. lokii, of which three were "mock specimens" (see Supplementary Table S2). The species delimitation analysis for S. contortum did not recover any well supported clustering schemes. The most common scheme in the posterior distribution of trees (four species, with one species each in the Arctic and Antarctic, and two species in the Gulf of Mexico) was only found in 3% of the trees. In all species delimitation analyses of N. lokii, the clustering scheme that received the highest support (49% in the analysis without outgroups, 71% in the analysis with outgroups and unlinked nuclear markers, and 68% in the analysis with outgroups and with linked nuclear markers), recovered N. lokii from the Arctic, the Barbados Trench and the Antarctic (including haplogroup A) as a single cluster. The second-best clustering scheme in all analyses (13% in the analysis without outgroups, 25% in the analysis with outgroups and unlinked nuclear markers, and 30% in the analysis with outgroups and with linked nuclear markers) recovered the same cluster as above, except haplogroup A, which was recovered as a separate cluster. All remaining clustering schemes were found in less than 2% of the trees.
In total 102 sequences were included in the phylogenetic analyses of Nicomache, all newly produced for this study (see Supplementary Table S2). Eight specimens were included from the Arctic, eleven specimens from the Barbados Trench and three specimens from the Antarctic (of which two were from haplogroup A). 18S was sequenced for one specimen from the Arctic, four specimens from the Barbados Trench and three specimens from the Antarctic, and showed no variation between these populations. The results were mainly the same for the phylogenetic analyses with the tree models for the nuclear markers 18S and 28S linked and unlinked. The analysis with linked 18S and 28S tree models converged faster (1 × 10 9 generations versus 3 × 10 9 generations), and gave somewhat higher node support, so only the tree from this analysis is presented here (Fig. 3; Figure S1 for the tree from the analysis with unlinked nuclear markers). The phylogenetic reconstruction recovered N. lokii from the Arctic, Barbados Trench and the Antarctic as a monophyletic group with maximum support (Fig. 3). N. minor and N. lumbricalis are recovered together with high support, but the position of N. quadrispinata relative to the other Nicomache species is unresolved (Fig. 3).

Discussion
The results presented here support that all the populations of Nicomache lokii are conspecific, extending the range of this species all the way from the Arctic to the Antarctic. Although species delimitation analyses did not provide any evidence in support of or against the hypothesis of conspecificity of the Sclerolinum contortum populations, the levels of genetic divergence between populations is low in all markers, both mitochondrial and nuclear, which corroborates the findings of Georgieva et al. 25 . We thereby conclude that the conspecificity of the populations of S. contortum is still supported. Although the geographic ranges of N. lokii and S. contortum are, to the best of our knowledge, the widest recorded of any species obligate to chemosynthesis-based ecosystems, wide ranges are not unusual among annelids from CBEs. A polynoid species from the Longqi vent field in the Indian Ocean is also found on the East Scotia Ridge (a range of ∼6000 km), and a phyllodocid species (Hesiolyra cf. bergi) is potentially shared between hydrothermal vents at Longqi, on the East Pacific Rise and the Mid Atlantic Ridge 49 . Two annelid species in the genus Archinome (Amphinomida) appear to be shared between CBEs in several ocean basins, but one of these (Archinome tethyana) has only been sequenced from the Mid-Atlantic Ridge and the other (Archinome jasoni) shows strong genetic structure between the West-Pacific and Indian Ocean sites 50 . Future efforts with genetic analyses of annelids from CBEs around the world is necessary to elucidate how common widespread distributions across several oceans are. In the case of N. lokii, it seems likely that the presence of two sympatric but very distinct mitochondrial lineages of N. lokii in the Antarctic, one of which has very similar haplotypes to those found in the Arctic and Barbados Trench, is the result of two independent colonization events. The genetic divergence between them could be the result of Haplogroup A having evolved independently in the Antarctic for a period of time before the arrival of the second lineage, possibly due to infrequent long-distance dispersal events such as those inferred to occur in the West-Pacific 51 . An alternative explanation is that Haplogroup A dispersed to the Antarctic from an as of yet unknown population of N. lokii that is isolated from the Arctic and Barbados Trench populations. When two divergent lineages come into secondary contact, they will, depending on how far along they are in the speciation process, either remain distinct or merge into one interbreeding population 52 . The presence of a shared 28S haplotype between these two lineages could be the result of either incomplete lineage sorting, or recent geneflow. Species delimitation results support that Haplogroup A in the Antarctic belongs to N. lokii, although this lineage shows levels of COI divergence usually considered to be in the interspecific range 53 . Species delimitation under the multispecies coalescent model assumes that all genetree discordance is caused by incomplete lineage sorting, and the analysis is known to be sensitive to violations of the model assumptions such as population structure or geneflow between species 54,55 . Simulation studies indicate that STACEY is more sensitive to the latter scenario, and may lump species if there has been recent geneflow between two distinct lineages due to secondary contact 54 . Thus, it is possible that the species delimitation results could have been confounded if the two distinct Antarctic lineages are a result of secondary contact between two lineages with incomplete reproductive barriers. To test these hypotheses and resolve whether Haplogroup A may constitute a distinct species from N. lokii would require a higher sampling of more variable nuclear markers.
The Antarctic populations of both N. lokii and Sclerolinum contortum show very low genetic diversity compared to the central Atlantic and Arctic populations. There are several possible explanations for this low diversity, such as recent founder events, population bottlenecks or selective sweeps. It is curious, however, that the same low diversity is evident in both species, and in N. lokii it is evident in both lineages, and across several sampling localities in the eastern Scotia Sea. A population genetic analysis of three other species (a kiwaid squat lobster, a peltospirid gastropod and a leptodrilid limpet) from the eastern Scotia Sea hydrothermal vents did not show a similar depressed genetic diversity 56 , which excludes the possibility of a local event that affected all of the vent fauna. Both S. contortum and N. lokii are able to reproduce asexually by fragmentation and regrowth 28,57 , which may explain low levels of genetic diversity in one locality, as the sampled specimens may be clonal individuals. This is a plausible explanation for the low genetic diversity in S. contortum, as all of the sequenced individuals in the present study is from the same sampling locality 25 .
The pairwise PhiPT values (Table 3) for both S. contortum and N. lokii could indicate greater isolation of the Antarctic population compared to the Arctic and central populations (GoM and Barbados Trench). Apart from the geographic distance, further isolation of the Antarctic populations could be caused by the Antarctic Circumpolar Current, which presents a barrier at around 50°S. The benthic fauna of the Antarctic shows high degrees of endemism, but the isolation is less pronounced in the deep sea than at shelf depths 58 . For N. lokii the high PhiPTs could be an effect of the very divergent Haplogroup A, and the conspecificity of this lineage should be further assessed before making any conclusions based on statistics including this lineage. It is also possible that the low diversity in the Antarctic populations of both N. lokii and S. contortum are inflating the pairwise PhiPT values in comparisons including the Antarctic populations. PhiPT is calculated as the proportion of the variance among populations, relative to the total variance (within population variance + between population variance). Thus, lower within population variation will lead to higher PhiPT values.
Although our results support that both S. contortum and N. lokii have very wide distributions ranging from the Arctic to the Antarctic, the lack of shared haplotypes in COI and high PhiPTs between all populations indicates that present-day levels of geneflow are low. However, the existence of unsampled populations between the currently known populations is highly likely, and sampling gaps can give a false impression of isolation 59 . Nicomache sp. is known from the hydrothermal vents of the Mid-Cayman Spreading Center 60 , and comparison of a COI sequence from that locality (GenBank accession number KJ566962) with our dataset shows that it is 1.5% different from the most similar haplotype of N. lokii (haplotype B from the Barbados Trench). This is within the range of COI-variation considered to be intraspecific for N. lokii in this study. The Pacific species Nicomache arwidsoni has been recorded from the Logatchev vent field on the mid-Atlantic ridge 61 , but given the geographic distribution of N. lokii presented here it is possible that this population could also be N. lokii. In addition, there are both unsampled and probably also unrecorded cold seeps along the continental margins of the Atlantic which could provide suitable habitats for N. lokii and S. contortum.
Although genetic variation was consistently lower in S. contortum than N. lokii, the haplotype networks of S. contortum showed a much stronger geographic structure, which could be interpreted as a higher degree of isolation. Two possible scenarios can be put forward to explain the differences in the degree of geographic structure between the two species: it could either be a result of a difference in dispersal ability between S. contortum and N. lokii, or of the more isolated placement of the GoM seeps compared to the Barbados Trench. Unfortunately, there is very little known about the reproductive biology of either species. Sclerolinum contortum has oval-shaped eggs of up to 430 μm in length and 110 μm in diameter 27 , while the eggs of N. lokii are disc-shaped and up to 240 μm in diameter. This size range is usually associated with lecithotrophic larvae, and this is the prevalent larval mode among the vestimentiferans, the sister group of Sclerolinum 9 . To disentangle the effects of geographic barriers and putative differences in dispersal capacity on the genetic connectivity it would be best to have samples of both species from the same localities, however this was not available for the present study. Given the geographic distance between the populations and lack of shared haplotypes, the low degree of geographic structure in the haplotype networks of N. lokii is peculiar. In the COI network, several of the Barbados Trench haplotypes were more similar to haplotypes from the Arctic, than to other haplotypes in the Barbados Trench. To evaluate possible explanations for this pattern would require a better knowledge of potential intermediate populations between the Arctic and the Barbados Trench. Possible asymmetries in geneflow and time of isolation could be inferred using Isolation with Migration models, but these models assume that there has been no geneflow between unsampled populations and the sampled populations under analysis 62 , which cannot reasonably be assumed for the current dataset.
Both N. lokii and S. contortum show an ability to inhabit a wide range of chemosynthesis-based ecosystems, including various kinds of hydrothermal vents and also cold seeps. Sclerolinum contortum is in addition able to inhabit sunken wood, which could provide stepping stones for dispersal between vents and seeps 63 . Nicomache lokii is able to occupy a wide depth range as well, with the sites sampled in the present study ranging from 1262 m at HMMV to 4930 m at Atalante East in the Barbados Trench. The low habitat-specificity of these two species may be part of the explanation for their wide distributions, as this increases the number of available habitat-patches. Despite the lack of habitat selectivity observed in both N. lokii and S. contortum, there may be within-habitat conditions that favor the establishment of both species together, or of each species individually. At Loki's Castle both species were found in the diffuse flow area of the vent field 30,64 , and N. lokii was also found on the walls of black-smoker chimneys, but there they were much smaller in size 28 . Studies of the habitat of S. contortum at the Hook Ridge vents in the Antarctic indicates it may be limited to low-temperature areas (∼20 °C) at which sulphide flux is also lower 31 , which is approximately the same temperature it occupies at Loki's Castle 23 . Sclerolinum contortum was found to be absent in parts of Hook Ridge influenced by fluid temperatures of up to 49 °C at which siliceous crusts had precipitated 31 , and higher temperatures combined with greater sulphide availability may be a possible explanation as to the absence of this species at the E2 and E9 Antarctic vent sites. Hydrothermal venting at Hook Ridge also appears to be largely ephemeral in nature and characterized by diffuse flow at low temperatures, which has been proposed as a reason for the absence of typical vent fauna 65 . Taxa better adapted to make use of resources below the sediment surface such as S. contortum may be favored at these vents, whereas N. lokii may be excluded from Hook Ridge for reasons of insufficient sulphide availability (and therefore insufficient food) near the sediment surface, and an inability to directly utilise sulphide at greater sediment depths. The presence of N. lokii on the Mid Cayman Spreading Center and putatively on the Mid-Atlantic Ridge, fits well with the hypothesis that this species is adapted to higher temperatures than S. contortum, which is not known from either of these vent sites. These observations support the view that environmental factors varying across different CBEs, such as fluid flux, may play a central role in determining the faunal composition of these habitats 66 .

Conclusions
The results presented here support that both Nicomache lokii and Sclerolinum contortum are distributed all the way from the Arctic to the Antarctic, giving them the widest recorded geographic range of any species obligate to chemosynthesis-based ecosystems, which is also supported by genetic data. In the Antarctic, there are two divergent mitochondrial lineages of N. lokii, which may be the result of two independent colonization events. However, the conspecificity of the most divergent lineage (Haplogroup A) to N. lokii should be reassessed with a higher number of more variable nuclear markers. S. contortum shows a higher degree of geographic structure in the genetic diversity than N. lokii, but whether this is due to intrinsic factors such as dispersal capacity or a reflection of different degrees of isolation of the sampled localities cannot be resolved with the present data. There are probably numerous unsampled populations of S. contortum and N. lokii, and future studies with more extensive sampling throughout the range of these species are needed to understand the genetic connectivity from the Arctic to the Antarctic. Observations on the environmental conditions inhabited by the two species indicate that N. lokii may be limited by flow intensity while S. contortum is probably able to tolerate lower flow conditions, but limited by temperature.