Polychaetes were collected from three locations within the Southern Ocean; the Scotia Arc containing six sampling sites and the Amundsen Sea and Weddell Sea both containing four sampling sites, Figure 1. All polychaetes were identified based on morphological characters, some individuals could not be identified to any currently described species and were considered new to science. These results are presented in Neal et al. (2017) and all voucher specimens were deposited in the Natural History Museum London, details can be found online using the Darwin Core Archive https://doi.org/10.5519/0068114. DNA was extracted and sequenced from nine polychaete morphospecies (211 specimens in all); for details on the phylogenetic analysis of these sequences including the identification of cryptic species see Brasier et al. (2016). Prior to DNA barcoding, a biogeographic distribution classification; cosmopolitan, circum-Antarctic or restricted, was allocated to each initial morphospecies analyzed (Table 1), where: “cosmopolitan” species are those that have been recorded throughout the majority of the world's oceans and both hemispheres; “cicrum-Antarctic” species are those that have been collected within different regions of the Southern Ocean and are considered to be widespread around the Antarctic and “restricted” species are those only recorded in one Antarctic region or location with no records within the Register of Antarctic Marine Species (RAMS) recorded or Basic Local Alignment Search Tool (BLAST) matches from other localities.

These definitions are coherent with the most recent biogeographic review of Antarctic polychaetes, Schüller and Ebbe (2014), who considered species circum-Antarctic if there were georeferenced RAMS records from at least the Weddell Sea, West Antarctic Peninsula (WAP) or the Scotia Arc, and the Ross Sea or Eastern Antarctica. The presumed distributions in Table 1 were defined using taxonomic databases including the Encyclopedia of Life (EOL) and RAMS. For undetermined morphospecies, their distribution was deciphered by their presence within samples collected in the each of the regions as well as any BLAST matches with georeferenced sequence including larval matches within the Ross Sea (Heimeier et al., 2010; Gallego et al., 2014).

The genetic diversity of the DNA barcodes generated from the nine morphospecies was investigated in Brasier et al. (2016). Phylogenetic and distance analyses found evidence for twelve cryptic species, increasing the number of species in this study from nine to seventeen. In Brasier et al. (2016) clades were assigned MB# to distinguish species, including cryptic species, e.g., Scalibregma sp. MB1 and Scalibregma sp. MB2. Cryptic species were referred to as either “Genus cf species MB#,” “Genus sp. MB#” or if no genus could be assigned “Family sp. MB#.” The use of “cf” or “sp” was associated with the type locality of the original species and likelihood of the Antarctic specimens being the “true” species, for more details see Brasier et al. (2016). Using haplotype networks we have analyzed the distribution of each genetically identified species, morphological and cryptic, listed in Table 1.

As the recovery of 16S sequences was greater than COI in all morphospecies examined, this marker was used to construct georeferenced haplotype networks of the sequenced specimens to visualize species distributions and speculate biogeographic patterns. An exception was Hesionidae sp. A for which cryptic species were revealed in COI but not 16S. If sequence matches within the GenBANK database were found during phylogenetic analysis these were also included in the network (Table 2). To avoid problems with gaps, all sequences of the same species were trimmed in Mesquite (Version 2.75) to the same length (Table 1) following MAFFT (for 16S) or MUSCLE (for COI) sequence alignment in Geneious (R7). GPS coordinates were assigned to each sequence for its given sample location, and networks constructed using statistical parsimony (Templeton et al., 1992) and the TCS programme (Clement et al., 2000) in PopART for editing in CorelDRAWX8. For the number of sequences per species by site used in the georeferenced haplotype networks, please see Supplementary Material. Haplotype networks with depth referenced sequences were also created using the five depth bins; < 500 m, 500–1,000 m, 1,001–1,500 m, 1,501–2,000 m and >2,000 m. The depth bins were chosen based on the discrete depth horizons sampled during collection expeditions (for details see Griffiths et al., 2008 and Linse et al., 2013). Most sequenced specimens were obtained from depths of 500, 1,000, 1,500, or 2,000 m, some specimens were collected from 100 and 200 m and were included in the < 500 m depth bin and those at depths greater than 2,000 m were also binned together.

To gain an insight into the potential pathways of larval dispersal between the sampled localities a particle tracking model was employed (ARIANE; Blanke and Raynaud, 1997) to the output of a NEMO (Madec, 2008) ocean general circulation model (OGCM). Five day mean ocean current data obtained from the National Oceanography Centre, U.K., NEMO 1/12° OGCM (Duchez et al., 2014), for the period 2000-2009, was used to calculate the 3D trajectories of passive particles released from 17 sites from four Antarctic locations, Table 3. Within each model grid cell co-located with a sampled site, eight evenly distributed particles were released at every model depth level and at 5 day intervals for the period 2000-2009. These particles were then tracked for 1 year. No vertical or horizontal dispersion was added to the particle motion and no buoyancy terms were included, so the subsequent particle pathway is purely an advective one. A 75 km search radius or “trap” at each sample site location for the full water column was used to determine whether foreign particles passed close by and were therefore potentially connected to another sample site.

Particle tracking analysis is relevant to polychaete reproductive modes because six of the eight families included in this study are considered to have pelagic larvae (Table 4), with five of these known to have feeding (plankotrophic) larvae. The reproductive traits of each species examined in this study were inferred from family level trait data available from the largest polychaete trait database Polytraits¹. Species level reproductive traits are rare, especially for those that are the hardest to collect and identify, e.g., deep sea and Antarctic species. We know from genetic studies that three of our species, Hesionidae sp. (MB1), Aglaophamus trissophyllus and Laonice weddellia have pelagic larvae in Antarctic waters (Heimeier et al., 2010; Gallego et al., 2014).

As discussed in our interpretation several aspects of larval biology were not incorporated into this model including larval mortality and behavior. These factors can lead to increased rentention of larvae within their source location and thus it is possible that our observations of the distance and direction of dispersal are overestimated (e.g., Cowen et al., 2006; Levin, 2006). Additionally we have not examined the seasonal differences in the release of larvae. Although some studies have found evidence for seasonal reproductive cycles in benthic fauna, different larval stages have been observed within the water column through the year (Bowden et al., 2009; Sewell and Jury, 2011). In this study particle tracking analyses are used purely to gain insight into how the oceanography of the Southern Ocean may theoretically affect population connectivity by its directionality and distance of dispersal. Furthermore, we do not know enough about Antarctic polychaete larvae to estimate these behavioral or biological model constraints.

The most common distribution observed was widespread Antarctic, with 71% of the species investigated recorded to have genetically similar specimens in at least two Antarctic locations. The most abundant were Laonice weddellia, Harmothoe fuligineum and Aricidea simplex, Figure 3. A. simplex was sequenced from all three Antarctic locations sampled. The distribution of L. weddellia was also extended to the Ross Sea by identical sequences to that of a larval specimen from the Ross Sea (Gallego et al., 2014). No specimens of H. fuligineum were identified from the Weddell Sea, but specimens were found within the Amundsen and Scotia Arc regions.

Most of the widespread Antarctic cryptic species were sequenced from two of the sampled locations, e.g., the Scotia Arc and Amundsen Sea [Glycera sp. (MB1) and (MB2), Scalibregma sp. (MB3)]. Aglaophamus sp. (MB2) was the only cryptic species sequenced from all three locations. Aglaophamus sp. (MB1) was only sequenced from specimens in the Scotia Arc but was recorded as widespread because sequenced specimens from the Scotia Arc matched larval sequences from the Ross Sea (Heimeier et al., 2010), Figure 4A. Further larval matches from the Ross Sea were found with Hesionidae sp. (MB2). In comparison to Hesionidae sp. (MB2) the distribution of Hesionidae sp. (MB1) was reduced with representation from the Scotia Arc (Southern Thule) and the Amundsen Sea only, Figure 4B.

Whilst both Scalibregma sp. (MB1) and Eurphrosinella cf cirratoformis (MB1) were collected from more than one Antarctic region, their cryptic counterparts Scalibregma sp. (MB2) and (MB3) and Eurphrosinella cf cirratoformis (MB2) were only sequenced from one location including; Southern Thule, Elephant Island and the Amundsen Sea respectively, Figure 2B. Scalibregma sp. (MB2) was only represented by a single individual but Scalibregma sp. (MB3) and Eurphrosinella cf cirratoformis (MB2) were represented by four and three individuals respectively. Although these clades appear more restricted compared to their cryptic counterparts they exist sympatrically within the same regions. The only cryptic species that appear to exist allopatrically are the three clades of the morphospecies Aricidea belgicae. Each species was sequenced from a single separate region, Aricidea cf belgicae (MB1) from the Amundsen Sea, Aricidea cf belgicae (MB2) from the Weddell Sea and Aricidea cf belgicae (MB3) from the Scotia Arc, Figure 3C.

The maximum depth range that could be recorded was 100–2,000 m. The greatest depth distribution recorded in this study was found in Laonice weddellia, which were collected at depths between 200 and 1,500 m. Aricidea simplex had a similar depth range but the deepest specimens were collected from 1,000 m, whilst Harmothoe fuligineum had the narrowest depth range for the non-cryptic species from 200 to 500 m. Figure 3.

Some patterns in the depth distributions of cryptic species were observed. For example Glycera sp. (MB1) was collected from 200 to 500 m and Glycera sp. (MB2) from 500 to 1,000 m depth, Figure 1A. Similar patterning was also observed between Hesionidae sp. (MB1) and Hesionidae sp. (MB2), Euphrosinella cf cirratoformis (MB1) and Euprhosinella cf cirratoformis (MB2) and, Aglaphamus cf trissophyllus (MB1) and Aglaophamus sp. (MB2), Figures 3, 4. Two of the cryptic species of Scalibregma, (MB1) and (MB3) exhibited the same depth distribution, 200–500 m, whilst the single representative of Scalibregma sp. (MB2) was only found at 1000 m, Figure 2B. For the allopatric cyprtic caldes of Aricidea belgicae, it was Aricidea cf belgicae (MB1) that exhibited the greatest depth distribution of 500–1,000 m, whilst Aricidea cf belgicae (MB2) and Aricidea cf belgicae (MB3) were only observed at 200 and 600 m respectively, Figure 3B.

Particles released within the Ross Sea shelf and slope sites were transported westward around the Antarctic continent in the counter current and thus away from the other sample locations, Figure 5. In contrast, particles released within the offshore site within the Ross Sea were carried eastward in the direction of the WAP but away from the continental slope in the ACC. Sites within the Amundsen Sea followed the same direction of dispersal to the west with those on the edge of the continental shelf (AS_BIO3 and AS_BIO6) transporting particles toward the Ross Sea sites, Figure 6 (AS_BIO3). Particles from release sites at the northern tip of the WAP (Elephant Island, Livingston Island and Powell Basin) were transported in both directions, but with greater dispersal westwards across the Scotia Sea reaching the Shag Rocks and South Georgia sites, Figure 7. Particles from Shag Rocks and South Georgia, were also transported east and northward being advected by the ACC. Particle release sites within the Weddell Sea were well connected by the Weddell Sea gyre however, the extent of dispersal varied between sites. For the inshore locations (WS_FTSE and WS_FSCE) particles were transported westward and up the WAP (Figure not shown). Whilst those released further offshore (WS_CS; Figure 8) reached further into the Scotia Sea and into the ACC, crossing though the Shag Rocks and South Georgia sites as well as those closer to the WAP: Powell Basin, Elephant Island and Livingston Island.

Statistical analysis in Figures 9, 10 show the percentage of particles passing through a 75 km radius or “trap” surrounding each of the other sites, and their mean time and depth of arrival. As expected, particles traveling the greatest distances were found in the surface layers as a result of the stronger advection (indicated by the lighter shaded squares in Figures 9, 10 and the lower percentages). In general, connectivity within regions was greater than between regions, for example 60% of particles from Elephant Island passed within 75 km of Livingston Island. However, this connectivity was unidirectional as only < 1% of particles released from Livingston Island were observed in the vicinity of Elephant Island. The level of inter- and intra-regional connectivity is highly dependent on the release site location in relation to local current pathways, and to a lesser degree the physical distance between the release sites. For example, there appears to be no connectivity between AS_BIO6 and AS_BIO4/5, yet the Ross Sea sites receive particles from all Amundsen Sea sites. Southern Thule was the least connected site receiving particles < 1% of particles from the Powell Basin site. Overall, the WS_CS had the highest connectivity to sites outside of its region with particles reaching all sites within the Scotia Arc.

Matching haplotypes between individuals of the same species from at least two of the sampled regions as well as matches with larval DNA from the Ross Sea (Heimeier et al., 2010; Gallego et al., 2014) suggests an abundance of widespread polychaete species. Our results suggest that the three morphospeices, Laonice weddellia, Harmothoe fuligineum, and Aricidea simplex, are widespread within the West Antarctic, with potentially circum-Antarctic distributions. Similar results were also found for the Dorvilleidae polychaete Ophryotrocha orensanzi by comparing COI data from west and east Antarctic populations (Paxton et al., 2016). Eight of the cryptic species sequenced in this study appear to be widespread occurring in multiple regions within the West Antarctic. The matching haplotypes between Antarctic regions may indicate continued dispersal and ongoing gene flow between these regions that maintain genetic similarity (Arango et al., 2011; Soler-Membrives et al., 2017). These results are supported by the model findings, documenting the potential for oceanographic currents to carry particles or “larvae” between Antarctic regions in both directions around the continent.

The majority of polychaete cryptic species studied exist sympatrically, covering the same or overlapping regions of the West Antarctic. How these species and their current distributions were established and maintained is most probably a result of several complex biophysical interactions over geological time. Their evolution and distribution is likely to be influenced by the relative roles of vicariance and dispersal (Clarke, 1992; Aronson et al., 2007; Waters, 2008; González-Wevar et al., 2011). Although many records of cryptic species in the Antarctic are restricted to certain locations, widespread cryptic species have been recorded in several taxa. Examples of which include: the crinoid Promachocrinus kerguelensis, where six genetically distinct phylogroups were considered circumpolar, sympatric and eurybathic (Hemery et al., 2012); the brittle star Ophionotus victoriae, and the amphipod Nymphon australe, where although some clades appear restricted others were considered widespread (Galaska et al., 2017; Soler-Membrives et al., 2017). Within the literature similar insights into larval dispersal have also been obtained from oceanographic observations in conjunction with genetic analysis. For example Matschiner et al. (2009) concluded that the lack of genetic structuring between populations of the notothenioid, Gobionotothen gibberifrons, throughout the Scotia Sea could be assigned to the passive transport of their larvae during the pelagic development phase by the ACC as indicated by surface drifter trajectories. The geographic distribution of marine taxa within the Southern Ocean reflects the species life-history traits (including bathymetric ranges, developmental modes and larval lifespans) and the influence of circum-Antarctic current systems (Young et al., 2015; González-Wevar et al., 2017; Moreau et al., 2017).

It is generally considered that the majority of cryptic species within Antarctic waters arose from physically separated populations over multiple glaciations (Allcock et al., 2001; Thatje et al., 2005; Galaska et al., 2017). Thus, the existence of sympatric species could suggest they evolved within the same area due to another method of isolation, for example differences in reproductive traits (Palumbi, 1994), responses to competition (Alizon et al., 2008) or predation (Wilson et al., 2013). An alternative explanation is that these species did evolve in geographic isolation and their widespread distribution was established after they evolved, possibly during post-glacial re-colonization. Given the limited size and mobility of these benthic polychaetes it is unlikely that adult specimens migrate between the three regions sampled. Instead, it is probable that genetic connectivity between adult populations is maintained by larval dispersal (Fraser et al., 2012).

The reproductive traits of polychaetes can vary on many taxonomic levels. For example within the family Polynoidae the brooding of eggs is considered to be relatively rare (Giangrande, 1997). However, Gambi et al. (2001) recorded brooding of eggs under dorsal elytra in three Antarctic Harmothoe species. Within Spionidae, different populations of Polydora species are known to have different PLDs and feeding modes, exhibiting adelphophagy, the production of both planktotophic and adelphophagic larvae (where unfertilized eggs are ingested by the developing larvae) (Blake, 1969; Mackay and Gibson, 1999). Other spionids such as Streblospio benedicti exhibit poecilogony, producing both planktotrophic and lecithotrophic larvae (Levin, 1984). Given the recorded variation in polychaete reproduction, the reproductive modes of the species studied may not reflect their generalized family level trait data (Table 4). However, many species may still have a larval dispersal phase (Blake and Arnofsky, 1999; Faulwetter et al., 2014). Even brooding species can still possess a dispersal phase as brooders can release pelagic planktotrophic larvae as indicated by egg size in Gambi et al. (2001). Thus, it is possible that the widespread distribution of polychaetes observed here, were established and are maintained by larval dispersal by circum-Antarctic oceanographic currents (Starmans et al., 1999).

Earlier studies of eurybathy suggested that the depth distributions of Antarctic polychaetes were comparable to that of European species. Thus, polychaetes did not confirm to the general “eurybathic” characteristics assigned to Antarctic taxa due to reduced physical changes with depth, which could influence physiology (Brey et al., 1996). In most oceans there is a noticeable change in faunal composition on the shelf break (Gage and Tyler, 1991). However, in Antarctica where the continental shelf is much deeper the change in species composition with depth occurs at about 2,000 m (Brandt et al., 2007). Sequenced specimens from the Amundsen Sea and Scotia Arc were collected from depths of 500–1,500 m and the additional samples from the Weddell Sea were collected between depths of 400–2,000 m. Compared to studies such as Brandt et al. (2007) the depth-range sampled may not be large enough to visualize any depth dependence in distribution patterns. Thus, for those species that appear eurybathic within this study may actually be a result of the limited depth range (up to 2,000 m). If deeper slope and abyssal communities were included observed depth related patterns may have been different. However, Neal et al. (2017) observed depth patterns within the range of this study. Neal et al. (2017) noted greatest similarity in polychaete community composition between 500 m stations on the inner and outer shelf of the Amundsen sea, when compared to the communities from the same station at 1,000 and 1,500 m depths.

In conjunction with their geographic distribution many cryptic species appear to exist sympatrically or have overlapping depth distributions. For example, Hesionidae sp. (MB1) was sampled from depths of 500 m or more, whereas Hesionidae sp. (MB2) was only collected at 500 m depth or shallower. It is possible that Hesionidae sp. (MB2) is more dominant at shallower depths, whilst Hesionidae sp. (MB1), is outcompeted. This interaction could then be reversed at the deeper sites and is potentially associated with functional differences between cryptic species. The true absence of species at certain depths, i.e., truly restricted species, are very difficult to confirm and in some cases are questioned by later studies. For example, Schüller (2011) described three cryptic clades of Glycera spp. from the Weddell Sea, two clades were thought to be restricted to 2,000 m. This is in contrast to this study, as one of the clades from Schüller (2011) matched the Glycera sp. (MB2) specimens that were sequenced from stations 500 m and shallower.

Previous studies uncovered restricted cryptic populations in species from several phyla. For example Held (2003) uncovered two cryptic clades of Ceratoserolis trilobitoides, one of which was found only on the WAP the other extending to the Weddell Sea. Furthermore, Linse et al. (2007) compared DNA sequences of the widely distributed bivalve Lissarca notorcadensis and uncovered four explicit haplotype groups within the Scotia Sea. Similar results have been recorded for the isopod Glyptonotus antarcticus (Held and Wägele, 2005) as well as the cephalopod genus Pareledone (Allcock et al., 2011). This high diversity has previously been assigned to the fragmented nature and limited accessibility of available habitats in this region favoring speciation by population fragmentation (Allcock et al., 2011), especially in species with limited dispersal capacities (Strugnell and Allcock, 2013). In comparison to pelagic species there is generally a higher level of genetic structuring in benthic invertebrates, in Rogers (2007) this was assigned to the lower dispersal capabilities of benthic species at both larval and adult life stages. Within this study some of the cryptic clades of presumed circum-Antarctic species were only found within a limited number of sampled stations within a single region. This was the case for Scalibregma sp. (MB2) and (MB3), and Euphrosinella cf cirratoformis (MB2) that were more restricted than their broadly distributed, potentially circum-Antarctic, sister cryptic clades. Most of the restricted clades listed above were limited to sites within the Scotia Arc, an area known for particularly high biodiversity within Antarctica (Linse et al., 2007; Allcock et al., 2011; Neal et al., 2017). However, we acknowledge that restricted species are difficult to define, as all these clades contained fewer individuals (< 4) this result could be an artifact of undersampling.

The three potential cryptic species of Aricidea cf. belgicae (MB1), (MB2), and (MB3) were each restricted to a single region; Scotia Arc (MB1), Amundsen Sea (MB2) and Weddell Sea (MB3). This is the only example within this study of allopatric cryptic species if their distribution evolved from geographically isolated populations and was maintained by limited dispersal capabilities. Again we do not know the reproductive mode of these clades and Paraonidae are known to undergo both direct and indirect development making it difficult to speculate their dispersal (Table 4). Furthermore, even broadcasting species may not be able to overcome geographic isolation (Galaska et al., 2017). As mentioned these “restricted” distsributions could again be an artifact of undersampling as clades (MB2) and (MB3) were only represented by two and one individuals respectively. Within undersampled environments it is difficult to determine the presence of both restricted and rare species (Smith et al., 2006). This is complicated further by the presence of cryptic diversity and thus the submission of DNA sequences to open access databases is extremely important in the future assessment of species biogeography.

Here we use the model results presented to speculate on the direction and potential distances of passive larval dispersal within the Southern Ocean from our sampled sites. Particles released in the Amundsen Sea were transported counter clockwise around the continent. This transport is likely to be due to the Antarctic coastal counter current. Topographically constrained by the Antarctic Slope Front (ASF) this westward current encircles the coastal margins of Antarctica (Thompson et al., 2009). The ASF is found consistently above or just offshore of the shelf break and the coastal counter current is found broadly over the continental shelf (Heywood et al., 1998). While the ACC is considered the dominant current in maintaining circumpolar connections, the coastal counter current has been shown to connect many regions of high krill density (Thorpe et al., 2007). In the Ross Sea, particles released on the shelf were tracked westward in the coastal counter current; whereas at the offshore location, the Ross Sea gyre and ACC advect the particles eastward. Releases on the Ross Sea slope appear to be influenced by both regimes, but with the majority of particles following the coastal counter current pathway. The ACC provides a pathway to connect regions in the opposite direction to the coastal counter current; particles released within the Scotia Arc transported eastward out into the Scotia Sea. This movement is likely to be driven by the ACC meandering as they are transported further off shelf (Young et al., 2015). Similar particle movements have been observed in Lagrangian tracking models used to estimate krill connectivity (Hofmann and Murphy, 2004; Piñones et al., 2013). Particle dispersal from the Weddell Sea sites is dictated by the cyclonic Weddell Sea gyre. Releases above the continental slope travel the farthest distance, with the majority of particles tracking eastwards once reaching the Scotia Arc. Sites closer to the coast sites within the Weddell Sea see a slower spreading toward the WAP, with particles advected both eastward toward the Scotia sites and westward in the coastal counter current.

The passive movement of particles around Antarctica within the PLD time frame indicates the possibility that larvae may be recruited into non-parent populations substantial distances from their origin. If larvae are able to settle, grow and reproduce in these locations this could maintain genetic connectivity and limit the potential for genetic differentiation between regions. Examples of source populations transporting larvae westward may include those within the Amundsen Sea supplying larvae to the Ross Sea and the Ross Sea shelf transporting larvae along the continental shelf. The Weddell Sea sites appear to be well connected and may provide larvae to sites along the WAP such as Livingston Island and Elephant Island, which in turn may supply larvae to northern sites within the Scotia Sea. However, it should be noted that given the observed depth ranges of these species (0–1,500 m) many of the particles advected over the deep open ocean may no longer be in suitable locations to establish populations, restricted by their physiology. If larvae are unable to disperse over unsuitable habitats this can prevent the degree of population connectivity (Rogers et al., 2006). This may be the case for dispersal from sites within the Scotia Arc, in which particles were transported north into the Scotia Sea. These sites included Shag Rocks and South Georgia but also Southern Thule from which no particles reached other locations. Thus, these sites may act as sink populations rather than source other regions.

The maintenance of widespread distributions by pelagic larval dispersal and recruitment may not applicable to all of the species, as we cannot confirm the presence of free-living larvae. Furthermore, Antarctic taxa in general are considered to lack free-living larval stages (Pearse et al., 1991, 2009) and the existence of circum-Antarctic brooding species is highly unlikely (Lörz et al., 2009). A potential explanation for the maintenance of genetic connectivity in brooding species or those lacking pelagic larvae includes the passive rafting of larvae or adults on floating substrate or ocean debris (Waters, 2008). Leese et al. (2010) suggested that this method maintained connected shallow water isolated populations of the isopod Septemserolis septemcarinata across the sub-Antarctic. Furthermore, direct evidence of rafting on kelp has been observed in the widespread sub-Antarctic brooding bivalve Gasimardia trapesina (Helmuth et al., 1994), the sea slug Onchidella marginata (Cumming et al., 2014) and two species of sub-Antarctic amphipods (Nikula et al., 2010). An additional factor that could influence population connectivity is anthropogenic transport (David and Loveday, 2017). We have limited knowledge of the extent or potential growth of this in the Southern Ocean however there is evidence of anthropogenic transport of non-native species into the Antarctic (Lee and Chown, 2007).

Given inherent difficulties of directly measuring dispersal when larvae are minute (~200 μm) compared to the potential scale of dispersal (~km) (Gilg and Hilbish, 2003) dispersal distance is more often estimated using coupled biophysical models. The model scenario used in this study was based limited biological traits data. Limited biological knowledge is considered to be the main challenge when attempting to predict and validate dispersal pathways and distance (Levin, 2006; Metaxas and Saunders, 2009; Hilário et al., 2015). In this study we tracked particles for a maximum of 1 year, however PLD can be highly variable, for example recorded PLD in polychaetes has ranged from 13 to 150 days for planktotrophic larvae and 1–25 days for lechitrophic larvae from California (Carson and Hentschel, 2006). Within the Southern Ocean region, the life histories of marine organisms are often much slower than similar temperate and tropical taxa (Pearse et al., 1991), for example the lifespan of Antarctic echinoderm larvae is thought to exceed 1 year at a temperature of −1.5°C (Shilling and Manahan, 1994; Marsh et al., 1999).

If the PLD of the species studied is less than 360 days, the sites reached by particles may be an overestimate of larval dispersal. Figure 9 presents the mean time of arrival of particles and can be used to estimate the distance of dispersal over shorter PLDs. With a shorter larval duration of 3–6 months (90–180 days in Figure 9) most particles will only reach locations within their region. For example most sites within the Weddell Sea would remain connected but particles would not reach the sites in the Scotia Arc. A shorter PLD would greatly increase particle retention within the site or region of release therefore reducing the potential connectivity between regions (e.g., Shanks, 2009; Faurby and Barber, 2012).

As well as PLD our knowledge of polychaete larval behavior is very limited. Matschiner et al. (2009) showed that many larvae are capable of active vertical migration, which could lead to the avoidance of advection and increased retention near their sight of dispersal (Swearer et al., 2002). Additionally, in our study, particles were only released from sites we sampled. The existence of populations in between our sites would be important if the PLD and dispersal distance is shorter than predicted. These in-between populations would provide additional sites for genetic mixing between populations and larval release, and so contribute to the maintenance of circum-Antarctic genetic connectivity.

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.