Distributional Patterns of Polychaetes Across the West Antarctic Based on DNA Barcoding and Particle Tracking Analyses

Brasier, Madeleine J.; Harle, James; Wiklund, Helena; Jeffreys, Rachel M.; Linse, Katrin; Ruhl, Henry A.; Glover, Adrian G.

doi:10.3389/fmars.2017.00356

ORIGINAL RESEARCH article

Front. Mar. Sci., 16 November 2017

Sec. Marine Evolutionary Biology, Biogeography and Species Diversity

Volume 4 - 2017 | https://doi.org/10.3389/fmars.2017.00356

This article is part of the Research Topic Antarctic Biology: Scale Matters View all 21 articles

Distributional Patterns of Polychaetes Across the West Antarctic Based on DNA Barcoding and Particle Tracking Analyses

$\r\nMadeleine J. Brasier,*$ Madeleine J. Brasier^1,2^*

Helena Wiklund²

¹School of Environmental Science, University of Liverpool, Liverpool, United Kingdom
²Life Sciences, Natural History Museum, London, United Kingdom
³National Oceanography Centre, University of Southampton, Southampton, United Kingdom
⁴BioSciences, British Antarctic Survey, Cambridge, United Kingdom

Recent genetic investigations have uncovered a high proportion of cryptic species within Antarctic polychaetes. It is likely that these evolved in isolation during periods of glaciation, and it is possible that cryptic populations would have remained geographically restricted from one another occupying different regions of Antarctica. By analysing the distributions of nine morphospecies, (six of which contained potential cryptic species), we find evidence for widespread distributions within the West Antarctic. Around 60% of the cryptic species exhibited sympatric distributions, and at least one cryptic clade was found to be widespread. Additional DNA barcodes from GenBank and morphological records extended the observed range of three species studied here, and indicate potential circum-Antarctic traits. Particle tracking analyses were used to model theoretical dispersal ranges of pelagic larvae. Data from these models suggest that the observed species distributions inferred from genetic similarity could have been established and maintained through the regional oceanographic currents, including the Antarctic Circumpolar Current (ACC) and its coastal counter current. Improved understanding of the distribution of Antarctic fauna is essential for predicting the impacts of environmental change and determining management strategies for the region.

Introduction

Enclosed by both the Antarctic Circumpolar Current (ACC) and frontal systems the Southern Ocean is often described as an isolated marine environment, Figure 1. These oceanographic features act as a physical barrier, which are thought to have prevented species movement into and out of the Southern Ocean. For this reason early predictions suggested that the majority of benthic fauna within Antarctic waters would be endemic to the Southern Ocean (Ekman, 1953; Hedgpeth, 1969). Endemism has since been observed in many major taxonomic groups based on species records from mostly morphological species identification (see reviews Dell, 1972; Arntz et al., 1997; Clarke and Johnston, 2003; Thorpe et al., 2007; De Broyer and Danis, 2011; Brandt et al., 2012; Kaiser et al., 2013). Current estimates of endemism, as reported in the recently published Southern Ocean Biogeographic Atlas, ranged from 50 to 97% between taxa (De Broyer et al., 2014).

FIGURE 1

Figure 1. The locations from which polychaetes were collected during the BIOPEARL and JR275 expeditions (black dots) as well as the locations of additonal GenBank sequences used in haplotype networks (gray dots). SR, Shag Rocks; SG, South Georgia; ST, Southern Thule; PB, Powel Basin; EI, Elephant Island; LI, Livingston Island; AS, Amundsen Sea; WS, Weddell Sea; RS, Ross Sea; RS_O, Ross Sea offshore. Position and direction of oceanographic currents including the Antarctic Circumpolar Current (ACC), East Wind Drift (counter current) (EWD) and the Weddell and Ross Sea Gyres. Position of the Polar Front (PF), the Southern Antarctic Circumpolar Current Front (SACCF), Subantarctic Front (SAF) and the Southern Subtropical Front (SSTF).

A further biogeographic pattern that has long been associated with Antarctic marine fauna is the circum-polarity or circum-Antarctic distributions (Arntz et al., 1997; Clarke and Johnston, 2003). These distributions could be a result of a combination of factors that could have a homogenizing effect on the faunal communities. The continuous coastline around the continent itself provides connectivity between the different seas around Antarctica. Furthermore, given the relatively uniform physical conditions across the continental shelf, individual settlement and survival is not restricted by their physiology, e.g., temperature tolerance (Arntz and Gallardo, 1994). Larval dispersal around the continent could be aided by oceanographic currents, including the ACC, its coastal counter current (also referred to as East Wind Drift), and the Weddell and Ross Sea gyres, Figure 1 (Orsi et al., 1993, 1995; Fahrbach et al., 1994; Linse et al., 2007).

Antarctic benthic fauna often have extended depth ranges and are considered to be eurybathic, i.e., capable of living in benthic habitats within both shallow and deep water. Some of the first suggestions for this were made by Dell (1972) and Knox and Lowry (1977) for various taxa including sponges, corals, polychaetes and molluscs. Broad depth distributions of Antarctic fauna are thought to be associated with the advance, and retreat of shelf ice during interglacial cycles (Brey et al., 1996). During periods of glacial expansion some “shelf fauna” may have been moved down slope into ice free habitats. In the following glacial retreat the now “slope fauna” could then recolonize shallower shelf areas thus evolving eurybathic distributions (Clarke and Johnston, 2003). This movement was possible due to similar physical conditions (e.g., temperature) on the shelf, slope and deep-sea floor, thus reducing the need for specific adaptations needed to survive in these environments (Clarke et al., 2009; Clarke and Crame, 2010).

Some of these generalized patterns have been challenged when it comes to the deep-sea benthos. Although the ACC and Polar Front may affect the movement of pelagic and shallow water species it may not be a barrier to the benthos (Clarke, 2003). Antarctic polychaetes seem to have some of the broadest geographical distribution ranges amongst the Antarctic benthic macrofaunal invertebrates (Schüller and Ebbe, 2007). During analysis of the ANtarctic benthic DEEP-sea biodiversity (ANDEEP) samples, it was noted that more than half of all polychaetes identified matched species found north of the Polar Front, 20% of which are also found in the northern hemisphere (Brandt et al., 2007). However, in recent works by Neal et al. (2017) depth was identified as the main factor structuring the polychaete communities within the Amundsen and Scotia Sea, contradicting the broad depth distributions often associated with Southern Ocean benthic marine fauna.

Large-scale sampling programs and species databases are providing valuable insight to the biogeography of many species, however, these analyses are predominantly based on morphological identifications. Furthermore, early identifications of Antarctic benthos reflect the species names on pre-existing monographs of, for example, European fauna. Given the abundance of undescribed deep-sea species, the number of specimens collected on large sampling programs and the potential for identifying features to be damaged on collection we use molecular taxonomy to validate the accuracy of biogeographic records. Schüller and Ebbe (2014), investigated all georeferenced Register of Antarctic Marine Species (RAMS) polychaetes within the Scientific Committee for Antarctic Research Marine Biodiversity Information Network (SCAR-MarBIN), concluding that while Southern Ocean polychaete taxonomy is improving with new species descriptions, there are still many “cosmopolitan” species which are most probably Antarctic species different to their Northern counterparts. This can be easily tested using DNA barcoding where comparative sequences exist (Brasier et al., 2016).

Our ability to define the geographic distribution, or biogeography, of species including how their distribution was established and is maintained, has also progressed with the development of DNA sequencing, phylogeography and population genetics (Riesgo et al., 2015). These methods allow us to visualize and, with sufficient sample numbers, calculate the level of population connectivity between known populations of species from different localities, which can be controlled by several interacting biological, physical and chemical factors. Since the application of genetics in Antarctic diversity and biogeographic studies many species presumed to be circum-Antarctic or cosmopolitan based on morphological analysis are actually comprised of more geographically restricted and separated cryptic clades. For example genetically distinct restricted populations of the isopod Betamorpha fusiformis (Raupach et al., 2007), the cephalopod Pareledone spp. (Allcock et al., 2011) and the crinoid Promachorinus kerguelensis (Wilson et al., 2007) have been identified by comparing mitochondrial DNA. In other cases, DNA barcoding has provided evidence for circum-Antarctic distributions, e.g., genetic homogeneity within benthic invertebrates has been recorded in the nemertean ribbon worm Parabolisa corrugatus (Thornhill et al., 2008), the two shrimp species Chorismus antarcticus and Nematocarcinus lanceopes (Raupach et al., 2010) and the pycnogonid Nymphon australe (Arango et al., 2011).

This study uses mitochondrial DNA barcoding data from one of the largest Antarctic polychaete barcoding projects to date (Brasier et al., 2016), to investigate the distributions of nine polychaete morphospecies previously considered to be cosmopolitan or circum-Antarctic. We combine our genetic analyses with particle tracking models to examine the directionality and distance of pelagic polychaete larvae and the potential to maintain genetic connectivity across widely distributed populations. Together with our knowledge of Southern Ocean glaciations and polychaete larval biology we aim to improve our understanding of the biological and physical factors influencing the distribution of polychaete species within the Southern Ocean. We hypothesize that although many morphospecies of polychaetes exhibit “characteristic” distributional patterns for Southern Ocean species including endemism, eurybathy and circum-Antarctic distributions (Griffiths et al., 2009; Brandt et al., 2012), their cryptic clades may be more restricted than their moprhospecies, i.e., cryptic clades are found within a single region of Antarctica, which may be an artifact of their evolution during glacial periods when populations were isolated from one another.

Methods

Specimen Selection and Presumed Distribution

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.

TABLE 1

Table 1. Species for which the geographic distribution will be investigated including their primary identification based on morphological analysis alone (in Neal et al., 2017) and their secondary identification from combined DNA and morphological analysis (in Brasier et al., 2016).

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).

DNA Sequencing and Analysis

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.

TABLE 2

Table 2. Additional sequences included in the haplotype networks, including the MB species they matched, the gene sequenced, specimen locality, GenBANk names and accession number.

Particle Tracking Model to Estimate Larval Dispersal

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.

TABLE 3

Table 3. The release sites used in particle tracking analysis by location with latitude and longitude position and the depth range in which particles were released.

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).

TABLE 4

Table 4. Reproductive traits of the 8 polychaete families containing species from which DNA barcodes were collected in this study.

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.

Results

The distribution patterns of seventeen genetically-determined polychaete species are based on the observed distributions from our haplotype networks and GenBank comparisons. In this study we found no evidence for cosmopolitan species. The genetic difference between our Antarctic specimens of Glycera capitata and Scalibregma inflatum, both considered to be cosmopolitan, with publically available sequence data from their Northern representatives [which were closer to their type localities, Greenland (G. capitata) and Norway (S. inflatum)], deemed the Antarctic clades to be cryptic species and so are examined for widespread Antarctic distributions, Figure 2.

FIGURE 2

Figure 2. Georeferenced (color) and depth binned (grayscale) haplotype networks for (A) Glycera spp. and (B) Scalibregma spp. Networks were constructed using 16S sequences, numbers indicate the number of nucleotide differences between haplotypes. Circle size reflects the number of individuals with that haplotype, black dots indicate missing haplotypes and n, the number of individuals for each potential species. Colored circles/segments relate to collection site where SR, Shag Rocks; SG, South Georgia; ST, Southern Thule; PB, Powel Basin; EI, Elephant Island; LI, Livingston Island; AS, Amundsen Sea; WS, Weddell Sea, RS, Ross Sea, grayscale for depth binned networks.

Observations of Widespread Species within West Antarctic

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.

FIGURE 3

Figure 3. Georeferenced (color) and depth binned (grayscale) haplotype networks for (A) Laonice weddellia, (B) Aricidea simplex, (C) Aricidea cf belgicae spp., (D) Harmothoe fuligineum, and (E) Euphrosinella cf cirratoformis spp. Networks were constructed using 16S sequences, numbers indicate the number of nucleotide differences between haplotypes. Circle size reflects the number of individuals with that haplotype, black dots indicate missing haplotypes and n, the number of individuals for each potential species. Colored circles/segments relate to collection site where SR, Shag Rocks; SG, South Georgia; ST, Southern Thule; PB, Powel Basin; EI, Elephant Island; LI, Livingston Island; AS, Amundsen Sea; WS, Weddell Sea, RS, Ross Sea; grayscale for depth binned networks.

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.

FIGURE 4

Figure 4. Georeferenced (color) and depth binned (grayscale) haplotype networks for (A) Aglaophamus spp. and (B) Hesionidae spp. Networks were constructed using 16S (Aglaophamus spp.) and COI (Hesionidae spp.) sequences, numbers indicate the number of nucleotide differences between haplotypes. Circle size reflects the number of individuals with that haplotype, black dots indicate missing haplotypes and n, the number of individuals for each potential species. Colored circles/segments relate to collection site where SR, Shag Rocks; SG, South Georgia; ST, Southern Thule; PB, Powel Basin; EI, Elephant Island; LI, Livingston Island; AS, Amundsen Sea; WS, Weddell Sea; RS, Ross Sea, grayscale for depth binned networks.

Observations of Restricted Species

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.

Depth Distributions

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.

Particle Tracking Analysis

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.

FIGURE 5

Figure 5. Binned particle trajectories released from the Ross Sea Shelf site (RS Sh; indicated by the white star) during the period 2000-2009. The mean 2°C isotherm at 100 m for this period is also plotted as a proxy indicator of the Polar Front (blue line). Values indicate the total number of particles passing through a hexbin over the sample period. For site abbreviations see Table 3.

FIGURE 6

Figure 6. Binned particle trajectories released from the Outer Amundsen Sea site (AS BIO3; indicated by the white star) during the period 2000-2009. The mean 2°C isotherm at 100 m for this period is also plotted as a proxy indicator of the Polar Front (blue line). Values indicate the total number of particles passing through a hexbin over the sample period. For site abbreviations see Table 3.

FIGURE 7

Figure 7. Binned particle trajectories released from the Livingston Island site (LI; indicated by the white star) during the period 2000-2009. The mean 2°C isotherm at 100 m for this period is also plotted as a proxy indicator of the Polar Front (blue line). Values indicate the total number of particles passing through a hexbin over the sample period. For site abbreviations see Table 3.

FIGURE 8

Figure 8. Binned particle trajectories released from the Weddell Sea Continental Slope site (WS CS; indicated by the white star) during the period 2000-2009. The mean 2°C isotherm at 100 m for this period is also plotted as a proxy indicator of the Polar Front (blue line). Values indicate the total number of particles passing through a hexbin over the sample period. For site abbreviations see Table 3.

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.

FIGURE 9

Figure 9. Site connectivity matrix. The values in the matrix indicate the percentage of particles that pass between two of the sample sites (origin: y-axis, destination: x-axis). The color shading indicates the mean transit time of the particles from point of origin to the site in question. The gray shading indicates that there is no connectivity (within 360 days) between two sites. For site abbreviations please see Table 3.

FIGURE 10

Figure 10. Site connectivity matrix. The values in the matrix indicate the percentage of particles that pass between two of the sample sites (origin: y-axis, destination: x-axis). The color shading indicates the mean depth at which the particles arrive at the site in question. The gray shading indicates that there is no connectivity (within 360 days) between two sites. For site abbreviations please see Table 3.

Discussion

In this study, nine of the seventeen genetically-identified species revealed different biogeographic patterns to their distributions inferred from morphology based identification records of the initial morphospecies and the RAMS database (Table 1). These nine species, have been described as cryptic species in Brasier et al. (2016), and their haplotype networks confirmed this displaying characteristics for cryptic speciation as described in Allcock and Strugnell (2012). In some cases, e.g., Aricidea cf belgicae (MB1), (MB2), and (MB3), more restricted distributions were observed but most remained widespread and potentially cicrum-Antarctic. Additionally, the networks showed the majority of cryptic species exhibited sympatric distributions, occupying the same geographic locations.

Widespread Antarctic Species

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).

Eurybathic Species

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.

Restricted Species

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.

The Use of Particle Tracking Analysis to Understand Antarctic Connectivity

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).

Limitations of Particle Tracking Models to Estimate Larval Dispersal of Antarctic Polychaetes

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.

Conclusions and Wider Implications

This study showed that the previously accepted biogeographic patterns of a third of the nine morphospecies examined should be questioned or re-described. Widespread distributions within the West Antarctic were recorded in 12 of the 17 species. These included 9 cryptic species existing sympatrically. The presence of widespread morphological and cryptic species is likely to be explained by their larval dispersal between populations as demonstrated by particle tracking models. The lack of consistency between the biogeography of cryptic species, some being widespread and potentially circum-Antarctic, whilst others are restricted, recorded within this study demonstrates the complexity of Southern Ocean biogeography (Brandão et al., 2010; Strugnell and Allcock, 2013; Chown et al., 2015). Fine scale differences in species distributions may be a result of variable life histories, habitat preferences, biological responses and ecological interactions within and between species through past and present physical conditions rather than a lack of transport connectivity from oceanographic currents. To fully appreciate why some species may be more dominant at certain depths, e.g., Hesionidae spp., or some widespread, whilst others are restricted, e.g., Scalibregma spp., would require investigations into their ecological traits to understand their functional differences. This is important because species with restricted distributions or limited dispersal capacities are often considered more vulnerable to extinction or less likely to recover from physical disturbances, as disturbed or removed populations may not be resupplied if they do not inhabit neighboring sites (Chown et al., 2015).

Our results presented here have valuable implications, improving our understanding of the drivers of biogeography and their implications for marine management under changing environmental conditions. The use of multiple data sets such as diversity, biogeographic and genetic data together with ocean physics model data, are valuable tools for the designation of effective marine management practices such as Marine Protected Areas (MPAs) and fishing restrictions (Robinson et al., 2017). The effectiveness of an MPA is, in part, reliant on the ability of species within the MPA to source external populations, thus mapping the likely dispersal pathways and distance of known species provides biological evidence for designation (e.g., Le Quesne and Codling, 2009; Planes et al., 2009). With increased benthic sampling and DNA barcoding these data can be used to assess the level of genetic connectivity between different polychaete populations, assess the abundance of rare and restricted species and provide further insight into the processes that determine species distributions.

Author Contributions

This paper includes the main finding of the second chapter of MB's PhD thesis. MB conducted all laboratory DNA data collection, led the analysis and wrote this manuscript. JH ran the particle tracking analyses, produced model figures and conducted statistical analyses on the particle data. HW significantly contributed to the acquisition, analysis and interpretation of the genetic data. RJ contributed to the conception of this project and manuscript editing. KL donated all of the specimens used in this project and assisted in the interpretation and discussion of the DNA data collected. HR contributed to hypothesis framing, model design making a substantial contribution to its application and contributed to the interpretation of results. AG led the polychaete sampling at sea of the BIOPEARLII samples, and project-managed the sorting and identification of all BIOPEARL polychaetes. AG also contributed to the conception of this project and manuscript editing. All authors contributed to the revised article and are accountable for all aspects of the work should the accuracy and integrity be questioned.

Funding

All DNA barcoding was supported by MB's PhD research training grant received from the University of Liverpool and the National Oceanography Centre. The expeditions JR144, JR179 (BIOPEARL I and II) and JR275 were part of the British Antarctic Survey core programmes “Global Science in the Antarctic Context” and “Polar Science for Planet Earth” funded by The Natural Environment Research Council.

Conflict of Interest Statement

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.

Acknowledgments

We would like to thank the crew and scientists who participated in the three research cruises that collected the polychaete specimens used in this study. We are particularly grateful to the BIOPEARL sample sorters that worked at sea, in the NHM Deep Sea Lab, and the BAS Lab: Rebekah Baker, David Barnes, Stefanie Kaiser, Ondine Cornubert, Adam Reed, Michael Mende, Wencke Krings, Moritz Stäbler and Chester Sands. Adrian Glover acknowledges the Natural Environment Research Council Collaborative Gearing Scheme for funding to participate in JR179 research cruise. Model data analyzed in this study was generated using the ARCHER UK National Supercomputing Service (http://www.archer.ac.uk). Particle tracking analysis was performed on the NERC funded JASMIN super-data-cluster facility (http://www.jasmin.ac.uk). Additional thanks are due to Lenka Neal for her help with specimen identification prior to DNA analyses. We are also grateful for the constructive comments provided by both the reviewers and editorial team.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmars.2017.00356/full#supplementary-material

Footnotes

1. ^PolyTraits Team. Data from: Polytraits: A database on biological traits of polychaetes. Hellenic Centre for Marine Research: Lifewatch, Greece (2017) http://polytraits.lifewatchgreece.eu

References

Alizon, S., Kucera, M., and Jansen, V. A. (2008). Competition between cryptic species explains variations in rates of lineage evolution. Proc. Natl. Acad. Sci. U.S.A. 105, 12382–12386. doi: 10.1073/pnas.0805039105

PubMed Abstract | CrossRef Full Text | Google Scholar

Allcock, A. L., and Strugnell, J. M. (2012). Southern Ocean diversity: new paradigms from molecular ecology. Trends Ecol. Evol. (Amst). 27, 520–528. doi: 10.1016/j.tree.2012.05.009

PubMed Abstract | CrossRef Full Text | Google Scholar

Allcock, A. L., Barratt, I., Eléaume, M., Linse, K., Norman, M. D., Smith, P. J., et al. (2011). Cryptic speciation and the circumpolarity debate: a case study on endemic Southern Ocean octopuses using the COI barcode of life. Deep-Sea Res. Part. II. Top. Stud. Oceanogr. 58, 242–249. doi: 10.1016/j.dsr2.2010.05.016