<?xml version="1.0" encoding="UTF-8" standalone="no"?>
<!DOCTYPE article PUBLIC "-//NLM//DTD Journal Publishing DTD v2.3 20070202//EN" "journalpublishing.dtd">
<article xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink" article-type="research-article">
<front>
<journal-meta>
<journal-id journal-id-type="publisher-id">Front. Ecol. Evol.</journal-id>
<journal-title>Frontiers in Ecology and Evolution</journal-title>
<abbrev-journal-title abbrev-type="pubmed">Front. Ecol. Evol.</abbrev-journal-title>
<issn pub-type="epub">2296-701X</issn>
<publisher>
<publisher-name>Frontiers Media S.A.</publisher-name>
</publisher>
</journal-meta>
<article-meta>
<article-id pub-id-type="doi">10.3389/fevo.2015.00063</article-id>
<article-categories>
<subj-group subj-group-type="heading">
<subject>Genetics</subject>
<subj-group>
<subject>Original Researche</subject>
</subj-group>
</subj-group>
</article-categories>
<title-group>
<article-title>Against the flow: evidence of multiple recent invasions of warmer continental shelf waters by a Southern Ocean brittle star</article-title>
</title-group>
<contrib-group>
<contrib contrib-type="author" corresp="yes">
<name><surname>Sands</surname> <given-names>Chester J.</given-names></name>
<xref ref-type="aff" rid="aff1"><sup>1</sup></xref>
<xref ref-type="author-notes" rid="fn001"><sup>&#x0002A;</sup></xref>
<uri xlink:href="http://loop.frontiersin.org/people/227126/overview"/>
</contrib>
<contrib contrib-type="author">
<name><surname>O&#x00027;Hara</surname> <given-names>Timothy D.</given-names></name>
<xref ref-type="aff" rid="aff2"><sup>2</sup></xref>
</contrib>
<contrib contrib-type="author">
<name><surname>Barnes</surname> <given-names>David K. A.</given-names></name>
<xref ref-type="aff" rid="aff1"><sup>1</sup></xref>
</contrib>
<contrib contrib-type="author">
<name><surname>Mart&#x000ED;n-Ledo</surname> <given-names>Rafael</given-names></name>
<xref ref-type="aff" rid="aff3"><sup>3</sup></xref>
<uri xlink:href="http://loop.frontiersin.org/people/234270/overview"/>
</contrib>
</contrib-group>
<aff id="aff1"><sup>1</sup><institution>Biodiversity, Evolution and Adaptation, Natural Environment Research Council, British Antarctic Survey</institution> <country>Cambridge, UK</country></aff>
<aff id="aff2"><sup>2</sup><institution>Museum Victoria</institution> <country>Melbourne, VIC, Australia</country></aff>
<aff id="aff3"><sup>3</sup><institution>&#x000C1;rea de Zoolog&#x000ED;a, Facultad de Ciencias, Universidad de Extremadura</institution> <country>Badajoz, Spain</country></aff>
<author-notes>
<fn fn-type="edited-by"><p>Edited by: Christian C. Figueroa, Universidad de Talca, Chile</p></fn>
<fn fn-type="edited-by"><p>Reviewed by: Florian Leese, Ruhr University Bochum, Germany; Pilar Haye Molina, Universidad Cat&#x000F3;lica del Norte, Chile</p></fn>
<fn fn-type="corresp" id="fn001"><p>&#x0002A;Correspondence: Chester J. Sands, British Antarctic Survey, High Cross, Madingley Road, Cambridge, Cambridgeshire, CB3 0ET, UK <email>cjsan&#x00040;bas.ac.uk</email></p></fn>
<fn fn-type="other" id="fn002"><p>This article was submitted to Evolutionary and Population Genetics, a section of the journal Frontiers in Ecology and Evolution</p></fn>
</author-notes>
<pub-date pub-type="epub">
<day>22</day>
<month>06</month>
<year>2015</year>
</pub-date>
<pub-date pub-type="collection">
<year>2015</year>
</pub-date>
<volume>3</volume>
<elocation-id>63</elocation-id>
<history>
<date date-type="received">
<day>01</day>
<month>04</month>
<year>2015</year>
</date>
<date date-type="accepted">
<day>03</day>
<month>06</month>
<year>2015</year>
</date>
</history>
<permissions>
<copyright-statement>Copyright &#x000A9; 2015 Sands, O&#x00027;Hara, Barnes and Mart&#x000ED;n-Ledo.</copyright-statement>
<copyright-year>2015</copyright-year>
<copyright-holder>Sands, O&#x00027;Hara, Barnes and Mart&#x000ED;n-Ledo</copyright-holder>
<license xlink:href="http://creativecommons.org/licenses/by/4.0/"><p>This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) or licensor are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.</p>
</license>
</permissions>
<abstract><p>The Southern Ocean is anomalously rich in benthos. This biodiversity is native, mostly endemic and perceived to be uniquely threatened from climate- and anthropogenically- mediated invasions. Major international scientific effort throughout the last decade has revealed more connectivity than expected between fauna north and south of the worlds strongest marine barrier&#x02014;the Polar Front (the strongest jet of the Antarctic Circumpolar Current). To date though, no research has demonstrated any radiations of marine taxa out from the Southern Ocean, except at abyssal depths (where conditions differ much less). Our phylogeographic investigation including mitochondrial sequences of 135 specimens of one of the most ubiquitous and abundant clades at high southern latitudes, the ophiuroids (brittlestars), shows that one of them, <italic>Ophiura lymani</italic>, has gone against the flow. Remarkably our phylogenetic and coalescent analyses suggest that <italic>O. lymani</italic> has successfully invaded the South American shelf from Antarctica at least three times, in recent (Pleistocene) radiation. Many previous studies have demonstrated links within clades across the Polar Front this is the first in which northwards directional movement of a shelf-restricted species is the only convincing explanation. Rapid, recent, regional warming is likely to facilitate multiple range shift invasions into the Southern Ocean, whereas movement of cold adapted fauna (considered highly stenothermal) out of the Antarctic to warmer shelves has, until now, seemed highly unlikely.</p></abstract>
<kwd-group>
<kwd>adaptation</kwd>
<kwd>climate change</kwd>
<kwd>cryptic species</kwd>
<kwd>Ophiuroidea</kwd>
<kwd>range expansion</kwd>
<kwd>radiation</kwd>
<kwd>phylogeography</kwd>
<kwd>Pleistocene</kwd>
</kwd-group>
<contract-sponsor id="cn001">Census of Antarctic Marine Life</contract-sponsor>
<contract-sponsor id="cn002">Antarctic Science Bursary</contract-sponsor>
<contract-sponsor id="cn003">BBSRC</contract-sponsor>
<contract-sponsor id="cn004">NERC</contract-sponsor>
<counts>
<fig-count count="7"/>
<table-count count="4"/>
<equation-count count="0"/>
<ref-count count="81"/>
<page-count count="13"/>
<word-count count="8777"/>
</counts>
</article-meta>
</front>
<body>
<sec sec-type="introduction" id="s1">
<title>Introduction</title>
<p>Marine animals inhabiting continental shelves can disperse widely over ecological time scales due to extensive shelf connectivity, current regimes, and reasonably consistent ocean chemistry (Paulay and Meyer, <xref ref-type="bibr" rid="B57">2006</xref>). Nevertheless, some local endemism occurs in most continental shelf regions (Sanford and Kelly, <xref ref-type="bibr" rid="B68">2011</xref>). Temperature is a strong predictor of biogeographic regions (Belanger et al., <xref ref-type="bibr" rid="B8">2012</xref>) and studies demonstrate that given a warming climate many marine shelf species track pole-wards and into deeper, cooler habitats (Norberg et al., <xref ref-type="bibr" rid="B54">2012</xref>; Pinsky et al., <xref ref-type="bibr" rid="B59">2013</xref>; Hiddink et al., <xref ref-type="bibr" rid="B36">2015</xref>). Shelf areas within the Southern Ocean are unusually isolated compared with those north of the Southern Ocean. They are surrounded by deep seas (&#x0003E;2000 m depth), the strongest oceanic current on the planet, the Antarctic Circumpolar Current (ACC) (Whitworth, <xref ref-type="bibr" rid="B80">1980</xref>; Orsi et al., <xref ref-type="bibr" rid="B55">1995</xref>), and encircled by the Polar Front, a narrow jet within the ACC, across which there are steep gradients in temperature and salinity. These isolating conditions, established in the Oligocene (Barker and Thomas, <xref ref-type="bibr" rid="B3">2004</xref>), have left a signature of high endemicity of the Southern Ocean shelf fauna (Dell, <xref ref-type="bibr" rid="B16">1972</xref>; Clarke and Crame, <xref ref-type="bibr" rid="B11">1989</xref>; Griffiths et al., <xref ref-type="bibr" rid="B32">2009</xref>).</p>
<p>Over the &#x0007E;24 million years since the opening of the Drake Passage, the Miocene Antarctic benthos has had to respond to the cooling conditions by either adapting <italic>in situ</italic> or range shifts into appropriate environmental conditions. This for instance could have resulted in species ranges being restricted to the northern continental shelves of South Georgia and South America, while Antarctic shelf populations went extinct (Clarke and Crame, <xref ref-type="bibr" rid="B11">1989</xref>). The onset of glaciation caused extinction in most, but not all, of the terrestrial flora and fauna (Convey et al., <xref ref-type="bibr" rid="B14">2008</xref>), but in the more constant marine environment, much of the fauna has adapted successfully and has become comparably rich (Clarke and Johnston, <xref ref-type="bibr" rid="B12">2003</xref>). There has long been debate regarding the origins of the Antarctic fauna and how much movement there has been of fauna into and out of the Antarctic realm over time since its isolation (Hedgpeth, <xref ref-type="bibr" rid="B34">1969</xref>; Dell, <xref ref-type="bibr" rid="B16">1972</xref>; Clarke and Crame, <xref ref-type="bibr" rid="B11">1989</xref>; Clarke et al., <xref ref-type="bibr" rid="B10">2005</xref>; Barnes et al., <xref ref-type="bibr" rid="B5">2006</xref>; Griffiths et al., <xref ref-type="bibr" rid="B32">2009</xref>). Although there is clearly biogeographic evidence of a South American influence in the Southern Ocean fauna (Griffiths et al., <xref ref-type="bibr" rid="B32">2009</xref>), phylogenetic studies suggest that octopus (Strugnell et al., <xref ref-type="bibr" rid="B70">2008</xref>), and pleurobranch gastropods (G&#x000F6;bbeler and Klussmann-Kolb, <xref ref-type="bibr" rid="B29">2010</xref>) have radiated out of the Antarctic, possibly assisted by the initiation of the formation of Antarctic bottom water, around the time of the Oligocene cooling. There is also molecular evidence suggesting that some Antarctic species and their northern lineages diverged much more recently in the Pliocene or Pliestocene (Page and Linse, <xref ref-type="bibr" rid="B56">2002</xref>; Hunter and Halanych, <xref ref-type="bibr" rid="B41">2008</xref>; Gonz&#x000E1;lez-Wevar et al., <xref ref-type="bibr" rid="B31">2010</xref>; Krabbe et al., <xref ref-type="bibr" rid="B44">2010</xref>; Poulin et al., <xref ref-type="bibr" rid="B61">2014</xref>), while several species are currently thought to be distributed on either side of the Polar Front indicating recurrent gene flow (see Sands et al., <xref ref-type="bibr" rid="B67">2013</xref> Supplementary Material for ophiuroid examples).</p>
<p>The recently observed increase in Southern Ocean surface temperatures (Meredith and King, <xref ref-type="bibr" rid="B50">2005</xref>; Whitehouse et al., <xref ref-type="bibr" rid="B79">2008</xref>; Turner et al., <xref ref-type="bibr" rid="B77">2014</xref>) is suggested to have serious implications for Southern Ocean biota. The fauna is widely considered to be cold adapted with no option of migration except to the deep sea (Peck, <xref ref-type="bibr" rid="B58">2005</xref>), which in itself requires specific adaptations for survival (Brown and Thatje, <xref ref-type="bibr" rid="B9">2015</xref>). Warming of the Southern Ocean is also likely to facilitate invasion of species from more northern regions. This is most likely to be observed first on the more northern shelf regions of the Southern Ocean, especially the large micro-continental shelf regions of South Georgia and Shag Rocks (Barnes et al., <xref ref-type="bibr" rid="B4">2009</xref>; Hogg et al., <xref ref-type="bibr" rid="B39">2011</xref>; Morley et al., <xref ref-type="bibr" rid="B52">2014</xref>) and the even larger but less studied waters of the Kerguelen Plateau. Despite being located south of the Polar Front, South Georgia is eastwards of the tip of South America and directly affected by the ACC flowing across from the South American Shelf. The benthic fauna of the South Georgia shelf has been described as &#x0201C;Antarctic&#x0201D; (Hogg et al., <xref ref-type="bibr" rid="B39">2011</xref>), although biogeographic work indicates that there is substantial South American influence (Griffiths et al., <xref ref-type="bibr" rid="B32">2009</xref>). There are several shelf-restricted brittle star (Ophiuroidea) species that are located on both the Patagonian and South Georgian shelves (Sands et al., <xref ref-type="bibr" rid="B67">2013</xref>). Given the strong ACC with the potential to promote gene flow we expect a pattern consistent with the South American influence hypothesis that predicts (1) decreasing genetic diversity from the Patagonian Shelf to Antarctic Shelf; (2) older lineages i.e., more basal clades in a phylogeny, will consist of Patagonian Shelf specimens and (3) tests for gene flow reveal asymmetrical migration in the direction of the ACC from the Patagonian shelf population to South Georgia. Using a fragment of the mitochondrial cytochrome c oxidase subunit 1 gene (CO1) amplified from collections of the brittle star <italic>Ophiura lymani</italic> and its sister species <italic>Ophiura carinifera</italic>, acquired throughout their ranges around the Southern Ocean, we tested these expectations under the hypothesis of South American origin against the alternative hypothesis of Antarctic origin.</p>
</sec>
<sec sec-type="materials and methods" id="s2">
<title>Materials and methods</title>
<p>Our approach was to first test for genetic structure within our samples using phylogenetic methods. We then conducted tests using coalescent techniques to best identify the processes that underlying the phylogenetic patterns identified.</p>
<p>Specimens of <italic>O. lymani</italic> and <italic>O. carinifera</italic> were collected on several international expeditions (JR 144, JR179, JR230, PS77, ANDEEP, ICEFISH, CEAMARC, and MacRidge 2) using Agassiz Trawls. Sampling locations are illustrated in Figure <xref ref-type="fig" rid="F1">1</xref>. Samples were sorted immediately on collection and stored in 96% ethanol. Individuals were identified morphologically to species level.</p>
<fig id="F1" position="float">
<label>Figure 1</label>
<caption><p><bold>Map of Antarctica and the Southern Ocean showing sampling locations of</bold> <italic><bold>Ophiura lymani</bold></italic> <bold>(purple) and</bold> <italic><bold>Ophiura carinifera</bold></italic> <bold>(blue) and</bold> <italic><bold>Ophiura rouchi</bold></italic> <bold>(green)</bold>. The approximate positions of the fronts defining the ACC are indicated with the Polar Front indicated in a black dashed line.</p></caption>
<graphic xlink:href="fevo-03-00063-g0001.tif"/>
</fig>
<p>Tissue samples taken from arm tips were used for DNA extraction. DNA extraction and sequencing of partial cytochrome c oxidase subunit 1 (CO1) was carried out by the Canadian Center for DNA Barcoding using the primers LCOech1aF1 and HCO2198 (Folmer et al., <xref ref-type="bibr" rid="B25">1994</xref>), with sequences and metadata stored and made publically available on the Barcode of Life Data Systems (Ratnasingham and Hebert, <xref ref-type="bibr" rid="B64">2007</xref>). Forward and reverse trace files were downloaded from the database and assembled and edited in CodonCode Aligner V5.0.1. Sequences were translated to ensure an open reading frame, and aligned using the MUSCLE (Edgar, <xref ref-type="bibr" rid="B22">2004</xref>) alignment process. We included sequences from <italic>Ophiura rouchi</italic> as the most reasonable outgroup and to stabilize more basal nodes (Felsenstein, <xref ref-type="bibr" rid="B24">1981</xref>) we included two further outgroups, <italic>Ophiocten dubium</italic> and <italic>Ophioperla koehleri</italic>.</p>
<p>Descriptive statistics produced using DnaSP v5.10.01 (Librado and Rozas, <xref ref-type="bibr" rid="B47">2009</xref>) and Arlequin v3.5 (Excoffier et al., <xref ref-type="bibr" rid="B23">2005</xref>). Departures from drift mutation equilibrium were tested by comparing estimates of nucleotide diversity with segregating sites&#x02014;Tajima&#x00027;s D (Tajima, <xref ref-type="bibr" rid="B72">1989</xref>)&#x02014;and by assessing observed vs. expected haplotype diversity&#x02014;Fu&#x00027;s Fs (Fu, <xref ref-type="bibr" rid="B27">1996</xref>). Both these statistics are affected by population size changes as well as selection. Significance was assessed from 1000 coalescent simulations in DnaSP. We conducted mismatch analysis (Rogers and Harpending, <xref ref-type="bibr" rid="B65">1992</xref>) to test whether the frequency distribution of pairwise differences conformed with the random (ragged) distribution expected under equilibrium or whether the distribution appeared non-random as expected under past population change. <italic>P</italic>-values were generated by conducting 1000 bootstrap pseudoreplicates in Arlequin. Within and among clade distances were calculated in MEGA v5.2.2 (Tamura et al., <xref ref-type="bibr" rid="B73">2011</xref>). Base compositional heterogeneity was examined using a match paired test in SeqVis v1.5 (Ho et al., <xref ref-type="bibr" rid="B37">2006</xref>).</p>
<p>To visualize the change of effective size through time in each population, Bayesian Skyline (Drummond et al., <xref ref-type="bibr" rid="B19">2005</xref>) and Bayesian Skyride plots (Minin et al., <xref ref-type="bibr" rid="B51">2008</xref>) were produced in BEAST v1.8.1(Drummond and Rambaut, <xref ref-type="bibr" rid="B20">2007</xref>) using coalescent Skyline and Gaussian Markov random field Bayesian Sky Ride tree priors respectively, and a mutation rate of 2.48 &#x000D7; 10 <sup>&#x02212;8</sup> (see below). Ten million generations of MCMC produced 10,000 trees. TRACER v1.6 (Rambaut et al., <xref ref-type="bibr" rid="B63">2014</xref>) was used to check diagnostics of convergence and effective sample sizes of parameters, as well as to analyse the trees and produce the plots.</p>
<p>To determine the relationships among <italic>O. lymani</italic> populations and the closely related High Antarctic species <italic>O. carinifera</italic>, we produced a Maximum Likelihood estimate of phylogeny using RAxML v8.0.26 (Stamatakis, <xref ref-type="bibr" rid="B69">2014</xref>) and a Bayesian estimate using the program MrBayes v3.2 (Huelsenbeck and Ronquist, <xref ref-type="bibr" rid="B40">2001</xref>; Ronquist and Huelsenbeck, <xref ref-type="bibr" rid="B66">2003</xref>). In order to account for violations of assumptions of stationarity and heterogeneity we used JModeltest v2 (Darriba et al., <xref ref-type="bibr" rid="B15">2012</xref>) and used AIC to determine that GTR&#x0002B;I&#x0002B;G was the most suitable substitution model. As proportion of invariant sites (I) is closely correlated with the gamma distribution (&#x0002B;G) we used the GTR&#x0002B;G option in RAxML, with 10 replicate runs on both a codon partitioned dataset and a non-partitioned dataset. We ran 1000 &#x0201C;thorough&#x0201D; (as opposed to the &#x0201C;rapid&#x0201D; option) bootstrap pseudoreplicates to estimate node support. A codon partitioned dataset was used in our Bayesian analysis to better account for variation in substitution rates across sites and lineages. Each partition was set to estimate all parameters (nst &#x0003D; 6) effectively allowing a GTR model with a gamma distribution. Phylogenies were estimated using four replicates of four heated chains over 10<sup>8</sup> generations with a 25% burnin that ensured trees were sampled only once split frequencies fell below 0.01. Convergence and effective sample size estimates of parameters were checked using TRACER v1.6. An alternative visualization of the relationships among haplotypes of <italic>O. lymani</italic> sampled from the Patagonian shelf and the South Georgia/Shag Rocks shelf areas was generated in the form of a statistical parsimony network (Templeton et al., <xref ref-type="bibr" rid="B75">1992</xref>) using the package TCS v1.21 (Clement et al., <xref ref-type="bibr" rid="B13">2000</xref>). From this reduced dataset we constructed a Bayesian phylogeny as described above to be used as the gene tree estimate in coalescent simulations.</p>
<p>Species delimitation was explored using popular analytic techniques. Automatic barcode gap discovery (Puillandre et al., <xref ref-type="bibr" rid="B62">2012</xref>) is a distance-based technique that explores a predicted &#x0201C;gap&#x0201D; in the frequencies of pairwise distances indicating the discrete change in intra- and interspecific evolutionary processes. Parameters <italic>P</italic><sub>min</sub> (lowest expected distance to include a barcode gap), <italic>P</italic><sub>max</sub> (maximum expected distance to include a barcode gap), X (relative width of the barcode gap), and distance correction (JC69, K2P, SimpleP) were adjusted to best encompass the range of variation in the data (<italic>P</italic><sub>min</sub> &#x0003D; 0.001, <italic>P</italic><sub>max</sub>, <italic>X</italic> &#x0003D; 0.1, K2P distance). The resulting number of species groups was estimated after the stabilization of recursive partitioning of the data. The general mixed Yule coalescent (GMYC) method is a phylogenetic technique (Pons et al., <xref ref-type="bibr" rid="B60">2006</xref>) that was implemented in the R package SPLITS (Fujisawa and Barraclough, <xref ref-type="bibr" rid="B28">2013</xref>) using an ultrametric tree produced in the software BEAST v1.8.1 (Drummond and Rambaut, <xref ref-type="bibr" rid="B20">2007</xref>) as described below.</p>
<p>Where results of non-monophyly between regions were observed we used a test of the cost of deep coalescence (cDC) events required to explain the data compared with that of distribution of coalescent trees simulated from the original gene tree. As the coalescent process is stochastic, any pattern inferred by phylogeny may be chance&#x02014;rather than process&#x02014;driven. In the software package Mesquite v3.01 (Maddison and Maddison, <xref ref-type="bibr" rid="B49">2014</xref>) we constructed a species tree specifying a topology representing two populations: Patagonian shelf and South Georgian shelf. The Bayesian gene tree estimated from the reduced dataset was forced into the species tree and the cDC recorded. The same gene tree was used as a starting point to simulate 10,000 coalescent trees, each one forced to fit the species tree and the number of deep coalescent events noted. The cDC from the gene tree was compared to a distribution of the cDC from the simulations.</p>
<p>The likelihood of the exchange of migrants between the Patagonian shelf and South Georgia shelf was tested in Migrate v3.6.4 (Beerli, <xref ref-type="bibr" rid="B6">2009</xref>) by examining the posterior distribution of migration estimates from an unrestricted model allowing full symmetrical migration, and then by ranking competing migration models using thermodynamic integration with a Bayes factor test (Beerli and Palczewski, <xref ref-type="bibr" rid="B7">2010</xref>). The models examined were symmetrical migration (PS&#x0003C;&#x02212;&#x0003E;SG), asymmetrical migration (PS-&#x0003E;SG), (PS&#x0003C;-SG), panmixia (PS &#x0003D; SG) and no migration (PS-x-SG). Prior to running these models, Migrate-n settings were optimized using a full model (symmetrical migration) to ensure priors were appropriate and runs would converge. Final run conditions were 10 replicate long chains each with four heated chains with temperatures 1, 1.5, 3, 10<sup>6</sup>. Each chain was run for 10<sup>6</sup> generations, sampling every 100 generations.</p>
<p>The timing of divergence events was estimated using a strict molecular clock in BEAST v1.8.1 (Drummond and Rambaut, <xref ref-type="bibr" rid="B20">2007</xref>). We allowed a GTR &#x0002B; G model on partitioned codon sites and a Yule process speciation prior. A lognormal prior was used for &#x0201C;clock rate&#x0201D; such that the median reflected the mutation rate per lineage per year (2.48 &#x000D7; 10 <sup>&#x02212;8</sup>) as determined by Naughton et al. (<xref ref-type="bibr" rid="B53">2014</xref>) and 95% distribution between 1.2 &#x000D7; 10 <sup>&#x02212;8</sup> and 8.7 &#x000D7; 10 <sup>&#x02212;8</sup> thus covering a broad range of uncertainty. Sixteen test runs of 10<sup>8</sup> generations were used to test and adjust priors and their effects on convergence, mixing and parameter estimates.</p>
<p>All data and meta-data used in this study has been made publically available on BOLD Systems (<ext-link ext-link-type="uri" xlink:href="http://www.boldsystems.org">www.boldsystems.org</ext-link>) under the project code DS-MS001 (DOI: dx.doi.org/10.5883/DS-MS001) and on GenBank with accession numbers KR861569-KR86164.</p>
</sec>
<sec sec-type="results" id="s3">
<title>Results</title>
<sec>
<title>The samples</title>
<p>A total of 69 <italic>O. lymani</italic> specimens were collected; 33 from the Patagonian shelf, 31 from the South Georgia/Shag Rocks shelf, three from South Sandwich Islands and two close to the remote Bouvetoya. Morphological examination of all specimens supported a single morpho-species.</p>
<p><italic>O. carinifera</italic>, the sister taxon to <italic>O. lymani</italic>, is clearly identifiable by strongly elevated and granulated disc plates and radial shields (flattened without granules in <italic>O. lymani</italic>), and narrow dorsal arm plates with a serrated keel (wide without serrated rest in <italic>O. lymani</italic>). We included 56 specimens sampled from around the Antarctic continental shelf&#x02014;from the Amundsen Sea, Bellingshausen Sea, Larsen shelf, South Orkney shelf and the south-east region of theWeddell Sea, east Antarctica near the Mertz Glacier and a specimen from the Macquarie Ridge. The Macquarie Ridge specimen (OPHAN310) was morphologically distinct from all others having a more flattened disc, similar to <italic>O. lymani</italic> and genetically divergent from all other <italic>O. carinifera</italic>, appearing to be the sole representative of a sister clade (Figure <xref ref-type="fig" rid="F2">2</xref>).</p>
<fig id="F2" position="float">
<label>Figure 2</label>
<caption><p><bold>Bayesian consensus phylogeny of partitioned CO1 sequences derived from</bold> <italic><bold>Ophiura lymani</bold></italic> <bold>from South Georgia shelf (red) Patagonian shelf (green), South Sandwich Islands (purple) and Bouvetoya (mustard)</bold>. <italic>Ophiura carinifera</italic> samples are blue, <italic>Ophiura rouchi</italic> samples are black. Posterior probabilities are given under the line by the associated adjacent node and bootstrap values are given above the line. No value was given if bootstrap was less than 50%. The identifiers on the tips are BOLD system accession numbers for the individual sequences.</p></caption>
<graphic xlink:href="fevo-03-00063-g0002.tif"/>
</fig>
</sec>
<sec>
<title>Molecular diversity</title>
<p>Our alignment of 132 sequences consisted of 658 bases, of which 138 were variable and 97 were parsimony informative. As expected with mitochondrial sequence there was a bias toward Adenosine and Thymine residues (GC &#x0003D; 39.2%).</p>
<p>Haplotype diversity (H<sub>D</sub>) was exceptionally high in all populations of <italic>O. lymani</italic> and <italic>O</italic>. <italic>carinifera</italic> (Table <xref ref-type="table" rid="T1">1</xref>). Tajima&#x00027;s D was negative, although not significantly different from zero in <italic>O. lymani</italic>. Significant negative values as seen in <italic>O. carinifera</italic> may indicate either purifying selection or population size change. Fu&#x00027;s Fs was also negative as to be expected given the high haplotype diversity. Fu&#x00027;s Fs was highly significant in the <italic>O. lymani</italic> populations from the Patagonian shelf and South Georgian shelf, both of which have very high haplotype diversity. Significant <italic>R</italic><sup>2</sup>-values also suggest population expansion, but only in <italic>O. carinifera</italic> and the South Georgia population of <italic>O. lymani</italic>. The frequency distribution of pairwise differences under a constant population is expected to be largely random, whereas the expectation under population growth is for a smooth curve. Although non-significant, the two <italic>O. lymani</italic> populations, particularly the South Georgia population, appeared to more closely fit the expectations of population expansion (Figure <xref ref-type="fig" rid="F3">3</xref>). Coalescent-based demographic reconstruction using Bayesian Skyline and Skyride plots indicated a population increase over the past few glacial cycles in all three populations (Figure <xref ref-type="fig" rid="F4">4</xref>). The Skyline plots show a sharp increase in population size starting 100&#x0002B; ka, whereas smoothing algorithms in the Skyride plots indicate a steady increase over time.</p>
<table-wrap position="float" id="T1">
<label>Table 1</label>
<caption><p><bold>Summary statistics of molecular diversity</bold>.</p></caption>
<table frame="hsides" rules="groups">
<thead>
<tr>
<th align="left"><bold>Population</bold></th>
<th align="center"><bold><italic>n</italic></bold></th>
<th align="center"><bold>Prob</bold></th>
<th align="center"><bold>No. haplotypes</bold></th>
<th align="center"><bold><italic>H</italic><sub>D</sub></bold></th>
<th align="center"><bold><italic>S</italic></bold></th>
<th align="center"><bold>&#x003C0;</bold></th>
<th align="center"><bold><italic>D<sub>T</sub></italic></bold></th>
<th align="center"><bold><italic>Fs</italic></bold></th>
<th align="center"><bold><italic>R</italic><sub>2</sub></bold></th>
<th align="center"><bold>Max. K</bold></th>
</tr>
</thead>
<tbody>
<tr>
<td align="left"><italic>O. lymani PS</italic></td>
<td align="center">33</td>
<td align="center">0.9412</td>
<td align="center">33</td>
<td align="center">1</td>
<td align="center">67</td>
<td align="center">0.018</td>
<td align="center">&#x02212;1.09</td>
<td align="center">&#x02212;22.17<xref ref-type="table-fn" rid="TN2"><sup>&#x0002A;&#x0002A;</sup></xref></td>
<td align="center">0.073169</td>
<td align="center">30</td>
</tr>
<tr>
<td align="left"><italic>O. lymani SG</italic></td>
<td align="center">31</td>
<td align="center">0.9375</td>
<td align="center">27</td>
<td align="center">0.989</td>
<td align="center">43</td>
<td align="center">0.010</td>
<td align="center">&#x02212;1.39</td>
<td align="center">&#x02212;15.62<xref ref-type="table-fn" rid="TN2"><sup>&#x0002A;&#x0002A;</sup></xref></td>
<td align="center">0.0637<xref ref-type="table-fn" rid="TN1"><sup>&#x0002A;</sup></xref></td>
<td align="center">18</td>
</tr>
<tr>
<td align="left"><italic>O. carinifera AS</italic></td>
<td align="center">55</td>
<td align="center">0.9643</td>
<td align="center">31</td>
<td align="center">0.926</td>
<td align="center">36</td>
<td align="center">0.006</td>
<td align="center">&#x02212;1.86<xref ref-type="table-fn" rid="TN1"><sup>&#x0002A;</sup></xref></td>
<td align="center">&#x02212;5.05</td>
<td align="center">0.0423<xref ref-type="table-fn" rid="TN2"><sup>&#x0002A;&#x0002A;</sup></xref></td>
<td align="center">12</td>
</tr>
<tr>
<td align="left"><italic>O. carinifera MR</italic></td>
<td align="center">1</td>
<td align="center">0</td>
<td align="center">1</td>
<td align="center">na</td>
<td align="center">na</td>
<td align="center">na</td>
<td align="center">na</td>
<td align="center">na</td>
<td align="center">na</td>
<td align="center">na</td>
</tr>
<tr>
<td align="left"><italic>O. lymani SSI</italic></td>
<td align="center">3</td>
<td align="center">0.5</td>
<td align="center">3</td>
<td align="center">1</td>
<td align="center">5</td>
<td align="center">0.005</td>
<td align="center">na</td>
<td align="center">na</td>
<td align="center">na</td>
<td align="center">4</td>
</tr>
<tr>
<td align="left"><italic>O. lymani BI</italic></td>
<td align="center">2</td>
<td align="center">0.325</td>
<td align="center">2</td>
<td align="center">1</td>
<td align="center">8</td>
<td align="center">0.012</td>
<td align="center">na</td>
<td align="center">na</td>
<td align="center">na</td>
<td align="center">8</td>
</tr>
</tbody>
</table>
<table-wrap-foot>
<fn id="TN1">
<label>&#x0002A;</label>
<p><italic>P &#x0003E; 0.05,</italic></p></fn>
<fn id="TN2">
<label>&#x0002A;&#x0002A;</label>
<p><italic>P &#x0003E; 0.01</italic>.</p></fn>
<p><italic>n, Number of individuals; Prob, Probablility of having captured the deepest coalescent event; No. haplotypes, Number of haplotypes; H<sub>D</sub>, Haplotype diversity; S, Number of segregating sites; &#x003C0;, Nucleotide diversity; D<sub>T</sub>, Tajima&#x00027;s D; Fs, Fu&#x00027;s Fs statistic; R<sub>2</sub>, RamosOnsins and Rozas&#x00027; R2 statistic; Max K, Maximum number of nucleotide differences between any two sequences within the population</italic>.</p>
</table-wrap-foot>
</table-wrap>
<fig id="F3" position="float">
<label>Figure 3</label>
<caption><p><bold>The expected (dotted line) and observed (solid line) distribution of pairwise frequencies between haplotypes under the expectation of constant size (A,C,E) and exponential growth (B,D,F) for</bold> <italic><bold>O. carinifera</bold></italic> <bold>(A,B)</bold>, <italic><bold>O. lymani</bold></italic> <bold>from the Patagonian Shelf (C,D) and</bold> <italic><bold>O. lymani</bold></italic> <bold>from the South Georgia shelf (E,F)</bold>.</p></caption>
<graphic xlink:href="fevo-03-00063-g0003.tif"/>
</fig>
<fig id="F4" position="float">
<label>Figure 4</label>
<caption><p><bold>Bayesian SkyRide (top) and SkyLine (bottom) plots indicating change in population size over time in</bold> <italic><bold>O. carinifera</bold></italic> <bold>(blue)</bold>, <italic><bold>O. lymani</bold></italic> <bold>from the Patagonian Shelf (green) and the South Georgia shelf (red)</bold>.</p></caption>
<graphic xlink:href="fevo-03-00063-g0004.tif"/>
</fig>
</sec>
<sec>
<title>Phylogenetic reconstruction</title>
<p>Phylogenetic reconstruction identified strongly supported clades of geographically defined populations of <italic>O. lymani</italic>, with the exception of five individuals sampled from the Patagonian shelf that grouped with the South Georgian shelf individuals, and a single Antarctic-wide clade including two lineages of <italic>O. carinifera</italic> (Figure <xref ref-type="fig" rid="F2">2</xref>). The low support of the node containing the South Sandwich Islands, South Georgia and Patagonian Shelf clades (pp. 0.62) persisted despite trialing various plausible evolutionary models and partitioning schemes. Low support for nodes should not always be interpreted as phylogenetic failure, but may rather reflect a genuine (hard) polytomy due to a radiation or series of speciation events between which there was insufficient time to accumulate mutational variation to permit accurate reconstruction (Hoelzer and Meinick, <xref ref-type="bibr" rid="B38">1994</xref>).</p>
<p>The statistical parsimony network highlighted five individuals sampled from the Patagonian shelf that either originated from the South Georgia clade, or are retained ancient lineages indicative of incomplete lineage sorting (Figure <xref ref-type="fig" rid="F5">5</xref>). In order to determine the likelihood of these lineages being retained, we determined the cost of deep coalescence (cDC) required to explain the The distribution of cDC of the simulated trees had a mean of 48 (41&#x02013;55, 95% CI). The cDC estimated from the Bayesian phylogeny was 4, well outside of this distribution, indicating that chance alone is very unlikely to explain the low number of deep coalescences required to produce this species tree from the gene tree.</p>
<fig id="F5" position="float">
<label>Figure 5</label>
<caption><p><bold>Statistical Parsimony network indicating sampling locations and genetic relationships among Patagonian (green) and South Georgian (red)</bold> <italic><bold>Ophiura lymani</bold></italic> <bold>haplotypes</bold>.</p></caption>
<graphic xlink:href="fevo-03-00063-g0005.tif"/>
</fig>
</sec>
<sec>
<title>Species delimitation</title>
<p>Genetic distances between populations were considerably greater than within-population distances, and generally greater than are expected for conspecific populations (Table <xref ref-type="table" rid="T2">2</xref>). Visual inspection suggests that clades representing <italic>O. lymani</italic> from Bouvet, the South Sandwich Islands and the Patagonian shelf/South Georgia are likely to be three distinct species, considering that the clearly morphologically distinct <italic>O. carinifera</italic> is grouped within the pan <italic>O. lymani</italic> clade. Depending on the relative gap width setting (X), ABGD identified these groupings, as well as OPHAN310 (<italic>O</italic>. <italic>carinifera</italic> from Macquarie Ridge) and DSOPH1855 (<italic>O. lymani</italic> from Patagonian shelf grouped into the South Georgia clade), implied to be two distinct species when X is set to 1. For larger values of X, ABGD failed to stabilize.</p>
<table-wrap position="float" id="T2">
<label>Table 2</label>
<caption><p><bold>Average pairwise distances between clades (below diagonal), associated standard error estimates (above diagonal) and within clades distances (along the diagonal)</bold>.</p></caption>
<table frame="hsides" rules="groups">
<thead>
<tr>
<th/>
<th align="left"><bold><italic>O. lymani BI</italic></bold></th>
<th align="center"><bold><italic>O. carinifera</italic></bold></th>
<th align="center"><bold><italic>O. lymani SSI</italic></bold></th>
<th align="center"><bold><italic>O. lymani PS</italic></bold></th>
<th align="center"><bold><italic>O. lymani SG</italic></bold></th>
</tr>
</thead>
<tbody>
<tr>
<td align="left"><italic>O. lymani BI</italic></td>
<td align="center"><bold>0.012</bold></td>
<td align="center">0.012</td>
<td align="center">0.010</td>
<td align="center">0.013</td>
<td align="center">0.012</td>
</tr>
<tr>
<td align="left"><italic>O. carinifera</italic></td>
<td align="center">0.051</td>
<td align="center"><bold>0.007</bold></td>
<td align="center">0.011</td>
<td align="center">0.011</td>
<td align="center">0.010</td>
</tr>
<tr>
<td align="left"><italic>O. lymani SSI</italic></td>
<td align="center">0.044</td>
<td align="center">0.047</td>
<td align="center"><bold>0.006</bold></td>
<td align="center">0.011</td>
<td align="center">0.010</td>
</tr>
<tr>
<td align="left"><italic>O. lymani PS</italic></td>
<td align="center">0.059</td>
<td align="center">0.053</td>
<td align="center">0.052</td>
<td align="center"><bold>0.013</bold></td>
<td align="center">0.007</td>
</tr>
<tr>
<td align="left"><italic>O. lymani SG</italic></td>
<td align="center">0.053</td>
<td align="center">0.046</td>
<td align="center">0.042</td>
<td align="center">0.029</td>
<td align="center"><bold>0.013</bold></td>
</tr>
</tbody>
</table>
</table-wrap>
<p>The phylogenetic method was less conservative: with the likelihood test finding the GMYC model a significant improvement over the null model (H<sub>GMYC</sub> 505.08 v H<sub>0</sub> 492.12, <italic>P</italic> &#x0003C; 0.0001). The test identified 14 entities (CI 9&#x02013;24) and 11 clusters (CI 6&#x02013;16).</p>
</sec>
<sec>
<title>Gene flow</title>
<p>Visual inspection of the phylogenies suggested distinct geographically distinct clades with the exception of five individuals from the Patagonian shelf that are more closely related to South Georgia haplotypes than other Patagonian haplotypes. As coalescent simulations suggest this pattern is non-random, we investigated migration between these regions.</p>
<p>The full model, allowing co-estimation of population size variation and free migration, indicates that migration is asymmetrical from South Georgia to the Patagonian shelf (Table <xref ref-type="table" rid="T3">3</xref>). The posterior distribution of migration rates has clear modes: at 0.3 for migration from Patagonia to South Georgia, and 27 for South Georgia to Patagonia (Figure <xref ref-type="fig" rid="F6">6</xref>). A comparative exploration of migration scenarios strongly supports asymmetrical migration scenarios over symmetrical migration, free migration (panmixia) or no migration (Table <xref ref-type="table" rid="T4">4</xref>). The scenario of migration from South Georgia to the Patagonian shelf was favored (probability of 0.59) and, although support for this model was only marginal compared with migration from the Patagonian shelf to South Georgia (probability of 0.41), it is hard to reconcile the pattern visualized in Figure <xref ref-type="fig" rid="F5">5</xref> with migration from the Patagonian shelf as no &#x0201C;Patagonian shelf&#x0201D;-like haplotypes were found around South Georgia. This seemingly marginal result lies at least partly with the incompatibility of the distributions being compared (Figure <xref ref-type="fig" rid="F6">6</xref>). Given that the posterior distribution must be positive, and the large amount of variation expected in a single marker analysis, even a &#x0201C;true&#x0201D; value of 0 will have a wide distribution and positive mean. In this case the distribution of 97% of the data was similar for both scenarios, producing a high mean (19.98) for the PS-&#x0003E;SG migration estimate (Table <xref ref-type="table" rid="T3">3</xref>) even though the mode (highest posterior estimate) was near 0, compared with SG-&#x0003E;PS where the mode (27) was much closer to the mean (33.1). Low confidence due to wide distribution of the posterior is largely down to the use of a single marker.</p>
<table-wrap position="float" id="T3">
<label>Table 3</label>
<caption><p><bold>Posterior distributions of migration estimates from the full migration model</bold>.</p></caption>
<table frame="hsides" rules="groups">
<thead>
<tr>
<th/>
<th align="left"><bold>2.5%</bold></th>
<th align="center"><bold>25%</bold></th>
<th align="center"><bold>Mode</bold></th>
<th align="center"><bold>Median</bold></th>
<th align="center"><bold>Mean</bold></th>
<th align="center"><bold>75%</bold></th>
<th align="center"><bold>97.5%</bold></th>
</tr>
</thead>
<tbody>
<tr>
<td align="left">SG-&#x0003E;PS</td>
<td align="center">0</td>
<td align="center">9.33</td>
<td align="center">27</td>
<td align="center">31</td>
<td align="center">33.1</td>
<td align="center">39.33</td>
<td align="center">74.67</td>
</tr>
<tr>
<td align="left">PS-&#x0003E;SG</td>
<td align="center">0</td>
<td align="center">0</td>
<td align="center">0.33</td>
<td align="center">13</td>
<td align="center">19.98</td>
<td align="center">12.67</td>
<td align="center">71.22</td>
</tr>
</tbody>
</table>
<table-wrap-foot>
<p><italic>The data is broken up into quantiles (25%, 75% of the distribution) and the values for the top and bottom 2.5% of the distribution is given. Of note is the large difference in the mode and the similarity of the extent of the distribution (97.5%) resulting in mean values skewed to the right</italic>.</p>
</table-wrap-foot>
</table-wrap>
<fig id="F6" position="float">
<label>Figure 6</label>
<caption><p><bold>The distribution of the posterior migration estimates from the analysis of the full migration model allowing migration between South Georgia to Patagonian shelf and Patagonian shelf to South Georgia</bold>.</p></caption>
<graphic xlink:href="fevo-03-00063-g0006.tif"/>
</fig>
<table-wrap position="float" id="T4">
<label>Table 4</label>
<caption><p><bold>Summary of the likelihoods of five competing migration models</bold>.</p></caption>
<table frame="hsides" rules="groups">
<thead>
<tr>
<th align="left"><bold>Model</bold></th>
<th align="center"><bold>Bezier lmL</bold></th>
<th align="center"><bold>Bayes factor</bold></th>
<th align="center"><bold>Choice</bold></th>
<th align="center"><bold>Model Prob</bold>.</th>
</tr>
</thead>
<tbody>
<tr>
<td align="left">FULL</td>
<td align="center">&#x02212;2491.92</td>
<td align="center">&#x02212;13.18</td>
<td align="center">3</td>
<td align="center">0.00</td>
</tr>
<tr>
<td align="left">PS-&#x0003E;SG</td>
<td align="center">&#x02212;2485.71</td>
<td align="center">&#x02212;0.76</td>
<td align="center">2</td>
<td align="center">0.41</td>
</tr>
<tr>
<td align="left">SG-&#x0003E;PS</td>
<td align="center">&#x02212;2485.33</td>
<td align="center">0</td>
<td align="center">1</td>
<td align="center">0.59</td>
</tr>
<tr>
<td align="left">PANMIXIA</td>
<td align="center">&#x02212;2504.85</td>
<td align="center">&#x02212;39.04</td>
<td align="center">4</td>
<td align="center">0.00</td>
</tr>
<tr>
<td align="left">NONE</td>
<td align="center">&#x02212;2585.86</td>
<td align="center">&#x02212;201.06</td>
<td align="center">5</td>
<td align="center">0.00</td>
</tr>
</tbody>
</table>
<table-wrap-foot>
<p><italic>Log likelihoods using thermodynamic integration method ranked with a Bayes factor (BF) is determined from the Bezier method</italic>.</p>
</table-wrap-foot>
</table-wrap>
</sec>
<sec>
<title>Timing of divergences</title>
<p>The root height of the tree was estimated as 8.3 million years ago (ma) (5.77, 11.07 95% HPD). The main radiation of the <italic>O. lymani</italic>/<italic>O. carinifera</italic> clades (nodes B&#x02013;D, Figure <xref ref-type="fig" rid="F7">7</xref>) was inferred to be around 2.5 ma. The most recent common ancestor of the Patagonian shelf and South Georgia shelf populations (node G) was estimated at around 1.5 ma. A similar time frame was estimated for the common ancestor of the Antarctic <italic>O. carinifera</italic> and its Macquarie Ridge relative (node E). The two groups of Patagonian shelf <italic>O. lymani</italic> that cluster within the South Georgia clade indicate that there were at least two separate migration events, one around 0.75 ma (node J) and the second around 0.22 ma (node L) although the time dependency of the molecular clock would affect these date estimates strongly and they should be considered to be more recent.</p>
<fig id="F7" position="float">
<label>Figure 7</label>
<caption><p><bold>Bayesian chronogram of partitioned CO1 sequences derived from</bold> <italic><bold>Ophiura lymani</bold></italic> <bold>from South Georgia shelf (red) Patagonian shelf (green), South Sandwich Islands (purple) and Bouvetoya (mustard)</bold>. <italic>Ophiura carinifera</italic> samples are blue, <italic>Ophiura rouchi</italic> samples are black. The distribution of uncertainty of node placement is indicated around each of the main nodes in light blue. The distribution of uncertainty around the nodes indicating migration events from South Georgia to Patagonian shelf are indicated in red. The scale bar is in millions of years.</p></caption>
<graphic xlink:href="fevo-03-00063-g0007.tif"/>
</fig>
</sec>
</sec>
<sec sec-type="discussion" id="s4">
<title>Discussion</title>
<p>The works of Hedgpeth (<xref ref-type="bibr" rid="B34">1969</xref>) and Griffiths et al. (<xref ref-type="bibr" rid="B32">2009</xref>) provide a biogeographic hypothesis in which the ACC strongly affect the distribution of Southern Ocean benthos by facilitating gene flow over vast geographic distances, and yet inhibiting gene flow particularly in and out of Antarctica. These processes have resulted in many species with circumpolar distributions, and assemblages that indicate a South American component, particularly at South Georgia. Our expectations under the South American input hypothesis were that there would be (1) appreciably more genetic diversity and (2) deeper coalescent times in the Patagonian Shelf populations, and (3) indication of asymmetrical movement along the direction of the ACC. However, our phylogeographic examination of <italic>O. lymani</italic> and <italic>O. carinifera</italic> demonstrated (1) very high diversity in all populations, (2) deeper coalescence times for the South Georgia population than any other, and (3) a mitochondrial signal of asymmetrical migration that appears to go <italic>against</italic> the ACC from South Georgia to the Patagonian shelf. Importantly, it appears that there are several independently evolving lineages of <italic>O. lymani</italic>, with the diverging Patagonian shelf lineage breaking out of the Antarctic realm and thereby rejecting the South American input hypothesis in this case.</p>
<sec>
<title>Crypsis</title>
<p>The populations of Bouvetoya and South Sandwich Islands are as genetically distinct from each other, and the Patagonian shelf and South Georgia populations, as each is also from <italic>O. carinifera</italic>, a clearly different morpho-species. Cryptic or unrecognized species (Janosik and Halanych, <xref ref-type="bibr" rid="B42">2010</xref>) are commonly identified by mitochondrial DNA studies of benthic invertebrates (For Antarctic examples see Linse et al., <xref ref-type="bibr" rid="B48">2007</xref>; Leese and Held, <xref ref-type="bibr" rid="B46">2008</xref>; Krabbe et al., <xref ref-type="bibr" rid="B44">2010</xref>; Arango et al., <xref ref-type="bibr" rid="B2">2011</xref>; Havermans et al., <xref ref-type="bibr" rid="B33">2011</xref>; Hemery et al., <xref ref-type="bibr" rid="B35">2012</xref>). There is considerable debate as to whether or not biological species should be delimited based solely on mitochondrial DNA sequence data (DeSalle et al., <xref ref-type="bibr" rid="B17">2005</xref>; Taylor and Harris, <xref ref-type="bibr" rid="B74">2012</xref>; Zhang et al., <xref ref-type="bibr" rid="B81">2013</xref>). In this case, based on clear genetic clades that correspond to discrete geographic areas, distinct morphology, and prior evidence that CO1 reliably distinguishes ophiuroid species (Ward et al., <xref ref-type="bibr" rid="B78">2008</xref>), we can assume that <italic>O. lymani</italic> could be conservatively divided into a Bouvetoya species, a South Sandwich Island species and a diverging Patagonian shelf/South Georgia species.</p>
</sec>
<sec>
<title>Against the flow</title>
<p>It appears counter intuitive that the most closely related populations are separated by what should be an intense barrier to gene flow. Similarly, despite currents that could facilitate movement from Patagonia to South Georgia, or South Georgia to the South Sandwich Islands or Bouvetoya, there is no evidence of movement in this direction or that currents themselves are responsible for population structuring. The only evidence of migration is from South Georgia to the Patagonian shelf, against the ACC and across a gradient between two water masses that strongly differ in temperature and salinity (Orsi et al., <xref ref-type="bibr" rid="B55">1995</xref>).</p>
<p>There are three possible mechanisms for dispersal: larval transport, active adult dispersal or passive adult dispersal. There is no published literature regarding the reproductive biology of <italic>O. lymani</italic> or <italic>O. carinifera</italic>. However, lack of obvious genetic structure of <italic>O. carinifera</italic> haplotypes around Antarctic continental shelves and <italic>O. lymani</italic> haplotypes between Shag Rocks and South Georgia shelf regions indicate moderate dispersal is possible. The narrow but deep gap between South Georgia and Shag Rock shelf fragments is known to inhibit gene flow in more mobile species (Allcock et al., <xref ref-type="bibr" rid="B1">1997</xref>; Strugnell et al., <xref ref-type="bibr" rid="B71">2012</xref>), and yet is observed in <italic>O. lymani</italic>, suggesting that <italic>O. lymani</italic> may have a more mobile life stage. Transport of larvae in ships ballast water is a possible mechanism, although the genetic distances of the haplotypes suggest the timing of the movement substantially predates shipping in the area.</p>
<p>Dispersal could occur by moving down slope and along abyssal sea floor. Such a deep sea link has been used to explain the close genetic connectivity between sub-Antarctic and Antarctic morphotypes of the echinoid <italic>Sterechinus</italic> (D&#x000ED;az et al., <xref ref-type="bibr" rid="B18">2011</xref>). The limits of adult <italic>O. lymani/carinifera</italic> dispersal are unknown. They are moderately large (20&#x0002B; mm disc size) and muscular, and certainly able to move across the sea floor (Sands, personal observation). Their depth range is unknown, largely due to the rarity of deep sampling events in the region. In the collections available to this study, <italic>O. lymani</italic> and <italic>O</italic>. <italic>carinifera</italic> were only recovered in the shallower trawls (194&#x02013;723 and 153&#x02013;556 m depth, respectively), and absent from deeper trawls (20 trawls between 1000 and 3213 m) indicating that they prefer shelf and upper slope depths. Therefore, while dispersal of adults remains possible, it appears unlikely.</p>
<p>Rafting, where individuals are attached to or entangled in kelp, pumice or other buoyant substrata and drift from one area to another, is a third possible mechanism. Drifting kelp is common in the Southern Ocean (Fraser et al., <xref ref-type="bibr" rid="B26">2009</xref>) and is a known dispersal mechanism for a variety of benthos (Barnes et al., <xref ref-type="bibr" rid="B5">2006</xref>; Leese et al., <xref ref-type="bibr" rid="B45">2010</xref>). However, rafting from South Georgia to Patagonian shelf would require either a reversal of the ACC or an almost complete circumpolar journey that, if sufficiently frequent to result in detectable gene flow, would be considerably higher in other closer shelf areas more convincingly in the path of the ACC and associated eddies, such as South Sandwich Islands and Bouvetoya.</p>
<p>Two other scenarios may explain the apparent migration pattern between South Georgia and Patagonia but without the need to invoke migration against the flow of the ACC. One is vicariance as a result of continental drift with retained ancestral lineages. Our results suggest that the most recent common ancestor of the clade (Node B, 2.8 ma) considerably post-dates the estimated opening of the North Scotia Ridge sufficient to allow the northward flow of the ACC around 17 Ma (Eagles, <xref ref-type="bibr" rid="B21">2010</xref>), and coalescent analyses indicate that retained ancestral lineages are an unlikely explanation. A second scenario is that both populations are supplied by an un-sampled source as in the case of the limpet <italic>Nacella magellanica</italic>, where samples from the Atlantic South American coast and Falkland Islands show two highly divergent lineages, but also more closely related lineages and one shared haplotype. Sampling the Pacific Patagonian shelf revealed that the shared haplotype between the Atlantic population and the Falkland Islands population was the most frequent haplotype from the Pacific region (Gonz&#x000E1;lez-Wevar et al., <xref ref-type="bibr" rid="B30">2012</xref>). We hope that further sampling opportunities will open up to test this possibility.</p>
</sec>
<sec>
<title>Out of antarctica</title>
<p>A clear pattern in our data is the Antarctic or sub-Antarctic origin of the <italic>O. lymani / O. carinifera</italic> clade. The change in morphology from <italic>O. lymani</italic> type to <italic>O. carinifera</italic> type occurred between 2.5 and 1.6 ma, and was maintained on the Antarctic continental shelf from around 0.8 ma. Movement out of the sub Antarctic occurred between 1.4 and 0.8 ma, with two subsequent migration events 0.74 and 0.22 ma. Poulin et al. (<xref ref-type="bibr" rid="B61">2014</xref>) summarized the divergences of several taxa across the Drake Passage and found that divergences generally occurred around 4&#x02013;5 ma, although the brittle star <italic>Astrotoma agassizi</italic> and the sea star <italic>Odontaster validus</italic> had more recent divergences of 1.5&#x02013;1 ma using a 3.5% mutation rate (Hunter and Halanych, <xref ref-type="bibr" rid="B41">2008</xref>; Janosik et al., <xref ref-type="bibr" rid="B43">2011</xref>) which, if standardized to the mutation rate used here, would be in line with our results. Although there is some evidence of benthic groups radiating out of the Antarctic during the Tertiary (Strugnell et al., <xref ref-type="bibr" rid="B70">2008</xref>; G&#x000F6;bbeler and Klussmann-Kolb, <xref ref-type="bibr" rid="B29">2010</xref>; Gonz&#x000E1;lez-Wevar et al., <xref ref-type="bibr" rid="B31">2010</xref>), and there are populations of some benthic species either side of the Polar Front (see Supplementary Tables in Sands et al., <xref ref-type="bibr" rid="B67">2013</xref>), to our knowledge this is the first evidence strongly suggesting Quaternary asymmetrical movement <italic>out</italic> of the Antarctic and into warmer waters. At a time when regional warming is causing a detectable shift in range distributions toward the poles (Thomas, <xref ref-type="bibr" rid="B76">2010</xref>), and there is an expectation that indigenous Antarctic cold-adapted species may be threatened by the warming marine environment (Peck, <xref ref-type="bibr" rid="B58">2005</xref>), the example of <italic>O. lymani</italic> shows that at least one species has had the ability to migrate and adapt to warmer temperature regimes over time.</p>
</sec>
<sec>
<title>Caveats and interpretation limits</title>
<p>Although evidence for crypsis, asymmetrical migration and radiation out of Antarctica appears compelling, there are caveats that apply to the current results and interpretations. The DNA sequence marker cytochrome c oxidase subunit 1 is widely used and heavily relied on in species identification and biodiversity studies. However, we emphasize that it is a single line of evidence and despite our best attempts at controlling for stochastic or confounding processes, without other gene sequences we cannot be overly confident that our gene trees reflect the species tree. Similarly, our mutation rates are best estimates, and the process they are derived from (the closing of the Isthmus of Panama) is in itself a best estimate, so the dates and error distributions provided here should be treated with great caution. The analyses here assume mutation drift equilibrium and selective neutrality, neither of which are safe assumptions in real populations. Finally, the sampling, although the best available, is still limited due to the enormous time and expense of carrying out surveys in the Southern Ocean and remote southern areas of other oceans. Perhaps the most significant sample related limitation for this study is the geopolitical boundaries on the South American continent, which have prevented wide sampling of the Patagonian shelf. However, with these caveats in mind, we can still interpret clear patterns of diversification within the Southern Ocean and what appears to be a recent movement from of lineages from South Georgia to the Patagonian Shelf.</p>
</sec>
</sec>
<sec sec-type="conclusions" id="s5">
<title>Conclusions</title>
<p>We have found evidence that the ophiuroid <italic>O. lymani</italic> has radiated in the Southern Ocean, with deep sea between isolated lineages likely to be the mechanism maintaining distinct populations or cryptic species. Our results indicate occasional migration over evolutionary time against the prevailing ACC and across the perceived barrier of the Polar Front. We provide evidence that <italic>O. lymani</italic> is a species complex of which <italic>O. carinifera</italic> is an element, with the South Georgia population highly divergent from other geographic populations in the Southern Ocean, and the likely source of an invasion into the cool temperate seas of South America. These results and hypotheses may cause a rethink in our general view of Antarctic biogeography, but will require further testing including better South American sampling, replication of markers, and extension to other taxa.</p>
</sec>
<sec>
<title>Author contributions</title>
<p>CS conceived the study, collected the majority of the specimens, acquired funding for the molecular sequencing, conducted the analyses and wrote the manuscript. TO provided additional specimens and sequences, helped shape the manuscript. DB assisted in the collections both financially and practically, contributed to the text and provided intellectual input throughout the project. RL identified the majority of the specimens, provided unique and relevant observations and contributed to the production of the manuscript.</p>
<sec>
<title>Conflict of interest statement</title>
<p>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.</p></sec>
</sec>
</body>
<back>
<ack>
<p>We would like to thank the Captains and crews of the RRS James Clarke Ross (JR144) and RV Polarstern (PS77) for their commitment to excellent science discovery, and the science leaders of the cruises that contributed to the collection of this material. Thanks to Elena Zaikina for assisting RL, Rachel Downey, Camille Moreau and Huw Griffiths for assistance in sorting material at BAS. Thanks also to Olivia Martin-Sanchez who assisted with the GIS data. Jennifer Jackson, Elise Biersma, Peter Convey and Will Goodall-Copestake provided useful discussion throughout the analysis and writing process. Funding for the sequencing service was supported by Census of Antarctic Marine Life. CS was supported by an Antarctic Science Bursary, a SCAR MarBIN grant and a SynTax award scheme supported by BBSRC and NERC.</p>
</ack>
<ref-list>
<title>References</title>
<ref id="B1">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Allcock</surname> <given-names>A. L.</given-names></name> <name><surname>Brierley</surname> <given-names>A. S.</given-names></name> <name><surname>Thorpe</surname> <given-names>J. P.</given-names></name> <name><surname>Rodhouse</surname> <given-names>P. G.</given-names></name></person-group> (<year>1997</year>). <article-title>Restricted gene flow and evolutionary divergence between geographically separated populations of the Antarctic octopus</article-title>. <source>Pareledone turqueti. Mar. Biol</source>. <volume>129</volume>, <fpage>97</fpage>&#x02013;<lpage>102</lpage>. <pub-id pub-id-type="doi">10.1007/s002270050150</pub-id></citation>
</ref>
<ref id="B2">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Arango</surname> <given-names>C. P.</given-names></name> <name><surname>Soler-Membrives</surname> <given-names>A.</given-names></name> <name><surname>Miller</surname> <given-names>K. J.</given-names></name></person-group> (<year>2011</year>). <article-title>Genetic differentiation in the circum&#x02014;Antarctic sea spider <italic>Nymphon australe</italic> (Pycnogonida; Nymphonidae)</article-title>. <source>Deep Sea Res. Part II Top. Stud. Oceanogr</source>. <volume>58</volume>, <fpage>212</fpage>&#x02013;<lpage>219</lpage>. <pub-id pub-id-type="doi">10.1016/j.dsr2.2010.05.019</pub-id></citation>
</ref>
<ref id="B3">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Barker</surname> <given-names>P. F.</given-names></name> <name><surname>Thomas</surname> <given-names>E.</given-names></name></person-group> (<year>2004</year>). <article-title>Origin, signature and palaeoclimatic influence of the Antarctic Circumpolar Current</article-title>. <source>Earth Sci. Rev</source>. <volume>66</volume>, <fpage>143</fpage>&#x02013;<lpage>162</lpage>. <pub-id pub-id-type="doi">10.1016/j.earscirev.2003.10.003</pub-id></citation>
</ref>
<ref id="B4">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Barnes</surname> <given-names>D. K. A.</given-names></name> <name><surname>Griffiths</surname> <given-names>H. J.</given-names></name> <name><surname>Kaiser</surname> <given-names>S.</given-names></name></person-group> (<year>2009</year>). <article-title>Geographic range shift responses to climate change by Antarctic benthos: where we should look</article-title>. <source>Mar. Ecol. Prog. Ser</source>. <volume>393</volume>, <fpage>13</fpage>&#x02013;<lpage>26</lpage>. <pub-id pub-id-type="doi">10.3354/meps08246</pub-id></citation>
</ref>
<ref id="B5">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Barnes</surname> <given-names>D. K. A.</given-names></name> <name><surname>Hodgson</surname> <given-names>D. A.</given-names></name> <name><surname>Convey</surname> <given-names>P.</given-names></name> <name><surname>Allen</surname> <given-names>C. S.</given-names></name> <name><surname>Clarke</surname> <given-names>A.</given-names></name></person-group> (<year>2006</year>). <article-title>Incursion and excursion of Antarctic biota: past, present and future</article-title>. <source>Glob. Ecol. Biogeogr</source>. <volume>15</volume>, <fpage>121</fpage>&#x02013;<lpage>142</lpage>. <pub-id pub-id-type="doi">10.1111/j.1466-822X.2006.00216.x</pub-id></citation>
</ref>
<ref id="B6">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Beerli</surname> <given-names>P.</given-names></name></person-group> (<year>2009</year>). <article-title>How to use MIGRATE or why are Markov chain Monte Carlo programs difficult to use</article-title>. <source>Popul. Genet. Anim. Conserv</source>. <volume>17</volume>, <fpage>42</fpage>&#x02013;<lpage>79</lpage>.</citation>
</ref>
<ref id="B7">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Beerli</surname> <given-names>P.</given-names></name> <name><surname>Palczewski</surname> <given-names>M.</given-names></name></person-group> (<year>2010</year>). <article-title>Unified framework to evaluate panmixia and migration direction among multiple sampling locations</article-title>. <source>Genetics</source> <volume>185</volume>, <fpage>313</fpage>&#x02013;<lpage>326</lpage>. <pub-id pub-id-type="doi">10.1534/genetics.109.112532</pub-id><pub-id pub-id-type="pmid">20176979</pub-id></citation>
</ref>
<ref id="B8">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Belanger</surname> <given-names>C. L.</given-names></name> <name><surname>Jablonski</surname> <given-names>D.</given-names></name> <name><surname>Roy</surname> <given-names>K.</given-names></name> <name><surname>Berke</surname> <given-names>S. K.</given-names></name> <name><surname>Krug</surname> <given-names>A. Z.</given-names></name> <name><surname>Valentine</surname> <given-names>J. W.</given-names></name></person-group> (<year>2012</year>). <article-title>Global environmental predictors of benthic marine biogeographic structure</article-title>. <source>Proc. Natl. Acad. Sci. U.S.A</source>. <volume>109</volume>, <fpage>14046</fpage>&#x02013;<lpage>14051</lpage>. <pub-id pub-id-type="doi">10.1073/pnas.1212381109</pub-id><pub-id pub-id-type="pmid">22904189</pub-id></citation>
</ref>
<ref id="B9">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Brown</surname> <given-names>A.</given-names></name> <name><surname>Thatje</surname> <given-names>S.</given-names></name></person-group> (<year>2015</year>). <article-title>The effects of changing climate on faunal depth distributions determine winners and losers</article-title>. <source>Glob. Change Biol</source>. <volume>21</volume>, <fpage>173</fpage>&#x02013;<lpage>180</lpage>. <pub-id pub-id-type="doi">10.1111/gcb.12680</pub-id><pub-id pub-id-type="pmid">25044552</pub-id></citation>
</ref>
<ref id="B10">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Clarke</surname> <given-names>A.</given-names></name> <name><surname>Barnes</surname> <given-names>D. K. A.</given-names></name> <name><surname>Hodgson</surname> <given-names>D. A.</given-names></name></person-group> (<year>2005</year>). <article-title>How isolated is Antarctica?</article-title> <source>Trends Ecol. Evol</source>. <volume>20</volume>, <fpage>1</fpage>&#x02013;<lpage>3</lpage>. <pub-id pub-id-type="doi">10.1016/j.tree.2004.10.004</pub-id><pub-id pub-id-type="pmid">16701330</pub-id></citation>
</ref>
<ref id="B11">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Clarke</surname> <given-names>A.</given-names></name> <name><surname>Crame</surname> <given-names>J. A.</given-names></name></person-group> (<year>1989</year>). <article-title>The origin of the Southern Ocean marine fauna</article-title>. <source>Geol. Soc. Lond. Spec. Publ</source>. <volume>47</volume>, <fpage>253</fpage>&#x02013;<lpage>268</lpage>. <pub-id pub-id-type="doi">10.1144/GSL.SP.1989.047.01.19</pub-id></citation>
</ref>
<ref id="B12">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Clarke</surname> <given-names>A.</given-names></name> <name><surname>Johnston</surname> <given-names>N. M.</given-names></name></person-group> (<year>2003</year>). <article-title>Antarctic marine benthic diversity</article-title>. <source>Oceanogr. Mar. Biol. Annu. Rev</source>. <volume>41</volume>, <fpage>47</fpage>&#x02013;<lpage>114</lpage>.</citation>
</ref>
<ref id="B13">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Clement</surname> <given-names>M.</given-names></name> <name><surname>Posada</surname> <given-names>D.</given-names></name> <name><surname>Crandall</surname> <given-names>K. A.</given-names></name></person-group> (<year>2000</year>). <article-title>TCS: a computer program to estimate gene genealogies</article-title>. <source>Mol. Ecol</source>. <volume>9</volume>, <fpage>1657</fpage>&#x02013;<lpage>1659</lpage>. <pub-id pub-id-type="doi">10.1046/j.1365-294x.2000.01020.x</pub-id><pub-id pub-id-type="pmid">11050560</pub-id></citation>
</ref>
<ref id="B14">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Convey</surname> <given-names>P.</given-names></name> <name><surname>Gibson</surname> <given-names>J. A. E.</given-names></name> <name><surname>Hillenbrand</surname> <given-names>C.-D.</given-names></name> <name><surname>Hodgson</surname> <given-names>D. A.</given-names></name> <name><surname>Pugh</surname> <given-names>P. J. A.</given-names></name> <name><surname>Smellie</surname> <given-names>J. L.</given-names></name> <etal/></person-group>. (<year>2008</year>). <article-title>Antarctic terrestrial life - challenging the history of the frozen continent?</article-title> <source>Biol. Rev</source>. <volume>83</volume>, <fpage>103</fpage>&#x02013;<lpage>117</lpage>. <pub-id pub-id-type="doi">10.1111/j.1469-185X.2008.00034.x</pub-id><pub-id pub-id-type="pmid">18429764</pub-id></citation>
</ref>
<ref id="B15">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Darriba</surname> <given-names>D.</given-names></name> <name><surname>Taboada</surname> <given-names>G. L.</given-names></name> <name><surname>Doallo</surname> <given-names>R.</given-names></name> <name><surname>Posada</surname> <given-names>D.</given-names></name></person-group> (<year>2012</year>). <article-title>jModelTest 2: more models, new heuristics and parallel computing</article-title>. <source>Nat. Methods</source> <volume>9</volume>, <fpage>772</fpage>&#x02013;<lpage>772</lpage>. <pub-id pub-id-type="doi">10.1038/nmeth.2109</pub-id><pub-id pub-id-type="pmid">22847109</pub-id></citation>
</ref>
<ref id="B16">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Dell</surname> <given-names>R. K.</given-names></name></person-group> (<year>1972</year>). <article-title>Antarctic benthos</article-title>. <source>Adv. Mar. Biol</source>. <volume>10</volume>, <fpage>1</fpage>&#x02013;<lpage>216</lpage>. <pub-id pub-id-type="doi">10.1016/s0065-2881(08)60416-2</pub-id></citation>
</ref>
<ref id="B17">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>DeSalle</surname> <given-names>R.</given-names></name> <name><surname>Egan</surname> <given-names>M. G.</given-names></name> <name><surname>Siddall</surname> <given-names>M.</given-names></name></person-group> (<year>2005</year>). <article-title>The unholy trinity: taxonomy, species delimitation and DNA barcoding</article-title>. <source>Philos. Trans. R. Soc. B Biol. Sci</source>. <volume>360</volume>, <fpage>1905</fpage>&#x02013;<lpage>1916</lpage>. <pub-id pub-id-type="doi">10.1098/rstb.2005.1722</pub-id></citation>
</ref>
<ref id="B18">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>D&#x000ED;az</surname> <given-names>A.</given-names></name> <name><surname>F&#x000E9;ral</surname> <given-names>J.-P.</given-names></name> <name><surname>David</surname> <given-names>B.</given-names></name> <name><surname>Sauc&#x000E8;de</surname> <given-names>T.</given-names></name> <name><surname>Poulin</surname> <given-names>E.</given-names></name></person-group> (<year>2011</year>). <article-title>Evolutionary pathways among shallow and deep-sea echinoids of the genus Sterechinus in the Southern Ocean</article-title>. <source>Deep Sea Res. Part II Top. Stud. Oceanogr</source>. <volume>58</volume>, <fpage>205</fpage>&#x02013;<lpage>211</lpage>. <pub-id pub-id-type="doi">10.1016/j.dsr2.2010.10.012</pub-id></citation>
</ref>
<ref id="B19">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Drummond</surname> <given-names>A. J.</given-names></name> <name><surname>Rambaut</surname> <given-names>A.</given-names></name> <name><surname>Shapiro</surname> <given-names>B.</given-names></name> <name><surname>Pybus</surname> <given-names>O. G.</given-names></name></person-group> (<year>2005</year>). <article-title>Bayesian coalescent inference of past population dynamics from molecular sequences</article-title>. <source>Mol. Biol. Evol</source>. <volume>22</volume>, <fpage>1185</fpage>&#x02013;<lpage>1192</lpage>. <pub-id pub-id-type="doi">10.1093/molbev/msi103</pub-id><pub-id pub-id-type="pmid">15703244</pub-id></citation>
</ref>
<ref id="B20">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Drummond</surname> <given-names>A.</given-names></name> <name><surname>Rambaut</surname> <given-names>A.</given-names></name></person-group> (<year>2007</year>). <article-title>BEAST: Bayesian evolutionary analysis by sampling trees</article-title>. <source>BMC Evol. Biol</source>. <volume>7</volume>:<fpage>214</fpage>. <pub-id pub-id-type="doi">10.1186/1471-2148-7-214</pub-id><pub-id pub-id-type="pmid">17996036</pub-id></citation>
</ref>
<ref id="B21">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Eagles</surname> <given-names>G.</given-names></name></person-group> (<year>2010</year>). <article-title>South Georgia and Gondwana&#x00027;s Pacific margin: lost in translation?</article-title> <source>J. South Am. Earth Sci</source>. <volume>30</volume>, <fpage>65</fpage>&#x02013;<lpage>70</lpage>. <pub-id pub-id-type="doi">10.1016/j.jsames.2010.04.004</pub-id></citation>
</ref>
<ref id="B22">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Edgar</surname> <given-names>R. C.</given-names></name></person-group> (<year>2004</year>). <article-title>MUSCLE: a multiple sequence alignment method with reduced time and space complexity</article-title>. <source>BMC Bioinformatics</source> <volume>5</volume>:<fpage>113</fpage>. <pub-id pub-id-type="doi">10.1186/1471-2105-5-113</pub-id><pub-id pub-id-type="pmid">15318951</pub-id></citation>
</ref>
<ref id="B23">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Excoffier</surname> <given-names>L.</given-names></name> <name><surname>Laval</surname> <given-names>L. G.</given-names></name> <name><surname>Schneider</surname> <given-names>S.</given-names></name></person-group> (<year>2005</year>). <article-title>Arlequin ver. 3.0: An integrated software package for population genetics data analysis</article-title>. <source>Evol. Bioinformatics Online</source> <volume>1</volume>, <fpage>47</fpage>&#x02013;<lpage>50</lpage>. <pub-id pub-id-type="pmid">19325852</pub-id></citation>
</ref>
<ref id="B24">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Felsenstein</surname> <given-names>J.</given-names></name></person-group> (<year>1981</year>). <article-title>Evolutionary trees from DNA sequences: a maximum likelihood approach</article-title>. <source>J. Mol. Evol</source>. <volume>17</volume>, <fpage>368</fpage>&#x02013;<lpage>376</lpage>. <pub-id pub-id-type="doi">10.1007/BF01734359</pub-id><pub-id pub-id-type="pmid">7288891</pub-id></citation>
</ref>
<ref id="B25">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Folmer</surname> <given-names>O.</given-names></name> <name><surname>Black</surname> <given-names>M.</given-names></name> <name><surname>Hoeh</surname> <given-names>W.</given-names></name> <name><surname>Lutz</surname> <given-names>R.</given-names></name> <name><surname>Vrijenhoek</surname> <given-names>R.</given-names></name></person-group> (<year>1994</year>). <article-title>DNA primers for amplification of mitochondrial cytochrome c oxidase subunit I from diverse metazoan invertebrates</article-title>. <source>Mol. Mar. Biol. Biotechnol</source>. <volume>3</volume>, <fpage>294</fpage>&#x02013;<lpage>299</lpage>. <pub-id pub-id-type="pmid">7881515</pub-id></citation>
</ref>
<ref id="B26">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Fraser</surname> <given-names>C. I.</given-names></name> <name><surname>Hay</surname> <given-names>C. H.</given-names></name> <name><surname>Spencer</surname> <given-names>H. G.</given-names></name> <name><surname>Waters</surname> <given-names>J. M.</given-names></name></person-group> (<year>2009</year>). <article-title>Genetic and morphological analyses of the southern bull kelp Durvillaea antarctica (Phaeophyceae: Durvillaeales) in New Zealand reveal cryptic species</article-title>. <source>J. Phycol</source>. <volume>45</volume>, <fpage>436</fpage>&#x02013;<lpage>443</lpage>. <pub-id pub-id-type="doi">10.1111/j.1529-8817.2009.00658.x</pub-id></citation>
</ref>
<ref id="B27">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Fu</surname> <given-names>Y. X.</given-names></name></person-group> (<year>1996</year>). <article-title>New statistical tests of neutrality for DNA samples from a population</article-title>. <source>Genetics</source> <volume>143</volume>, <fpage>557</fpage>&#x02013;<lpage>570</lpage>. <pub-id pub-id-type="pmid">8722804</pub-id></citation>
</ref>
<ref id="B28">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Fujisawa</surname> <given-names>T.</given-names></name> <name><surname>Barraclough</surname> <given-names>T. G.</given-names></name></person-group> (<year>2013</year>). <article-title>Delimiting species using single-locus data and the generalized mixed Yule Coalescent Approach: a revised method and evaluation on simulated data sets</article-title>. <source>Syst. Biol</source>. <volume>62</volume>, <fpage>707</fpage>&#x02013;<lpage>724</lpage>. <pub-id pub-id-type="doi">10.1093/sysbio/syt033</pub-id><pub-id pub-id-type="pmid">23681854</pub-id></citation>
</ref>
<ref id="B29">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>G&#x000F6;bbeler</surname> <given-names>K.</given-names></name> <name><surname>Klussmann-Kolb</surname> <given-names>A.</given-names></name></person-group> (<year>2010</year>). <article-title>Out of Antarctica? &#x02013; New insights into the phylogeny and biogeography of the Pleurobranchomorpha (Mollusca, Gastropoda)</article-title>. <source>Mol. Phylogenet. Evol</source>. <volume>55</volume>, <fpage>996</fpage>&#x02013;<lpage>1007</lpage>. <pub-id pub-id-type="doi">10.1016/j.ympev.2009.11.027</pub-id><pub-id pub-id-type="pmid">19995612</pub-id></citation>
</ref>
<ref id="B30">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Gonz&#x000E1;lez-Wevar</surname> <given-names>C. A.</given-names></name> <name><surname>H&#x000FC;ne</surname> <given-names>M.</given-names></name> <name><surname>Ca&#x000F1;ete</surname> <given-names>J. I.</given-names></name> <name><surname>Mansilla</surname> <given-names>A.</given-names></name> <name><surname>Nakano</surname> <given-names>T.</given-names></name> <name><surname>Poulin</surname> <given-names>E.</given-names></name></person-group> (<year>2012</year>). <article-title>Towards a model of postglacial biogeography in shallow marine species along the Patagonian Province: lessons from the limpet <italic>Nacella magellanica</italic> (Gmelin, 1791)</article-title>. <source>BMC Evol. Biol</source>. <volume>12</volume>:<fpage>139</fpage>. <pub-id pub-id-type="doi">10.1186/1471-2148-12-139</pub-id><pub-id pub-id-type="pmid">22871029</pub-id></citation>
</ref>
<ref id="B31">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Gonz&#x000E1;lez-Wevar</surname> <given-names>C. A.</given-names></name> <name><surname>Nakano</surname> <given-names>T.</given-names></name> <name><surname>Ca&#x000F1;ete</surname> <given-names>J. I.</given-names></name> <name><surname>Poulin</surname> <given-names>E.</given-names></name></person-group> (<year>2010</year>). <article-title>Molecular phylogeny and historical biogeography of Nacella (Patellogastropoda: Nacellidae) in the Southern Ocean</article-title>. <source>Mol. Phylogenet. Evol</source>. <volume>56</volume>, <fpage>115</fpage>&#x02013;<lpage>124</lpage>. <pub-id pub-id-type="doi">10.1016/j.ympev.2010.02.001</pub-id><pub-id pub-id-type="pmid">20139020</pub-id></citation>
</ref>
<ref id="B32">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Griffiths</surname> <given-names>H. J.</given-names></name> <name><surname>Barnes</surname> <given-names>D. K. A.</given-names></name> <name><surname>Linse</surname> <given-names>K.</given-names></name></person-group> (<year>2009</year>). <article-title>Towards a generalized biogeography of the Southern Ocean benthos</article-title>. <source>J. Biogeogr</source>. <volume>36</volume>, <fpage>162</fpage>&#x02013;<lpage>177</lpage>. <pub-id pub-id-type="doi">10.1111/j.1365-2699.2008.01979.x</pub-id></citation>
</ref>
<ref id="B33">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Havermans</surname> <given-names>C.</given-names></name> <name><surname>Nagy</surname> <given-names>Z. T.</given-names></name> <name><surname>Sonet</surname> <given-names>G.</given-names></name> <name><surname>De Broyer</surname> <given-names>C.</given-names></name> <name><surname>Martin</surname> <given-names>P.</given-names></name></person-group> (<year>2011</year>). <article-title>DNA barcoding reveals new insights into the diversity of Antarctic species of Orchomene sensu lato (Crustacea: Amphipoda: Lysianassoidea)</article-title>. <source>Deep Sea Res. Part II Top. Stud. Oceanogr</source>. <volume>58</volume>, <fpage>230</fpage>&#x02013;<lpage>241</lpage>. <pub-id pub-id-type="doi">10.1016/j.dsr2.2010.09.028</pub-id></citation>
</ref>
<ref id="B34">
<citation citation-type="book"><person-group person-group-type="author"><name><surname>Hedgpeth</surname> <given-names>J. W.</given-names></name></person-group> (<year>1969</year>). <article-title>Introduction to Antarctic zoogeography</article-title>, in <source>Distribution of Selected Groups of Marine Invertebrates in Waters South of 35&#x000B0; S</source>. Antarctic Map Folio Series, eds <person-group person-group-type="editor"><name><surname>Bushnell</surname> <given-names>V. C.</given-names></name> <name><surname>Hedgpeth</surname> <given-names>J. W.</given-names></name></person-group> (<publisher-loc>New York, NY</publisher-loc>: <publisher-name>American Geographical Society</publisher-name>), <fpage>1</fpage>&#x02013;<lpage>29</lpage>.</citation>
</ref>
<ref id="B35">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Hemery</surname> <given-names>L. G.</given-names></name> <name><surname>El&#x000E9;aume</surname> <given-names>M.</given-names></name> <name><surname>Roussel</surname> <given-names>V.</given-names></name> <name><surname>Am&#x000E9;ziane</surname> <given-names>N.</given-names></name> <name><surname>Gallut</surname> <given-names>C.</given-names></name> <name><surname>Steinke</surname> <given-names>D.</given-names></name> <etal/></person-group>. (<year>2012</year>). <article-title>Comprehensive sampling reveals circumpolarity and sympatry in seven mitochondrial lineages of the Southern Ocean crinoid species <italic>Promachocrinus kerguelensis</italic> (Echinodermata)</article-title>. <source>Mol. Ecol</source>. <volume>21</volume>, <fpage>2502</fpage>&#x02013;<lpage>2518</lpage>. <pub-id pub-id-type="doi">10.1111/j.1365-294X.2012.05512.x</pub-id><pub-id pub-id-type="pmid">22469156</pub-id></citation>
</ref>
<ref id="B36">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Hiddink</surname> <given-names>J. G.</given-names></name> <name><surname>Burrows</surname> <given-names>M. T.</given-names></name> <name><surname>Garc&#x000ED;a Molinos</surname> <given-names>J.</given-names></name></person-group> (<year>2015</year>). <article-title>Temperature tracking by North Sea benthic invertebrates in response to climate change</article-title>. <source>Glob. Change Biol</source>. <volume>21</volume>, <fpage>117</fpage>&#x02013;<lpage>129</lpage>. <pub-id pub-id-type="doi">10.1111/gcb.12726</pub-id><pub-id pub-id-type="pmid">25179407</pub-id></citation>
</ref>
<ref id="B37">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Ho</surname> <given-names>J. W. K.</given-names></name> <name><surname>Adams</surname> <given-names>C. E.</given-names></name> <name><surname>Lew</surname> <given-names>J. B.</given-names></name> <name><surname>Matthews</surname> <given-names>T. J.</given-names></name> <name><surname>Ng</surname> <given-names>C. C.</given-names></name> <name><surname>Shahabi-Sirjani</surname> <given-names>A.</given-names></name> <etal/></person-group>. (<year>2006</year>). <article-title>SeqVis: visualization of compositional heterogeneity in large alignments of nucleotides</article-title>. <source>Bioinformatics</source> <volume>22</volume>, <fpage>2162</fpage>&#x02013;<lpage>2163</lpage>. <pub-id pub-id-type="doi">10.1093/bioinformatics/btl283</pub-id><pub-id pub-id-type="pmid">16766557</pub-id></citation>
</ref>
<ref id="B38">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Hoelzer</surname> <given-names>G. A.</given-names></name> <name><surname>Meinick</surname> <given-names>D. J.</given-names></name></person-group> (<year>1994</year>). <article-title>Patterns of speciation and limits to phylogenetic resolution</article-title>. <source>Trends Ecol. Evol</source>. <volume>9</volume>, <fpage>104</fpage>&#x02013;<lpage>107</lpage>. <pub-id pub-id-type="doi">10.1016/0169-5347(94)90207-0</pub-id><pub-id pub-id-type="pmid">21236789</pub-id></citation>
</ref>
<ref id="B39">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Hogg</surname> <given-names>O. T.</given-names></name> <name><surname>Barnes</surname> <given-names>D. K. A.</given-names></name> <name><surname>Griffiths</surname> <given-names>H. J.</given-names></name></person-group> (<year>2011</year>). <article-title>Highly diverse, poorly studied and uniquely threatened by climate change: an assessment of marine biodiversity on South Georgia&#x00027;s continental shelf</article-title>. <source>PLoS ONE</source> <volume>6</volume>:<fpage>e19795</fpage>. <pub-id pub-id-type="doi">10.1371/journal.pone.0019795</pub-id><pub-id pub-id-type="pmid">21647236</pub-id></citation>
</ref>
<ref id="B40">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Huelsenbeck</surname> <given-names>J. P.</given-names></name> <name><surname>Ronquist</surname> <given-names>F.</given-names></name></person-group> (<year>2001</year>). <article-title>MRBAYES: Bayesian inference of phylogenetic trees</article-title>. <source>Bioinformatics</source> <volume>17</volume>, <fpage>754</fpage>&#x02013;<lpage>755</lpage>. <pub-id pub-id-type="doi">10.1093/bioinformatics/17.8.754</pub-id><pub-id pub-id-type="pmid">11524383</pub-id></citation>
</ref>
<ref id="B41">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Hunter</surname> <given-names>R. L.</given-names></name> <name><surname>Halanych</surname> <given-names>K. M.</given-names></name></person-group> (<year>2008</year>). <article-title>Evaluating connectivity in the brooding brittle star <italic>Astrotoma agassizii</italic> across the Drake Passage in the Southern Ocean</article-title>. <source>J. Hered</source>. <volume>99</volume>, <fpage>137</fpage>&#x02013;<lpage>148</lpage>. <pub-id pub-id-type="doi">10.1093/jhered/esm119</pub-id><pub-id pub-id-type="pmid">18245798</pub-id></citation>
</ref>
<ref id="B42">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Janosik</surname> <given-names>A. M.</given-names></name> <name><surname>Halanych</surname> <given-names>K. M.</given-names></name></person-group> (<year>2010</year>). <article-title>Unrecognized Antarctic biodiversity: a case study of the genus <italic>Odontaster</italic> (Odontasteridae; Asteroidea)</article-title>. <source>Integr. Comp. Biol</source>. <volume>50</volume>, <fpage>981</fpage>&#x02013;<lpage>992</lpage>. <pub-id pub-id-type="doi">10.1093/icb/icq119</pub-id><pub-id pub-id-type="pmid">21558254</pub-id></citation>
</ref>
<ref id="B43">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Janosik</surname> <given-names>A. M.</given-names></name> <name><surname>Mahon</surname> <given-names>A. R.</given-names></name> <name><surname>Halanych</surname> <given-names>K. M.</given-names></name></person-group> (<year>2011</year>). <article-title>Evolutionary history of Southern Ocean <italic>Odontaster</italic> sea star species (Odontasteridae; Asteroidea)</article-title>. <source>Polar Biol</source>. <volume>34</volume>, <fpage>575</fpage>&#x02013;<lpage>586</lpage>. <pub-id pub-id-type="doi">10.1007/s00300-010-0916-7</pub-id></citation>
</ref>
<ref id="B44">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Krabbe</surname> <given-names>K.</given-names></name> <name><surname>Leese</surname> <given-names>F.</given-names></name> <name><surname>Mayer</surname> <given-names>C.</given-names></name> <name><surname>Tollrian</surname> <given-names>R.</given-names></name> <name><surname>Held</surname> <given-names>C.</given-names></name></person-group> (<year>2010</year>). <article-title>Cryptic mitochondrial lineages in the widespread pycnogonid Colossendeis megalonyx Hoek, 1881 from Antarctic and Subantarctic waters</article-title>. <source>Polar Biol</source>. <volume>33</volume>, <fpage>281</fpage>&#x02013;<lpage>292</lpage>. <pub-id pub-id-type="doi">10.1007/s00300-009-0703-5</pub-id></citation>
</ref>
<ref id="B45">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Leese</surname> <given-names>F.</given-names></name> <name><surname>Agrawal</surname> <given-names>S.</given-names></name> <name><surname>Held</surname> <given-names>C.</given-names></name></person-group> (<year>2010</year>). <article-title>Long-distance island hopping without dispersal stages: transportation across major zoogeographic barriers in a Southern Ocean isopod</article-title>. <source>Naturwissenschaften</source> <volume>97</volume>, <fpage>583</fpage>&#x02013;<lpage>594</lpage>. <pub-id pub-id-type="doi">10.1007/s00114-010-0674-y</pub-id><pub-id pub-id-type="pmid">20454771</pub-id></citation>
</ref>
<ref id="B46">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Leese</surname> <given-names>F.</given-names></name> <name><surname>Held</surname> <given-names>C.</given-names></name></person-group> (<year>2008</year>). <article-title>Identification and characterization of microsatellites from the Antarctic isopod Ceratoserolis trilobitoides: nuclear evidence for cryptic species</article-title>. <source>Conserv. Genet</source>. <volume>9</volume>, <fpage>1369</fpage>&#x02013;<lpage>1372</lpage>. <pub-id pub-id-type="doi">10.1007/s10592-007-9491-z</pub-id></citation>
</ref>
<ref id="B47">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Librado</surname> <given-names>P.</given-names></name> <name><surname>Rozas</surname> <given-names>J.</given-names></name></person-group> (<year>2009</year>). <article-title>DnaSP v5: a software for comprehensive analysis of DNA polymorphism data</article-title>. <source>Bioinformatics</source> <volume>25</volume>, <fpage>1451</fpage>&#x02013;<lpage>1452</lpage>. <pub-id pub-id-type="doi">10.1093/bioinformatics/btp187</pub-id><pub-id pub-id-type="pmid">19346325</pub-id></citation>
</ref>
<ref id="B48">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Linse</surname> <given-names>K.</given-names></name> <name><surname>Cope</surname> <given-names>T.</given-names></name> <name><surname>L&#x000F6;rz</surname> <given-names>A.-N.</given-names></name> <name><surname>Sands</surname> <given-names>C.</given-names></name></person-group> (<year>2007</year>). <article-title>Is the Scotia Sea a centre of Antarctic marine diversification? Some evidence of cryptic speciation in the circum-Antarctic bivalve <italic>Lissarca notorcadensis</italic> (Arcoidea: Philobryidae)</article-title>. <source>Polar Biol</source>. <volume>30</volume>, <fpage>1059</fpage>&#x02013;<lpage>1068</lpage>. <pub-id pub-id-type="doi">10.1007/s00300-007-0265-3</pub-id></citation>
</ref>
<ref id="B49">
<citation citation-type="web"><person-group person-group-type="author"><name><surname>Maddison</surname> <given-names>W. P.</given-names></name> <name><surname>Maddison</surname> <given-names>D. R.</given-names></name></person-group> (<year>2014</year>). <source>Mesquite: A Modular System for Evolutionary Analysis</source>. Version 3.01. Available online at: <ext-link ext-link-type="uri" xlink:href="http://mesquiteproject.org">mesquiteproject.org</ext-link></citation>
</ref>
<ref id="B50">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Meredith</surname> <given-names>M. P.</given-names></name> <name><surname>King</surname> <given-names>J. C.</given-names></name></person-group> (<year>2005</year>). <article-title>Rapid climate change in the ocean west of the Antarctic Peninsula during the second half of the 20th century</article-title>. <source>Geophys. Res. Lett</source>. <volume>32</volume>, <fpage>L19604</fpage>. <pub-id pub-id-type="doi">10.1029/2005GL024042</pub-id></citation>
</ref>
<ref id="B51">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Minin</surname> <given-names>V. N.</given-names></name> <name><surname>Bloomquist</surname> <given-names>E. W.</given-names></name> <name><surname>Suchard</surname> <given-names>M. A.</given-names></name></person-group> (<year>2008</year>). <article-title>Smooth skyride through a rough skyline: Bayesian coalescent-based inference of population dynamics</article-title>. <source>Mol. Biol. Evol</source>. <volume>25</volume>, <fpage>1459</fpage>&#x02013;<lpage>1471</lpage>. <pub-id pub-id-type="doi">10.1093/molbev/msn090</pub-id><pub-id pub-id-type="pmid">18408232</pub-id></citation>
</ref>
<ref id="B52">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Morley</surname> <given-names>S. A.</given-names></name> <name><surname>Belchier</surname> <given-names>M.</given-names></name> <name><surname>Sands</surname> <given-names>C.</given-names></name> <name><surname>Barnes</surname> <given-names>D. K. A.</given-names></name> <name><surname>Peck</surname> <given-names>L. S.</given-names></name></person-group> (<year>2014</year>). <article-title>Geographic isolation and physiological mechanisms underpinning species distributions at the range limit hotspot of South Georgia</article-title>. <source>Rev. Fish Biol. Fish</source>. <volume>24</volume>, <fpage>485</fpage>&#x02013;<lpage>492</lpage>. <pub-id pub-id-type="doi">10.1007/s11160-013-9308-8</pub-id></citation>
</ref>
<ref id="B53">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Naughton</surname> <given-names>K. M.</given-names></name> <name><surname>O&#x00027;Hara</surname> <given-names>T. D.</given-names></name> <name><surname>Appleton</surname> <given-names>B.</given-names></name> <name><surname>Cisternas</surname> <given-names>P. A.</given-names></name></person-group> (<year>2014</year>). <article-title>Antitropical distributions and species delimitation in a group of ophiocomid brittle stars (Echinodermata: Ophiuroidea: Ophiocomidae)</article-title>. <source>Mol. Phylogenet. Evol</source>. <volume>78</volume>, <fpage>232</fpage>&#x02013;<lpage>244</lpage>. <pub-id pub-id-type="doi">10.1016/j.ympev.2014.05.020</pub-id><pub-id pub-id-type="pmid">24875252</pub-id></citation>
</ref>
<ref id="B54">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Norberg</surname> <given-names>J.</given-names></name> <name><surname>Urban</surname> <given-names>M. C.</given-names></name> <name><surname>Vellend</surname> <given-names>M.</given-names></name> <name><surname>Klausmeier</surname> <given-names>C. A.</given-names></name> <name><surname>Loeuille</surname> <given-names>N.</given-names></name></person-group> (<year>2012</year>). <article-title>Eco-evolutionary responses of biodiversity to climate change</article-title>. <source>Nat. Clim. Change</source> <volume>2</volume>, <fpage>747</fpage>&#x02013;<lpage>751</lpage>. <pub-id pub-id-type="doi">10.1038/nclimate1588</pub-id></citation>
</ref>
<ref id="B55">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Orsi</surname> <given-names>A. H.</given-names></name> <name><surname>Whitworth</surname> <given-names>T.</given-names> <suffix>III</suffix></name> <name><surname>Nowlin</surname> <given-names>W. D.</given-names> <suffix>Jr.</suffix></name></person-group> (<year>1995</year>). <article-title>On the meridional extent and fronts of the Antarctic Circumpolar Current</article-title>. <source>Deep Sea Res. Part Oceanogr. Res. Pap</source>. <volume>42</volume>, <fpage>641</fpage>&#x02013;<lpage>673</lpage>. <pub-id pub-id-type="doi">10.1016/0967-0637(95)00021-W</pub-id></citation>
</ref>
<ref id="B56">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Page</surname> <given-names>T.</given-names></name> <name><surname>Linse</surname> <given-names>K.</given-names></name></person-group> (<year>2002</year>). <article-title>More evidence of speciation and dispersal across the Antarctic Polar Front through molecular systematics of Southern Ocean Limatula (Bivalvia: Limidae)</article-title>. <source>Polar Biol</source>. <volume>25</volume>, <fpage>818</fpage>&#x02013;<lpage>826</lpage>. <pub-id pub-id-type="doi">10.1007/s00300-002-0414-7</pub-id></citation>
</ref>
<ref id="B57">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Paulay</surname> <given-names>G.</given-names></name> <name><surname>Meyer</surname> <given-names>C.</given-names></name></person-group> (<year>2006</year>). <article-title>Dispersal and divergence across the greatest ocean region: do larvae matter?</article-title> <source>Integr. Comp. Biol</source>. <volume>46</volume>, <fpage>269</fpage>&#x02013;<lpage>281</lpage>. <pub-id pub-id-type="doi">10.1093/icb/icj027</pub-id><pub-id pub-id-type="pmid">21672741</pub-id></citation>
</ref>
<ref id="B58">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Peck</surname> <given-names>L.</given-names></name></person-group> (<year>2005</year>). <article-title>Prospects for surviving climate change in Antarctic aquatic species</article-title>. <source>Front. Zool</source>. <volume>2</volume>:<fpage>9</fpage>. <pub-id pub-id-type="doi">10.1186/1742-9994-2-9</pub-id><pub-id pub-id-type="pmid">15938749</pub-id></citation>
</ref>
<ref id="B59">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Pinsky</surname> <given-names>M. L.</given-names></name> <name><surname>Worm</surname> <given-names>B.</given-names></name> <name><surname>Fogarty</surname> <given-names>M. J.</given-names></name> <name><surname>Sarmiento</surname> <given-names>J. L.</given-names></name> <name><surname>Levin</surname> <given-names>S. A.</given-names></name></person-group> (<year>2013</year>). <article-title>Marine taxa track local climate velocities</article-title>. <source>Science</source> <volume>341</volume>, <fpage>1239</fpage>&#x02013;<lpage>1242</lpage>. <pub-id pub-id-type="doi">10.1126/science.1239352</pub-id><pub-id pub-id-type="pmid">24031017</pub-id></citation>
</ref>
<ref id="B60">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Pons</surname> <given-names>J.</given-names></name> <name><surname>Barraclough</surname> <given-names>T. G.</given-names></name> <name><surname>Gomez-Zurita</surname> <given-names>J.</given-names></name> <name><surname>Cardoso</surname> <given-names>A.</given-names></name> <name><surname>Duran</surname> <given-names>D. P.</given-names></name> <name><surname>Hazell</surname> <given-names>S.</given-names></name> <etal/></person-group>. (<year>2006</year>). <article-title>Sequence-based species delimitation for the DNA taxonomy of undescribed insects</article-title>. <source>Syst. Biol</source>. <volume>55</volume>, <fpage>595</fpage>&#x02013;<lpage>609</lpage>. <pub-id pub-id-type="doi">10.1080/10635150600852011</pub-id><pub-id pub-id-type="pmid">16967577</pub-id></citation>
</ref>
<ref id="B61">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Poulin</surname> <given-names>E.</given-names></name> <name><surname>Gonz&#x000E1;lez-Wevar</surname> <given-names>C.</given-names></name> <name><surname>D&#x000ED;az</surname> <given-names>A.</given-names></name> <name><surname>G&#x000E9;rard</surname> <given-names>K.</given-names></name> <name><surname>H&#x000FC;ne</surname> <given-names>M.</given-names></name></person-group> (<year>2014</year>). <article-title>Divergence between Antarctic and South American marine invertebrates: what molecular biology tells us about Scotia Arc geodynamics and the intensification of the Antarctic Circumpolar Current</article-title>. <source>Glob. Planet. Change</source> <volume>123</volume>, <fpage>392</fpage>&#x02013;<lpage>399</lpage>. <pub-id pub-id-type="doi">10.1016/j.gloplacha.2014.07.017</pub-id></citation>
</ref>
<ref id="B62">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Puillandre</surname> <given-names>N.</given-names></name> <name><surname>Lambert</surname> <given-names>A.</given-names></name> <name><surname>Brouillet</surname> <given-names>S.</given-names></name> <name><surname>Achaz</surname> <given-names>G.</given-names></name></person-group> (<year>2012</year>). <article-title>ABGD, Automatic Barcode Gap Discovery for primary species delimitation</article-title>. <source>Mol. Ecol</source>. <volume>21</volume>, <fpage>1864</fpage>&#x02013;<lpage>1877</lpage>. <pub-id pub-id-type="doi">10.1111/j.1365-294X.2011.05239.x</pub-id><pub-id pub-id-type="pmid">21883587</pub-id></citation>
</ref>
<ref id="B63">
<citation citation-type="web"><person-group person-group-type="author"><name><surname>Rambaut</surname> <given-names>A.</given-names></name> <name><surname>Suchard</surname> <given-names>M. A.</given-names></name> <name><surname>Xie</surname> <given-names>D.</given-names></name> <name><surname>Drummond</surname> <given-names>A.</given-names></name></person-group> (<year>2014</year>). <source>Tracer v1.6</source>. Available online at: <ext-link ext-link-type="uri" xlink:href="http://beast.bio.ed.ac.uk/Tracer">http://beast.bio.ed.ac.uk/Tracer</ext-link></citation>
</ref>
<ref id="B64">
<citation citation-type="web"><person-group person-group-type="author"><name><surname>Ratnasingham</surname> <given-names>S.</given-names></name> <name><surname>Hebert</surname> <given-names>P. D. N.</given-names></name></person-group> (<year>2007</year>). <article-title>Bold: the barcode of life data system (<ext-link ext-link-type="uri" xlink:href="http://www.barcodinglife.org">http://www.barcodinglife.org</ext-link>)</article-title>. <source>Mol. Ecol. Notes</source> <volume>7</volume>, <fpage>355</fpage>&#x02013;<lpage>364</lpage>. <pub-id pub-id-type="doi">10.1111/j.1471-8286.2007.01678.x</pub-id></citation>
</ref>
<ref id="B65">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Rogers</surname> <given-names>A. R.</given-names></name> <name><surname>Harpending</surname> <given-names>H.</given-names></name></person-group> (<year>1992</year>). <article-title>Population growth makes waves in the distribution of pairwise genetic differences</article-title>. <source>Mol. Biol. Evol</source>. <volume>9</volume>, <fpage>552</fpage>&#x02013;<lpage>569</lpage>. <pub-id pub-id-type="pmid">1316531</pub-id></citation>
</ref>
<ref id="B66">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Ronquist</surname> <given-names>F.</given-names></name> <name><surname>Huelsenbeck</surname> <given-names>J. P.</given-names></name></person-group> (<year>2003</year>). <article-title>MrBayes 3: Bayesian phylogenetic inference under mixed models</article-title>. <source>Bioinformatics</source> <volume>19</volume>, <fpage>1572</fpage>&#x02013;<lpage>1574</lpage>. <pub-id pub-id-type="doi">10.1093/bioinformatics/btg180</pub-id><pub-id pub-id-type="pmid">12912839</pub-id></citation>
</ref>
<ref id="B67">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Sands</surname> <given-names>C. J.</given-names></name> <name><surname>Griffiths</surname> <given-names>H. J.</given-names></name> <name><surname>Downey</surname> <given-names>R. V.</given-names></name> <name><surname>Barnes</surname> <given-names>D. K. A.</given-names></name> <name><surname>Linse</surname> <given-names>K.</given-names></name> <name><surname>Mart&#x000ED;n-Ledo</surname> <given-names>R.</given-names></name></person-group> (<year>2013</year>). <article-title>Observations of the ophiuroids from the West Antarctic sector of the Southern Ocean</article-title>. <source>Antarct. Sci</source>. <volume>25</volume>, <fpage>3</fpage>&#x02013;<lpage>10</lpage>. <pub-id pub-id-type="doi">10.1017/S0954102012000612</pub-id></citation>
</ref>
<ref id="B68">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Sanford</surname> <given-names>E.</given-names></name> <name><surname>Kelly</surname> <given-names>M. W.</given-names></name></person-group> (<year>2011</year>). <article-title>Local adaptation in marine invertebrates</article-title>. <source>Annu. Rev. Mar. Sci</source>. <volume>3</volume>, <fpage>509</fpage>&#x02013;<lpage>535</lpage>. <pub-id pub-id-type="doi">10.1146/annurev-marine-120709-142756</pub-id><pub-id pub-id-type="pmid">21329215</pub-id></citation>
</ref>
<ref id="B69">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Stamatakis</surname> <given-names>A.</given-names></name></person-group> (<year>2014</year>). <article-title>RAxML version 8: a tool for phylogenetic analysis and post-analysis of large phylogenies</article-title>. <source>Bioinformatics</source> <volume>30</volume>, <fpage>1312</fpage>&#x02013;<lpage>1313</lpage>. <pub-id pub-id-type="doi">10.1093/bioinformatics/btu033</pub-id><pub-id pub-id-type="pmid">24451623</pub-id></citation>
</ref>
<ref id="B70">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Strugnell</surname> <given-names>J. M.</given-names></name> <name><surname>Rogers</surname> <given-names>A. D.</given-names></name> <name><surname>Prodohl</surname> <given-names>P. A.</given-names></name> <name><surname>Collins</surname> <given-names>M. A.</given-names></name> <name><surname>Allcock</surname> <given-names>A. L.</given-names></name></person-group> (<year>2008</year>). <article-title>The thermohaline expressway: the Southern Ocean as a centre of origin for deep-sea octopuses</article-title>. <source>Cladistics</source> <volume>24</volume>, <fpage>853</fpage>&#x02013;<lpage>860</lpage>. <pub-id pub-id-type="doi">10.1111/j.1096-0031.2008.00234.x</pub-id></citation>
</ref>
<ref id="B71">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Strugnell</surname> <given-names>J. M.</given-names></name> <name><surname>Watts</surname> <given-names>P. C.</given-names></name> <name><surname>Smith</surname> <given-names>P. J.</given-names></name> <name><surname>Allcock</surname> <given-names>A. L.</given-names></name></person-group> (<year>2012</year>). <article-title>Persistent genetic signatures of historic climatic events in an Antarctic octopus</article-title>. <source>Mol. Ecol</source>. <volume>21</volume>, <fpage>2775</fpage>&#x02013;<lpage>2787</lpage>. <pub-id pub-id-type="doi">10.1111/j.1365-294X.2012.05572.x</pub-id><pub-id pub-id-type="pmid">22494503</pub-id></citation>
</ref>
<ref id="B72">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Tajima</surname> <given-names>F.</given-names></name></person-group> (<year>1989</year>). <article-title>Statistical method for testing the neutral mutation hypothesis by DNA polymorphism</article-title>. <source>Genetics</source> <volume>123</volume>, <fpage>585</fpage>&#x02013;<lpage>595</lpage>. <pub-id pub-id-type="pmid">2513255</pub-id></citation>
</ref>
<ref id="B73">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Tamura</surname> <given-names>K.</given-names></name> <name><surname>Peterson</surname> <given-names>D.</given-names></name> <name><surname>Peterson</surname> <given-names>N.</given-names></name> <name><surname>Stecher</surname> <given-names>G.</given-names></name> <name><surname>Nei</surname> <given-names>M.</given-names></name> <name><surname>Kumar</surname> <given-names>S.</given-names></name></person-group> (<year>2011</year>). <article-title>MEGA5: molecular evolutionary genetics analysis using maximum likelihood, evolutionary distance, and maximum parsimony methods</article-title>. <source>Mol. Biol. Evol</source>. <volume>28</volume>, <fpage>2731</fpage>&#x02013;<lpage>2739</lpage>. <pub-id pub-id-type="doi">10.1093/molbev/msr121</pub-id><pub-id pub-id-type="pmid">21546353</pub-id></citation>
</ref>
<ref id="B74">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Taylor</surname> <given-names>H. R.</given-names></name> <name><surname>Harris</surname> <given-names>W. E.</given-names></name></person-group> (<year>2012</year>). <article-title>An emergent science on the brink of irrelevance: a review of the past 8 years of DNA barcoding</article-title>. <source>Mol. Ecol. Resour</source>. <volume>12</volume>, <fpage>377</fpage>&#x02013;<lpage>388</lpage>. <pub-id pub-id-type="doi">10.1111/j.1755-0998.2012.03119.x</pub-id><pub-id pub-id-type="pmid">22356472</pub-id></citation>
</ref>
<ref id="B75">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Templeton</surname> <given-names>A. R.</given-names></name> <name><surname>Crandall</surname> <given-names>K. A.</given-names></name> <name><surname>Sing</surname> <given-names>C. F.</given-names></name></person-group> (<year>1992</year>). <article-title>A cladistic analysis of phenotypic associations with haplotypes inferred from restriction endonucleaase mapping and DNA sequence data. III. Cladogram estimation</article-title>. <source>Genetics</source> <volume>132</volume>, <fpage>619</fpage>&#x02013;<lpage>633</lpage>. <pub-id pub-id-type="pmid">1385266</pub-id></citation>
</ref>
<ref id="B76">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Thomas</surname> <given-names>C. D.</given-names></name></person-group> (<year>2010</year>). <article-title>Climate, climate change and range boundaries</article-title>. <source>Divers. Distrib</source>. <volume>16</volume>, <fpage>488</fpage>&#x02013;<lpage>495</lpage>. <pub-id pub-id-type="doi">10.1111/j.1472-4642.2010.00642.x</pub-id></citation>
</ref>
<ref id="B77">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Turner</surname> <given-names>J.</given-names></name> <name><surname>Barrand</surname> <given-names>N. E.</given-names></name> <name><surname>Bracegirdle</surname> <given-names>T. J.</given-names></name> <name><surname>Convey</surname> <given-names>P.</given-names></name> <name><surname>Hodgson</surname> <given-names>D. A.</given-names></name> <name><surname>Jarvis</surname> <given-names>M.</given-names></name> <etal/></person-group>. (<year>2014</year>). <article-title>Antarctic climate change and the environment: an update</article-title>. <source>Polar Rec</source>. <volume>50</volume>, <fpage>237</fpage>&#x02013;<lpage>259</lpage>. <pub-id pub-id-type="doi">10.1017/S0032247413000296</pub-id></citation>
</ref>
<ref id="B78">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Ward</surname> <given-names>R. D.</given-names></name> <name><surname>Holms</surname> <given-names>B. H.</given-names></name> <name><surname>O&#x00027;Hara</surname> <given-names>T. D.</given-names></name></person-group> (<year>2008</year>). <article-title>DNA barcoding discriminates echinoderm species</article-title>. <source>Mol. Ecol. Resour</source>. <volume>8</volume>, <fpage>1202</fpage>&#x02013;<lpage>1211</lpage>. <pub-id pub-id-type="doi">10.1111/j.1755-0998.2008.02332.x</pub-id><pub-id pub-id-type="pmid">21586007</pub-id></citation>
</ref>
<ref id="B79">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Whitehouse</surname> <given-names>M. J.</given-names></name> <name><surname>Meredith</surname> <given-names>M. P.</given-names></name> <name><surname>Rothery</surname> <given-names>P.</given-names></name> <name><surname>Atkinson</surname> <given-names>A.</given-names></name> <name><surname>Ward</surname> <given-names>P.</given-names></name> <name><surname>Korb</surname> <given-names>R. E.</given-names></name></person-group> (<year>2008</year>). <article-title>Rapid warming of the ocean around South Georgia, Southern Ocean, during the 20th century: forcings, characteristics and implications for lower trophic levels</article-title>. <source>Deep Sea Res. Part Oceanogr. Res. Pap</source>. <volume>55</volume>, <fpage>1218</fpage>&#x02013;<lpage>1228</lpage>. <pub-id pub-id-type="doi">10.1016/j.dsr.2008.06.002</pub-id></citation>
</ref>
<ref id="B80">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Whitworth</surname> <given-names>T.</given-names> <suffix>III.</suffix></name></person-group> (<year>1980</year>). <article-title>Zonation and geostrophic flow of the Antarctic circumpolar current at Drake Passage</article-title>. <source>Deep Sea Res. Part Oceanogr. Res. Pap</source>. <volume>27</volume>, <fpage>497</fpage>&#x02013;<lpage>507</lpage>. <pub-id pub-id-type="doi">10.1016/0198-0149(80)90036-9</pub-id></citation>
</ref>
<ref id="B81">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Zhang</surname> <given-names>J.</given-names></name> <name><surname>Kapli</surname> <given-names>P.</given-names></name> <name><surname>Pavlidis</surname> <given-names>P.</given-names></name> <name><surname>Stamatakis</surname> <given-names>A.</given-names></name></person-group> (<year>2013</year>). <article-title>A general species delimitation method with applications to phylogenetic placements</article-title>. <source>Bioinformatics</source> <volume>29</volume>, <fpage>2869</fpage>&#x02013;<lpage>2876</lpage>. <pub-id pub-id-type="doi">10.1093/bioinformatics/btt499</pub-id><pub-id pub-id-type="pmid">23990417</pub-id></citation>
</ref>
</ref-list>
</back>
</article>
