Phylogeography of a pan-Atlantic abyssal protobranch bivalve: implications for evolution in the Deep Atlantic

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.

Form and function of the mantle edge in Protobranchia (Mollusca: Bivalvia)

We analyzed, by optical and transmission electron microscopy, the morphology and function of the mantle edge, including the formation of the periostracum, of ten species of protobranchs. Five species from the order Nuculida, four species from the order Nuculanida and one species from the order Solemyida were studied. A second outer fold, which seems to function as a template for the internal marginal crenulations of the valves, is present in the crenulated species of Nucula. The minute non-crenulated Ennucula aegeensis shows the glandular basal cells displaced toward the periostracal groove, resembling a minute additional fold between the outer and middle folds. Intense secretion of glycocalyx, together with active uptake of particles, have been observed in the inner epithelium of the middle mantle fold and the whole epithelium of the inner mantle fold in all the studied species. Contrary to the rest of the bivalves, all the protobranchs analyzed have two basal cells involved in the formation of the external nanometric pellicle of the periostracum, a character that would support the monophyly of protobranchs. A three-layered pattern is the general rule for the periostracum in protobranchs, like for other bivalves. The presence of pouches of translucent layer inside the tanned dark layer under periostracal folds is characteristic of the species with a folded periostracum; its function is unclear but could give flexibility to the periostracum. The non-nacreous internal shell layer and the presence of translucent pouches under periostracal folds in Sarepta speciosa resemble those found in nuculanids. However, the free periostracum is rather similar to those of N. hanleyi and E. aegeensis, with a continuous vesicular layer. All the latter supports the inclusion of Sarepta in the order Nuculanida but could indicate either a basal lineage or that the translucent vesicular layer is an adaptive trait.

A global phylogeny of the deep-sea gastropod family Scaphandridae (Heterobranchia: Cephalaspidea): redefinition and generic classification

We present the most comprehensive phylogeny of a globally distributed deep-sea group of gastropods published to date including over 80% of the recognized diversity of the family Scaphandridae. The definition and taxonomic composition of the Scaphandridae has been hampered by the lack of a sound phylogenetic framework and definition of synapomorphic traits. We used a combination of molecular phylogenetics (Bayesian Inference and Maximum Likelihood) based on five gene markers (cytochrome c oxidase subunit I, 12S rRNA, 16S rRNA, 18S rRNA, and 28S rRNA) and morpho-anatomical characters to redefine the Scaphandridae and its genera. A new classification is proposed with the three genera Nipponoscaphander, Sabatia, and Scaphander. Main differences between genera lie on the shells (shape, parietal callus, spire) and male reproductive system (prostate). The species Hamineobulla kawamurai is reassigned to the closely related family Eoscaphandridae, currently defined mostly based on pleisiomorphic traits. Biogeographically the genus Nipponoscaphander is restricted to the Indo-West Pacific; Sabatia is mostly circumscribed to the Indo-West Pacific, but has one lineage present in the north Atlantic Ocean. Polyphyly across ocean realms prevails in the specious and globally distributed genus Scaphander with multiple speciation events between Indo-Pacific and Atlantic lineages but also with several episodes of cladogenesis within realms. Two rare cases of species with a broad distribution spanning the Indo-West Pacific and Atlantic realms are confirmed (S. meridionalis and S. nobilis).

Exploring the diversity of the deep sea—four new species of the amphipod genus Oedicerina described using morphological and molecular methods

Collections of the amphipod genus Oedicerina were obtained during six expeditions devoted to the study of deep-sea environments of the Pacific Ocean. The material revealed four species new to science. Two species (Oedicerina henricisp. nov. and Oedicerina teresae sp. nov.) were found at abyssal depths of the central eastern Pacific in the Clarion-Clipperton Zone; one species (Oedicerina claudei sp. nov.) was recovered in the Sea of Okhotsk (north-west Pacific), and one (Oedicerina lesci sp. nov.) in the abyss adjacent to the Kuril-Kamchatka Trench (KKT). The four new species differ from each other and known species by the shapes of the rostrum, coxae 1 and 4, basis of pereopod 7, armatures of pereonite 7, pleonites and urosomites. An identification key for all known species is provided. The study of the cytochrome c oxidase subunit I gene of the four new species and Oedicerina ingolfi collected in the North Atlantic confirmed their genetic distinction. However, small intraspecific variation within each of the studied species was observed. In the case of the new species occurring across the KKT, the same haplotype was found on both sides of the trench, providing evidence that the trench does not constitute an insurmountable barrier for population connectivity. None of the species have so far been found on both sides of the Pacific.

An approach using ddRADseq and machine learning for understanding speciation in Antarctic Antarctophilinidae gastropods

Sampling impediments and paucity of suitable material for molecular analyses have precluded the study of speciation and radiation of deep-sea species in Antarctica. We analyzed barcodes together with genome-wide single nucleotide polymorphisms obtained from double digestion restriction site-associated DNA sequencing (ddRADseq) for species in the family Antarctophilinidae. We also reevaluated the fossil record associated with this taxon to provide further insights into the origin of the group. Novel approaches to identify distinctive genetic lineages, including unsupervised machine learning variational autoencoder plots, were used to establish species hypothesis frameworks. In this sense, three undescribed species and a complex of cryptic species were identified, suggesting allopatric speciation connected to geographic or bathymetric isolation. We further observed that the shallow waters around the Scotia Arc and on the continental shelf in the Weddell Sea present high endemism and diversity. In contrast, likely due to the glacial pressure during the Cenozoic, a deep-sea group with fewer species emerged expanding over great areas in the South-Atlantic Antarctic Ridge. Our study agrees on how diachronic paleoclimatic and current environmental factors shaped Antarctic communities both at the shallow and deep-sea levels, promoting Antarctica as the center of origin for numerous taxa such as gastropod mollusks.

High diversity and pan-oceanic distribution of deep-sea polychaetes: Prionospio and Aurospio (Annelida: Spionidae) in the Atlantic and Pacific Ocean

Prionospio Malmgren 1867 and Aurospio Maciolek 1981 (Annelida: Spionidae) are polychaete genera commonly found in the deep sea. Both genera belong to the Prionospio complex, whose members are known to have limited distinguishing characters. Morphological identification of specimens from the deep sea is challenging, as fragmentation and other damages are common during sampling. These issues impede investigations into the distribution patterns of these genera in the deep sea. In this study, we employ two molecular markers (16S rRNA and 18S) to study the diversity and the distribution patterns of Prionospio and Aurospio from the tropical North Atlantic, the Puerto Rico Trench and the central Pacific. Based on different molecular analyses (Automated Barcode Gap Discovery, GMYC, pairwise genetic distances, phylogenetics, haplotype networks), we were able to identify and differentiate 21 lineages (three lineages composed solely of GenBank entries) that represent putative species. Seven of these lineages exhibited pan-oceanic distributions (occurring in the Atlantic as well as the Pacific) in some cases even sharing identical 16S rRNA haplotypes in both oceans. Even the lineages found to be restricted to one of the oceans were distributed over large regional scales as for example across the Mid-Atlantic Ridge from the Caribbean to the eastern Atlantic (> 3389 km). Our results suggest that members of Prionospio and Aurospio may have the potential to disperse across large geographic distances, largely unaffected by topographic barriers and possibly even between oceans. Their high dispersal capacities are probably explained by their free-swimming long-lived planktonic larvae.

Appropriate Assignment of Fossil Calibration Information Minimizes the Difference between Phylogenetic and Pedigree Mutation Rates in Humans

Studies that measured mutation rates in human populations using pedigrees have reported values that differ significantly from rates estimated from the phylogenetic comparison of humans and chimpanzees. Consequently, exchanges between mutation rate values across different timescales lead to conflicting divergence time estimates. It has been argued that this variation of mutation rate estimates across hominoid evolution is in part caused by incorrect assignment of calibration information to the mean coalescent time among loci, instead of the true genetic isolation (speciation) time between humans and chimpanzees. In this study, we investigated the feasibility of estimating the human pedigree mutation rate using phylogenetic data from the genomes of great apes. We found that, when calibration information was correctly assigned to the human⁻chimpanzee speciation time (and not to the coalescent time), estimates of phylogenetic mutation rates were statistically equivalent to the estimates previously reported using studies of human pedigrees. We conclude that, within the range of biologically realistic ancestral generation times, part of the difference between whole-genome phylogenetic and pedigree mutation rates is due to inappropriate assignment of fossil calibration information to the mean coalescent time instead of the speciation time. Although our results focus on the human⁻chimpanzee divergence, our findings are general, and relevant to the inference of the timescale of the tree of life.

Incomplete lineage sorting impacts the inference of macroevolutionary regimes from molecular phylogenies when concatenation is employed: An analysis based on Cetacea

Interest in methods that estimate speciation and extinction rates from molecular phylogenies has increased over the last decade. The application of such methods requires reliable estimates of tree topology and node ages, which are frequently obtained using standard phylogenetic inference combining concatenated loci and molecular dating. However, this practice disregards population-level processes that generate gene tree/species tree discordance. We evaluated the impact of employing concatenation and coalescent-based phylogeny inference in recovering the correct macroevolutionary regime using simulated data based on the well-established diversification rate shift of delphinids in Cetacea. We found that under scenarios of strong incomplete lineage sorting, macroevolutionary analysis of phylogenies inferred by concatenating loci failed to recover the delphinid diversification shift, while the coalescent-based tree consistently retrieved the correct rate regime. We suggest that ignoring microevolutionary processes reduces the power of methods that estimate macroevolutionary regimes from molecular data.

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.

Deep-Sea Phylogeographic Structure Shaped by Paleoenvironmental Changes and Ongoing Ocean Currents Around the Sea of Japan in a Crangonid Shrimp, Argis lar

The deep-sea crangonid shrimp, Argis lar, is a highly abundant species from the northern Pacific Ocean. We investigated its phylogeographic and demographic structure across the species' extensive range, using mitochondrial DNA sequence variation to evaluate the impact of deep-sea paleoenvironmental dynamics in the Sea of Japan on population histories. The haplotype network detected three distinct lineages with allopatric isolation, which roughly corresponded to the Sea of Japan (Lineage A), the northwestern Pacific off the Japanese Archipelago (Lineage B), and the Bering Sea/Gulf of Alaska (Lineage C). Lineage A showed relatively low haplotype and nucleotide diversity, a significantly negative value of Tajima's D, and a star-shaped network, suggesting that anoxic bottom-water in the Sea of Japan over the last glacial period may have brought about a reduction in the Sea of Japan population. Furthermore, unexpectedly, the distributions of Lineage A and B were closely related to the pathways of the two ocean currents, especially along the Sanriku Coast. This result indicated that A. lar could disperse across shallow straits through the ocean current, despite their deep-sea adult habitat. Bayesian inference of divergence time revealed that A. lar separated into three lineages approximately 1 million years before present (BP) in the Pleistocene, and then had been influenced by deep-sea paleoenvironmental change in the Sea of Japan during the last glacial period, followed by a more recent larval dispersal with the ocean current since ca. 6 kilo years BP.

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.

Population genetic structure in Sabatieria (Nematoda) reveals intermediary gene flow and admixture between distant cold seeps from the Mediterranean Sea

BACKGROUND

There is a general lack of information on the dispersal and genetic structuring for populations of small-sized deep-water taxa, including free-living nematodes which inhabit and dominate the seafloor sediments. This is also true for unique and scattered deep-sea habitats such as cold seeps. Given the limited dispersal capacity of marine nematodes, genetic differentiation between such geographically isolated habitat patches is expected to be high. Against this background, we examined genetic variation in both mitochondrial (COI) and nuclear (18S and 28S ribosomal) DNA markers of 333 individuals of the genus Sabatieria, abundantly present in reduced cold-seep sediments. Samples originated from four Eastern Mediterranean cold seeps, separated by hundreds of kilometers, and one seep in the Southeast Atlantic.

RESULTS

Individuals from the Mediterranean and Atlantic were divided into two separate but closely-related species clades. Within the Eastern Mediterranean, all specimens belonged to a single species, but with a strong population genetic structure (ΦST = 0.149). The haplotype network of COI contained 19 haplotypes with the most abundant haplotype (52% of the specimens) shared between all four seeps. The number of private haplotypes was high (15), but the number of mutations between haplotypes was low (1-8). These results indicate intermediary gene flow among the Mediterranean Sabatieria populations with no evidence of long-term barriers to gene flow.

CONCLUSIONS

The presence of shared haplotypes and multiple admixture events indicate that Sabatieria populations from disjunct cold seeps are not completely isolated, with gene flow most likely facilitated through water current transportation of individuals and/or eggs. Genetic structure and molecular diversity indices are comparable to those of epiphytic shallow-water marine nematodes, while no evidence of sympatric cryptic species was found for the cold-seep Sabatieria.

Increasing the Depth of Current Understanding: Sensitivity Testing of Deep-Sea Larval Dispersal Models for Ecologists

Larval dispersal is an important ecological process of great interest to conservation and the establishment of marine protected areas. Increasing numbers of studies are turning to biophysical models to simulate dispersal patterns, including in the deep-sea, but for many ecologists unassisted by a physical oceanographer, a model can present as a black box. Sensitivity testing offers a means to test the models’ abilities and limitations and is a starting point for all modelling efforts. The aim of this study is to illustrate a sensitivity testing process for the unassisted ecologist, through a deep-sea case study example, and demonstrate how sensitivity testing can be used to determine optimal model settings, assess model adequacy, and inform ecological interpretation of model outputs. Five input parameters are tested (timestep of particle simulator (TS), horizontal (HS) and vertical separation (VS) of release points, release frequency (RF), and temporal range (TR) of simulations) using a commonly employed pairing of models. The procedures used are relevant to all marine larval dispersal models. It is shown how the results of these tests can inform the future set up and interpretation of ecological studies in this area. For example, an optimal arrangement of release locations spanning a release area could be deduced; the increased depth range spanned in deep-sea studies may necessitate the stratification of dispersal simulations with different numbers of release locations at different depths; no fewer than 52 releases per year should be used unless biologically informed; three years of simulations chosen based on climatic extremes may provide results with 90% similarity to five years of simulation; and this model setup is not appropriate for simulating rare dispersal events. A step-by-step process, summarising advice on the sensitivity testing procedure, is provided to inform all future unassisted ecologists looking to run a larval dispersal simulation.

The role of ontogeny in physiological tolerance: decreasing hydrostatic pressure tolerance with development in the northern stone crab Lithodes maja

Extant deep-sea invertebrate fauna represent both ancient and recent invasions from shallow-water habitats. Hydrostatic pressure may present a significant physiological challenge to organisms seeking to colonize deeper waters or migrate ontogenetically. Pressure may be a key factor contributing to bottlenecks in the radiation of taxa and potentially drive speciation. Here, we assess shifts in the tolerance of hydrostatic pressure through early ontogeny of the northern stone crab Lithodes maja, which occupies a depth range of 4-790 m in the North Atlantic. The zoea I, megalopa and crab I stages were exposed to hydrostatic pressures up to 30.0 MPa (equivalent of 3000 m depth), and the relative fold change of genes putatively coding for the N-methyl-D-aspartate receptor-regulated protein 1 (narg gene), two heat-shock protein 70 kDa (HSP70) isoforms and mitochondrial Citrate Synthase (CS gene) were measured. This study finds a significant increase in the relative expression of the CS and hsp70a genes with increased hydrostatic pressure in the zoea I stage, and an increase in the relative expression of all genes with increased hydrostatic pressure in the megalopa and crab I stages. Transcriptional responses are corroborated by patterns in respiratory rates in response to hydrostatic pressure in all stages. These results suggest a decrease in the acute high-pressure tolerance limit as ontogeny advances, as reflected by a shift in the hydrostatic pressure at which significant differences are observed.

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.

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.

Investigating the Bivalve Tree of Life – an exemplar-based approach combining molecular and novel morphological characters

To re-evaluate the relationships of the major bivalve lineages, we amassed detailed morpho-anatomical, ultrastructural and molecular sequence data for a targeted selection of exemplar bivalves spanning the phylogenetic diversity of the class. We included molecular data for 103 bivalve species (up to five markers) and also analysed a subset of taxa with four additional nuclear protein-encoding genes. Novel as well as historically employed morphological characters were explored, and we systematically disassembled widely used descriptors such as gill and stomach ‘types’. Phylogenetic analyses, conducted using parsimony direct optimisation and probabilistic methods on static alignments (maximum likelihood and Bayesian inference) of the molecular data, both alone and in combination with morphological characters, offer a robust test of bivalve relationships. A calibrated phylogeny also provided insights into the tempo of bivalve evolution. Finally, an analysis of the informativeness of morphological characters showed that sperm ultrastructure characters are among the best morphological features to diagnose bivalve clades, followed by characters of the shell, including its microstructure. Our study found support for monophyly of most broadly recognised higher bivalve taxa, although support was not uniform for Protobranchia. However, monophyly of the bivalves with protobranchiate gills was the best-supported hypothesis with incremental morphological and/or molecular sequence data. Autobranchia, Pteriomorphia, Heteroconchia, Palaeoheterodonta, Archiheterodonta, Euheterodonta, Anomalodesmata and Imparidentia new clade ( = Euheterodonta excluding Anomalodesmata) were recovered across analyses, irrespective of data treatment or analytical framework. Another clade supported by our analyses but not formally recognised in the literature includes Palaeoheterodonta and Archiheterodonta, which emerged under multiple analytical conditions. The origin and diversification of each of these major clades is Cambrian or Ordovician, except for Archiheterodonta, which diverged from Palaeoheterodonta during the Cambrian, but diversified during the Mesozoic. Although the radiation of some lineages was shifted towards the Palaeozoic (Pteriomorphia, Anomalodesmata), or presented a gap between origin and diversification (Archiheterodonta, Unionida), Imparidentia showed steady diversification through the Palaeozoic and Mesozoic. Finally, a classification system with six major monophyletic lineages is proposed to comprise modern Bivalvia: Protobranchia, Pteriomorphia, Palaeoheterodonta, Archiheterodonta, Anomalodesmata and Imparidentia.

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.

Into the deep: a phylogenetic approach to the bivalve subclass Protobranchia

A molecular phylogeny of Protobranchia, the subclass of bivalve mollusks sister to the remaining Bivalvia, has long proven elusive, because many constituent lineages are deep-sea endemics, which creates methodological challenges for collecting and preserving genetic material. We obtained 74 representatives of all 12 extant protobranch families and investigated the internal phylogeny of this group using sequence data from five molecular loci (16S rRNA, 18S rRNA, 28S rRNA, cytochrome c oxidase subunit I, and histone H3). Model-based and dynamic homology parsimony approaches to phylogenetic reconstruction unanimously supported four major clades of Protobranchia, irrespective of treatment of hypervariable regions in the nuclear ribosomal genes 18S rRNA and 28S rRNA. These four clades correspond to the superfamilies Nuculoidea (excluding Sareptidae), Nuculanoidea (including Sareptidae), Solemyoidea, and Manzanelloidea. Salient aspects of the phylogeny include (1) support for the placement of the family Sareptidae with Nuculanoidea; (2) the non-monophyly of the order Solemyida (Solemyidae+Nucinellidae); (3) and the non-monophyly of most nuculoid and nuculanoid genera and families. In light of this first family-level phylogeny of Protobranchia, we present a revised classification of the group. Estimation of divergence times in concert with analyses of diversification rates demonstrate the signature of the end-Permian mass extinction in the phylogeny of extant protobranchs.