Genetic and morphological divergences in the cosmopolitan deep-sea amphipod Eurythenes gryllus reveal a diverse abyss and a bipolar species

Multi-ocean distribution of a brooding predator in the abyssal benthos

How far are species distributed on the abyssal plains? Spanning from 3000 to 6000 m below sea level, abyssal plains cover three-quarters of the ocean floor and are the largest but also least explored habitat on Earth. The question of vertical and horizontal distribution is central to understanding biogeographic and population genetic processes within species inhabiting the deep-sea benthos. Amphipod crustaceans are an important and dominant taxon in this ecosystem. As they are brooders, their dispersal capacities are more limited compared to species with free-swimming larvae, and with the exception of a few scavenging species deep-sea amphipods are restricted to a single ocean. Based on an integrative taxonomic approach (morphology, COI, 16S and 18S) we demonstrate the occurrence of a predatory amphipod species, Rhachotropis abyssalis, in three oceans: the Antarctic Ross Sea, the Northwest Pacific and the North Atlantic; regions more than 20,000 km apart. Although such extensive geographic distributions may represent a rare exception for brooding predators, these findings might also be no exception at all, but a reflection of the rare sampling and rare taxonomic investigation of invertebrate predators in the deep-sea. Our findings highlight our abysmal state of knowledge regarding biodiversity and biogeography on abyssal plains.

Morphology, nuclear SNPs and mate selection reveal that COI barcoding overestimates species diversity in a Mediterranean freshwater amphipod by an order of magnitude

DNA sequence information has revealed many morphologically cryptic species worldwide. For animals, DNA-based assessments of species diversity usually rely on the mitochondrial cytochrome c oxidase subunit I (COI) gene. However, a growing amount of evidence indicate that mitochondrial markers alone can lead to misleading species diversity estimates due to mito-nuclear discordance. Therefore, reports of putative species based solely on mitochondrial DNA should be verified by other methods, especially in cases where COI sequences are identical for different morphospecies or where divergence within the same morphospecies is high. Freshwater amphipods are particularly interesting in this context because numerous putative cryptic species have been reported. Here, we investigated the species status of the numerous mitochondrial molecular operational taxonomic units (MOTUs) found within Echinogammarus sicilianus. We used an integrative approach combining DNA barcoding with mate selection observations, detailed morphometrics and genome-wide double digest restriction site-associated DNA sequencing (ddRAD-seq). Within a relatively small sampling area, we detected twelve COI MOTUs (divergence = 1.8-20.3%), co-occurring in syntopy at two-thirds of the investigated sites. We found that pair formation was random and there was extensive nuclear gene flow among the ten MOTUs co-occurring within the same river stretch. The four most common MOTUs were also indistinguishable with respect to functional morphology. Therefore, the evidence best fits the hypothesis of a single, yet genetically diverse, species within the main river system. The only two MOTUs sampled outside the focal area were genetically distinct at the nuclear level and may represent distinct species. Our study reveals that COI-based species delimitation can significantly overestimate species diversity, highlighting the importance of integrative taxonomy for species validation, especially in hyperdiverse complexes with syntopically occurring mitochondrial MOTUs.

Large effective population size masks population genetic structure in Hirondellea amphipods within the deepest marine ecosystem, the Mariana Trench

The examination of genetic structure in the deep-ocean hadal zone has focused on divergence between tectonic trenches to understand how environment and geography may drive species divergence and promote endemism. There has been little attempt to examine localized genetic structure within trenches, partly because of logistical challenges associated with sampling at an appropriate scale, and the large effective population sizes of species that can be sampled adequately may mask underlying genetic structure. Here we examine genetic structure in the superabundant amphipod Hirondellea gigas in the Mariana Trench at depths of 8126-10,545 m. RAD sequencing was used to identify 3182 loci containing 43,408 single nucleotide polymorphisms (SNPs) across individuals after stringent pruning of loci to prevent paralogous multicopy genomic regions being erroneously merged. Principal components analysis of SNP genotypes resolved no genetic structure between sampling locations, consistent with a signature of panmixia. However, discriminant analysis of principal components identified divergence between all sites driven by 301 outlier SNPs in 169 loci and significantly associated with latitude and depth. Functional annotation of loci identified differences between singleton loci used in analysis and paralogous loci pruned from the data set and also between outlier and nonoutlier loci, all consistent with hypotheses explaining the role of transposable elements driving genome dynamics. This study challenges the traditional perspective that highly abundant amphipods within a trench form a single panmictic population. We discuss the findings in relation to eco-evolutionary and ontogenetic processes operating in the deep sea, and highlight key challenges associated with population genetic analysis in nonmodel systems with inherent large effective population sizes and genomes.

Barriers to gene flow in the deepest ocean ecosystems: Evidence from global population genomics of a cosmopolitan amphipod

The deepest marine ecosystem, the hadal zone, hosts endemic biodiversity resulting from geographic isolation and environmental selection pressures. However, the pan-ocean distribution of some fauna challenges the concept that the hadal zone is a series of isolated island-like habitats. Whether this remains true at the population genomic level is untested. We investigated phylogeographic patterns of the amphipod, Bathycallisoma schellenbergi, from 12 hadal features across the Pacific, Atlantic, Indian, and Southern oceans and analyzed genome-wide single-nucleotide polymorphism markers and two mitochondrial regions. Despite a cosmopolitan distribution, populations were highly restricted to individual features with only limited gene flow between topographically connected features. This lack of connectivity suggests that populations are on separate evolutionary trajectories, with evidence of potential cryptic speciation at the Atacama Trench. Together, this global study demonstrates that the shallower ocean floor separating hadal features poses strong barriers to dispersal, driving genetic isolation and creating pockets of diversity to conserve.

Trace Elemental Analysis of the Exoskeleton, Leg Muscle, and Gut of Three Hadal Amphipods

Hadal trenches are the deepest areas worldwide. Amphipods are considered a key factor in hadal ecosystems because of their important impacts on the hadal environment. Amphipods have benthic habits, and therefore, serve as good metal biomonitors. However, little is known about the hadal amphipod metal accumulations. In the present study, Alicella gigantea, Hirondellea gigas, and Scopelocheirus schellenbergi were sampled from the New Britain Trench (8824m, 7.02S 149.16E), Mariana Trench (10,839m, 11.38N 142.42E), and Marceau Trench (6690m, 1.42N 148.74E) in the West Pacific Ocean, respectively. The elemental concentrations of the three hadal amphipods were subsequently investigated. Nine trace elements (V, Cr, Mn, Co, Ni, Se, Mo, Ag, and Cd) of three tissues (exoskeleton, leg muscle, and gut) of the hadal amphipods were detected by using inductively coupled plasma mass spectrometry (ICP-MS) method. The concentrations of Cr, Cd, and Mn were comparably higher among those nine examined elements. The greatest accumulations of the elements Cr, Ag, and V in the exoskeleton and leg muscle were observed in H. gigas, and elements Mn, Co, and Se showed the highest accumulations in the gut in H. gigas among the three hadal amphipods. In addition, comparisons of the leg muscle trace element accumulation between the hadal amphipods and non-abyssal and shallow water decapoda and amphipoda species showed that the hadal amphipods possessed comparably higher concentrations of the trace elements Cd, Co, Mo, Ag, and V. This finding suggested a bottom-up effect of food availability and indicated the effects of human activities within the hadal environments. This study reveals the trace element bio-accumulation of three hadal amphipods, and suggests that deep-sea amphipods are potential indicator species for trace element bioavailability in the deep-sea environment.

Macroecology of Southern Ocean benthic Ostracoda (Crustacea) from the continental margin and abyss

The Southern Ocean (SO) is the least understood environment on earth, but anthropogenic impacts related to climate change, fishing activities and tourism are already well-established in the region. Herein, we investigate biodiversity patterns in the Atlantic Sector of the SO by investigating a considerable number of samples collected from a wide depth range (89–6224 m). Three topics are addressed: (1) which taxa occur in the SO; (2) if there are distinct shelf, slope and abyssal faunas; and (3) which abiotic parameters explain the ostracod assemblages. We also tested putative ecological preferences of ostracod genera, widely used in palaeoceanographical reconstruction. For this, almost 5000 living ostracods from 51 epibenthic sledge samples were identified to 140 species. Statistically significant differences were observed among the faunas from the shelf, slope and abyssal faunas, indicating that the shelf fauna survived in refugia during the glacial periods of the Quaternary. Our analyses showed that the main variables influencing ostracod distribution is depth, followed by nitrate and phosphate. Temperature, salinity and oxygen content do not seem to be as important as previously proposed. Putative ecological preferences of ostracod genera, for example, low oxygen for Cytherella, could not be confirmed, suggesting that this topic requires further investigation.

Recent speciation and hybridization in Icelandic deep-sea isopods: An integrative approach using genomics and proteomics

The crustacean marine isopod species Haploniscus bicuspis (Sars, 1877) shows circum-Icelandic distribution in a wide range of environmental conditions and along well-known geographic barriers, such as the Greenland-Iceland-Faroe (GIF) Ridge. We wanted to explore population genetics, phylogeography and cryptic speciation as well as investigate whether previously described, but unaccepted subspecies have any merit. Using the same set of specimens, we combined mitochondrial COI sequences, thousands of nuclear loci (ddRAD), and proteomic profiles, plus selected morphological characters using confocal laser scanning microscopy (CLSM). Five divergent genetic lineages were identified by COI and ddRAD, two south and three north of the GIF Ridge. Assignment of populations to the three northern lineages varied and detailed analyses revealed hybridization and gene flow between them, suggesting a single northern species with a complex phylogeographic history. No apparent hybridization was observed among lineages south of the GIF Ridge, inferring the existence of two more species. Differences in proteomic profiles between the three putative species were minimal, implying an ongoing or recent speciation process. Population differentiation was high, even among closely associated populations, and higher in mitochondrial COI than nuclear ddRAD loci. Gene flow is apparently male-biased, leading to hybrid zones and instances of complete exchange of the local nuclear genome through immigrating males. This study did not confirm the existence of subspecies defined by male characters, which probably instead refer to different male developmental stages.

Eurythenes atacamensis sp. nov. (Crustacea: Amphipoda) exhibits ontogenetic vertical stratification across abyssal and hadal depths in the Atacama Trench, eastern South Pacific Ocean

Eurythenes S.I. Smith in Scudder, 1882 (Crustacea: Amphipoda) are prevalent scavengers of the benthopelagic community from bathyal to hadal depths. While a well-studied genus, molecular systematic studies have uncovered cryptic speciation and multiple undescribed lineages. Here, we apply an integrative taxonomic approach and describe the tenth species, Eurythenes atacamensis sp. nov., based on specimens from the 2018 Atacamex and RV Sonne SO261 Expeditions to the southern sector of the Peru-Chile Trench, the Atacama Trench (24-⁠21°S). Eurythenes atacamensis sp. nov. is a large species, max. observed length 83.2 mm, possesses diagnostic features, including a short gnathopod 1 palm and a chelate gnathopod 2 palm, and a distinct genetic lineage based on a 16S rRNA and COI phylogeny. This species is a dominant bait-attending fauna with an extensive bathymetric range, spanning from 4974 to 8081 m. The RV Sonne SO261 specimens were recovered along a 10-station transect from abyssal to hadal depths and further examined for demographic and bathymetric-related patterns. Ontogenetic vertical stratification was evident across the trench axis, with only juveniles present at abyssal depths (4974-6025 m). Total length-depth analysis revealed that the size of females was unrelated to depth, whereas juveniles followed a sigmoidal relationship with a step-up in size at depths >7200 m. Thus, these bathymetric trends suggest that juveniles and females employ differing ecological strategies in subduction trench environments. This study highlights that even dominant and ecologically important species are still being discovered within the abyssal and hadal environments. Continued systematic expeditions will lead to an improved understanding of the eco-evolutionary drivers of speciation in the world's largest ecosystem.

Phylogeography and genetic diversity of the microbivalve Kidderia subquadrata, reveals new data from West Antarctic Peninsula

It is well established that Antarctic biodiversity has been strongly influenced by rapid climatic fluctuations during the Quaternary. Marine invertebrates from Antarctica constitute an interesting lens through which to study the impacts of the last glacial periods as glaciation impacted the distribution and intraspecific genetic variation of these animals. However, the impact on the spatial genetic distribution and historical demography of local processes in areas adjacent to the West Antarctic Peninsula (WAP) is less clear. Here we present new genetic information on the bivalve Kidderia subquadrata, a small mollusk that inhabits intertidal rocky island ecosystems throughout the WAP. Using a phylogeographical approach, we examined the spatial patterns of genetic diversity in this brooder species to test the hypothesis of strong genetic structure in incubating organisms and the hypothesis of glacial refugia in organisms with limited dispersion. We found evidence of strong genetic structure among populations of the WAP and a recent expansion in the South Shetland Islands. Our findings are concordant with the predictions that incubating organisms, abundant in Antarctica, present a strong genetic structure among their populations and also support the hypothesis of glacial refugia in organisms with limited dispersion. The effect of the coastal current pattern in the WAP is suggested as a driver to the local spatial dynamics of the genetic diversity distribution. Although genetic information about this microbivalve is still scarce, the knowledge reported here has increased our understanding of the evolutionary patterns of this organism that is endemic to the Southern Ocean.

Species boundaries and phylogeographic patterns in new species of Nannoniscus (Janiroidea: Nannoniscidae) from the equatorial Pacific nodule province inferred from mtDNA and morphology

Spatial patterns of genetic variation (based on COI and 16S mtDNA) for morphologically similar species in the isopod genus Nannoniscus G.O. Sars. 1870 were examined that occur broadly across the Clarion Clipperton Fracture Zone (CCZ). Samples were obtained from five different licence areas as well as an Area of Particular Environmental Interest (APEI-6) with sites located at various distances (a few to several hundred kilometres) from one another. Applying three different species delimitation (SD) methods (sGMYC, mPTP and ABGD) of the molecular data, we could distinguish between four and 12 different molecular taxonomic operational units (MOTUs). Morphological analyses could confirm five distinct phenotypic clades that represent species new to science and are described here: Nannoniscus brenkei sp. nov., Nannoniscus hilario sp. nov., Nannoniscus magdae sp. nov., Nannoniscus menoti sp. nov. and Nannoniscus pedro sp. nov. Despite the assumed limited dispersal capacity of Nannoniscus species, we found haplotypes of two species to be geographically widespread (up to > 1400 km apart), as opposed to several divergent clades occurring in close vicinity or even sympatry. Geographic distance appeared to explain the phylogeographic structure of Nannoniscus species to some extent, although oceanographic features and level of environmental heterogeneity were probably equally important.

Integrative taxonomy of giant crested Eusirus in the Southern Ocean, including the description of a new species (Crustacea: Amphipoda: Eusiridae)

Among Antarctic amphipods of the genus Eusirus, a highly distinctive clade of giant species is characterized by a dorsal, blade-shaped tooth on pereionites 5–7 and pleonites 1–3. This lineage, herein named ‘crested Eusirus’, includes two potential species complexes, the Eusirus perdentatus and Eusirus giganteus complexes, in addition to the more distinctive Eusirus propeperdentatus. Molecular phylogenies and statistical parsimony networks (COI, CytB and ITS2) of crested Eusirus are herein reconstructed. This study aims to formally revise species diversity within crested Eusirus by applying several species delimitation methods (Bayesian implementation of the Poisson tree processes model, general mixed Yule coalescent, multi-rate Poisson tree processes and automatic barcode gap discovery) on the resulting phylogenies. In addition, results from the DNA-based methods are benchmarked against a detailed morphological analysis of all available specimens of the E. perdentatus complex. Our results indicate that species diversity of crested Eusirus is underestimated. Overall, DNA-based methods suggest that the E. perdentatus complex is composed of three putative species and that the E. giganteus complex includes four or five putative species. The morphological analysis of available specimens from the E. perdentatus complex corroborates molecular results by identifying two differentiable species, the genuine E. perdentatus and a new species, herein described as Eusirus pontomedon sp. nov.

Continental-scale patterns of hyper-cryptic diversity within the freshwater model taxon Gammarus fossarum (Crustacea, Amphipoda)

Traditional morphological diagnoses of taxonomic status remain widely used while an increasing number of studies show that one morphospecies might hide cryptic diversity, i.e. lineages with unexpectedly high molecular divergence. This hidden diversity can reach even tens of lineages, i.e. hyper cryptic diversity. Even well-studied model-organisms may exhibit overlooked cryptic diversity. Such is the case of the freshwater crustacean amphipod model taxon Gammarus fossarum. It is extensively used in both applied and basic types of research, including biodiversity assessments, ecotoxicology and evolutionary ecology. Based on COI barcodes of 4926 individuals from 498 sampling sites in 19 European countries, the present paper shows (1) hyper cryptic diversity, ranging from 84 to 152 Molecular Operational Taxonomic Units, (2) ancient diversification starting already 26 Mya in the Oligocene, and (3) high level of lineage syntopy. Even if hyper cryptic diversity was already documented in G. fossarum, the present study increases its extent fourfold, providing a first continental-scale insight into its geographical distribution and establishes several diversification hotspots, notably south-eastern and central Europe. The challenges of recording hyper cryptic diversity in the future are also discussed.

Just the once will not hurt: DNA suggests species lumping over two oceans in deep-sea snails (Cryptogemma)

The practice of species delimitation using molecular data commonly leads to the revealing of species complexes and an increase in the number of delimited species. In a few instances, however, DNA-based taxonomy has led to lumping together of previously described species. Here, we delimit species in the genus Cryptogemma (Gastropoda: Conoidea: Turridae), a group of deep-sea snails with a wide geographical distribution, primarily by using the mitochondrial COI gene. Three approaches of species delimitation (ABGD, mPTP and GMYC) were applied to define species partitions. All approaches resulted in eight species. According to previous taxonomic studies and shell morphology, 23 available names potentially apply to the eight Cryptogemma species that were recognized herein. Shell morphometrics, radular characters and geographical and bathymetric distributions were used to link type specimens to these delimited species. In all, 23 of these available names are here attributed to seven species, resulting in 16 synonymizations, and one species is described as new: Cryptogemma powelli sp. nov. We discuss the possible reasons underlying the apparent overdescription of species within Cryptogemma, which is shown here to constitute a rare case of DNA-based species lumping in the hyper-diversified superfamily Conoidea.

The complete mitochondrial genome of the deep-sea amphipod Eurythenes magellanicus (Crustacea: Amphipoda: Lysianassidae)

The complete mitochondrial genome of the deep sea amphipod Eurythenes magellanicus was determined in this paper. This molecular was 14,988 bp in length, and contained the typical 13 protein coding genes (PCGs), 22 transfer RNAs (tRNAs), two ribosomal RNAs (rRNAs) and one control region (CR). The gene order of E. magellanicus was identical to that from E. maldoror, a deep sea amphipod inhabiting in a deeper habitat than E. magellanicus. A maximum-likelihood tree based on the 13 PCGs from 25 amphipods indicated that E. magellanicus and E. maldoror were closely related and the origin of deep sea amphipods was not monophyletic.

Integrative systematics and ecology of a new deep-sea family of tanaidacean crustaceans

A new family of paratanaoidean Tanaidacea – Paranarthrurellidae fam. nov. – is erected to accommodate two genera without family classification (Paratanaoidea incertae sedis), namely Armatognathia Kudinova-Pasternak, 1987 and Paranarthrurella Lang, 1971. Seven new species of Paranarthrurella and two of Armatognathia are described from material taken in different deep-sea areas of the Atlantic and Pacific oceans. The type species of Paranarthrurella — P. caudata (Kudinova-Pasternak, 1965) — is redescribed based on the paratype. The genus Cheliasetosatanais Larsen and Araújo-Silva, 2014 originally classified within Colletteidae is synonymised with Paranarthrurella, and Arthrura shiinoi Kudinova-Pasternak, 1973 is transferred to Armatognathia. Amended diagnoses of Armatognathia and Paranarthrurella genera are given. Choosing characters for distinguishing and defining both genera was supported by Principal Component Analysis. Designation of the new family is supported by molecular phylogenetic analysis of COI and 18S datasets. The distribution of all species currently included in the new family was visualised and their bathymetric distribution analysed.

Deep ocean seascape and Pseudotanaidae (Crustacea: Tanaidacea) diversity at the Clarion-Clipperton Fracture Zone

Understanding the diversity and spatial distribution of benthic species is fundamental to properly assess the impact of deep sea mining. Tanaidacea provide an exceptional opportunity for assessing spatial patterns in the deep-sea, given their low mobility and limited dispersal potential. The diversity and distribution of pseudotanaid species is characterized here for the Clarion and Clipperton Fractures Zone (CCZ), which is the most extensive deposit field of metallic nodules. Samples were taken from the Belgian, German and French license areas, but also from the APEI 3 (Area of Particular Environmental Interest 3) of the Interoceanmetal consortium associates. The combination of morphological and genetic data uncovered one new pseudotanaid genus (Beksitanais n. gen.) and 14 new species of Pseudotanais (2 of them virtual taxa). Moreover, our results suggest that spatial structuring of pseudotanaid diversity is correlated with deep-sea features, particularly the presence of fractures and seamount chains crossing the CCZ. The presence of geographical barriers delimiting species distributions has important implications for the establishment of protected areas, and the APEI3 protected area contains only one third of the total pseudotanaid species in CCZ. The specimen collection studied here is extremely valuable and represents an important first step in characterizing the diversity and distribution of pseudotanaids within the Tropical Eastern Pacific.

The complete mitochondrial genome of the largest amphipod, Alicella gigantea: Insight into its phylogenetic relationships and deep sea adaptive characters

Alicella gigantea (Alicelloidae) is a scavenger with the largest body size among amphipods. It is a participant in the foodweb of deepsea ecosystem and distributed with vast bathymetric and geographic ranges. In this study, the mitochondrial genome of A. gigantea was completely assembled and characterized. The complete sequence has a total length of 16,851 bp, comprising the usual eukaryotic components, with 13 protein-coding genes (PCGs), 2 ribosomal RNA genes (rRNAs), 22 transfer RNA genes (tRNAs), and 2 noncoding control regions (CRs). The gene rearrangement and reverse nucleotide strand bias of its mitochondrial genome are similar to those observed in the deepsea amphipod Eurythenes maldoror (Eurytheneidae), but different from the characters of Halice sp. MT-2017 (Dexaminoidea), an inhabitant of a deeper environment. Phylogenetic analysis indicates that A. gigantea occupies the basal branch of deepsea species-E. maldoror and Hirondellea gigas. This phylogeny supports the hypothesis that the evolution of hadal amphipods has undergone a transition from the abyssal depth. Compared to 41 available shallow water equivalents, the four accessible mitochondrial genomes from the deep sea, including the one produced in this study, show significantly fewer charged amino acids in the 13 PCGs, which suggests an adaption to the deepsea environment.

Genetic Diversity and Population Structure Analysis of Three Deep-Sea Amphipod Species from Geographically Isolated Hadal Trenches in the Pacific Ocean

Amphipods of the superfamily Lysianassoidea that inhabit the hadal zone ( > 6000 m) have large bathymetric ranges and play a key role in deep ocean ecosystems. The endemism of these amphipod species makes them a good model for investigating potent natural selection and restricted dispersal in deep ocean trenches. Here, we describe genetic diversity and intraspecific population differentiation among three amphipod species from four Pacific trenches based on a mtDNA concatenated dataset (CO Ι and 16S rRNA genes) from 150 amphipod individuals. All amphipod populations had low genetic diversity, as indicated by haplotype and nucleotide diversity values. Population geographic relationship analysis of two Alicella gigantea populations revealed no genetic differentiation between these two localities (pairwise genetic differentiation coefficient = 0.00032, gene flow = 784.58), and the major variation (99.97%) was derived from variation within the populations. Historical demographic events were investigated using Tajima's D and Fu's F neutrality tests and analysis of mismatch distribution. Consistent results provided strong evidence to support the premise that demographic expansion occurred only for the Mariana population of Hirondellea gigas, possibly within the last 2.1-3.4 million years. These findings suggest that the formation of amphipod population structure might be the result of multiple factors including high hydrostatic pressure, food distribution, trench topographic forcing and potential ecological interactions.

The use of spatially explicit genetic variation data from four deep-sea sponges to inform the protection of Vulnerable Marine Ecosystems

The United Nations General Assembly has called for greater protection of the world’s deep-sea species and of features such as Vulnerable Marine Ecosystems (VMEs). Sponges are important components of VMEs and information about their spatially explicit genetic diversity can inform management decisions concerning the placement of protected areas. We employed a spatially explicit hierarchical testing framework to examine genetic variation amongst archived samples of four deep-sea sponges in the New Zealand region. For Poecillastra laminaris Sollas 1886, significant mitochondrial (COI, Cytb) and nuclear DNA (microsatellite) genetic differences were observed between provinces, amongst north-central-south regions and amongst geomorphic features. For Penares sp. no significant structure was detected (COI, 12S) across the same areas. For both Neoaulaxinia persicum Kelly, 2007 (COI, 12S) and Pleroma menoui Lévi & Lévi 1983 (COI) there was no evidence of genetic differentiation within their northern only regional distributions. Of 10 separate species-by-marker tests for isolation-by-distance and isolation-by-depth, only the isolation-by-depth test for N. persicum for COI was significant. The use of archived samples highlights how historical material may be used to support national and international management decisions. The results are discussed in the broader context of existing marine protected areas, and possible future design of spatial management measures for protecting VMEs in the New Zealand region.

Mapping knowledge gaps in marine diversity reveals a latitudinal gradient of missing species richness

A reliable description of any spatial pattern in species richness requires accurate knowledge about species geographical distribution. However, sampling bias may generate artefactual absences within species range and compromise our capacity to describe biodiversity patterns. Here, we analysed the spatial distribution of 35,000 marine species (varying from copepods to sharks) to identify missing occurrences (gaps) across their latitudinal range. We find a latitudinal gradient of species absence peaking near the equator, a pattern observed in both shallow and deep waters. The tropical gap in species distribution seems a consequence of reduced sampling effort at low latitudes. Overall, our results suggest that spatial gaps in species distribution are the main cause of the bimodal pattern of marine diversity. Therefore, only increasing sampling effort at low latitudes will reveal if the absence of species in the tropics, and the consequent dip in species richness, are artefacts of sampling bias or a natural phenomenon.

Accurate understanding of species biogeographic patterns is contingent upon adequate sampling effort across space. Here, the authors analyse the distribution records for 35,000 marine species, highlighting data gaps caused by undersampling in the tropics.

A new predator connecting the abyssal with the hadal in the Kuril-Kamchatka Trench, NW Pacific

The bathyal to hadal deep sea of north-west Pacific Ocean was recently intensively sampled during four international expeditions (KuramBio I and II, SoJaBio and SokhoBio). A large amphipod, Rhachotropis saskia n. sp., was sampled in the Kuril-Kamchatka Trench and increases the number of described hadal species of that area to eight. A detailed description of the new species is provided, including illustrations, scanning-microscope images and molecular analysis. This predatory species was sampled at both continental and ocean abyssal margins of the Kuril-Kamchatka Trench as well as at hadal depths of the trench. The wide bathymetric distribution of the new species over more than 3,000 m is confirmed by molecular analysis, indicating that the Kuril Kamchatka Trench is not a distribution barrier for this species. However, the molecular analysis indicated the presence of isolation by distance of the populations of the studied taxon.

Species identification and connectivity of marine amphipods in Canada's three oceans

Monitoring the distribution of marine biodiversity is a crucial step to better assess the impacts of global changes. Arctic marine fauna is dominated by amphipods in terms of abundance and biomass. These peracarids are an important marine order of crustaceans but the number of species found in the different Canadian oceans is currently unknown. Furthermore, most species are difficult to identify due to poor taxonomic descriptions and morphological convergence. The aim of this study was to assess the species diversity of marine amphipods in the three Canadian oceans using DNA barcoding. To do so, we produced a database of DNA barcodes of amphipods from the three Canadian Oceans publicly available from the BOLD website to which we added 310 new sequences from the Canadian Arctic Archipelago. We first delimited amphipod species based on barcode gap detection techniques and tree based method (bPTP) and then compared the composition of amphipods among the three oceans in order to assess the influence of past transarctic exchanges on Arctic diversity. Our analysis of 2309 sequences which represent more than 250 provisional species revealed a high connectivity between the Atlantic and Arctic Oceans. Our results also suggest that a single threshold to delimitate species is not suitable for amphipods. This study highlights the challenges involved in species delimitation and the need to obtain complete barcoding inventories in marine invertebrates.

Cryptic species in a well-known habitat: applying taxonomics to the amphipod genus Epimeria (Crustacea, Peracarida)

Taxonomy plays a central role in biological sciences. It provides a communication system for scientists as it aims to enable correct identification of the studied organisms. As a consequence, species descriptions should seek to include as much available information as possible at species level to follow an integrative concept of 'taxonomics'. Here, we describe the cryptic species Epimeria frankei sp. nov. from the North Sea, and also redescribe its sister species, Epimeria cornigera. The morphological information obtained is substantiated by DNA barcodes and complete nuclear 18S rRNA gene sequences. In addition, we provide, for the first time, full mitochondrial genome data as part of a metazoan species description for a holotype, as well as the neotype. This study represents the first successful implementation of the recently proposed concept of taxonomics, using data from high-throughput technologies for integrative taxonomic studies, allowing the highest level of confidence for both biodiversity and ecological research.

Polysaccharide hydrolase of the hadal zone amphipods Hirondellea gigas

Hirondellea species are common inhabitants in the hadal region deeper than 7,000 m. We found that Hirondellea gigas thrived in the Challenger Deep possessed polysaccharide hydrolases as digestive enzymes. To obtain various enzymes of other H. gigas, we captured amphipods from the Japan Trench, and Izu-Ogasawara (Bonin) Trench. A phylogenetic analysis based on the cytochrome oxidase I gene showed close relationships among amphipods, despite the geographic distance between the localities. However, several differences in enzymatic properties were observed in these H. gigas specimens. We also carried out RNA sequencing of H. gigas from the Izu-Ogasawara Trench. The cellulase gene of H. gigas was highly homologous to cellobiohydrolase of Glucosyl Hydrolase family 7 (GH7). On the other hand, enzymatic properties of H. gigas's cellulase were different from those of typical GH7 cellobiohydrolase. Thus, these results indicate that hadal-zone amphipod can be good candidates as the new enzyme resource.

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

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

Amphipod family distributions around Iceland

Amphipod crustaceans were collected at all 55 stations sampled with an epibenthic sledge during two IceAGE expeditions (Icelandic marine Animals: Genetics and Ecology) in 2011 and 2013. In total, 34 amphipod families and three superfamilies were recorded in the samples. Distribution maps are presented for each taxon along with a summary of the regional taxonomy for the group. Statistical analyses based on presence/absence data revealed a pattern of family distributions that correlated with sampling depth. Clustering according to the geographic location of the stations (northernmost North Atlantic Sea and Arctic Ocean) can also be observed. IceAGE data for the Amphilochidae and Oedicerotidae were analysed on species level; in case of the Amphilochidae they were compared to the findings from a previous Icelandic benthic survey, BIOICE (Benthic Invertebrates of Icelandic waters), which also identified a high abundance of amphipod fauna.

A genetic fingerprint of Amphipoda from Icelandic waters - the baseline for further biodiversity and biogeography studies

Amphipods constitute an abundant part of Icelandic deep-sea zoobenthos yet knowledge of the diversity of this fauna, particularly at the molecular level, is scarce. The present work aims to use molecular methods to investigate genetic variation of the Amphipoda sampled during two IceAGE collecting expeditions. The mitochondrial cytochrome oxidase subunit 1 (COI) of 167 individuals originally assigned to 75 morphospecies was analysed. These targeted morhospecies were readily identifiable by experts using light microscopy and representative of families where there is current ongoing taxonomic research. The study resulted in 81 Barcode Identity Numbers (BINs) (of which >90% were published for the first time), while Automatic Barcode Gap Discovery revealed the existence of 78 to 83 Molecular Operational Taxonomic Units (MOTUs). Six nominal species (Rhachotropis helleri, Arrhis phyllonyx, Deflexilodes tenuirostratus, Paroediceros propinquus, Metopa boeckii, Astyra abyssi) appeared to have a molecular variation higher than the 0.03 threshold of both p-distance and K2P usually used for amphipod species delineation. Conversely, two Oedicerotidae regarded as separate morphospecies clustered together with divergences in the order of intraspecific variation. The incongruence between the BINs associated with presently identified species and the publicly available data of the same taxa was observed in case of Paramphithoe hystrix and Amphilochus manudens. The findings from this research project highlight the necessity of supporting molecular studies with thorough morphology species analyses.

Invertebrate population genetics across Earth's largest habitat: The deep-sea floor

Despite the deep sea being the largest habitat on Earth, there are just 77 population genetic studies of invertebrates (115 species) inhabiting non-chemosynthetic ecosystems on the deep-sea floor (below 200 m depth). We review and synthesize the results of these papers. Studies reveal levels of genetic diversity comparable to shallow-water species. Generally, populations at similar depths were well connected over 100s-1,000s km, but studies that sampled across depth ranges reveal population structure at much smaller scales (100s-1,000s m) consistent with isolation by adaptation across environmental gradients, or the existence of physical barriers to connectivity with depth. Few studies were ocean-wide (under 4%), and 48% were Atlantic-focused. There is strong emphasis on megafauna and commercial species with research into meiofauna, "ecosystem engineers" and other ecologically important species lacking. Only nine papers account for ~50% of the planet's surface (depths below 3,500 m). Just two species were studied below 5,000 m, a quarter of Earth's seafloor. Most studies used single-locus mitochondrial genes revealing a common pattern of non-neutrality, consistent with demographic instability or selective sweeps; similar to deep-sea hydrothermal vent fauna. The absence of a clear difference between vent and non-vent could signify that demographic instability is common in the deep sea, or that selective sweeps render single-locus mitochondrial studies demographically uninformative. The number of population genetics studies to date is miniscule in relation to the size of the deep sea. The paucity of studies constrains meta-analyses where broad inferences about deep-sea ecology could be made.

Distinct genetic differentiation and species diversification within two marine nematodes with different habitat preference in Antarctic sediments

Background

Dispersal ability, population genetic structure and species divergence in marine nematodes are still poorly understood, especially in remote areas such as the Southern Ocean. We investigated genetic differentiation of species and populations of the free-living endobenthic nematode genera Sabatieria and Desmodora using nuclear 18S rDNA, internal transcribed spacer (ITS) rDNA, and mitochondrial cytochrome oxidase I (COI) gene sequences. Specimens were collected at continental shelf depths (200–500 m) near the Antarctic Peninsula, Scotia Arc and eastern side of the Weddell Sea. The two nematode genera co-occurred at all sampled locations, but with different vertical distribution in the sediment. A combination of phylogenetic (GMYC, Bayesian Inference, Maximum Likelihood) and population genetic (AMOVA) analyses were used for species delimitation and assessment of gene flow between sampling locations.

Results

Sequence analyses resulted in the delimitation of four divergent species lineages in Sabatieria, two of which could not be discriminated morphologically and most likely constitute cryptic species. Two species were recognised in Desmodora, one of which showed large intraspecific morphological variation. Both genera comprised species that were restricted to one side of the Weddell Sea and species that were widely spread across it. Population genetic structuring was highly significant and more pronounced in the deeper sediment-dwelling Sabatieria species, which are generally less prone to resuspension and passive dispersal in the water column than surface Desmodora species.

Conclusions

Our results indicate that gene flow is restricted at large geographic distance in the Southern Ocean, which casts doubt on the efficiency of the Weddell gyre and Antarctic Circumpolar Current in facilitating circum-Antarctic nematode species distributions. We also show that genetic structuring and cryptic speciation can be very different in nematode species isolated from the same geographic area, but with different habitat preferences (surface versus deeper sediment layers).

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-017-0968-1) contains supplementary material, which is available to authorized users.

Pelagic and benthic communities of the Antarctic ecosystem of Potter Cove: Genomics and ecological implications

Molecular technologies are more frequently applied in Antarctic ecosystem research and the growing amount of sequence-based information available in databases adds a new dimension to understanding the response of Antarctic organisms and communities to environmental change. We apply molecular techniques, including fingerprinting, and amplicon and metagenome sequencing, to understand biodiversity and phylogeography to resolve adaptive processes in an Antarctic coastal ecosystem from microbial to macrobenthic organisms and communities. Interpretation of the molecular data is not only achieved by their combination with classical methods (pigment analyses or microscopy), but furthermore by combining molecular with environmental data (e.g., sediment characteristics, biogeochemistry or oceanography) in space and over time. The studies form part of a long-term ecosystem investigation in Potter Cove on King-George Island, Antarctica, in which we follow the effects of rapid retreat of the local glacier on the cove ecosystem. We formulate and encourage new approaches to integrate molecular tools into Antarctic ecosystem research, environmental conservation actions, and polar ocean observatories.

Sensitivity of the deep-sea amphipod Eurythenes gryllus to chemically dispersed oil

In the context of an oil spill accident and the following oil spill response, much attention is given to the use of dispersants. Dispersants are used to disperse an oil slick from the sea surface into the water column generating a cloud of dispersed oil droplets. The main consequence is an increasing of the sea water-oil interface which induces an increase of the oil biodegradation. Hence, the use of dispersants can be effective in preventing oiling of sensitive coastal environments. Also, in case of an oil blowout from the seabed, subsea injection of dispersants may offer some benefits compared to containment and recovery of the oil or in situ burning operation at the sea surface. However, biological effects of dispersed oil are poorly understood for deep-sea species. Most effects studies on dispersed oil and also other oil-related compounds have been focusing on more shallow water species. This is the first approach to assess the sensitivity of a macro-benthic deep-sea organism to dispersed oil. This paper describes a toxicity test which was performed on the macro-benthic deep-sea amphipod (Eurythenes gryllus) to determine the concentration causing lethality to 50% of test individuals (LC50) after an exposure to dispersed Brut Arabian Light (BAL) oil. The LC50 (24 h) was 101 and 24 mg L(-1) after 72 h and 12 mg L(-1) at 96 h. Based on EPA scale of toxicity categories to aquatic organisms, an LC50 (96 h) of 12 mg L(-1) indicates that the dispersed oil was slightly to moderately toxic to E. gryllus. As an attempt to compare our results to others, a literature study was performed. Due to limited amount of data available for dispersed oil and amphipods, information on other crustacean species and other oil-related compounds was also collected. Only one study on dispersed oil and amphipods was found, the LC50 value in this study was similar to the LC50 value of E. gryllus in our study. Since toxicity data are important input to risk assessment and net environmental benefit analyses, and since such data are generally lacking on deep-sea species, the data set produced in this study is of interest to the industry, stakeholders, environmental management, and ecotoxicologists. However, studies including more deep-sea species covering different functional groups are needed to evaluate the sensitivity of the deep-sea compartments to dispersed oil relative to other environmental compartments.

Isolation and Characterization of Microsatellite DNA Markers in the Deep-Sea Amphipod Paralicella tenuipes by Illumina MiSeq Sequencing

Here, we describe the development of 16 polymorphic microsatellite markers using an Illumina MiSeq sequencing approach in the deep-sea amphipod Paralicella tenuipes A total of 25 577 844 DNA sequences were filtered for microsatellite motifs of which 197 873 sequences were identified. From these sequences, 64 had sufficient flanking regions for primer design and 16 of these loci were polymorphic. Between 5 and 30 alleles were detected per locus, with an average of 13.63 alleles per locus, across a total of 120 individuals from 5 separate deep sea trenches from the Pacific Ocean. For the 16 loci, observed and expected heterozygosity values ranged from 0.116 to 0.414 and 0.422 to 0.820, respectively, with one locus displaying significant deviation from Hardy-Weinberg equilibrium. The microsatellite loci that have been isolated and described here are the first molecular markers developed for deep sea amphipods and will be invaluable for elucidating the genetic population structure and the extent of connectivity between deep ocean trenches.

Abyssal fauna of the UK-1 polymetallic nodule exploration claim, Clarion-Clipperton Zone, central Pacific Ocean: Echinodermata

We present data from a DNA taxonomy register of the abyssal benthic Echinodermata collected as part of the Abyssal Baseline (ABYSSLINE) environmental survey cruise ‘AB01’ to the UK Seabed Resources Ltd (UKSRL) polymetallic-nodule exploration claim ‘UK-1’ in the eastern Clarion-Clipperton Zone (CCZ), central Pacific Ocean abyssal plain. Morphological and genetic data are presented for 17 species (4 Asteroidea, 4 Crinoidea, 2 Holothuroidea and 7 Ophiuroidea) identified by a combination of morphological and genetic data. No taxa matched previously published genetic sequences, but 8 taxa could be assigned to previously-described species based on morphology, although here we have used a precautionary approach in taxon assignments to avoid over-estimating species ranges. The Clarion-Clipperton Zone is a region undergoing intense exploration for potential deep-sea mineral extraction. We present these data to facilitate future taxonomic and environmental impact study by making both data and voucher materials available through curated and accessible biological collections.

The Antarctic Circumpolar Current as a diversification trigger for deep-sea octocorals

BACKGROUND

Antarctica is surrounded by the Antarctic Circumpolar Current (ACC), the largest and strongest current in the world. Despite its potential importance for shaping biogeographical patterns, the distribution and connectivity of deep-sea populations across the ACC remain poorly understood. In this study we conducted the first assessment of phylogeographical patterns in deep-sea octocorals in the South Pacific and Southern Ocean, specifically a group of closely related bottlebrush octocorals (Primnoidae: Tokoprymno and Thourella), as a test case to study the effect of the ACC on the population structure of brooding species. We assessed the degree to which the ACC constitutes a barrier to gene flow between northern and southern populations and whether the onset of diversification of these corals coincides with the origin of the ACC (Oligocene-Miocene boundary).

RESULTS

Based on DNA sequences of two nuclear genes from 80 individuals and a combination of phylogeographic model-testing approaches we found a phylogenetic break corresponding to the spatial occurrence of the ACC. We also found significant genetic structure among our four regional populations. However, we uncovered shared haplotypes among certain population pairs, suggesting long-distance, asymmetrical migration. Our divergence time analyses indicated that the separation of amphi-ACC populations took place during the Middle Miocene around 12.6 million years ago, i.e., after the formation of the ACC.

CONCLUSION

We suggest that the ACC constitutes a semi-permeable barrier to these deep-sea octocorals capable of separating and structuring populations, while allowing short periods of gene flow. The fluctuations in latitudinal positioning of the ACC during the Miocene likely contributed to the diversification of these octocorals. Additionally, we provide evidence that the populations from each of our four sampling regions could actually constitute different species.

A chemosynthetic weed: the tubeworm Sclerolinum contortum is a bipolar, cosmopolitan species

BACKGROUND

Sclerolinum (Annelida: Siboglinidae) is a genus of small, wiry deep-sea tubeworms that depend on an endosymbiosis with chemosynthetic bacteria for their nutrition, notable for their ability to colonise a multitude of reducing environments. Since the early 2000s, a Sclerolinum population has been known to inhabit sediment-hosted hydrothermal vents within the Bransfield Strait, Southern Ocean, and whilst remaining undescribed, it has been suggested to play an important ecological role in this ecosystem. Here, we show that the Southern Ocean Sclerolinum population is not a new species, but more remarkably in fact belongs to the species S. contortum, first described from an Arctic mud volcano located nearly 16,000 km away.

RESULTS

Our new data coupled with existing genetic studies extend the range of this species across both polar oceans and the Gulf of Mexico. Our analyses show that the populations of this species are structured on a regional scale, with greater genetic differentiation occurring between rather than within populations. Further details of the external morphology and tube structure of S. contortum are revealed through confocal and SEM imaging, and the ecology of this worm is discussed.

CONCLUSIONS

These results shed further insight into the plasticity and adaptability of this siboglinid group to a range of reducing conditions, and into the levels of gene flow that occur between populations of the same species over a global extent.

Speciation in the dark: diversification and biogeography of the deep-sea gastropod genus Scaphander in the Atlantic Ocean

AIM

The aim of this work was to improve understanding about the mode, geography and tempo of diversification in deep-sea organisms, using a time-calibrated molecular phylogeny of the heterobranch gastropod genus Scaphander.

LOCATION

Atlantic and Indo-West Pacific (IWP) oceans.

METHODS

Two mitochondrial gene markers (COI and 16S) and one nuclear ribosomal gene (28S) from six Atlantic species of Scaphander, and four IWP species were used to generate a multilocus phylogenetic hypothesis using uncorrelated relaxed-clock Bayesian methods implemented in beast and calibrated with the first occurrence of Scaphander in the fossil record (58.7-55.8 Ma).

RESULTS

Two main clades were supported: clade A, with sister relationships between species and subclades from the Atlantic and IWP; and clade B, with two western Atlantic sister species. Our estimates indicate that the two earliest divergences in clade A occurred between the middle Eocene and late Miocene and the most recent speciation occurred within the middle Miocene to Pleistocene. The divergence between the two western Atlantic species in clade B was estimated at late Oligocene-Pliocene.

MAIN CONCLUSIONS

The prevailing mode of speciation in Scaphander was allopatric, but one possible case of sympatric speciation was detected between two western Atlantic species. Sister relationships between IWP and Atlantic lineages suggest the occurrence both of vicariance events caused by the closure of the Tethyan Seaway and of dispersal between the two ocean basins, probably around South Africa during episodic disruptions of the deep-sea regional current system caused by glacial-interglacial cycles. Cladogenetic estimates do not support comparatively older diversification of deep-sea faunas, but corroborate the hypothesis of a pulse of diversification centred in the Oligocene and Miocene epochs. Amphi-Atlantic species were found to occur at deeper depths (bathyal-abyssal) and we hypothesize that trans-Atlantic connectivity is maintained by dispersal between neighbouring reproductive populations inhabiting the abyssal sea floor and by dispersal across the shelf and slope of Arctic and sub-Arctic regions.

A reverse taxonomic approach to assess macrofaunal distribution patterns in abyssal Pacific polymetallic nodule fields

Heightened interest in the exploitation of deep seafloor minerals is raising questions on the consequences for the resident fauna. Assessing species ranges and determination of processes underlying current species distributions are prerequisites to conservation planning and predicting faunal responses to changing environmental conditions. The abyssal central Pacific nodule belt, located between the Clarion and Clipperton Fracture Zones (CCZ), is an area prospected for mining of polymetallic nodules. We examined variations in genetic diversity and broad-scale connectivity of isopods and polychaetes across the CCZ. Faunal assemblages were studied from two mining claims (the eastern German and French license areas) located 1300 km apart and influenced by different productivity regimes. Using a reverse taxonomy approach based on DNA barcoding, we tested to what extent distance and large-scale changes in environmental parameters lead to differentiation in two macrofaunal taxa exhibiting different functions and life-history patterns. A fragment of the mitochondrial gene Cytochrome Oxidase Subunit 1 (COI) was analyzed. At a 97% threshold the molecular operational taxonomic units (MOTUs) corresponded well to morphological species. Molecular analyses indicated high local and regional diversity mostly because of large numbers of singletons in the samples. Consequently, variation in composition of genotypic clusters between sites was exceedingly large partly due to paucity of deep-sea sampling and faunal patchiness. A higher proportion of wide-ranging species in polychaetes was contrasted with mostly restricted distributions in isopods. Remarkably, several cryptic lineages appeared to be sympatric and occurred in taxa with putatively good dispersal abilities, whereas some brooding lineages revealed broad distributions across the CCZ. Geographic distance could explain variation in faunal connectivity between regions and sites to some extent, while assumed dispersal capabilities were not as important.

A hierarchical classification of benthic biodiversity and assessment of protected areas in the Southern Ocean

An international effort is underway to establish a representative system of marine protected areas (MPAs) in the Southern Ocean to help provide for the long-term conservation of marine biodiversity in the region. Important to this undertaking is knowledge of the distribution of benthic assemblages. Here, our aim is to identify the areas where benthic marine assemblages are likely to differ from each other in the Southern Ocean including near-shore Antarctica. We achieve this by using a hierarchical spatial classification of ecoregions, bathomes and environmental types. Ecoregions are defined according to available data on biogeographic patterns and environmental drivers on dispersal. Bathomes are identified according to depth strata defined by species distributions. Environmental types are uniquely classified according to the geomorphic features found within the bathomes in each ecoregion. We identified 23 ecoregions and nine bathomes. From a set of 28 types of geomorphic features of the seabed, 562 unique environmental types were classified for the Southern Ocean. We applied the environmental types as surrogates of different assemblages of biodiversity to assess the representativeness of existing MPAs. We found that 12 ecoregions are not represented in MPAs and that no ecoregion has their full range of environmental types represented in MPAs. Current MPA planning processes, if implemented, will substantially increase the representation of environmental types particularly within 8 ecoregions. To meet internationally agreed conservation goals, additional MPAs will be needed. To assist with this process, we identified 107 spatially restricted environmental types, which should be considered for inclusion in future MPAs. Detailed supplementary data including a spatial dataset are provided.