Evidence for a role of bathymetry and emergence in speciation in the genus Glycera (Glyceridae, Polychaeta) from the deep Eastern Weddell Sea

A latitudinal gradient of deep-sea invasions for marine fishes

Although the tropics harbor the greatest species richness globally, recent work has demonstrated that, for many taxa, speciation rates are faster at higher latitudes. Here, we explore lability in oceanic depth as a potential mechanism for this pattern in the most biodiverse vertebrates - fishes. We demonstrate that clades with the highest speciation rates also diversify more rapidly along the depth gradient, drawing a fundamental link between evolutionary and ecological processes on a global scale. Crucially, these same clades also inhabit higher latitudes, creating a prevailing latitudinal gradient of deep-sea invasions concentrated in poleward regions. We interpret these findings in the light of classic ecological theory, unifying the latitudinal variation of oceanic features and the physiological tolerances of the species living there. This work advances the understanding of how niche lability sculpts global patterns of species distributions and underscores the vulnerability of polar ecosystems to changing environmental conditions.

On the Diversity of Phyllodocida (Annelida: Errantia), with a Focus on Glyceridae, Goniadidae, Nephtyidae, Polynoidae, Sphaerodoridae, Syllidae, and the Holoplanktonic Families

Phyllodocida is a clade of errantiate annelids characterized by having ventral sensory palps, anterior enlarged cirri, axial muscular proboscis, compound chaetae (if present) with a single ligament, and of lacking dorsolateral folds. Members of most families date back to the Carboniferous, although the earliest fossil was dated from the Devonian. Phyllodocida holds 27 well-established and morphologically homogenous clades ranked as families, gathering more than 4600 currently accepted nominal species. Among them, Syllidae and Polynoidae are the most specious polychaete groups. Species of Phyllodocida are mainly found in the marine benthos, although a few inhabit freshwater, terrestrial and planktonic environments, and occur from intertidal to deep waters in all oceans. In this review, we (1) explore the current knowledge on species diversity trends (based on traditional species concept and molecular data), phylogeny, ecology, and geographic distribution for the whole group, (2) try to identify the main knowledge gaps, and (3) focus on selected families: Alciopidae, Goniadidae, Glyceridae, Iospilidae, Lopadorrhynchidae, Polynoidae, Pontodoridae, Nephtyidae, Sphaerodoridae, Syllidae, Tomopteridae, Typhloscolecidae, and Yndolaciidae. The highest species richness is concentrated in European, North American, and Australian continental shelves (reflecting a strong sampling bias). While most data come from shallow coastal and surface environments most world oceans are clearly under-studied. The overall trends indicate that new descriptions are constantly added through time and that less than 10% of the known species have molecular barcode information available.

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.

Hidden diversity in Antarctica: Molecular and morphological evidence of two different species within one of the most conspicuous ascidian species

The Southern Ocean is one of the most isolated marine ecosystems, characterized by high levels of endemism, diversity, and biomass. Ascidians are among the dominant groups in Antarctic benthic assemblages; thus, recording the evolutionary patterns of this group is crucial to improve our current understanding of the assembly of this polar ocean. We studied the genetic variation within Cnemidocarpa verrucosa sensu lato, one of the most widely distributed abundant and studied ascidian species in Antarctica. Using a mitochondrial and a nuclear gene (COI and 18S), the phylogeography of fifteen populations distributed along the West Antarctic Peninsula and Burdwood Bank/MPA Namuncurá (South American shelf) was characterized, where the distribution of the genetic distance suggested the existence of, at least, two species within nominal C. verrucosa. When reevaluating morphological traits to distinguish between genetically defined species, the presence of a basal disk in one of the genotypes could be a diagnostic morphological trait to differentiate the species. These results are surprising due to the large research that has been carried out with the conspicuous C. verrucosa with no differentiation between species. Furthermore, it provides important tools to distinguish species in the field and laboratory. But also, these results give new insights into patterns of differentiation between closely related species that are distributed in sympatry, where the permeability of species boundaries still needs to be well understood.

The complete mitochondrial genome of Glycera chirori Izuka (Annelida: Polychaeta): an evidence of conservativeness between gene arrangement and phylogenesis in Glycera

The complete mitochondrial genome of Glycera chirori Izuka (Annelida: Polychaeta) was presented, which is a circular molecule of 15,930 bp nucleotides. It encodes 37 genes, including 13 PCGs, 22 tRNAs, and two rRNAs. The length of non-coding regions is 1428 bp, and the longest one (1346 bp) is speculated as the control region, which is located between trnA and trnL₂ and is longer than most species in Glycera. The complete mitogenome of G. chirori Izuka consists of 31.2% A, 23.6% C, 12.9% G, and 32.2% T, which has T vs. A skew (−0.02) and C vs. G skew (−0.29), respectively. Phylogenetic analysis indicates the classification status of G. chirori Izuka and the relationship with other species in Glycera, which is closer with Glycera unicornis and Glycera fallax (bootstrap = 100). By comparisons, the gene arrangement of G. chirori Izuka and other seven species in Glycera are identical and they also cluster together in phylogenetic tree with higher support rate, which indicates the conservativeness between gene arrangement and phylogenetic analysis in Glycera. In conclusion, the complete mitochondrial genome of G. chirori Izuka can provide supportive data for further molecular and evolutionary analysis of Glycera.

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.

Closely related octopus species show different spatial genetic structures in response to the Antarctic seascape

Determining whether comparable processes drive genetic divergence among marine species is relevant to molecular ecologists and managers alike. Sympatric species with similar life histories might be expected to show comparable patterns of genetic differentiation and a consistent influence of environmental factors in shaping divergence. We used microsatellite loci to quantify genetic differentiation across the Scotia Arc in three species of closely related benthic octopods, Pareledone turqueti, P. charcoti, and Adelieledone polymorpha. The relative importance of environmental factors (latitude, longitude, depth, and temperature) in shaping genetic structure was investigated when significant spatial genetic structure was uncovered. Isolated populations of P. turqueti and A. polymorpha at these species' range margins were genetically different to samples close to mainland Antarctica; however, these species showed different genetic structures at a regional scale. Samples of P. turqueti from the Antarctic Peninsula, Elephant Island, and Signy Island were genetically different, and this divergence was associated primarily with sample collection depth. By contrast, weak or nonsignificant spatial genetic structure was evident across the Antarctic Peninsula, Elephant Island, and Signy Island region for A. polymorpha, and slight associations between population divergence and temperature or depth (and/or longitude) were detected. Pareledone charcoti has a limited geographic range, but exhibited no genetic differentiation between samples from a small region of the Scotia Arc (Elephant Island and the Antarctic Peninsula). Thus, closely related species with similar life history strategies can display contrasting patterns of genetic differentiation depending on spatial scale; moreover, depth may drive genetic divergence in Southern Ocean benthos.

Genetic variation and local differences in Pacific cod Gadus macrocephalus around Japan

The population structure of the Pacific cod Gadus macrocephalus was examined using 15 microsatellite loci and mitochondrial DNA (ND2 region). In total, 274 individuals were sampled from 16 locations around Japan to estimate the level of genetic differentiation and effective population size (N_e ). Pairwise F_ST , analysis of molecular variance and Bayesian clustering analysis suggested the presence of two genetically distinct groups in waters around Japan, with a higher N_e value in the eastern group than in the western group. A possible factor that restricts gene flow between groups may be related to the water temperature differences in the south-western part of the Sea of Japan, where the Tsushima Warm Current flows around the area inhabited by the western group, which may limit migration between the west and east.

DNA barcoding uncovers cryptic diversity in 50% of deep-sea Antarctic polychaetes

The Antarctic marine environment is a diverse ecosystem currently experiencing some of the fastest rates of climatic change. The documentation and management of these changes requires accurate estimates of species diversity. Recently, there has been an increased recognition of the abundance and importance of cryptic species, i.e. those that are morphologically identical but genetically distinct. This article presents the largest genetic investigation into the prevalence of cryptic polychaete species within the deep Antarctic benthos to date. We uncover cryptic diversity in 50% of the 15 morphospecies targeted through the comparison of mitochondrial DNA sequences, as well as 10 previously overlooked morphospecies, increasing the total species richness in the sample by 233%. Our ability to describe universal rules for the detection of cryptic species within polychaetes, or normalization to expected number of species based on genetic data is prevented by taxon-specific differences in phylogenetic outputs and genetic variation between and within potential cryptic species. These data provide the foundation for biogeographic and functional analysis that will provide insight into the drivers of species diversity and its role in ecosystem function.

Trilobodrilus itoi sp. nov., with a Re-Description of T. nipponicus (Annelida: Dinophilidae) and a Molecular Phylogeny of the Genus

The marine interstitial annelid Trilobodrilus itoi sp. nov., the sixth member of the genus, is described on the basis of specimens collected intertidally at Ishikari Beach, Hokkaido, Japan; this is the second species in the genus described from the Pacific Rim. In addition, T. nipponicus Uchida and Okuda, 1943 is re-described based on fresh topotypic material from Akkeshi, Hokkaido, Japan. From both species, we determined sequences of the nuclear 18S and 28S rRNA genes, and the mitochondrial cytochrome c oxidase subunit I (COI) gene. Molecular phylogenetic trees based on concatenated sequences of the three genes showed that T. itoi and T. nipponicus form a clade, which was the sister group to a clade containing the two European congeners T. axi Westheide, 1967 and T. heideri Remane, 1925. The Kimura two-parameter distance for COI was 22.5-22.7% between T. itoi and T. nipponicus, comparable with interspecific values in other polychaete genera. We assessed the taxonomic utility of epidermal inclusions and found that the known six species can be classified into three groups.

The Utility of Genome Skimming for Phylogenomic Analyses as Demonstrated for Glycerid Relationships (Annelida, Glyceridae)

Glyceridae (Annelida) are a group of venomous annelids distributed worldwide from intertidal to abyssal depths. To trace the evolutionary history and complexity of glycerid venom cocktails, a solid backbone phylogeny of this group is essential. We therefore aimed to reconstruct the phylogenetic relationships of these annelids using Illumina sequencing technology. We constructed whole-genome shotgun libraries for 19 glycerid specimens and 1 outgroup species (Glycinde armigera). The chosen target genes comprise 13 mitochondrial proteins, 2 ribosomal mitochondrial genes, and 4 nuclear loci (18SrRNA, 28SrRNA, ITS1, and ITS2). Based on partitioned maximum likelihood as well as Bayesian analyses of the resulting supermatrix, we were finally able to resolve a robust glycerid phylogeny and identified three clades comprising the majority of taxa. Furthermore, we detected group II introns inside the cox1 gene of two analyzed glycerid specimens, with two different insertions in one of these species. Moreover, we generated reduced data sets comprising 10 million, 4 million, and 1 million reads from the original data sets to test the influence of the sequencing depth on assembling complete mitochondrial genomes from low coverage genome data. We estimated the coverage of mitochondrial genome sequences in each data set size by mapping the filtered Illumina reads against the respective mitochondrial contigs. By comparing the contig coverage calculated in all data set sizes, we got a hint for the scalability of our genome skimming approach. This allows estimating more precisely the number of reads that are at least necessary to reconstruct complete mitochondrial genomes in Glyceridae and probably non-model organisms in general.

Regional differentiation and extensive hybridization between mitochondrial clades of the Southern Ocean giant sea spider Colossendeis megalonyx

Assessing the enormous diversity of Southern Ocean benthic species and their evolutionary histories is a central task in the era of global climate change. Based on mitochondrial markers, it was recently suggested that the circumpolar giant sea spider Colossendeis megalonyx comprises a complex of at least six cryptic species with mostly small and non-overlapping distribution ranges. Here, we expand the sampling to include over 500 mitochondrial COI sequences of specimens from around the Antarctic. Using multiple species delimitation approaches, the number of distinct mitochondrial OTUs increased from six to 15-20 with our larger dataset. In contrast to earlier studies, many of these clades show almost circumpolar distributions. Additionally, analysis of the nuclear internal transcribed spacer region for a subset of these specimens showed incongruence between nuclear and mitochondrial results. These mito-nuclear discordances suggest that several of the divergent mitochondrial lineages can hybridize and should not be interpreted as cryptic species. Our results suggest survival of C. megalonyx during Pleistocene glaciations in multiple refugia, some of them probably located on the Antarctic shelf, and emphasize the importance of multi-gene datasets to detect the presence of cryptic species, rather than their inference based on mitochondrial data alone.

Does a global DNA barcoding gap exist in Annelida?

Accurate identification of unknown specimens by means of DNA barcoding is contingent on the presence of a DNA barcoding gap, among other factors, as its absence may result in dubious specimen identifications - false negatives or positives. Whereas the utility of DNA barcoding would be greatly reduced in the absence of a distinct and sufficiently sized barcoding gap, the limits of intraspecific and interspecific distances are seldom thoroughly inspected across comprehensive sampling. The present study aims to illuminate this aspect of barcoding in a comprehensive manner for the animal phylum Annelida. All cytochrome c oxidase subunit I sequences (cox1 gene; the chosen region for zoological DNA barcoding) present in GenBank for Annelida, as well as for "Polychaeta", "Oligochaeta", and Hirudinea separately, were downloaded and curated for length, coverage and potential contaminations. The final datasets consisted of 9782 (Annelida), 5545 ("Polychaeta"), 3639 ("Oligochaeta"), and 598 (Hirudinea) cox1 sequences and these were either (i) used as is in an automated global barcoding gap detection analysis or (ii) further analyzed for genetic distances, separated into bins containing intraspecific and interspecific comparisons and plotted in a graph to visualize any potential global barcoding gap. Over 70 million pairwise genetic comparisons were made and results suggest that although there is a tendency towards separation, no distinct or sufficiently sized global barcoding gap exists in either of the datasets rendering future barcoding efforts at risk of erroneous specimen identifications (but local barcoding gaps may still exist allowing for the identification of specimens at lower taxonomic ranks). This seems to be especially true for earthworm taxa, which account for fully 35% of the total number of interspecific comparisons that show 0% divergence.

Explaining bathymetric diversity patterns in marine benthic invertebrates and demersal fishes: physiological contributions to adaptation of life at depth

Bathymetric biodiversity patterns of marine benthic invertebrates and demersal fishes have been identified in the extant fauna of the deep continental margins. Depth zonation is widespread and evident through a transition between shelf and slope fauna from the shelf break to 1000 m, and a transition between slope and abyssal fauna from 2000 to 3000 m; these transitions are characterised by high species turnover. A unimodal pattern of diversity with depth peaks between 1000 and 3000 m, despite the relatively low area represented by these depths. Zonation is thought to result from the colonisation of the deep sea by shallow-water organisms following multiple mass extinction events throughout the Phanerozoic. The effects of low temperature and high pressure act across hierarchical levels of biological organisation and appear sufficient to limit the distributions of such shallow-water species. Hydrostatic pressures of bathyal depths have consistently been identified experimentally as the maximum tolerated by shallow-water and upper bathyal benthic invertebrates at in situ temperatures, and adaptation appears required for passage to deeper water in both benthic invertebrates and demersal fishes. Together, this suggests that a hyperbaric and thermal physiological bottleneck at bathyal depths contributes to bathymetric zonation. The peak of the unimodal diversity-depth pattern typically occurs at these depths even though the area represented by these depths is relatively low. Although it is recognised that, over long evolutionary time scales, shallow-water diversity patterns are driven by speciation, little consideration has been given to the potential implications for species distribution patterns with depth. Molecular and morphological evidence indicates that cool bathyal waters are the primary site of adaptive radiation in the deep sea, and we hypothesise that bathymetric variation in speciation rates could drive the unimodal diversity-depth pattern over time. Thermal effects on metabolic-rate-dependent mutation and on generation times have been proposed to drive differences in speciation rates, which result in modern latitudinal biodiversity patterns over time. Clearly, this thermal mechanism alone cannot explain bathymetric patterns since temperature generally decreases with depth. We hypothesise that demonstrated physiological effects of high hydrostatic pressure and low temperature at bathyal depths, acting on shallow-water taxa invading the deep sea, may invoke a stress-evolution mechanism by increasing mutagenic activity in germ cells, by inactivating canalisation during embryonic or larval development, by releasing hidden variation or mutagenic activity, or by activating or releasing transposable elements in larvae or adults. In this scenario, increased variation at a physiological bottleneck at bathyal depths results in elevated speciation rate. Adaptation that increases tolerance to high hydrostatic pressure and low temperature allows colonisation of abyssal depths and reduces the stress-evolution response, consequently returning speciation of deeper taxa to the background rate. Over time this mechanism could contribute to the unimodal diversity-depth pattern.

Population differentiation and species formation in the deep sea: the potential role of environmental gradients and depth

Ecological speciation probably plays a more prominent role in diversification than previously thought, particularly in marine ecosystems where dispersal potential is great and where few obvious barriers to gene flow exist. This may be especially true in the deep sea where allopatric speciation seems insufficient to account for the rich and largely endemic fauna. Ecologically driven population differentiation and speciation are likely to be most prevalent along environmental gradients, such as those attending changes in depth. We quantified patterns of genetic variation along a depth gradient (1600-3800m) in the western North Atlantic for a protobranch bivalve (Nuculaatacellana) to test for population divergence. Multilocus analyses indicated a sharp discontinuity across a narrow depth range, with extremely low gene flow inferred between shallow and deep populations for thousands of generations. Phylogeographical discordance occurred between nuclear and mitochondrial loci as might be expected during the early stages of species formation. Because the geographic distance between divergent populations is small and no obvious dispersal barriers exist in this region, we suggest the divergence might reflect ecologically driven selection mediated by environmental correlates of the depth gradient. As inferred for numerous shallow-water species, environmental gradients that parallel changes in depth may play a key role in the genesis and adaptive radiation of the deep-water fauna.

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

Eurythenes gryllus is one of the most widespread amphipod species, occurring in every ocean with a depth range covering the bathyal, abyssal and hadal zones. Previous studies, however, indicated the existence of several genetically and morphologically divergent lineages, questioning the assumption of its cosmopolitan and eurybathic distribution. For the first time, its genetic diversity was explored at the global scale (Arctic, Atlantic, Pacific and Southern oceans) by analyzing nuclear (28S rDNA) and mitochondrial (COI, 16S rDNA) sequence data using various species delimitation methods in a phylogeographic context. Nine putative species-level clades were identified within E. gryllus. A clear distinction was observed between samples collected at bathyal versus abyssal depths, with a genetic break occurring around 3,000 m. Two bathyal and two abyssal lineages showed a widespread distribution, while five other abyssal lineages each seemed to be restricted to a single ocean basin. The observed higher diversity in the abyss compared to the bathyal zone stands in contrast to the depth-differentiation hypothesis. Our results indicate that, despite the more uniform environment of the abyss and its presumed lack of obvious isolating barriers, abyssal populations might be more likely to show population differentiation and undergo speciation events than previously assumed. Potential factors influencing species' origins and distributions, such as hydrostatic pressure, are discussed. In addition, morphological findings coincided with the molecular clades. Of all specimens available for examination, those of the bipolar bathyal clade seemed the most similar to the 'true' E. gryllus. We present the first molecular evidence for a bipolar distribution in a macro-benthic deep-sea organism.

Niche divergence by deep-sea octocorals in the genus Callogorgia across the continental slope of the Gulf of Mexico

Environmental variables that are correlated with depth have been suggested to be among the major forces underlying speciation in the deep sea. This study incorporated phylogenetics and ecological niche models (ENM) to examine whether congeneric species of Callogorgia (Octocorallia: Primnoidae) occupy different ecological niches across the continental slope of the Gulf of Mexico (GoM) and whether this niche divergence could be important in the evolution of these closely related species. Callogorgia americana americana, Callogorgia americana delta and Callogorgia gracilis were documented at 13 sites in the GoM (250-1000 m) from specimen collections and extensive video observations. On a first order, these species were separated by depth, with C. gracilis occurring at the shallowest sites, C. a. americana at mid-depths and C. a. delta at the deepest sites. Callogorgia a. delta was associated with areas of increased seep activity, whereas C. gracilis and C. a. americana were associated with narrow, yet warmer, temperature ranges and did not occur near cold seeps. ENM background and identity tests revealed little to no overlap in ecological niches between species. Temporal calibration of the phylogeny revealed the formation of the Isthmus of Panama was a vicariance event that may explain some of the patterns of speciation within this genus. These results elucidate the potential mechanisms for speciation in the deep sea, emphasizing both bathymetric speciation and vicariance events in the evolution of a genus across multiple regions.