Register      Login
Invertebrate Systematics Invertebrate Systematics Society
Systematics, phylogeny and biogeography
RESEARCH ARTICLE (Open Access)

One Antarctic slug to confuse them all: the underestimated diversity of Doris kerguelenensis

Paige J. Maroni https://orcid.org/0000-0002-1974-3977 A B * , Bill J. Baker C , Amy L. Moran D , H. Arthur Woods E , Conxita Avila F , Glenn J. Johnstone G , Jonathan S. Stark G , Kevin M. Kocot H , Susanne Lockhart I , Thomas Saucède J , Greg W. Rouse https://orcid.org/0000-0001-9036-9263 K and Nerida G. Wilson https://orcid.org/0000-0002-0784-0200 A B K L
+ Author Affiliations
- Author Affiliations

A The University of Western Australia, School of Biological Sciences (M092), 35 Stirling Highway, Crawley, WA 6009, Australia.

B Western Australian Museum, Research & Collections, 49 Kew Street, Welshpool, WA 6106, Australia.

C University of South Florida, Department of Chemistry (CHE 205), 4202 East Fowler Street, Tampa, FL 33620, USA.

D University of Hawai’i at Mānoa, School of Life Sciences (EDM 216), 2538 McCarthy Mall, Honolulu, HI 96822, USA.

E University of Montana, Division of Biological Sciences (HS 104), 32 Campus Drive, Missoula, MN 59812, USA.

F University of Barcelona, Department of Evolutionary Biology, Ecology and Environmental Sciences & Biodiversity Research Institute (IRBio), Diagonal Avenida 643, E-08028 Barcelona, Catalonia, Spain.

G Australian Antarctic Division, Environmental Protection Program, 203 Channel Highway, Kingston, Tas. 7050, Australia.

H University of Alabama, Department of Biological Sciences and Alabama Museum of Natural History, 500 Hackberry Lane, Tuscaloosa, AL 35487, USA.

I California Academy of Sciences, Department of Invertebrate Zoology and Geology, 55 Music Concourse Drive, San Francisco, CA 94118, USA.

J Université Bourgogne Franche-Comté, Biogéosciences, 6 Gabriel Boulevard, UMR 6282 CNRS, F-21000 Dijon, France.

K University of California San Diego, Scripps Institution of Oceanography, 9500 Gilman Drive, La Jolla, CA 92093, USA.

L Securing Antarctica’s Environmental Future, Western Australian Museum, 49 Kew Street, Welshpool, WA 6106, Australia.


Handling Editor: Ana Riesgo

Invertebrate Systematics 36(5) 419-435 https://doi.org/10.1071/IS21073
Submitted: 22 October 2021  Accepted: 24 March 2022   Published: 1 June 2022

© 2022 The Author(s) (or their employer(s)). Published by CSIRO Publishing. This is an open access article distributed under the Creative Commons Attribution-NonCommercial-NoDerivatives 4.0 International License (CC BY-NC-ND)

Abstract

The Antarctic marine environment, although rich in life, is predicted to experience rapid and significant effects from climate change. Despite a revolution in the approaches used to document biodiversity, less than one percent of Antarctic marine invertebrates are represented by DNA barcodes and we are at risk of losing biodiversity before discovery. The ease of sequencing mitochondrial DNA barcodes has promoted this relatively ‘universal’ species identification system across most metazoan phyla and barcode datasets are currently readily used for exploring questions of species-level taxonomy. Here we present the most well-sampled phylogeny of the direct-developing, Southern Ocean nudibranch mollusc, Doris kerguelenensis to date. This study sampled over 1000 new Doris kerguelenensis specimens spanning the Southern Ocean and sequenced the mitochondrial COI gene. Results of a maximum likelihood phylogeny and multiple subsequent species delimitation analyses identified 27 new species in this complex (now 59 in total). Using rarefaction techniques, we infer more species are yet to be discovered. Some species were only collected from southern South America or the sub-Antarctic islands, while at least four species were found spanning the Polar Front. This is contrary to dispersal predictions for species without a larval stage such as Doris kerguelenensis. Our work demonstrates the value of increasing geographic scope in sampling and highlights what could be lost given the current global biodiversity crisis.

Keywords: allopatry, Antarctica, Antarctic marine biodiversity, cryptic species, cytochrome oxidase I, direct development, mtDNA, nudibranch mollusc, phylogeny, refugia, species delimitation.

Introduction

Biodiversity and natural ecosystems provide humankind with significant economic benefits (e.g. Wallace 1997; McClintock and Baker 2001; Pertierra et al. 2021) by way of indirect essential services (e.g. maintaining water cycles – McNeely et al. 1990; CO2 emission mitigation – Domke et al. 2020). Despite all known benefits, the rapid destruction of the world’s most diverse ecosystems has led most experts to conclude that the earth’s biological diversity is in danger (Singh 2002). Although the total number of extant species is currently estimated to be between 5.3 million and 1 trillion, only 1.85 million species have been formally described (Mora et al. 2011; Locey and Lennon 2016). We also know that our oceans house a large array of marine species (146 969 accepted species; Ocean Biodiversity Information System, see https://obis.org/), with another 1.4–1.6 million species hypothesised to exist, and are awaiting discovery and description (Bouchet 2006). Concurrent with the evolution of human societies (c. 11 000 years ago) in the form of advanced infrastructure, farming and transport, there has been continual biodiversity loss, both of species, and wider ecosystem integrity and functionality (e.g. Rogers-Bennett and Catton 2019; Prates and Perez 2021). This loss is underestimated due to taxonomic uncertainty and data gaps. Throughout geological time, at least five putative mass extinctions have occurred (Barnosky et al. 2011). Unlike past natural events, however, that spanned millions of years, the present anthropogenically driven mass extinction is likely to occur over only hundreds of years (as few as 200 years) (see: Ceballos et al. 2015). This loss of biodiversity has significant effects on ecosystem function that has implications for human populations and these may still be vastly underestimated.

Despite being the most isolated continent on Earth, Antarctica has not escaped the negative effects of human activity (Vaughan et al. 2003; Aronson et al. 2011; Chown et al. 2015; Stephens 2018). Such impacts include pollution, overfishing, increased melting of ice and the introduction of invasive or alien species (Tin et al. 2009; Aronson et al. 2011; Stark et al. 2019; Avila et al. 2020; De Castro-Fernández et al. 2021). Globally, the Antarctic region is known to act as a significant carbon sink and surrounding ecosystems are experiencing temperature increases due to rising global atmospheric carbon dioxide concentrations (Ito et al. 2010). Parts of the Antarctic Peninsula, including the west Antarctic Peninsula for example, are particularly vulnerable to these human-induced ecological disasters, as these are already experiencing the greatest increases in mean annual atmospheric temperatures on Earth (Chapman and Walsh 2007; Clarke et al. 2007; Ingels et al. 2012; Torre et al. 2017). Failure to address these global threats will not only result in the degradation of Antarctic marine ecosystems but will potentially lead to enormous losses of global biodiversity.

Technical advances in molecular phylogenetics over the past several decades have resulted in the development of many more sensitive rapid tools for detecting new species (Féral 2002; Goetze 2003; Baird et al. 2011). The ease of sequencing mitochondrial DNA barcodes has promoted a relatively ‘universal’ species identification system across most metazoan phyla and now barcode datasets are more readily used for exploring questions of species-level taxonomy (Hebert et al. 2003; Puillandre et al. 2012; Taberlet et al. 2012; Eberle et al. 2020). In current phylogenetic studies the types of molecular information extracted for rapid assessments of animal diversity have evolved from interpreting single-locus (e.g. Hebert et al. 2003) to whole-genome datasets (e.g. Jensen et al. 2021). This expansion of available sequences, facilitated by the automation of sequencing and decrease in sequencing costs, has accelerated phylogenetic studies and whole-genome-based research (e.g. Johnson et al. 2008; Layton et al. 2018; Cai et al. 2019). This is particularly significant when barcoding cryptic and pseudocryptic species (organisms that are first discerned by non-morphological methods, e.g. Brasier et al. 2016; Matsuda and Gosliner 2018; Tyagi et al. 2019). When divergent interspecific traits are not morphologically obvious, traditional taxonomic methods fail to detect speciation events and biodiversity will remain under-reported (Knowlton 1993; Baird et al. 2011).

The Antarctic continental shelf is one marine realm that has recently revealed apparent high levels of cryptic species (e.g. Linse et al. 2007; Wilson et al. 2009; Baird et al. 2011; Brasier et al. 2016), specifically in organisms with poor dispersal abilities (Griffiths 2010; De Broyer and Danis 2011; Grant et al. 2011; Neusser et al. 2011; Brandt et al. 2012). The growth and decay of ice sheets during Antarctica’s history has been one of the most significant disturbances acting at an ecosystem level across the shelf and underpins one of the fundamental frameworks for understanding Antarctic diversity, the ‘Antarctic Biodiversity Pump’ hypothesis (Clarke and Crame 1989, 2010; Gutt and Starmans 2002; Thatje et al. 2005). This hypothesis proposes that overall, increased speciation rates are the result of Milankovitch-cycle driven glacial oscillations that drive cryptic and allopatric speciation by gene flow inhibition and speciation events (Clarke and Crame 1989; Crame 1997; Griffiths 2010). A substantial amount of Antarctic species diversity is hypothesised to be the result of species flocks that have been identified as monophyletic taxa displaying high levels of endemic species that are ecologically diverse and abundant in relation to the surrounding habitat (Ribbink 1984; Eastman and McCune 2000; Lecointre et al. 2013; Chenuil et al. 2018). During these glacial oscillations and ensuing periodic habitat disruptions, ice-free refugia are thought to have sustained reduced populations of once widely distributed species (e.g. Smith et al. 2010; Lau et al. 2020), allowing many to diverge during this time.

One Southern Ocean marine invertebrate that is emerging as a case study in cryptic speciation is Doris kerguelenensis (Bergh, 1884). Recent studies on this nominal species have highlighted either new genetic lineages (Wilson et al. 2009, 2013) or new chemical compounds (Iken et al. 2002; Maschek et al. 2012; Avila 2020). Doris kerguelenensis is a direct-developing sea slug with limited dispersal potential and long generation times (Hain and Arnaud 1992; Moles et al. 2017a). Belonging to the Dorididea, this nudibranch is a simultaneous hermaphrodite that feeds exclusively on sponges and synthesises secondary metabolites de novo (Maschek et al. 2012). Wägele (1990) reviewed and ultimately synonymised ten Southern Ocean dorid nudibranch species within the single, morphologically variable species Austrodoris kerguelenensis (Bergh, 1884). That work also designated another two species as nomina dubia due to misplaced holotypes and inadequate descriptions. The genus Austrodoris Odhner, 1926 was subsequently revised and synonymised along with five other cryptobranch dorid nudibranch genera into Doris Linnaeus, 1758 (Valdes 2001). Wilson et al. (2009) examined the mitochondrial protein-coding gene Cytochrome Oxidase I (COI) and revealed 29 putative lineages within this nominal species. These lineages were subsequently corroborated with nuclear and metabolomic trait data to infer biological species (Wilson et al. 2013); therefore, this sets a sound basis for using COI clades as a proxy for species in this study. Interestingly, three new species were recovered by resampling in the same geographic regions (Wilson et al. 2013), hinting at further undetected diversity. In this study, we explored the identity of over 1000 new Doris kerguelenensis specimens using mitochondrial DNA sequence data, with the aim of (i) creating an expanded phylogenetic hypothesis for the group, (ii) testing for additional cryptic species within this species complex and (iii) exploring the distributional patterns of these species.


Methods

Specimen collection and preservation

In this study, 1275 individuals of Doris kerguelenensis were included from 146 sites from depths between the intertidal and 798 m (Fig. 1, Supplementary Table S1). Our specimens were collected during various Antarctic field expeditions using a Blake trawl, Smith-McIntyre grab, Agassiz trawl, wire dredge, epibenthic sled or hand collected by SCUBA diving. Samples were collected from various locations in the Southern Ocean between 2006 and 2018 (Fig. 1). Nineteen geographical regions were defined a priori (Table 1). Owing to differing depths, distances between regions, coastal currents and ocean circulation patterns (Smith et al. 1999) we separated the Antarctic Peninsula region into four regions: (i) Palmer Archipelago, (ii) Bransfield Strait, (iii) South Shetland Islands and (iv) Elephant Island. Tissue subsamples were taken from specimens preserved in 96–100% ethanol or frozen. Sequenced specimens were housed in the Western Australian Museum (WAM), Scripps Institution of Oceanography, Benthic Invertebrates Collection (SIO-BIC), Yale Peabody Museum (YPM), California Academy of Sciences (CAS), National Museum of Natural History, Smithsonian Institution (USNM), University of Barcelona and the Baker Laboratory, University of South Florida. Data for all specimens sequenced in this study are available in Supplementary Table S1 including GenBank accession numbers for outgroup taxa.


Fig. 1.  Map of Antarctica showing the sites of all Doris kerguelenensis samples sequenced for the mitochondrial DNA gene Cytochrome Oxidase I (COI). The Antarctic Polar Front is denoted by the solid white line with two adjustments proposed by Park et al. (2014) and Sokolov and Rintoul (2009) (hashed white lines depicting the APF moving south of Kerguelen Island). Colours indicate geographic sampling regions. Base map generated through Quantarctica (ver. 3.2, see https://www.npolar.no/quantarctica/; Matsuoka et al. 2021).
Click to zoom


Table 1.  Table of COI diversity of Doris kerguelenensis individuals and species by region.
Click to zoom

In addition to the ingroup specimens, a range of outgroup sequences were selected from GenBank (Wollscheid-Lengeling et al. 2001; Shields 2009; Pola and Gosliner 2010; Jung et al. 2014; Palomar et al. 2014; Hulett et al. 2015; Mahguib and Valdés 2015; Goodheart et al. 2018) (see Supplementary Table S1) along with other available sequences for Doris kerguelenensis (Wilson et al. 2009, 2013). This was done as a robust sister group for Dorididae, across the Nudibranchia phylogeny is not yet certain (e.g. Wollscheid-Lengeling et al. 2001; Shields 2009; Mahguib and Valdés 2015; Korshunova et al. 2020); however, the tree was ultimately rooted with Prodoris clavigera.

DNA extraction, amplification and DNA sequencing

Total genomic DNA was extracted from 1077 ethanol-fixed or frozen samples using a DNeasy Blood & Tissue Kit (Qiagen) according to the manufacturer’s instructions (excluding the optional repeat of a final elution step). For analysis, an additional 198 previously published sequences were also incorporated (see Wilson et al. 2009, 2013). The same primers were used for PCR and sequencing (LCO1490/HCO2198, Folmer et al. 1994; or modified degenerate versions of the former, jgLCO1490/jgHCO2198, Geller et al. 2013). These primers amplified a fragment of the cytochrome oxidase I gene (COI). Each PCR reaction included: 16.8 μL of molecular grade (deionised) water, 5.0 μL 5 × MyTaq PCR buffer (1 × = 1 mM of dNTPs, 3 mM of MgC12) (Bioline), 0.8 μL of forward primer, 0.8 μL of reverse primer, 0.2 μL of Platinum Taq polymerase (Invitrogen), with 1.5 μL of DNA template added. The thermocycling protocol for COI using primers LCO1490/HCO2198 followed the procedure: 5 min. at 95°C, 7 cycles 30 s at 95°C, 30 s at 45°C, 1 min. at 72°C followed by 35 cycles of 30 s at 95°C, 30 s at 50°C, 1 min. at 72°C with a final extension time of 10 min. with a temperature of 72°C. The thermocycling protocol using primers jgLCO1490/jgHCO2198 was as follows: 3 min. at 95°C, 8 cycles of 30 s at 95°C, 30 s at 50°C, 45 s at 72°C followed by 32 cycles of 30 s at 95°C, 30 s at 48°C and 45 s at 72°C with a final extension time of 5 min at 72°C. Amplicons were screened on E-gels (Invitrogen) and the positive reactions, determined through Bio-Rad Image Laboratory software, were outsourced to the Australian Genome Research Facility (Perth) for enzyme purification and Sanger sequencing using an Applied Biosystems 3730 capillary sequencer. All sequences were assembled and edited in Geneious Prime 2020.1 (Kearse et al. 2012). In cases where amplicon bands were deemed weak or if amplicons failed to sequence, experiments were run to improve sequencing success either by increasing DNA or DNA dilutions (1:10, 1:20, 1:50 and 1:100).

Phylogenetic reconstruction and primary species hypotheses (PSH)

COI sequences from Doris kerguelenensis and outgroups were aligned with the MAFFT plugin for Geneious (ver. 1.4.0, see https://www.geneious.com/plugins/mafft-plugin/; Katoh and Standley 2013) using default settings. The sequences were used to create a single gene, maximum likelihood (ML) phylogeny using IQ-TREE (ver. 1.6.12, see http://www.iqtree.org/; Nguyen et al. 2015; Trifinopoulos et al. 2016), herein referred to as the Primary Species Hypotheses (PSH) (Fig. 2, Supplementary Fig. S1) that will be further tested with delimitation methods. The –m TEST (Kalyaanamoorthy et al. 2017) option in IQ-TREE identified TN + F + I + G4 as the best-fit model chosen according to Bayesian Information Criterion (BIC). Nodal support was assessed with 1000 ultrafast bootstrap replicates (Hoang et al. 2018). This phylogenetic tree was used to make semi-arbitrary decisions about inferred species clades, much in the same way that taxonomists intuitvely assign names based on morphological variation (Dayrat 2005). Some of the ‘arbitrary’ species-level clades had already been corroborated with nuclear and trait data (see Wilson et al. 2013).


Fig. 2.  Maximum likelihood (ML) phylogeny of the unique COI haplotypes dataset, including 320 Doris sequences and 21 outgroups (not shown) using a TN + F + I + G4 best-fit model. Nodes with ultrafast bootstrap support values of 95 or higher have been denoted by a circular node shape. Triangles represent collapsed clades. Boxes represent the Primary Species Hypothesis (PSH) followed by the partition results for species delimitation analyses (ABGD P2, ABGD P6, ASAP 1st [all input files except ASAP K80 uncorrected], ASAP 1st K80 uncorrected, TCS and mPTP). Grey boxes indicate species with samples collected from both sides of the Antarctic Polar Front (southern South America and Antarctica).
F2

In previous publications assessing diversity in D. kerguelenensis, species-level clades were numbered (Wilson et al. 2009, 2013) from 1 to 32. These original numbers have been incorporated into this new work and we have continued numbering new clades from 33 onwards, with one exception. The original clade 12, originally consisting of a single specimen (USNM1120712) has now been synonymised within clade 11. The newly named clade 12 in this study comprises a previously unknown clade.

Pairwise distances (uncorrected p-distances) or Tamura and Nei (1993) (TN93) corrected genetic distances were both calculated from the mtDNA sequences using the ‘dist.dna’ function from the R.cran (ver. 4.2.0, R Foundation for Statistical Computing, Vienna, Austria, see https://cran.r-project.org/) package Ape (ver. 5.6, see https://cran.r-project.org/web/packages/ape/index.html; Paradis and Schliep 2019) with either ‘raw’ or ‘TN93’ selected as the evolutionary model (Supplementary Fig. S2). The gamma shape α and base frequencies were also parameters specified for the distance calculations. Intra- and interspecific (PSH clades) p-distances and TN93 corrected genetic distances were also calculated in MEGA X (ver. 10.2.6, see https://www.megasoftware.net/; Kumar et al. 2018) (Supplementary Table S2).

Species delimitation

The process of delimitation was carried out in a series of steps. We tested our PSH through a range of species delimitation methods that included the Multi-rate Poisson Tree Process analysis ((m)PTP, Kapli et al. 2017), statistical parsimony networks (TCS Java program, ver. 1.21, see https://bio.tools/tcs; Clement et al. 2000; Hart and Sunday 2007), the Automatic Barcode Gap Definition (ABGD, Puillandre et al. 2012) and Assemble Species by Automatic Partitioning (ASAP, Puillandre et al. 2021). All delimitation analyses were run without outgroups included in the input files.

Multi-rate PTP is designed to better accommodate sampling and population specific characteristics of a broad range of datasets as this incorporates different levels of intraspecific genetic diversity (based on specific branching events or species-specific sampling) (Kapli et al. 2017). This analysis has been shown to outperform PTP (Zhang et al. 2013) in yielding more accurate delimitations with respect to taxonomy (i.e. identifies more taxonomically accepted species).

TCS (Clement et al. 2000) bins sequences into haplotypes and calculates the frequencies of the haplotypes in the sample (Clement et al. 2000). Relationships among haplotype networks were explored using this method using ‘parsimnet’, a function in the R.cran package Haplotypes (ver. 1.1.2, C. Aktas, see https://cran.r-project.org/web/packages/haplotypes/index.html) that finds the most parsimonious networks and is an implementation of the TCS methods proposed by Templeton et al. (1992). TCS (Clement et al. 2000) has been shown to be useful for species delimitation as the parsimony connection limit (95%) successfully delimits the same number of subnetworks as taxa (see Hart and Sunday 2007).

ABGD and ASAP are both ascending hierarchal clustering programs that merge sequences into groups defined as partitions. The partitions are defined by probability and barcode gap width however, unlike ABGD that is solely based on pairwise distances (Puillandre et al. 2012), there is no need for an a priori defined P (when P is the prior maximum divergence of intraspecific diversity) within ASAP’s input functions. ASAP produces a scoring system that ranks the partitions based on a combination of both probability and gap width (‘asap-score’) (Puillandre et al. 2021). In ABGD and ASAP analyses, there are two available substitution models, Jukes and Cantor (1969) (JC69) and Kimura (1980) (K80). However, the evolutionary model selected by IQ-TREE (a variant of Tamura and Nei 1993) is not an available model option when partitioning sequences in these programs, so Tamura and Nei (1993) corrected genetic distances were imported into these delimitation analyses for comparison. We used ABGD alongside ASAP because although ASAP has been developed to replace ABGD, comparisons between the two programs will allow for direct comparision to results from previously published literature. The default parameters employed within ABGD were Pmin = 0.001, Pmax = 0.10, 10 steps, X = 1.5, Nb bins of 20. Only the initial partitions were examined from the ABGD algorithm outputs (following Puillandre et al. 2021). In ASAP, the default parameters were also used.

When assessing ABGD and ASAP partitions, we examined a range of partitioning definitions. This was done to encapsulate all inter- and intraspecific variability and accommodate for over-splitting or coalescing of recently radiated taxa, rather than only selecting a single partition to represent species boundaries. ABGD initially divides the data into groups based on a range of prior intraspecific divergences and a statistically inferred barcode gap, and subsequently recursively applies the same procedure to the initial groups (Puillandre et al. 2012). For this reason, the initial partitions generally represent species or operational taxonomic units (OTUs) whereas the recursive partitions generally reflect population level diversity. Here we only examined the initial partition outputs and only considered ABGD partition 2 (P = 0.001; to highlight the intraspecific genetic variability) and partition 6 (P = 0.01), as defining species based on less than 1% sequence divergence is not appropriate. Partition 1 was omitted from reports as the primary partition assumes that a single gap can be defined for the entire dataset however the gap distances are highly likely to differ between groups within a dataset (Puillandre et al. 2012). We considered these two partitions from both substitution model outputs (JC69 and K80) (Supplementary Fig. S3). ASAP-scores are the average of the P-value and relative barcode gap width of the partitions. In terms of the ASAP partition outputs, three input files were examined. The original, uncorrected alignment file was entered into the webserver (https://bioinfo.mnhn.fr/abi/public/asap/asapweb.html#) and both substitution models were applied to the file (=uncorrected p-distances). The Tamura and Nei (1993) corrected MEGA X CSV genetic distance file was imported and tested against both JC69 and K80. A MEGA X CSV TN93 corrected distance file with the gamma shape α (0.692) and base frequencies (A = 0.245, C = 0.174, G = 0.189 and t = 0.391) specified, was also tested with both substitution models. ASAP 1st and ASAP 2nd partitions were considered due to the low (and thus better supported) ASAP scores of between 5.00 and 8.00 for all input files and substitution models (see Supplementary Fig. S4).

Species diversity estimates

Considering the large percentage of species clades consisting of single specimens (Fig. 3), we assessed whether further diversity remained hidden due to undersampling, by generating species rarefaction–extrapolation (R–E) curves and the associated 95% confidence intervals, using iNEXT (ver. 2.0.2, see https://cran.r-project.org/web/packages/iNEXT/index.html; Hsieh et al. 2016). Through iNEXT (interpolation and extrapolation), an R.cran package (Hsieh et al. 2016), our diversity estimates used abundance-based data (Gotelli and Colwell 2011) that tally the abundance of each species (SSH results were utilised as input) across pre-determined geographic regions (Table 1). A sample-size R–E curve was produced to determine whether the rate of discovery of new species slowed down or reached saturation with an increase in sample size. To avoid discarding data, a sample-size R–E curve that rarefies to smaller sample sizes or extrapolates to larger sample sizes (i.e. plots with respect to sample size) was utilised (Colwell et al. 2012; Hsieh et al. 2016). We conducted both whole dataset and region-based R–E curves that are represented by the means of repeated resampling (1000 bootstrap replicates). Examining the slope of these curves allowed us to assess the impact of undersampling on estimating diversity.


Fig. 3.  Frequency plot depicting the number of samples present within species of Doris kerguelenensis. For example, 14 species are represented only by singletons.
F3


Results

Species delimitation

A total of 1275 COI sequences (433 haplotypes) was generated and the final alignment constituted 658 base pairs (bp) in length. The primary species hypotheses (PSH) generated with a maximum likelihood (ML) phylogeny recovered 59 species-level clades within the monophyletic Doris kerguelenensis species complex, 27 more than previously known. The ML phylogeny generated through IQ-TREE showed low bootstrap support at interior nodes but mostly high support for terminal clusters (ultrafast bootstrap support >95–100), herein referred to as species. The only two species-level clades that were not supported by high bootstrap values were species 15 and 59 (bootstrap support values 69 and 58 respectively) (Fig. 2, Supplementary Fig. S1). The proportion of invariable sites was 0.475, the gamma shape α distribution parameter was 0.692 and the base frequencies were A = 0.245, C = 0.174, G = 0.189, t = 0.391 (unequal base frequencies). Rates for the six substitution types estimated from the dataset were AC = 1.0000, AG = 10.6386, AT = 1.0000, CG = 1.0000, CT = 7.1609 and GT = 1.0000 (unequal transition rates, equal transversion rates and unequal purine and pyrimidine rates).

To assess the robustness of the PSH, results were compared to those from four delimitation methods (mPTP, TCS, ABGD and ASAP) (Fig. 2, 4, Supplementary Table S1, Supplementary Fig. S1–S4). The Multi-rate Poisson Tree Process (mPTP) yielded 62 taxonomic units and TCS recovered 56 groups. Species delimitation based on genetic distance using ABGD analysis detected between 49 and 52 groups respectively (initial partitions 6 (P = 0.01) and 2 (P = 0.001)) using JC69 and K80 models (Supplementary Fig. S3). All ABGD results were identical between both substitution models (JC69 and K80, TS/TV = 2.0) (see Supplementary Fig. S3). The ASAP species hypothesis partition, ranked by best ASAP score (ASAP 1st) delimited 59 (JC69) and 61 (K80) species, and ASAP 2nd delimited 57 (JC69) and 63 (K80) species. When the TN93 corrected genetic distance files were incorporated into analyses and a model subsequently applied, 59 species were delimited (for both JC69 or K80) in ASAP 1st, and 48 or 60 (JC69 and K80 respectively) in ASAP 2nd. When the additional gamma shape α and base frequency parameters were defined, ASAP 1st and 2nd (JC69) delimited 59 and 60 species respectively, and ASAP 1st and 2nd (K80) delimited 59 and 48 species (see all ASAP results: Supplementary Fig. S4).


Fig. 4.  Comparison of species delimitation results for the Secondary Species Hypotheses (SSH). Results are ordered as shown in Fig. 2 (phylogenetic tree). ASAP 1st (lowest score, stippled for five of the six ASAP 1st results) was chosen as the SSH for this dataset. The Primary Species Hypotheses (PSH) also represented 59 species.
Click to zoom

The species delimitation results for mPTP, TCS, ABGD (initial partitions 2 and 6) and ASAP (ASAP 1st and 2nd (all input file structures)); are plotted against the ML phylogeny (Fig. 2, 4, Supplementary Fig. S1). Overall, for the Secondary Species Hypotheses (SSH), we chose to accept 59 species (Fig. 4), as this number was recovered by the best scoring ASAP partition in five of six conditions. This number of partitions was always recovered by ASAP 1st when the input alignment was adjusted with corrected distances, and for gamma shape and base frequencies.

Intraspecific colour variation (and to a lesser extent, dorsal tubercle exaggeration) was assessed to determine whether colour could be utilised as a taxonomically distinguishing characteristic. Live specimen photos (Supplementary Fig. S5) showed substantial colour variation within an example species (24) and was therefore deemed an uninformative trait for taxonomic purposes.

Sampling effort and geographic distributions

Most clades contained between 1 and 40 samples (mean = 21) and overall, 14 species were represented by single specimens (Fig. 3). The slope of the R–E curve approached saturation (an asymptotic shape) at around 68 species (Fig. 5), indicating that our sampling coverage was relatively high (accurately sampled enough of the landscape to capture species diversity) and we have therefore approached a realistic estimate for the number of species in this complex. The R–E curves reflecting sampling effort for specific locations such as the Palmer Archipelago (PAL), Elephant Island (EL) (Fig. 6a), Shag Rocks (SR) and Herdman Bank (HB) (Fig. 6b) reached an asymptote. However, the opposite result was detected for regions such as Discovery Bank (DB), the Weddell Sea (WS), Davis Station (DS), Casey Station and others.


Fig. 5.  Estimation of Doris kerguelenensis species diversity based on COI abundance data. Rarefaction curve for observed samples and the 95% upper and lower confidence intervals based on the abundance-based rarefaction curve.
F5


Fig. 6.  Rarefaction curves were constructed through a single analysis that was subsequently split for visual ease. (a) Estimation of Doris kerguelenensis species diversity based on COI abundance data separated by geographic sampling region (Burdwood Bank, Bransfield Strait, Elephant Island, Palmer Archipelago, South Orkney Island and South Shetland Island) for regions that include 95 individual specimens or more. Rarefaction curve for observed samples and the 95% upper and lower confidence intervals are based on an abundance-based rarefaction curve. (b) Estimation of Doris kerguelenensis species diversity based on COI abundance data separated by geographic sampling region (Bouvet Island, Casey Station, Discovery Bank, Davis Station, Falkland Island, Herdman Bank, Kerguelen Plateau, Ross Sea, South Georgia, Mount Siple, Shag Rocks, South Sandwich Islands and Weddell Sea) for regions that include less than 95 individual specimens. Rarefaction curves for observed samples within locations and the 95% upper and lower confidence intervals based on the abundance-based rarefaction curve.
Click to zoom

One of the most well-sampled species, species 29 (n = 280), has a circum-Antarctic distribution (collected from over 11 000 km) and contains samples collected from Prydz Bay, the Ross Sea and the Antarctic Peninsula. An additional six species have very large distributions (defined here as >2000 km) and data show that at least four species distributions span the Polar Front (clades 14, 24, 26 and 42) (grey highlighted boxes within the PSH, Fig. 3). Thirteen clades, including the aforementioned four species (1, 2, 12, 28, 35, 44, 47, 48 and 50) have been collected above the Polar Front (Falkland Islands, Burdwood Bank, Bouvet Island and Kerguelen Plateau) and three of these species have an Antarctic continental shelf sister-species (1 with 10 + 55, 2 with 29 and 44 with 30, Supplementary Table S1). The deepest depth at which a specimen was collected was 798 m.

Many species occurred in sympatry. Twenty-one species were collected from the Bransfield Strait, nineteen on the South Orkney microcontinent and thirteen from the Burdwood Bank, Palmer Archipelago and South Shetland Islands (Table 1). Overall, species richness remains a direct reflection of sampling intensity within these results, e.g. seven ‘well sampled’ regions (each with over 40 samples) had more than 10 species within that corresponding location (Table 1).


Discussion

Doris kerguelenensis was long thought to consist of a single, widespread species (since Wägele 1990). With the application of molecular data, multiple species delimitation algorithms and increased sampling, the species complex is currently recognised to comprise at least 59 geographically overlapping species. Here, we show: (i) 27 previously undiscovered species, (ii) further evidence that more species remain to be discovered within previously sampled regions of the Southern Ocean and (iii) new insights into the distributions of these animals.

Even more new species delimited

The primary species hypotheses (PSH) outlined 59 species within the Doris kerguelenensis complex, of which 27 were new species-level clades. Subsequent species delimitation analyses showed highly congruent results (48 – 63 total range), with the best supported ASAP scores converging on 59 species. Our results corroborated the D. kerguelenensis clades that were previously detected through mitochondrial, nuclear and metabolomic datasets (Wilson et al. 2009, 2013), with one exception (the original clade 12, Wilson et al. 2009).

As expected, species delimitation methods based on different underlying assumptions rendered varying species hypotheses. Overall, ASAP 2nd (K80) showed the highest evidence for over-splitting (63 species), relative to the clades recognised in the PSH and all other partition analyses. We suspected that ASAP 2nd (K80) may have been over-splitting clades as mPTP also suggested there were 62 species and mPTP performed poorly when the number of species was small or consisted of a small sample size (Puillandre et al. 2021). Also, many of these new D. kerguelenensis species comprised one or only a few individuals and this may explain the similar result of 62 species delimited through this method. Our TCS analyses recovered 56 species and this was ultimately identical to the PSH with the single exception that clades 15 + 31 + 46 + 60 were merged into one. The ABGD initial partitions delimited between 49 and 52 species (JC69 and K80). This outcome was generally expected as the initial partitions of ABGD results will delimit a more conservative number of taxa, while recursive have been suggested to represent population level diversity rather than interspecific diversity (Puillandre et al. 2012). Overall, ASAP's 2nd partitions delimited between 48 and 63 species. Finally, ASAP 1st partition results (from all input file structures, except K80) all converged on 59 species and this was accepted as the most well-supported delimitation estimation.

Despite some colour variations between specimens collected (ranging from white, yellow, orange, pink, including white with pink colouration around the mouth; N. G. Wilson, pers. obs. and Supplementary Fig. S5), most were morphologically similar. This means that gill or body colour, or tubercle structure could not be used to visually differentiate species. Additionally, when preserved, the colour of specimens disappears and the animals become a uniform off-white colour (first noted in Wägele 1990). Specimens also ranged considerably in length, from ~1 to 20 cm. Morphological studies mention that external traits such as the size and density of the dorsal papillae should only be tentatively relied on for species determination, because these change in appearance when preserved (Wägele 1990; Cattaneo-Vietti 1991). Wägele (1990) highlighted considerable morphological variation in the nervous system, the digestive system, and the colour and shape of the radula. None of these variations were consistent within those species concepts, therefore all animals examined were considered to be one species. Morphology was, therefore, not further considered in species delimitation efforts throughout this study. A phylogenetic framework does offer an opportunity to re-examine morpho-anatomical systems for additional characters, to correlate traits with molecular clades. However, given the relatively uninformative radular and reproductive systems, finding enough informative characters that could act as synapomorphies to resolve more than 59 species is unlikely (Fišer et al. 2018). Also, although synonyms proposed by Wägele (1990) may represent genetically distinguishable units, these are mostly preserved in ways that are not conducive to connecting specimens to these numbered species clades. Possibilities may however, still lie in sequencing holotypes and fixing those names to clades.

Still more species to be discovered

To date, expansion in either sampling or geographic scope (Wilson et al. 2009, 2013) has increased the number of species known in this complex. Here, we assessed the impact of potential undersampling by extrapolating known species records and found that the true number of cryptic species within this complex may be even greater, as the R–E curves illustrate undersampling of genetic diversity and helped identify geographic regions that require more investigation. The R–E curve results did not reach an asymptotic shape, but ~90% of D. kerguelenensis species expected in the sampled areas appeared to be recovered for this study. At a smaller scale, the Palmer Archipelago R–E curve reaches an asymptote at ~12 species and 250 sampling units. Other well-sampled locations that are approaching an asymptote include the Bransfield Strait and Elephant Island. Even locations such as Shag Rocks, Herdman Bank and the South Sandwich Islands, that consist of 50 specimens or less, still approach asymptotic levels, and likely reflect lower diversity in these areas, perhaps due to size or oceanic isolation (separated from other shallow subsea land masses by deep sea or abyssal plains).

This is interesting compared to regions that are well sampled (<100 specimens collected) but that do not approach an asymptotic shape within the R–E curves, such as Burdwood Bank, the South Orkney Islands or the South Shetland Islands. This likely indicates that high levels of D. kerguelenensis species occur in these three regions. Owing to the location (and ancient supercontinent link), the Burdwood Bank plays an important role in diverting circumpolar ocean flow (Fraysse et al. 2018), and is associated with the uplift of nutrient rich bottom water that in turn supports an abundant production of phytoplankton and subsequently rich levels of biodiversity (Schejter et al. 2020). The South Orkney Islands are bordered by two current regimes (the ACC and Weddell Sea Gyre) and have been recorded to host approximately one-fifth (Barnes et al. 2009; Brasier et al. 2018) of all benthic marine species recorded for the entire Southern Ocean (as estimated by Clarke and Johnston 2003). The South Shetland Islands are considered a transitional ecosystem from which new species are still regularly recorded despite being documented as one of the most well-sampled regions across the Antarctic (Barnes et al. 2008). Given the indications that further Doris species diversity is to be found in these areas, these are likely to also be very rich in general benthic marine fauna and should be considered for future biodiversity work (e.g. Barnes et al. 2009). For D. kerguelenensis, areas that should be of interest for future exploration (due to the low sample numbers but high levels of species diversity detected) include Herdman Bank, the Falkland Islands, Bouvet Island, Kerguelen Plateau, Prydz Bay, and the areas surrounding Davis and Casey stations. Totally unexplored locations for this species also include the south-west of South America where the Humboldt Current comes into contact with this coast upon breaking from the northern ACC.

Large distributions for a directly developing slug

Some of the cryptic species in the D. kerguelenensis complex documented here spanned thousands of kilometres. The Southern Ocean, while connected by major currents, still contains numerous potential dispersal barriers for benthic invertebrates, even for those that do have planktonic larvae. Barriers can include temperature and salinity changes, plankton availability, potentially uninhabitable deep sea areas, and the strength and direction of currents and gyres. These large-scale distributions therefore raise many questions regarding the mechanisms of dispersal in D. kerguelenensis.

In aquatic ecosystems, connectivity can be related to dispersal ability between populations (Jablonski 1991; Palumbi 1994; Hellberg et al. 2002; Shanks et al. 2003; Cowen and Sponaugle 2009; González-Wevar et al. 2021). Broadcast spawners and species with pelagic larval development should exhibit higher levels of connectivity and less genetic structure when compared to benthic or direct developing organisms (Ronce 2007; Gillespie et al. 2012). However, there are examples that contradict this prediction (e.g. Marko 2004; Weersing and Toonen 2009; Mercier et al. 2013; Segovia et al. 2017). There are also many benthic Southern Ocean species with large distributions and benthic development (typically showing small-scale and geographically structured distributions) that have been subject to phylogeographic studies. These include gastropods (Nikula et al. 2011a; Cumming et al. 2014; González-Wevar et al. 2021), chitons (Nikula et al. 2011b) and crustaceans (Nikula et al. 2010).

Dayton et al. (1970) proposed several dispersal mechanisms for benthic organisms including (i) adults rafting on detached benthic organisms, (ii) rafting on mobile organisms, either as egg masses or adults, or (iii) anchor ice removing organisms from the benthos and depositing these elsewhere. Alternatively, egg masses could be laid on sessile organisms that could be subsequently dislodged (see Wilson et al. 2009). All of these aforementioned examples are based on dispersal by forms of rafting and rafting appears to be unlikely to be the major method of dispersal for adult individuals in the D. kerguelenensis complex due to weak adhesive ability of the foot (N. G. Wilson, pers. obs.). Doris kerguelenensis species are direct developers (embryonic period recorded up to 21 months; Hain 1989; Moles et al. 2017a) that feed on sedentary organisms such as demosponges and hexactinellid sponges (reviewed by McDonald and Nybakken 1997; Iken et al. 2002). Although there is a small possibility that some species in the D. kerguelenensis complex have planktonic larvae, there is very little evidence in support of this. If rafting is the method of dispersal, this most likely involves egg masses or juveniles.

Current geological evidence shows that small, transient, ice-free areas with reduced levels of primary production did persist, even through glacial maxima across the Southern Ocean continental shelves (including during the Last Glacial Maxima, Poulin et al. 2002; Thatje et al. 2005; Convey et al. 2009; Pearse et al. 2009). Direct development appears to have experienced strong positive selection within these refugia due to low food conditions, long developmental times and speciation driven by allopatry (Poulin et al. 2002; Pearse et al. 2009; Lau et al. 2020). As ice-free areas on the continental shelf could have enabled in situ survival of benthic fauna, populations may have experienced demographic bottlenecks within these refugia instead of complete eradication (Allcock and Strugnell 2012). If secondary contact of populations occurred before reproductive isolation and post glacial expansion, species would demonstrate widespread distributions (Allcock and Strugnell 2012).

In addition to the exentsive geographic distributions, these nudibranchs have been widely reported from depths of up to 1550 m (Iken et al. 2002). The morphologically identified sample from the deepest depth was collected from Halley Bay, in 1998 from 1549 m (Avila et al. 1999). This sample is stored at the University of Barcelona and was fixed in formalin (C. Avila, pers. obs.) and was therefore not used in this study. Here we provide the deepest barcoded record of D. kerguelenensis (798 m). Wägele (1987) similarly reports specimens from 788 m from the Weddell Sea. Although one D. kerguelenensis specimen was reported from New Caledonia from a depth of 680 m (Valdes 2001), given the large geographic gap between this single specimen and all other collected records of this species complex, this specimen likely represents a different taxon. The general absence of D. kerguelenensis from the deep sea areas surrounding Antarctica is worth noting; throughout the ANDEEP expeditions (I, II and III: Brandt et al. 2004, 2007), no D. kerguelenensis were collected from abyssal areas. The ANDEEP project set out to understand the abyssal diversity and faunal exchange between South America and Antarctica (Brandt et al. 2004). This question remains unanswered for this group of organisms and how these benthic organisms disperse across the abyssal plains of the Southern Ocean is currently unknown.

Evidence for dispersal across the Polar Front

Despite all records of D. kerguelenensis being restricted to the continental shelf (see above), some species within this complex have clearly dispersed long distances, including across the APF. These trans-APF species appear in our COI tree in four separate places and nine species were restricted to the Southern American continental shelf. The Antarctic Polar Front represents a strong geographic and oceanographic divide that across time, has split evolutionary lineages between the Antarctic and northern oceans (Page and Linse 2002; Lee et al. 2004; Hunter and Halanych 2008; Thornhill et al. 2008; Wilson et al. 2009, 2013; Krabbe et al. 2010). Very few benthic marine taxa have been recorded to span this region except a sea star (Moore et al. 2018), a brittle star (Galaska et al. 2017), sea spiders (e.g. Linse et al. 2006; Munilla and Membrives 2009; Dietz et al. 2019), an isopod species (Leese et al. 2010) and a tritoniid nudibranch (Moles et al. 2021). Most of these species have a dispersive larval stage except for the isopod that has had long-distance dispersal linked to rafting. Additionally, the only directly developing nudibranch example listed is Tritonia vorax, a sub-Antarctic species that has been recorded from the Southern South American continental shelf and South Georgia in the Scotia Arc (Moles et al. 2021). All of these individuals or populations, unless already in deep water, also overcame the temperature gradient at the Antarctic Polar Front (APF) (3–4°C) that differentiates the SO from more northerly oceans (see Park et al. 2014). The APF is one of the world’s strongest open ocean barriers that is known to isolate the Southern Ocean from warmer waters at lower latitudes and was formed after the opening of a deep-water passage (the Drake Passage) between South America and Antarctica c. 35 Ma (Livermore et al. 2004; Pfuhl and McCave 2005; Barker et al. 2007). The Antarctic Circumpolar Current (ACC) and major Antarctic gyres (in the Weddell and Ross Seas), however, can act as large-scale dispersal vectors that either jet eastward around Antarctica and meet the southernmost region of Southern America (ACC) or rotate clockwise (gyres) due to the interaction between the APF and ACC. In principle, these currents may have driven the dispersal of D. kerguelenensis (Barker and Thomas 2004), however this mechanism remains unclear.


Future directions

Antarctica has long been associated with the diversification and speciation of many benthic marine taxa (Briggs 2003; Rogers 2007; Strugnell et al. 2008), including heterobranch gastropods (Wägele et al. 2008; Martynov and Schrödl 2009; Moles et al. 2017b). In support of this, we found that there is a multitude of undescribed species within the Antarctic species complex D. kerguelenensis (Bergh, 1884) (59 detected and more than 68 estimated). Although the monophyly of the complex is well supported, the interspecific relationships remain unresolved. Genetic barcoding of the complex, that is characterised by intra- and interspecific morphological similarities, appears effective in delimiting species and highlights how molecular phylogenetics can uncover cryptic diversity. To gain phylogenetic resolution within the complex, we now need large quantities of sequence data from across the genome. Future work could incorporate transcriptomes, exon capture or ultra-conserved elements (UCEs). Several recent studies have shown the power of these methods for marine invertebrate taxa (e.g. Horowitz et al. 2020; Layton et al. 2020; Richards et al. 2020). Better knowledge of this species complex will help us to understand the evolutionary dynamics of benthic taxa in Antarctica and provide a phylogenetic scaffold for tracing the evolution of secondary metabolites (chemical defence compounds) through time and across species.

Also, very few Antarctic marine studies have sampled widely enough to assess whether benthic Antarctic invertebrates are truly circumpolar (but see Allcock et al. 2011; Arango et al. 2011; Hemery et al. 2012; Moore et al. 2018; Moles et al. 2021). Most studies are restricted by the logistical challenges that occur when sampling such a vast and remote ecosystem. Future research should have a strategic focus on exploring distributions and understudied regions, specifically the Kerguelen Plateau and East Antarctica as these are likely to contain additional undiscovered diversity. Understanding the extent of diversity is especially important because marine systems worldwide are being affected by climate change, pollution, biological invasions and a host of other anthropogenic disturbances. If no immediate action is taken to address the increasing loss to biodiversity, Antarctic marine systems, along with other unique marine environments spanning the globe, will be at risk of degradation and homogenisation.


Data availability

Supplementary data for this article can be accessed in the Supplementary material. COI sequences are available through GenBank (https://www.ncbi.nlm.nih.gov/): ON419127-ON419135.


Conflicts of interest

The authors declare that Dr Nerida Wilson and Dr Greg Rouse are both Editors of Invertebrate Systematics, but did not at any stage have editor-level access to this manuscript while in peer review, as is the standard practice when handling manuscripts submitted by an editor to this journal. Invertebrate Systematics encourages editors to publish in the journal and are kept totally separate from the decision-making processes for manuscripts. The authors have no further conflicts of interest to declare.


Declaration of funding

The work was supported by the Antarctic Circumnavigation Expedition (carried out by the Swiss Polar Institute, supported by the ACE Foundation and Ferring Pharmaceuticals), the French Polar Institute and LTSER ZATA (#1044), the US National Science Foundation (PLR-1341485, ANT-0551969 to A. L. Moran, ANT-0440577 to H. A. Woods, and DEB-1846174 to K. M. Kocot), BLUEBIO (CTM2016–78901/ANT), the Society of Australian Systematic Biologists (SASB), the University of Western Australia Oceans Institute (UWA-OI) Robson and Robertson award, the Malacological Society of Australasia (MSA) and the Antarctic Science Foundation (ASF). This work was also supported by ARC SRIEAS Grant SR200100005 Securing Antarctica’s Environmental Future.


Supplementary material

Supplementary material is available online.



Acknowledgements

We thank staff and students at these institutions for helping make samples available: Western Australian Museum (WAM), the Baker Laboratory; University of South Florida (USF), Scripps Institution of Oceanography (SIO-BIC), the Australian Antarctic Division (AAD), the Alfred Wegener Institute (AWI), the Yale Peabody Museum (YPM), the California Academy of Sciences (CAS) and the National Museum of Natural History, Smithsonian Institution (USNM) for making samples available. In particular, we appreciate the efforts of Charlotte Seid (SIO-BIC), Sarah Dietrick and Nicole Avalon (USF), Elie Poulin (UCHILE), Javier Naretto (Costa Humboldt) and Christoph Held (AWI) in supporting this work. We are grateful to Mia Hillyer, Alex Hickling and Linette Umbrello (all WAM) for extraction and sequencing assistance, and also to Linette Umbrello for analytical advice and assistance. We thank all expedition participants for help with collecting samples. Finally, we thank Dr Ana Riesgo, the Invertebrate Systematics publishing team and the two anonymous reviewers for their insight and advice.


References

Allcock, AL, and Strugnell, JM (2012). Southern Ocean diversity: new paradigms from molecular ecology. Trends in Ecology & Evolution 27, 520–528.
Southern Ocean diversity: new paradigms from molecular ecology.Crossref | GoogleScholarGoogle Scholar |

Allcock, AL, Barratt, I, Eléaume, M, Linse, K, Norman, MD, Smith, PJ, Steinke, D, Stevens, DW, and Strugnell, JM (2011). Cryptic speciation and the circumpolarity debate: a case study on endemic Southern Ocean octopuses using the COI barcode of life. Deep-sea Research – II. Topical Studies in Oceanography 58, 242–249.
Cryptic speciation and the circumpolarity debate: a case study on endemic Southern Ocean octopuses using the COI barcode of life.Crossref | GoogleScholarGoogle Scholar |

Arango, CP, Soler-Membrives, A, and Miller, KJ (2011). Genetic differentiation in the circum-Antarctic sea spider Nymphon australe (Pycnogonida; Nymphonidae). Deep-sea Research – II. Topical Studies in Oceanography 58, 212–219.
Genetic differentiation in the circum-Antarctic sea spider Nymphon australe (Pycnogonida; Nymphonidae).Crossref | GoogleScholarGoogle Scholar |

Aronson, RB, Thatje, S, McClintock, JB, and Hughes, KA (2011). Anthropogenic impacts on marine ecosystems in Antarctica. Annals of the New York Academy of Sciences 1223, 82–107.
Anthropogenic impacts on marine ecosystems in Antarctica.Crossref | GoogleScholarGoogle Scholar | 21449967PubMed |

Avila, C (2020). Terpenoids in marine heterobranch molluscs. Marine Drugs 18, 162.
Terpenoids in marine heterobranch molluscs.Crossref | GoogleScholarGoogle Scholar |

Avila, C, Iken, K, and Beyer, K (1999). Chemical ecology of opisthobranchs and related species.  Ber Polarforsch 301, 184–189.

Avila, C, Angulo-Preckler, C, Martín-Martín, RP, Figuerola, B, Griffiths, HJ, and Waller, CL (2020). Invasive marine species discovered on non-native kelp rafts in the warmest Antarctic island. Scientific Reports 10, 1639.
Invasive marine species discovered on non-native kelp rafts in the warmest Antarctic island.Crossref | GoogleScholarGoogle Scholar | 32005904PubMed |

Baird, HP, Miller, KJ, and Stark, JS (2011). Evidence of hidden biodiversity, ongoing speciation and diverse patterns of genetic structure in giant Antarctic amphipods. Molecular Ecology 20, 3439–3454.
Evidence of hidden biodiversity, ongoing speciation and diverse patterns of genetic structure in giant Antarctic amphipods.Crossref | GoogleScholarGoogle Scholar | 21733028PubMed |

Barker, PF, and Thomas, E (2004). Origin, signature and palaeoclimatic influence of the Antarctic Circumpolar Current. Earth-Science Reviews 66, 143–162.
Origin, signature and palaeoclimatic influence of the Antarctic Circumpolar Current.Crossref | GoogleScholarGoogle Scholar |

Barker, PF, Filippelli, GM, Florindo, F, Martin, EE, and Scher, HD (2007). Onset and role of the Antarctic Circumpolar Current. Deep-sea Research – II. Topical Studies in Oceanography 54, 2388–2398.
Onset and role of the Antarctic Circumpolar Current.Crossref | GoogleScholarGoogle Scholar |

Barnes, DKA, Linse, K, Enderlein, P, Smale, D, Fraser, KPP, and Brown, M (2008). Marine richness and gradients at Deception Island, Antarctica. Antarctic Science 20, 271–280.
Marine richness and gradients at Deception Island, Antarctica.Crossref | GoogleScholarGoogle Scholar |

Barnes, DKA, Kaiser, S, Griffiths, HJ, and Linse, K (2009). Marine, intertidal, freshwater and terrestrial biodiversity of an isolated polar archipelago. Journal of Biogeography 36, 756–769.
Marine, intertidal, freshwater and terrestrial biodiversity of an isolated polar archipelago.Crossref | GoogleScholarGoogle Scholar |

Barnosky, AD, Matzke, N, Tomiya, S, Wogan, GOU, Swartz, B, Quental, TB, Marshall, C, McGuire, JL, Lindsey, EL, Maguire, KC, Mersey, B, and Ferrer, EA (2011). Has the Earth’s sixth mass extinction already arrived? Nature 471, 51–57.
Has the Earth’s sixth mass extinction already arrived?Crossref | GoogleScholarGoogle Scholar | 21368823PubMed |

Bouchet P (2006) The magnitude of marine biodiversity. In ‘The Exploration of Marine Biodiversity: Scientific and Technological Challenges’. (Ed. CM Duarte) pp. 31–64. (Fundación BBVA)

Brandt, A, De Broyer, C, Gooday, AJ, Hilbig, B, and Thomson, MRA (2004). Introduction to ANDEEP (Antarctic benthic deep-sea biodiversity: colonization history and recent community patterns) – a tribute to Howard L. Sanders. Deep-sea Research – II. Topical Studies in Oceanography 51, 1457–1465.
Introduction to ANDEEP (Antarctic benthic deep-sea biodiversity: colonization history and recent community patterns) – a tribute to Howard L. Sanders.Crossref | GoogleScholarGoogle Scholar |

Brandt, A, Ebbe, B, and Gooday, AJ (2007). Introduction to ANDEEP, summary and outlook. Deep-sea Research – II. Topical Studies in Oceanography 54, 1645–1651.
Introduction to ANDEEP, summary and outlook.Crossref | GoogleScholarGoogle Scholar |

Brandt, A, Błażewicz-Paszkowycz, M, Bamber, R, Mühlenhardt-Siegel, U, Malyutina, M, Kaiser, S, De Broyer, C, and Havermans, C (2012). Are there widespread peracarid species in the deep sea (Crustacea: Malacostraca)? Polish Polar Research 33, 139–162.
Are there widespread peracarid species in the deep sea (Crustacea: Malacostraca)?Crossref | GoogleScholarGoogle Scholar |

Brasier, MJ, Wiklund, H, Neal, L, Jeffreys, R, Linse, K, Ruhl, H, and Glover, AG (2016). DNA barcoding uncovers cryptic diversity in 50% of deep-sea Antarctic polychaetes. Royal Society Open Science 3, 160432.
DNA barcoding uncovers cryptic diversity in 50% of deep-sea Antarctic polychaetes.Crossref | GoogleScholarGoogle Scholar | 28018624PubMed |

Brasier, MJ, Grant, SM, Trathan, PN, Allcock, L, Ashford, O, Blagbrough, H, Brandt, A, Danis, B, Downey, R, Eléaume, MP, Enderlein, P, Ghiglione, C, Hogg, O, Linse, K, Mackenzie, M, Moreau, C, Robinson, LF, Rodriguez, E, Spiridonov, V, Tate, A, Taylor, M, Waller, C, Wiklund, H, and Griffiths, HJ (2018). Benthic biodiversity in the South Orkney Islands Southern Shelf Marine Protected Area. Biodiversity 19, 5–19.
Benthic biodiversity in the South Orkney Islands Southern Shelf Marine Protected Area.Crossref | GoogleScholarGoogle Scholar |

Briggs, JC (2003). Marine centres of origin as evolutionary engines. Journal of Biogeography 30, 1–18.
Marine centres of origin as evolutionary engines.Crossref | GoogleScholarGoogle Scholar |

Cai, H, Li, Q, Fang, X, Li, J, Curtis, NE, Altenburger, A, Shibata, T, Feng, M, Maeda, T, Schwartz, JA, Shigenobu, S, Lundholm, N, Nishiyama, T, Yang, H, Hasebe, M, Li, S, Pierce, SK, and Wang, J (2019). A draft genome assembly of the solar-powered sea slug Elysia chlorotica. Scientific Data 6, 190022.
A draft genome assembly of the solar-powered sea slug Elysia chlorotica.Crossref | GoogleScholarGoogle Scholar | 30778257PubMed |

Cattaneo-Vietti, R (1991). Nudibranch molluscs from the Ross Sea, Antarctica. The Journal of Molluscan Studies 57, 223–228.
Nudibranch molluscs from the Ross Sea, Antarctica.Crossref | GoogleScholarGoogle Scholar |

Ceballos, G, Ehrlich, PR, Barnosky, AD, García, A, Pringle, RM, and Palmer, TM (2015). Accelerated modern human–induced species losses: entering the sixth mass extinction. Science Advances 1, e1400253.
Accelerated modern human–induced species losses: entering the sixth mass extinction.Crossref | GoogleScholarGoogle Scholar | 26601195PubMed |

Chapman, WL, and Walsh, JE (2007). A synthesis of Antarctic temperatures. Journal of Climate 20, 4096–4117.
A synthesis of Antarctic temperatures.Crossref | GoogleScholarGoogle Scholar |

Chenuil, A, Saucède, T, Hemery, LG, Eléaume, M, Féral, J-P, Améziane, N, David, B, Lecointre, G, and Havermans, C (2018). Understanding processes at the origin of species flocks with a focus on the marine Antarctic fauna. Biological Reviews of the Cambridge Philosophical Society 93, 481–504.
Understanding processes at the origin of species flocks with a focus on the marine Antarctic fauna.Crossref | GoogleScholarGoogle Scholar | 28799256PubMed |

Chown, SL, Clarke, A, Fraser, CI, Cary, SC, Moon, KL, and McGeoch, MA (2015). The changing form of Antarctic biodiversity. Nature 522, 431–438.
The changing form of Antarctic biodiversity.Crossref | GoogleScholarGoogle Scholar | 26108852PubMed |

Clarke, A, and Crame, JA (1989). The origin of the Southern Ocean marine fauna. Geological Society of London, Special Publications 47, 253–268.
The origin of the Southern Ocean marine fauna.Crossref | GoogleScholarGoogle Scholar |

Clarke, A, and Crame, JA (2010). Evolutionary dynamics at high latitudes: speciation and extinction in polar marine faunas. Philosophical Transactions of the Royal Society of London – B. Biological Sciences 365, 3655–3666.
Evolutionary dynamics at high latitudes: speciation and extinction in polar marine faunas.Crossref | GoogleScholarGoogle Scholar | 20980314PubMed |

Clarke, A, and Johnston, NM (2003). Antarctic marine benthic diversity. Oceanography and Marine Biology - an Annual Review 41, 47–114.

Clarke, A, Murphy, EJ, Meredith, MP, King, JC, Peck, LS, Barnes, DKA, and Smith, RC (2007). Climate change and the marine ecosystem of the western Antarctic Peninsula. Philosophical Transactions of the Royal Society of London – B. Biological Sciences 362, 149–166.
Climate change and the marine ecosystem of the western Antarctic Peninsula.Crossref | GoogleScholarGoogle Scholar | 17405211PubMed |

Clement, M, Posada, D, and Crandall, KA (2000). TCS: a computer program to estimate gene genealogies. Molecular Ecology 9, 1657–1659.
TCS: a computer program to estimate gene genealogies.Crossref | GoogleScholarGoogle Scholar | 11050560PubMed |

Colwell, RK, Chao, A, Gotelli, NJ, Lin, S-Y, Mao, CX, Chazdon, RL, and Longino, JT (2012). Models and estimators linking individual-based and sample-based rarefaction, extrapolation and comparison of assemblages. Journal of Plant Ecology 5, 3–21.
Models and estimators linking individual-based and sample-based rarefaction, extrapolation and comparison of assemblages.Crossref | GoogleScholarGoogle Scholar |

Convey, P, Stevens, MI, Hodgson, DA, Smellie, JL, Hillenbrand, C-D, Barnes, DKA, Clarke, A, Pugh, PJA, Linse, K, and Cary, SC (2009). Exploring biological constraints on the glacial history of Antarctica. Quaternary Science Reviews 28, 3035–3048.
Exploring biological constraints on the glacial history of Antarctica.Crossref | GoogleScholarGoogle Scholar |

Cowen, RK, and Sponaugle, S (2009). Larval dispersal and marine population connectivity. Annual Review of Marine Science 1, 443–466.
Larval dispersal and marine population connectivity.Crossref | GoogleScholarGoogle Scholar | 21141044PubMed |

Crame, JA (1997). An evolutionary framework for the polar regions. Journal of Biogeography 24, 1–9.
An evolutionary framework for the polar regions.Crossref | GoogleScholarGoogle Scholar |

Cumming, RA, Nikula, R, Spencer, HG, and Waters, JM (2014). Transoceanic genetic similarities of kelp-associated sea slug populations: long-distance dispersal via rafting? Journal of Biogeography 41, 2357–2370.
Transoceanic genetic similarities of kelp-associated sea slug populations: long-distance dispersal via rafting?Crossref | GoogleScholarGoogle Scholar |

Dayrat, B (2005). Towards integrative taxonomy. Biological Journal of the Linnean Society. Linnean Society of London 85, 407–415.
Towards integrative taxonomy.Crossref | GoogleScholarGoogle Scholar |

Dayton, P, Robilliard, GA, and Paine, RT (1970). Benthic faunal zonation as a result of anchor ice at McMurdo Sound, Antartica. Antarctic. Ecology 1, 244–258.

De Broyer, C, and Danis, B (2011). How many species in the Southern Ocean? Towards a dynamic inventory of the Antarctic marine species. Deep-sea Research – II. Topical Studies in Oceanography 58, 5–17.
How many species in the Southern Ocean? Towards a dynamic inventory of the Antarctic marine species.Crossref | GoogleScholarGoogle Scholar |

De Castro-Fernández, P, Cardona, L, and Avila, C (2021). Distribution of trace elements in benthic infralittoral organisms from the western Antarctic Peninsula reveals no latitudinal gradient of pollution. Scientific Reports 11, 16266.
Distribution of trace elements in benthic infralittoral organisms from the western Antarctic Peninsula reveals no latitudinal gradient of pollution.Crossref | GoogleScholarGoogle Scholar | 34381092PubMed |

Dietz, L, Dömel, JS, Leese, F, Mahon, AR, and Mayer, C (2019). Phylogenomics of the longitarsal Colossendeidae: the evolutionary history of an Antarctic sea spider radiation. Molecular Phylogenetics and Evolution 136, 206–214.
Phylogenomics of the longitarsal Colossendeidae: the evolutionary history of an Antarctic sea spider radiation.Crossref | GoogleScholarGoogle Scholar | 31002869PubMed |

Domke, GM, Oswalt, SN, Walters, BF, and Morin, RS (2020). Tree planting has the potential to increase carbon sequestration capacity of forests in the United States. Proceedings of the National Academy of Sciences of the United States of America 117, 24649–24651.
Tree planting has the potential to increase carbon sequestration capacity of forests in the United States.Crossref | GoogleScholarGoogle Scholar | 32958649PubMed |

Eastman, JT, and McCune, AR (2000). Fishes on the Antarctic continental shelf: evolution of amarine species flock? Journal of Fish Biology 57, 84–102.
Fishes on the Antarctic continental shelf: evolution of amarine species flock?Crossref | GoogleScholarGoogle Scholar |

Eberle, J, Ahrens, D, Mayer, C, Niehuis, O, and Misof, B (2020). A plea for standardized nuclear markers in metazoan DNA taxonomy. Trends in Ecology & Evolution 35, 336–345.
A plea for standardized nuclear markers in metazoan DNA taxonomy.Crossref | GoogleScholarGoogle Scholar |

Féral, J-P (2002). How useful are the genetic markers in attempts to understand and manage marine biodiversity? Journal of Experimental Marine Biology and Ecology 268, 121–145.
How useful are the genetic markers in attempts to understand and manage marine biodiversity?Crossref | GoogleScholarGoogle Scholar |

Fišer, C, Robinson, CT, and Malard, F (2018). Cryptic species as a window into the paradigm shift of the species concept. Molecular Ecology 27, 613–635.
Cryptic species as a window into the paradigm shift of the species concept.Crossref | GoogleScholarGoogle Scholar | 29334414PubMed |

Folmer, O, Black, M, Wr, H, Lutz, R, and Vrijenhoek, R (1994). DNA primers for amplification of mitochondrial cytochrome c oxidase subunit I from diverse metazoan invertebrates. Molecular Marine Biology and Biotechnology 3, 294–299.
| 7881515PubMed |

Fraysse, C, Calcagno, J, and Pérez, AF (2018). Asteroidea of the southern tip of South America, including Namuncurá Marine Protected Area at Burdwood Bank and Tierra del Fuego Province, Argentina. Polar Biology 41, 2423–2433.
Asteroidea of the southern tip of South America, including Namuncurá Marine Protected Area at Burdwood Bank and Tierra del Fuego Province, Argentina.Crossref | GoogleScholarGoogle Scholar |

Galaska, MP, Sands, CJ, Santos, SR, Mahon, AR, and Halanych, KM (2017). Crossing the divide: admixture across the Antarctic polar front revealed by the brittle star Astrotoma agassizii. The Biological Bulletin 232, 198–211.
Crossing the divide: admixture across the Antarctic polar front revealed by the brittle star Astrotoma agassizii.Crossref | GoogleScholarGoogle Scholar | 28898598PubMed |

Geller, J, Meyer, C, Parker, M, and Hawk, H (2013). Redesign of PCR primers for mitochondrial cytochrome c oxidase subunit I for marine invertebrates and application in all-taxa biotic surveys. Molecular Ecology Resources 13, 851–861.
Redesign of PCR primers for mitochondrial cytochrome c oxidase subunit I for marine invertebrates and application in all-taxa biotic surveys.Crossref | GoogleScholarGoogle Scholar | 23848937PubMed |

Gillespie, RG, Baldwin, BG, Waters, JM, Fraser, CI, Nikula, R, and Roderick, GK (2012). Long-distance dispersal: a framework for hypothesis testing. Trends in Ecology & Evolution 27, 47–56.
Long-distance dispersal: a framework for hypothesis testing.Crossref | GoogleScholarGoogle Scholar |

Goetze, E (2003). Cryptic speciation on the high seas; global phylogenetics of the copepod family Eucalanidae. Proceedings of the Royal Society of London – B. Biological Sciences 270, 2321–2331.
Cryptic speciation on the high seas; global phylogenetics of the copepod family Eucalanidae.Crossref | GoogleScholarGoogle Scholar |

González-Wevar, CA, Segovia, NI, Rosenfeld, S, Noll, D, Maturana, CS, Hüne, M, Naretto, J, Gérard, K, Díaz, A, Spencer, HG, Saucède, T, Féral, J-P, Morley, SA, Brickle, P, Wilson, NG, and Poulin, E (2021). Contrasting biogeographical patterns in Margarella (Gastropoda: Calliostomatidae: Margarellinae) across the Antarctic Polar Front. Molecular Phylogenetics and Evolution 156, 107039.
Contrasting biogeographical patterns in Margarella (Gastropoda: Calliostomatidae: Margarellinae) across the Antarctic Polar Front.Crossref | GoogleScholarGoogle Scholar | 33310059PubMed |

Goodheart, JA, Bleidißel, S, Schillo, D, Strong, EE, Ayres, DL, Preisfeld, A, Collins, AG, Cummings, MP, and Wägele, H (2018). Comparative morphology and evolution of the cnidosac in Cladobranchia (Gastropoda: Heterobranchia: Nudibranchia). Frontiers in Zoology 15, 43.
Comparative morphology and evolution of the cnidosac in Cladobranchia (Gastropoda: Heterobranchia: Nudibranchia).Crossref | GoogleScholarGoogle Scholar | 30473719PubMed |

Gotelli NJ, Colwell RK (2011) Estimating species richness. In ‘Biological diversity: frontiers in measurement and assessment’. (Eds AE Magurran, BJ McGill) pp. 39–54. (Oxford University Press)

Grant, RA, Griffiths, HJ, Steinke, D, Wadley, V, and Linse, K (2011). Antarctic DNA barcoding; a drop in the ocean? Polar Biology 34, 775–780.
Antarctic DNA barcoding; a drop in the ocean?Crossref | GoogleScholarGoogle Scholar |

Griffiths, HJ (2010). Antarctic Marine Biodiversity – What Do We Know About the Distribution of Life in the Southern Ocean? PLoS One 5, e11683.
Antarctic Marine Biodiversity – What Do We Know About the Distribution of Life in the Southern Ocean?Crossref | GoogleScholarGoogle Scholar | 20689841PubMed |

Gutt J, Starmans A (2002) Quantification of iceberg impact and benthic recolonisation patterns in the Weddell Sea (Antarctica). In ‘Ecological Studies in the Antarctic Sea Ice Zone: Results of EASIZ Midterm Symposium’. (Eds WE Arntz, A Clarke) pp. 210–214. (Springer)
| Crossref |

Hain S (1989) Beiträge zur Biologie der beschalten Mollusken (Kl. Gastropoda & Bivalvia) des Weddellmeeres, Antarktis. University of Bremen, Bremen, Germany.

Hain S, Arnaud PM (1992) Notes on the reproduction of high-Antarctic molluscs from the Weddell Sea. In ‘Weddell Sea Ecology’. (Ed. G Hempel) pp. 303–312. (Springer) 10.1007/978-3-642-77595-6_35

Hart, MW, and Sunday, J (2007). Things fall apart: biological species form unconnected parsimony networks. Biology Letters 3, 509–512.
Things fall apart: biological species form unconnected parsimony networks.Crossref | GoogleScholarGoogle Scholar | 17650475PubMed |

Hebert, PDN, Cywinska, A, Ball, SL, and deWaard, JR (2003). Biological identifications through DNA barcodes. Proceedings of the Royal Society of London – B. Biological Sciences 270, 313–321.
Biological identifications through DNA barcodes.Crossref | GoogleScholarGoogle Scholar |

Hellberg, ME, Burton, RS, Neigel, JE, and Palumbi, SR (2002). Genetic assessment of connectivity among marine populations. Bulletin of Marine Science 70, 273–290.

Hemery, LG, Eléaume, M, Roussel, V, Améziane, N, Gallut, C, Steinke, D, Cruaud, C, Couloux, A, and Wilson, NG (2012). Comprehensive sampling reveals circumpolarity and sympatry in seven mitochondrial lineages of the Southern Ocean crinoid species Promachocrinus kerguelensis (Echinodermata). Molecular Ecology 21, 2502–2518.
Comprehensive sampling reveals circumpolarity and sympatry in seven mitochondrial lineages of the Southern Ocean crinoid species Promachocrinus kerguelensis (Echinodermata).Crossref | GoogleScholarGoogle Scholar | 22469156PubMed |

Hoang, DT, Chernomor, O, von Haeseler, A, Minh, BQ, and Vinh, LS (2018). UFBoot2: improving the ultrafast bootstrap approximation. Molecular Biology and Evolution 35, 518–522.
UFBoot2: improving the ultrafast bootstrap approximation.Crossref | GoogleScholarGoogle Scholar | 29077904PubMed |

Horowitz, J, Brugler, MR, Bridge, TCL, and Cowman, PF (2020). Morphological and molecular description of a new genus and species of black coral (Cnidaria: Anthozoa: Hexacorallia: Antipatharia: Antipathidae: Blastopathes) from Papua New Guinea. Zootaxa 4821, 553–569.
Morphological and molecular description of a new genus and species of black coral (Cnidaria: Anthozoa: Hexacorallia: Antipatharia: Antipathidae: Blastopathes) from Papua New Guinea.Crossref | GoogleScholarGoogle Scholar |

Hsieh, TC, Ma, KH, and Chao, A (2016). iNEXT: an R package for rarefaction and extrapolation of species diversity (Hill numbers). Methods in Ecology and Evolution 7, 1451–1456.
iNEXT: an R package for rarefaction and extrapolation of species diversity (Hill numbers).Crossref | GoogleScholarGoogle Scholar |

Hulett, RE, Mahguib, J, Gosliner, TM, and Valdés, Á (2015). Molecular evaluation of the phylogenetic position of the enigmatic species Trivettea papalotla (Bertsch) (Mollusca : Nudibranchia). Invertebrate Systematics 29, 215–222.
Molecular evaluation of the phylogenetic position of the enigmatic species Trivettea papalotla (Bertsch) (Mollusca : Nudibranchia).Crossref | GoogleScholarGoogle Scholar |

Hunter, RL, and Halanych, KM (2008). Evaluating connectivity in the brooding brittle star Astrotoma agassizii across the Drake Passage in the Southern Ocean. The Journal of Heredity 99, 137–148.
Evaluating connectivity in the brooding brittle star Astrotoma agassizii across the Drake Passage in the Southern Ocean.Crossref | GoogleScholarGoogle Scholar | 18245798PubMed |

Iken, K, Avila, C, Fontana, A, and Gavagnin, M (2002). Chemical ecology and origin of defensive compounds in the Antarctic nudibranch Austrodoris kerguelenensis (Opisthobranchia: Gastropoda). Marine Biology 141, 101–109.
Chemical ecology and origin of defensive compounds in the Antarctic nudibranch Austrodoris kerguelenensis (Opisthobranchia: Gastropoda).Crossref | GoogleScholarGoogle Scholar |

Ingels, J, Vanreusel, A, Brandt, A, Catarino, AI, David, B, Ridder, CD, Dubois, P, Gooday, AJ, Martin, P, Pasotti, F, and Robert, H (2012). Possible effects of global environmental changes on Antarctic benthos: a synthesis across five major taxa. Ecology and Evolution 2, 453–485.
Possible effects of global environmental changes on Antarctic benthos: a synthesis across five major taxa.Crossref | GoogleScholarGoogle Scholar | 22423336PubMed |

Ito, T, Woloszyn, M, and Mazloff, M (2010). Anthropogenic carbon dioxide transport in the Southern Ocean driven by Ekman flow. Nature 463, 80–83.
Anthropogenic carbon dioxide transport in the Southern Ocean driven by Ekman flow.Crossref | GoogleScholarGoogle Scholar | 20054394PubMed |

Jablonski, D (1991). Extinctions: a paleontological perspective. Science 253, 754–757.
| 17835491PubMed |

Jensen, A, Lillie, M, Bergström, K, Larsson, P, and Höglund, J (2021). Whole genome sequencing reveals high differentiation, low levels of genetic diversity and short runs of homozygosity among Swedish wels catfish. Heredity 127, 79–91.
Whole genome sequencing reveals high differentiation, low levels of genetic diversity and short runs of homozygosity among Swedish wels catfish.Crossref | GoogleScholarGoogle Scholar | 33963302PubMed |

Johnson, SB, Warén, A, and Vrijenhoek, RC (2008). DNA barcoding of Lepetodrilus limpets reveals cryptic species. Journal of Shellfish Research 27, 43–51.
DNA barcoding of Lepetodrilus limpets reveals cryptic species.Crossref | GoogleScholarGoogle Scholar |

Jukes TH, Cantor CR (1969) Evolution of protein molecules. In ‘Mammalian Protein Metabolism’. (Ed. HN Munro) pp. 21–132. (Academic Press: New York, NY, USA)

Jung, D, Lee, J, and Kim, C-B (2014). A Report on five new records of nudibranch molluscs from Korea. Animal Systematics, Evolution and Diversity 30, 124–131.
A Report on five new records of nudibranch molluscs from Korea.Crossref | GoogleScholarGoogle Scholar |

Kalyaanamoorthy, S, Minh, BQ, Wong, TKF, von Haeseler, A, and Jermiin, LS (2017). ModelFinder: fast model selection for accurate phylogenetic estimates. Nature Methods 14, 587–589.
ModelFinder: fast model selection for accurate phylogenetic estimates.Crossref | GoogleScholarGoogle Scholar | 28481363PubMed |

Kapli, P, Lutteropp, S, Zhang, J, Kobert, K, Pavlidis, P, Stamatakis, A, and Flouri, T (2017). Multi-rate Poisson tree processes for single-locus species delimitation under maximum likelihood and Markov chain Monte Carlo. Bioinformatics 33, 1630–1638.
Multi-rate Poisson tree processes for single-locus species delimitation under maximum likelihood and Markov chain Monte Carlo.Crossref | GoogleScholarGoogle Scholar | 28108445PubMed |

Katoh, K, and Standley, DM (2013). MAFFT multiple sequence alignment software version 7: improvements in performance and usability. Molecular Biology and Evolution 30, 772–780.
MAFFT multiple sequence alignment software version 7: improvements in performance and usability.Crossref | GoogleScholarGoogle Scholar | 23329690PubMed |

Kearse, M, Moir, R, Wilson, A, Stones-Havas, S, Cheung, M, Sturrock, S, Buxton, S, Cooper, A, Markowitz, S, Duran, C, Thierer, T, Ashton, B, Meintjes, P, and Drummond, A (2012). Geneious Basic: an integrated and extendable desktop software platform for the organization and analysis of sequence data. Bioinformatics 28, 1647–1649.
Geneious Basic: an integrated and extendable desktop software platform for the organization and analysis of sequence data.Crossref | GoogleScholarGoogle Scholar | 22543367PubMed |

Kimura, M (1980). A simple method for estimating evolutionary rates of base substitutions through comparative studies of nucleotide sequences. Journal of Molecular Evolution 16, 111–120.
A simple method for estimating evolutionary rates of base substitutions through comparative studies of nucleotide sequences.Crossref | GoogleScholarGoogle Scholar | 7463489PubMed |

Knowlton, N (1993). Sibling species in the sea. Annual Review of Ecology and Systematics 24, 189–216.
Sibling species in the sea.Crossref | GoogleScholarGoogle Scholar |

Korshunova, T, Malmberg, K, Prkić, J, Petani, A, Fletcher, K, Lundin, K, and Martynov, A (2020). Fine-scale species delimitation: speciation in process and periodic patterns in nudibranch diversity. ZooKeys 917, 15–50.
Fine-scale species delimitation: speciation in process and periodic patterns in nudibranch diversity.Crossref | GoogleScholarGoogle Scholar | 32206016PubMed |

Krabbe, K, Leese, F, Mayer, C, Tollrian, R, and Held, C (2010). Cryptic mitochondrial lineages in the widespread pycnogonid Colossendeis megalonyx Hoek, 1881 from Antarctic and Subantarctic waters. Polar Biology 33, 281–292.
Cryptic mitochondrial lineages in the widespread pycnogonid Colossendeis megalonyx Hoek, 1881 from Antarctic and Subantarctic waters.Crossref | GoogleScholarGoogle Scholar |

Kumar, S, Stecher, G, Li, M, Knyaz, C, and Tamura, K (2018). MEGA X: molecular evolutionary genetics analysis across computing platforms. Molecular Biology and Evolution 35, 1547–1549.
MEGA X: molecular evolutionary genetics analysis across computing platforms.Crossref | GoogleScholarGoogle Scholar | 29722887PubMed |

Lau, SCY, Wilson, NG, Silva, CNS, and Strugnell, JM (2020). Detecting glacial refugia in the Southern Ocean. Ecography 43, 1639–1656.
Detecting glacial refugia in the Southern Ocean.Crossref | GoogleScholarGoogle Scholar |

Layton, KKS, Gosliner, TM, and Wilson, NG (2018). Flexible colour patterns obscure identification and mimicry in Indo-Pacific Chromodoris nudibranchs (Gastropoda: Chromodorididae). Molecular Phylogenetics and Evolution 124, 27–36.
Flexible colour patterns obscure identification and mimicry in Indo-Pacific Chromodoris nudibranchs (Gastropoda: Chromodorididae).Crossref | GoogleScholarGoogle Scholar | 29476907PubMed |

Layton, KKS, Carvajal, JI, and Wilson, NG (2020). Mimicry and mitonuclear discordance in nudibranchs: new insights from exon capture phylogenomics. Ecology and Evolution 10, 11966–11982.
Mimicry and mitonuclear discordance in nudibranchs: new insights from exon capture phylogenomics.Crossref | GoogleScholarGoogle Scholar | 33209263PubMed |

Lecointre, G, Améziane, N, Boisselier, M-C, Bonillo, C, Busson, F, Causse, R, Chenuil, A, Couloux, A, Coutanceau, J-P, Cruaud, C, d’Acoz d’Udekem, C, Ridder, CD, Denys, G, Dettaï, A, Duhamel, G, Eléaume, M, Féral, J-P, Gallut, C, Havermans, C, Held, C, Hemery, L, Lautrédou, A-C, Martin, P, Ozouf-Costaz, C, Pierrat, B, Pruvost, P, Puillandre, N, Samadi, S, Saucède, T, Schubart, C, and David, B (2013). Is the species flock concept operational? The Antarctic Shelf case. PLoS One 8, e68787.
Is the species flock concept operational? The Antarctic Shelf case.Crossref | GoogleScholarGoogle Scholar | 23936311PubMed |

Lee, Y-H, Song, M, Lee, S, Leon, R, Godoy, SO, and Canete, I (2004). Molecular phylogeny and divergence time of the Antarctic sea urchin (Sterechinus neumayeri) in relation to the South American sea urchins. Antarctic Science 16, 29–36.
Molecular phylogeny and divergence time of the Antarctic sea urchin (Sterechinus neumayeri) in relation to the South American sea urchins.Crossref | GoogleScholarGoogle Scholar |

Leese, F, Agrawal, S, and Held, C (2010). Long-distance island hopping without dispersal stages: transportation across major zoogeographic barriers in a Southern Ocean isopod. Naturwissenschaften 97, 583–594.
Long-distance island hopping without dispersal stages: transportation across major zoogeographic barriers in a Southern Ocean isopod.Crossref | GoogleScholarGoogle Scholar | 20454771PubMed |

Linse, K, Griffiths, HJ, Barnes, DKA, and Clarke, A (2006). Biodiversity and biogeography of Antarctic and sub-Antarctic mollusca. Deep-sea Research – II. Topical Studies in Oceanography 53, 985–1008.
Biodiversity and biogeography of Antarctic and sub-Antarctic mollusca.Crossref | GoogleScholarGoogle Scholar |

Linse, K, Cope, T, Lörz, A-N, and Sands, C (2007). Is the Scotia Sea a centre of Antarctic marine diversification? Some evidence of cryptic speciation in the circum-Antarctic bivalve Lissarca notorcadensis (Arcoidea: Philobryidae). Polar Biology 30, 1059–1068.
Is the Scotia Sea a centre of Antarctic marine diversification? Some evidence of cryptic speciation in the circum-Antarctic bivalve Lissarca notorcadensis (Arcoidea: Philobryidae).Crossref | GoogleScholarGoogle Scholar |

Livermore, R, Eagles, G, Morris, P, and Maldonado, A (2004). Shackleton Fracture Zone: No barrier to early circumpolar ocean circulation. Geology 32, 797–800.
Shackleton Fracture Zone: No barrier to early circumpolar ocean circulation.Crossref | GoogleScholarGoogle Scholar |

Locey, KJ, and Lennon, JT (2016). Scaling laws predict global microbial diversity. Proceedings of the National Academy of Sciences of the United States of America 113, 5970–5975.
Scaling laws predict global microbial diversity.Crossref | GoogleScholarGoogle Scholar | 27140646PubMed |

Mahguib, J, and Valdés, Á (2015). Molecular investigation of the phylogenetic position of the polar nudibranch Doridoxa (Mollusca, Gastropoda, Heterobranchia). Polar Biology 38, 1369–1377.
Molecular investigation of the phylogenetic position of the polar nudibranch Doridoxa (Mollusca, Gastropoda, Heterobranchia).Crossref | GoogleScholarGoogle Scholar |

Marko, PB (2004). ‘What’s larvae got to do with it?’ Disparate patterns of post-glacial population structure in two benthic marine gastropods with identical dispersal potential. Molecular Ecology 13, 597–611.
‘What’s larvae got to do with it?’ Disparate patterns of post-glacial population structure in two benthic marine gastropods with identical dispersal potential.Crossref | GoogleScholarGoogle Scholar | 14871364PubMed |

Martynov, AV, and Schrödl, M (2009). The new Arctic side-gilled sea slug genus Boreoberthella (Gastropoda, Opisthobranchia): pleurobranchoidean systematics and evolution revisited. Polar Biology 32, 53–70.
The new Arctic side-gilled sea slug genus Boreoberthella (Gastropoda, Opisthobranchia): pleurobranchoidean systematics and evolution revisited.Crossref | GoogleScholarGoogle Scholar |

Maschek, JA, Mevers, E, Diyabalanage, T, Chen, L, Ren, Y, McClintock, JB, Amsler, CD, Wu, J, and Baker, BJ (2012). Palmadorin chemodiversity from the Antarctic nudibranch Austrodoris kerguelenensis and inhibition of Jak2/STAT5-dependent HEL leukemia cells. Tetrahedron 68, 9095–9104.
Palmadorin chemodiversity from the Antarctic nudibranch Austrodoris kerguelenensis and inhibition of Jak2/STAT5-dependent HEL leukemia cells.Crossref | GoogleScholarGoogle Scholar |

Matsuda, SB, and Gosliner, TM (2018). Molecular phylogeny of Glossodoris (Ehrenberg, 1831) nudibranchs and related genera reveals cryptic and pseudocryptic species complexes. Cladistics 34, 41–56.
Molecular phylogeny of Glossodoris (Ehrenberg, 1831) nudibranchs and related genera reveals cryptic and pseudocryptic species complexes.Crossref | GoogleScholarGoogle Scholar | 34641635PubMed |

Matsuoka, K, Skoglund, A, Roth, G, de Pomereu, J, Griffiths, H, Headland, R, Herried, B, Katsumata, K, Le Brocq, A, Licht, K, Morgan, F, Neff, PD, Ritz, C, Scheinert, M, Tamura, T, Van de Putte, A, van den Broeke, M, von Deschwanden, A, Deschamps-Berger, C, Van Liefferinge, B, Tronstad, S, and Melvær, Y (2021). Quantarctica, an integrated mapping environment for Antarctica, the Southern Ocean, and sub-Antarctic islands. Environmental Modelling & Software 140, 105015.
Quantarctica, an integrated mapping environment for Antarctica, the Southern Ocean, and sub-Antarctic islands.Crossref | GoogleScholarGoogle Scholar |

McClintock JB, Baker BJ (2001) ‘Marine Chemical Ecology.’ (CRC Press)

McDonald, GR, and Nybakken, JW (1997). List of the worldwide food habits of nudibranchs. The Veliger 40, 1–426.

McNeely JA, Miller KR, Reid WV, Mittermeier RA, Werner TB (1990) ‘Conserving the world’s biological diversity.’ (IUCN: Gland, Switzerland; WRI, Conservation International, WWF-US, WorldBank: Washington, DC, USA)

Mercier, A, Sewell, MA, and Hamel, J-F (2013). Pelagic propagule duration and developmental mode: reassessment of a fading link. Global Ecology and Biogeography 22, 517–530.
Pelagic propagule duration and developmental mode: reassessment of a fading link.Crossref | GoogleScholarGoogle Scholar |

Moles, J, Wägele, H, Cutignano, A, Fontana, A, Ballesteros, M, and Avila, C (2017a). Giant embryos and hatchlings of Antarctic nudibranchs (Mollusca: Gastropoda: Heterobranchia). Marine Biology 164, 114.
Giant embryos and hatchlings of Antarctic nudibranchs (Mollusca: Gastropoda: Heterobranchia).Crossref | GoogleScholarGoogle Scholar |

Moles, J, Wägele, H, Schrödl, M, and Avila, C (2017b). A new Antarctic heterobranch clade is sister to all other Cephalaspidea (Mollusca: Gastropoda). Zoologica Scripta 46, 127–137.
A new Antarctic heterobranch clade is sister to all other Cephalaspidea (Mollusca: Gastropoda).Crossref | GoogleScholarGoogle Scholar |

Moles, J, Berning, MI, Hooker, Y, Padula, V, Wilson, NG, and Schrödl, M (2021). Due south: the evolutionary history of Sub-Antarctic and Antarctic Tritoniidae nudibranchs. Molecular Phylogenetics and Evolution 162, 107209.
Due south: the evolutionary history of Sub-Antarctic and Antarctic Tritoniidae nudibranchs.Crossref | GoogleScholarGoogle Scholar | 34044161PubMed |

Moore, JM, Carvajal, JI, Rouse, GW, and Wilson, NG (2018). The Antarctic Circumpolar Current isolates and connects: structured circumpolarity in the sea star Glabraster antarctica. Ecology and Evolution 8, 10621–10633.
The Antarctic Circumpolar Current isolates and connects: structured circumpolarity in the sea star Glabraster antarctica.Crossref | GoogleScholarGoogle Scholar | 30464833PubMed |

Mora, C, Tittensor, DP, Adl, S, Simpson, AGB, and Worm, B (2011). How many species are there on earth and in the ocean? PLoS Biology 9, e1001127.
How many species are there on earth and in the ocean?Crossref | GoogleScholarGoogle Scholar | 21886479PubMed |

Munilla, T, and Membrives, AS (2009). Check-list of the pycnogonids from Antarctic and sub-Antarctic waters: zoogeographic implications. Antarctic Science 21, 99–111.
Check-list of the pycnogonids from Antarctic and sub-Antarctic waters: zoogeographic implications.Crossref | GoogleScholarGoogle Scholar |

Neusser, TP, Jörger, KM, and Schrödl, M (2011). Cryptic species in tropic sands – interactive 3D anatomy, molecular phylogeny and evolution of meiofaunal Pseudunelidae (Gastropoda, Acochlidia). PLoS One 6, e23313.
Cryptic species in tropic sands – interactive 3D anatomy, molecular phylogeny and evolution of meiofaunal Pseudunelidae (Gastropoda, Acochlidia).Crossref | GoogleScholarGoogle Scholar | 21912592PubMed |

Nguyen, L-T, Schmidt, HA, von Haeseler, A, and Minh, BQ (2015). IQ-TREE: a fast and effective stochastic algorithm for estimating maximum-likelihood phylogenies. Molecular Biology and Evolution 32, 268–274.
IQ-TREE: a fast and effective stochastic algorithm for estimating maximum-likelihood phylogenies.Crossref | GoogleScholarGoogle Scholar | 25371430PubMed |

Nikula, R, Fraser, CI, Spencer, HG, and Waters, JM (2010). Circumpolar dispersal by rafting in two subantarctic kelp-dwelling crustaceans. Marine Ecology Progress Series 405, 221–230.
Circumpolar dispersal by rafting in two subantarctic kelp-dwelling crustaceans.Crossref | GoogleScholarGoogle Scholar |

Nikula, R, Spencer, HG, and Waters, JM (2011a). Comparison of population-genetic structuring in congeneric kelp- versus rock-associated snails: a test of a dispersal-by-rafting hypothesis. Ecology and Evolution 1, 169–180.
Comparison of population-genetic structuring in congeneric kelp- versus rock-associated snails: a test of a dispersal-by-rafting hypothesis.Crossref | GoogleScholarGoogle Scholar | 22393493PubMed |

Nikula, R, Spencer, HG, and Waters, JM (2011b). Evolutionary consequences of microhabitat: population-genetic structuring in kelp- vs. rock-associated chitons. Molecular Ecology 20, 4915–4924.
Evolutionary consequences of microhabitat: population-genetic structuring in kelp- vs. rock-associated chitons.Crossref | GoogleScholarGoogle Scholar | 22026515PubMed |

Page, TJ, and Linse, K (2002). More evidence of speciation and dispersal across the Antarctic Polar Front through molecular systematics of Southern Ocean Limatula (Bivalvia: Limidae). Polar Biology 25, 818–826.
More evidence of speciation and dispersal across the Antarctic Polar Front through molecular systematics of Southern Ocean Limatula (Bivalvia: Limidae).Crossref | GoogleScholarGoogle Scholar |

Palomar, G, Pola, M, and Garcia-Vazquez, E (2014). First molecular phylogeny of the subfamily Polycerinae (Mollusca, Nudibranchia, Polyceridae). Helgoland Marine Research 68, 143–153.
First molecular phylogeny of the subfamily Polycerinae (Mollusca, Nudibranchia, Polyceridae).Crossref | GoogleScholarGoogle Scholar |

Palumbi, SR (1994). Genetic divergence, reproductive isolation, and marine speciation. Annual Review of Ecology and Systematics 25, 547–572.
Genetic divergence, reproductive isolation, and marine speciation.Crossref | GoogleScholarGoogle Scholar |

Paradis, E, and Schliep, K (2019). ape 5.0: an environment for modern phylogenetics and evolutionary analyses in R. Bioinformatics 35, 526–528.
ape 5.0: an environment for modern phylogenetics and evolutionary analyses in R.Crossref | GoogleScholarGoogle Scholar | 30016406PubMed |

Park, Y-H, Durand, I, Kestenare, E, Rougier, G, Zhou, M, d’Ovidio, F, Cotte, C, and Lee, J-H (2014). Polar Front around the Kerguelen Islands: an up-to-date determination and associated circulation of surface/subsurface waters.  Journal of Geophysical Research: Oceans 119, 6575–6592.
Polar Front around the Kerguelen Islands: an up-to-date determination and associated circulation of surface/subsurface waters. Crossref | GoogleScholarGoogle Scholar |

Pearse JS, Mooi R, Lockhart SJ, Brandt A (2009) Brooding and species diversity in the Southern Ocean: selection for brooders or speciation within brooding clades? In ‘Smithsonian at the Poles: Contributions to International Polar Year Science’. (Eds I Krupnik, MA Lang, SE Miller) pp. 181–196. (Smithsonian Institution Scholarly Press).
| Crossref |

Pertierra, LR, Santos-Martin, F, Hughes, KA, Avila, C, Caceres, JO, De Filippo, D, Gonzalez, S, Grant, SM, Lynch, H, Marina-Montes, C, Quesada, A, Tejedo, P, Tin, T, and Benayas, J (2021). Ecosystem services in Antarctica: global assessment of the current state, future challenges and managing opportunities. Ecosystem Services 49, 101299.
Ecosystem services in Antarctica: global assessment of the current state, future challenges and managing opportunities.Crossref | GoogleScholarGoogle Scholar |

Pfuhl, HA, and McCave, IN (2005). Evidence for late Oligocene establishment of the Antarctic Circumpolar Current. Earth and Planetary Science Letters 235, 715–728.
Evidence for late Oligocene establishment of the Antarctic Circumpolar Current.Crossref | GoogleScholarGoogle Scholar |

Pola, M, and Gosliner, TM (2010). The first molecular phylogeny of cladobranchian opisthobranchs (Mollusca, Gastropoda, Nudibranchia). Molecular Phylogenetics and Evolution 56, 931–941.
The first molecular phylogeny of cladobranchian opisthobranchs (Mollusca, Gastropoda, Nudibranchia).Crossref | GoogleScholarGoogle Scholar | 20460158PubMed |

Poulin, E, Palma, AT, and Féral, J-P (2002). Evolutionary versus ecological success in Antarctic benthic invertebrates. Trends in Ecology & Evolution 17, 218–222.
Evolutionary versus ecological success in Antarctic benthic invertebrates.Crossref | GoogleScholarGoogle Scholar |

Prates, L, and Perez, SI (2021). Late Pleistocene South American megafaunal extinctions associated with rise of Fishtail points and human population. Nature Communications 12, 2175.
Late Pleistocene South American megafaunal extinctions associated with rise of Fishtail points and human population.Crossref | GoogleScholarGoogle Scholar | 33846353PubMed |

Puillandre, N, Lambert, A, Brouillet, S, and Achaz, G (2012). ABGD, automatic barcode gap discovery for primary species delimitation. Molecular Ecology 21, 1864–1877.
ABGD, automatic barcode gap discovery for primary species delimitation.Crossref | GoogleScholarGoogle Scholar | 21883587PubMed |

Puillandre, N, Brouillet, S, and Achaz, G (2021). ASAP: assemble species by automatic partitioning. Molecular Ecology Resources 21, 609–620.
ASAP: assemble species by automatic partitioning.Crossref | GoogleScholarGoogle Scholar | 33058550PubMed |

Ribbink AJ (1984) Is the species flock concept tenable? In ‘Evolution of Fish Species Flocks’. (Eds AA Echelle, I Kornfield) pp. 21–25. (University of Maine at Orono Press: Orono, ME, USA)

Richards, ZT, Carvajal, JI, Wallace, CC, and Wilson, NG (2020). Phylotranscriptomics confirms Alveopora is sister to Montipora within the family Acroporidae. Marine Genomics 50, 100703.
Phylotranscriptomics confirms Alveopora is sister to Montipora within the family Acroporidae.Crossref | GoogleScholarGoogle Scholar | 31466869PubMed |

Rogers, AD (2007). Evolution and biodiversity of Antarctic organisms: a molecular perspective. Philosophical Transactions of the Royal Society of London – B. Biological Sciences 362, 2191–2214.
Evolution and biodiversity of Antarctic organisms: a molecular perspective.Crossref | GoogleScholarGoogle Scholar | 17553774PubMed |

Rogers-Bennett, L, and Catton, CA (2019). Marine heat wave and multiple stressors tip bull kelp forest to sea urchin barrens. Scientific Reports 9, 15050.
Marine heat wave and multiple stressors tip bull kelp forest to sea urchin barrens.Crossref | GoogleScholarGoogle Scholar | 31636286PubMed |

Ronce, O (2007). How does it feel to be like a rolling stone? Ten questions about dispersal evolution. Annual Review of Ecology, Evolution, and Systematics 38, 231–253.
How does it feel to be like a rolling stone? Ten questions about dispersal evolution.Crossref | GoogleScholarGoogle Scholar |

Schejter, L, Genzano, G, Gaitán, E, Perez, CD, and Bremec, CS (2020). Benthic communities in the Southwest Atlantic Ocean: conservation value of animal forests at the Burdwood Bank slope. Aquatic Conservation 30, 426–439.
Benthic communities in the Southwest Atlantic Ocean: conservation value of animal forests at the Burdwood Bank slope.Crossref | GoogleScholarGoogle Scholar |

Segovia, NI, Gallardo-Escárate, C, Poulin, E, and Haye, PA (2017). Lineage divergence, local adaptation across a biogeographic break, and artificial transport, shape the genetic structure in the ascidian Pyura chilensis. Scientific Reports 7, 44559.
Lineage divergence, local adaptation across a biogeographic break, and artificial transport, shape the genetic structure in the ascidian Pyura chilensis.Crossref | GoogleScholarGoogle Scholar | 28300177PubMed |

Shanks, AL, Grantham, BA, and Carr, MH (2003). Propagule dispersal distance and the size and spacing of marine reserves. Ecological Applications 13, 159–169.
Propagule dispersal distance and the size and spacing of marine reserves.Crossref | GoogleScholarGoogle Scholar |

Shields C (2009) Nudibranchs of the Ross Sea, Antarctica: phylogeny, diversity, and divergence. MSc Thesis, Clemson University, SC, USA. Available at https://tigerprints.clemson.edu/all_theses/637/

Singh, JS (2002). The biodiversity crisis: a multifaceted review. Current Science 82, 638–647.

Smith, JA, Hillenbrand, C-D, Pudsey, CJ, Allen, CS, and Graham, AGC (2010). The presence of polynyas in the Weddell Sea during the last glacial period with implications for the reconstruction of sea-ice limits and ice sheet history. Earth and Planetary Science Letters 296, 287–298.
The presence of polynyas in the Weddell Sea during the last glacial period with implications for the reconstruction of sea-ice limits and ice sheet history.Crossref | GoogleScholarGoogle Scholar |

Smith, DA, Hofmann, EE, Klinck, JM, and Lascara, CM (1999). Hydrography and circulation of the West Antarctic Peninsula Continental Shelf. Deep Sea Research – I. Oceanographic Research Papers 46, 925–949.
Hydrography and circulation of the West Antarctic Peninsula Continental Shelf.Crossref | GoogleScholarGoogle Scholar |

Sokolov, S, and Rintoul, SR (2009). Circumpolar structure and distribution of the Antarctic Circumpolar Current fronts: 1. Mean circumpolar paths. Journal of Geophysical Research: Oceans 114, C11018.
Circumpolar structure and distribution of the Antarctic Circumpolar Current fronts: 1. Mean circumpolar paths.Crossref | GoogleScholarGoogle Scholar |

Stark JS, Raymond T, Deppeler SL, Morrison AK (2019) Chapter 1 - Antarctic Seas. In ‘World Seas: An Environmental Evaluation’, 2nd edn. (Ed. C Sheppard) pp. 1–44, (Academic Press).
| Crossref |

Stephens, T (2018). The Antarctic Treaty system and the Anthropocene. The Polar Journal 8, 29–43.
The Antarctic Treaty system and the Anthropocene.Crossref | GoogleScholarGoogle Scholar |

Strugnell, JM, Rogers, AD, Prodöhl, PA, Collins, MA, and Allcock, AL (2008). The thermohaline expressway: the Southern Ocean as a centre of origin for deep-sea octopuses. Cladistics 24, 853–860.
The thermohaline expressway: the Southern Ocean as a centre of origin for deep-sea octopuses.Crossref | GoogleScholarGoogle Scholar | 34892888PubMed |

Taberlet, P, Coissac, E, Pompanon, F, Brochmann, C, and Willerslev, E (2012). Towards next-generation biodiversity assessment using DNA metabarcoding. Molecular Ecology 21, 2045–2050.
Towards next-generation biodiversity assessment using DNA metabarcoding.Crossref | GoogleScholarGoogle Scholar | 22486824PubMed |

Tamura, K, and Nei, M (1993). Estimation of the number of nucleotide substitutions in the control region of mitochondrial DNA in humans and chimpanzees. Molecular Biology and Evolution 10, 512–526.
Estimation of the number of nucleotide substitutions in the control region of mitochondrial DNA in humans and chimpanzees.Crossref | GoogleScholarGoogle Scholar | 8336541PubMed |

Templeton, AR, Crandall, KA, and Sing, CF (1992). A cladistic analysis of phenotypic associations with haplotypes inferred from restriction endonuclease mapping and DNA sequence data. III. Cladogram estimation. Genetics 132, 619–633.
A cladistic analysis of phenotypic associations with haplotypes inferred from restriction endonuclease mapping and DNA sequence data. III. Cladogram estimation.Crossref | GoogleScholarGoogle Scholar | 1385266PubMed |

Thatje, S, Hillenbrand, C-D, and Larter, R (2005). On the origin of Antarctic marine benthic community structure. Trends in Ecology & Evolution 20, 534–540.
On the origin of Antarctic marine benthic community structure.Crossref | GoogleScholarGoogle Scholar |

Thornhill, DJ, Mahon, AR, Norenburg, JL, and Halanych, KM (2008). Open-ocean barriers to dispersal: a test case with the Antarctic Polar Front and the ribbon worm Parborlasia corrugatus (Nemertea: Lineidae). Molecular Ecology 17, 5104–5117.
Open-ocean barriers to dispersal: a test case with the Antarctic Polar Front and the ribbon worm Parborlasia corrugatus (Nemertea: Lineidae).Crossref | GoogleScholarGoogle Scholar | 18992005PubMed |

Tin, T, Fleming, ZL, Hughes, KA, Ainley, DG, Convey, P, Moreno, CA, Pfeiffer, S, Scott, J, and Snape, I (2009). Impacts of local human activities on the Antarctic environment. Antarctic Science 21, 3–33.
Impacts of local human activities on the Antarctic environment.Crossref | GoogleScholarGoogle Scholar |

Torre, L, Tabares, PCC, Momo, F, Meyer, JFCA, and Sahade, R (2017). Climate change effects on Antarctic benthos: a spatially explicit model approach. Climatic Change 141, 733–746.
Climate change effects on Antarctic benthos: a spatially explicit model approach.Crossref | GoogleScholarGoogle Scholar |

Trifinopoulos, J, Nguyen, L-T, von Haeseler, A, and Minh, BQ (2016). W-IQ-TREE: a fast online phylogenetic tool for maximum likelihood analysis. Nucleic Acids Research 44, W232–W235.
W-IQ-TREE: a fast online phylogenetic tool for maximum likelihood analysis.Crossref | GoogleScholarGoogle Scholar | 27084950PubMed |

Tyagi, K, Kumar, V, Kundu, S, Pakrashi, A, Prasad, P, Caleb, JTD, and Chandra, K (2019). Identification of Indian spiders through DNA barcoding: cryptic species and species complex. Scientific Reports 9, 14033.
Identification of Indian spiders through DNA barcoding: cryptic species and species complex.Crossref | GoogleScholarGoogle Scholar | 31575965PubMed |

Valdes, A (2001). Deep-sea cryptobranch dorid nudibranchs (Mollusca, Opisthobranchia) from the tropical West Pacific, with descriptions of two new genera and eighteen new species. Malacologia 43, 237–311.

Vaughan, DG, Marshall, GJ, Connolley, WM, Parkinson, C, Mulvaney, R, Hodgson, DA, King, JC, Pudsey, CJ, and Turner, J (2003). Recent rapid regional climate warming on the Antarctic Peninsula. Climatic Change 60, 243–274.
Recent rapid regional climate warming on the Antarctic Peninsula.Crossref | GoogleScholarGoogle Scholar |

Wägele, H (1987). The distribution of some Antarctic nudibranchs (Opisthobranchia). The Journal of Molluscan Studies 53, 179–188.
The distribution of some Antarctic nudibranchs (Opisthobranchia).Crossref | GoogleScholarGoogle Scholar |

Wägele, H (1990). Revision of the genus Austrodoris Odhner, 1926 (Gastropoda, Opisthobranchia). The Journal of Molluscan Studies 56, 163–180.
Revision of the genus Austrodoris Odhner, 1926 (Gastropoda, Opisthobranchia).Crossref | GoogleScholarGoogle Scholar |

Wägele H, Klussmann-Kolb A, Vonnemann V, Medina M (2008) Heterobranchia I: the Opisthobranchia. In ‘Phylogeny and Evolution of the Mollusca’. (Ed. W Ponder) pp. 385–408. (California Scholarship Online).
| Crossref |

Wallace, RW (1997). Drugs from the sea: harvesting the results of aeons of chemical evolution. Molecular Medicine Today 3, 291–295.
Drugs from the sea: harvesting the results of aeons of chemical evolution.Crossref | GoogleScholarGoogle Scholar | 9257296PubMed |

Weersing, K, and Toonen, RJ (2009). Population genetics, larval dispersal, and connectivity in marine systems. Marine Ecology Progress Series 393, 1–12.
Population genetics, larval dispersal, and connectivity in marine systems.Crossref | GoogleScholarGoogle Scholar |

Wilson, NG, Schrödl, M, and Halanych, KM (2009). Ocean barriers and glaciation: evidence for explosive radiation of mitochondrial lineages in the Antarctic sea slug Doris kerguelenensis (Mollusca, Nudibranchia). Molecular Ecology 18, 965–984.
Ocean barriers and glaciation: evidence for explosive radiation of mitochondrial lineages in the Antarctic sea slug Doris kerguelenensis (Mollusca, Nudibranchia).Crossref | GoogleScholarGoogle Scholar | 19207248PubMed |

Wilson, NG, Maschek, JA, and Baker, BJ (2013). A species flock driven by predation? Secondary metabolites support diversification of slugs in Antarctica. PLoS One 8, e80277.
A species flock driven by predation? Secondary metabolites support diversification of slugs in Antarctica.Crossref | GoogleScholarGoogle Scholar | 24303002PubMed |

Wollscheid-Lengeling, E, Boore, J, Brown, W, and Wägele, H (2001). The phylogeny of Nudibranchia (Opisthobranchia, Gastropoda, Mollusca) reconstructed by three molecular markers. Organisms, Diversity & Evolution 1, 241–256.
The phylogeny of Nudibranchia (Opisthobranchia, Gastropoda, Mollusca) reconstructed by three molecular markers.Crossref | GoogleScholarGoogle Scholar |

Zhang, J, Kapli, P, Pavlidis, P, and Stamatakis, A (2013). A general species delimitation method with applications to phylogenetic placements. Bioinformatics 29, 2869–2876.
A general species delimitation method with applications to phylogenetic placements.Crossref | GoogleScholarGoogle Scholar | 23990417PubMed |