Distinct genetic differentiation and species diversification within two marine nematodes with different habitat preference in Antarctic sediments
© The Author(s). 2017
Received: 9 December 2016
Accepted: 12 May 2017
Published: 30 May 2017
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.
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.
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).
KeywordsAntarctica Continental shelf Cryptic species Desmodora Dispersal Nematoda Phylogeny Population genetics Sabatieria
Marine nematodes are the most abundant metazoan inhabitants of seafloor sediments and estimates of total marine species numbers (including parasites) are believed to exceed 50,000 . Yet most of this diversity remains undescribed due to the difficult and time-consuming taxonomy, and logistically challenging sampling in remote (e.g., deep-sea) environments [2, 3]. To date, the number of described nematode species in the marine environment is ca. 12,000 (of which 6900 are free-living; ), which obviously covers only a limited fraction of total estimates [4, 5]. As a consequence, accurate characterisation of species diversity and biogeographic distributions for this highly abundant phylum is currently lacking and the study of macroecological patterns is inevitably limited to genus-level data. Additionally, the observation of extensive cryptic species diversity in species with different life history traits [6–11] further hampers correct estimation of global and local species diversity. Globally distributed nematode species may in fact constitute a series of cryptic species with more restricted geographic distribution for which morphological differences are not readily observable . Coexistence of such cryptic nematode species at local scales may then partly be driven by differential ecological tolerances, preferences for abiotic factors and/or resource differentiation [12–14]. A profound understanding of species-specific preferences and life history traits (e.g., habitat preference, dispersal ability), in combination with knowledge on physical drivers of connectivity among marine populations (e.g., hydrodynamic forces, habitat characteristics) is thus imperative in the study of nematode species distribution patterns across various spatial scales and habitats.
The marine environment presents few obvious barriers to gene flow, and this has led to predictions of little genetic structure of marine species over large spatial scales , and speciation being mainly driven by broad-scale allopatric processes (e.g., [16, 17]). Alternatively, speciation in the sea can proceed without absolute barriers to gene flow along ecological boundaries, when divergent selection overwhelms the homogenising effect of gene flow (sympatric speciation; [18, 19]). In the Southern Ocean, genetic exchange between locations around the continent may be facilitated by the eastward Antarctic Circumpolar Current (ACC) and westward Antarctic Coastal Current (ACoC) systems, as well as the Weddell gyre [20–22]. Several Antarctic marine benthic invertebrates indeed have circum-Antarctic and eurybathic distributions (e.g., [23, 24]), reflecting a long history of isolation from other water masses and migrations up and down the slope during glacial cycles [22, 25]. However, evidence from DNA markers showed that their populations present substantial genetic differentiation and may be isolated over smaller spatial scales and depth ranges than previously thought [25, 26].
In light of current knowledge on cryptic speciation, cosmopolitan distribution and genetic structure in nematodes we expected to find 1/cryptic nematode species and strong genetic structuring in view of the large geographic distances between locations; 2/increased population genetic structure with increasing geographic distance (cf. isolation-by-distance principle IBD; [39, 40]), given the presumed limited dispersal capacity for nematodes (see also ); 3/ stronger population genetic structuring in Sabatieria than in Desmodora based on its preference for deeper sediment depths, assuming that surface dwellers have higher dispersal probability than species that occur deeper in the sediment.
Nematode collection, isolation and vouchering
Overview of the five sampling locations, sampled specimens and number of sequences for species of both genera
total per location
total per location
114 | 10 | 2
25 | 4 | -
8 | 5 | -
- | - | -
147 | 19 | 2
17 | - | 9
- | - | -
17 | - | 9
8 | - | -
25 | 3 | -
19 | 4 | -
- | - | -
52 | 7 | -
5 | - | 8
- | - | 12
5 | - | 20
27 | 4 | -
1 | - | -
8 | 1 | -
- | - | -
36 | 5 | -
- | - | -
- | - | -
- | - | -
4 | 1 | -
1 | - | 1
- | - | -
2 | - | -
7 | 1 | 1
1 | - | 1
- | - | -
1 | - | 1
46 | 5 | -
16 |2 | 11
- | - | -
22 | 3 | 2
84 | 10 | 13
2 | - | 7
- | - | -
2 | - | 7
199 | 20 | 2
68 | 9 | 12
35 | 10 | -
24 | 3 | 2
326 | 42 | 16
25 | - | 25
- | - | 12
25 | - | 37
DESS samples were carefully screened under a stereomicroscope (50 × magnification) and individuals from both targeted genera were handpicked with a fine needle and washed in three separate dishes with sterile distilled water to remove all remaining DESS compounds. Individuals were mounted on a temporary microscopic slide in a drop of distilled water and identified under a Leica DLMS compound microscope (1000 × magnification). During this ‘vouchering’ process, each specimen was assigned to a certain morphological group based on conspicuous body features, which were photographed at different magnifications. For Sabatieria, we distinguished three morphological groups, with differences in tail shape, number of amphid turns and male copulatory organs (see Additional file 1: Table S1.1). For Desmodora at least three distinct morphological groups (cf. D. campbelli, D. sp.A/B and D. sp.D of ; Additional file 1: Table S1.1) were recognised based on body length, position and length of somatic setae, male precloacal supplements and spicule apparatus, and presence of lateral body lines. After the vouchering process (5–10 min per specimen), each nematode was transferred individually into a microcentrifuge tube containing 20 μL Worm Lysis Buffer (WLB: 50 mM KCl, 10 mM Tris–HCl pH 8.3, 2.5 mM MgCl2, 0.45% NP40, 0.45% Tween 20; ), and stored at −20 °C.
DNA extraction, amplification and sequencing
Proteinase K (1 μL; 10 mg mL−1) was added to the WLB-stored specimens for digestion after which samples were incubated at 65 °C for 1 h, followed by 10 min at 95 °C. They were centrifuged for 1 min at 14,000 rpm prior to usage of the DNA. Three markers were amplified by polymerase chain reaction (PCR): the nuclear ribosomal DNA (rDNA) Internal Transcribed Spacer (ITS) region, part of the mitochondrial cytochrome c oxidase subunit 1 (COI) gene, and for a subset of Sabatieria specimens, part of the nuclear small subunit (18S) rDNA. Details on the amplification process can be found in Additional files 2, 3, and 4 of the Supplementary Information. Sequences can be found in GenBank under accession numbers LT577954 – LT578168.
DNA sequence alignments
Electropherograms of the COI, ITS and 18S sequences were analysed and assembled with lasergene® Version 7.1.0 (DNASTAR, Madison, WI) and trimmed to remove primer ends. Sequence length after trimming varied between 307 and 313 bp for COI of Sabatieria, 647–662 for COI of Desmodora, 647–662 for ITS of Sabatieria, 591–599 for ITS of Desmodora, and 845–864 bp for 18S of Sabatieria. Sequences were aligned for the two genera and each gene separately using clustalw v2 with default gap opening/extension costs of 15/6.66 in mega v6.0 [46, 47]. COI sequences were translated to amino acids using genetic translation Table 5 (invertebrate mitochondrial) to assign the correct reading frame and to verify that no indels or stop codons were present in the alignment. For each alignment, the best fit substitution model was selected in jModelTest [48, 49], using the Bayesian Information Criterion (BIC) (Additional file 1: Table S1.3).
The different alignments were analysed using different tree construction algorithms to inspect for the presence of concordant terminal clades among different markers. Maximum likelihood (ML) trees (bootstrap replication = 1000) were generated with RAxML v8.2.4 . Ultrametric trees were produced using beauti v1.8.2 and beast v1.8.2 (Bayesian Evolutionary Analysis Sampling Trees; ) under different substitution models (Additional file 1: Table S1.3), lognormal relaxed clock model, and coalescent tree prior. A Markov Chain Monte Carlo analysis was run for 10 million generations, of which every 1000th generation was sampled, resulting in 10,000 Bayesian trees. Convergence of runs was checked in Tracer v1.6 , after which the first 5000 trees were discarded as burn-in, while the last 5000 trees were used to construct a consensus tree in TreeAnnotator v1.8.2 (beast package) and define posterior probabilities. Resulting consensus trees for all markers were visualised in FigTree (http://tree.bio.ed.ac.uk/software/figtree/) and used in further analysis. ML and beast analyses were run on the xsede server of the cipres Science Gateway v3.3 (https://www.phylo.org; ).
DNA-based species delimitation
To test whether sequence datasets constituted a single or multiple species, a General Mixed Yule Coalescent (GMYC) model approach was applied . Using the ultrametric gene tree as input, the GMYC algorithm compares two alternative models: i) a single coalescence model that assumes a single species, and ii) a model that combines a coalescent model of intraspecific branching with a Yule model for interspecific branching, thus assuming multiple species. The location of the switch (threshold T) from speciation to coalescence nodes is then fitted on the tree, resulting in an estimation of species diversity. Species delimitation under a single-threshold GMYC model was assessed in R  using packages ‘ape’  and ‘splits’ . Lineages-Through-Time (LTT) plots marking the position of threshold T on a relative timescale were constructed in R.
The presence of species-level lineages in sequence variation was also assessed by means of statistical parsimony . tcs v1.21 software  partitioned the data into independent haplotype networks (gaps = missing data), connected by changes that are non-homoplastic with a probability of 95%. Final TCS haplotype networks [58, 60] were built using the PopART software (http://popart.otago.ac.nz), which only takes unambiguous sites into account.
We relied on a conservative consensus approach towards reconciling the results of the different species delimitation methods to maximise the reliability of species boundaries. More specifically, we recognised species clades that 1/received high nodal support (at least 75% bootstrap support in the ML tree), 2/showed compatible patterns based on statistical parsimony and GMYC analyses, 3/formed concordant clades in the trees inferred from nuclear and mitochondrial markers and/or expressed different morphological characteristics. Mean inter- and intraspecific differences (using a K2P (+ G) correction; Additional file 1: Table S1.3) were calculated in mega v6.0.
Population genetic analyses were performed on ITS for Sabatieria and COI for Desmodora species as these were the most complete datasets (see later). Single-level Analysis of Molecular Variance (AMOVA; 1000 permutations, 0.05 significance level) was carried out in Arlequin v220.127.116.11  to calculate fixation index Φ st . The different sampling sites were considered the different populations. Only species clades consisting of more than two populations with more than five individuals each were included in population genetic analyses. Standard measures of genetic variation within populations, such as nucleotide diversity (π; ) and gene diversity (h; [63, 64]) were also assessed in Arlequin. Intra-population and pairwise inter-population divergences were calculated where appropriate, using pairwise deletion of gaps and K2P-corrected distances (based on jModelTest results, Additional file 1: Table S1.3). Finally, isolation by distance (IBD) was assessed through Mantel testing in ibdws v3.23  based on DNA sequences (ignoring gaps; between-population distance Φ st; between-sequence distance K2P) and 1000 randomisations.
DNA-based species delimitation
Statistical parsimony analysis collapsed the 326 ITS sequences into 95 haplotypes (sequence divergence based on K2P distances = 0.2–26%) and 7 separate haplotype networks (Ia/b, Ic, IIa, IIb, IIc, III and IV; connection limit = 95% or 11 mutations), all corresponding to clades or sub-clades of the Bayesian tree (Fig. 2). The GMYC model gave a significantly better fit for the ITS data (likelihood ratio = 20.6; P < 0.001) than did the null model assuming uniform branching rates. The position of the threshold time T, marking the transition from between- to within-species rate of lineage branching, was estimated at −0.004 on a relative timescale (Fig. 2 upper left), and coincided with 21 species clusters. The confidence interval for the estimated number of species ranged from 8 to 26. As opposed to ITS, the GMYC model was insignificant when applied to 18S and COI data (P > 0.1), possibly as a consequence of the low number of sequences available.
Mean intra- and interspecific genetic divergence for Sabatieria
ITS Sabatieria (n = 326; 679 bp)
1.40 ± 0.28
11.09 ± 1.24
3.73 ± 0.50
15.16 ± 1.59
20.71 ± 1.93
1.26 ± 0.15
14.92 ± 1.60
18.56 ± 1.73
19.86 ± 1.90
0.22 ± 0.08
18S Sabatieria (n = 42; 864 bp)
0.15 ± 0.05
0.24 ± 0.08
0.24 ± 0.11
1.57 ± 0.41
1.70 ± 0.41
0.22 ± 0.09
1.48 ± 0.33
1.58 ± 0.34
2.87 ± 0.52
1.13 ± 0.28
COI Sabatieria (n = 16; 313 bp)
0.00 ± 0.00
25.20 ± 3.24
1.49 ± 0.38
37.78 ± 4.37
37.09 ± 4.14
0.64 ± 0.45
Single-level AMOVA results for each Sabatieria species based on ITS sequence data
Source of variation
Pairwise Φ st values between populations of the different Sabatieria species based on ITS sequence data
Species I (n = 200)
Species II (n = 66)
Species III (n = 35)
Species III and IV were restricted to one side. Species III occurred at the western side of the Weddell Sea and consisted of three populations (SG, SO & KG) for which genetic structuring was significant, but considerably lower than for species I and II for the same populations on this side of the Weddell Sea (AMOVA Φ st species III = 0.178, P < 0.001; Table 3; Φ st species I & II = 0.589–0.599, P < 0.001; results not shown). Within-population variation for species III (~ 82%) exceeded that between populations (17.8%). Genetic differences were non-significant between locations SG and SO (Table 4), which also shared one haplotype (Additional file 1: Table S1.4). Average K2P distances between these populations were also clearly lower than for the other two species (Additional file 1: Table S1.5). Species IV was restricted to the two locations at the eastern Weddell Sea, and comprised 11 haplotypes. Within-population divergence was comparable or even larger than between-population differences, which were non-significant (Additional file 1: Table S1.5).
Despite the observation that main differences between populations of species were situated between different sides of the Weddell Sea (hence, at a large spatial scale), genetic divergence did not consistently decrease with increasing geographic distance (IBD r-values for species I, II and III were non-significant; P > 0.05; Additional file 1: Table S1.6).
DNA-based species delimitation
Mean intra- and interspecific genetic divergence for COI of Desmodora
COI Desmodora (n = 37; 662 bp)
1.76 ± 0.25
23.44 ± 2.08
1.59 ± 0.25
Single-level AMOVA main and pairwise results for Desmodora species I based on COI sequence data
Source of variation
Pairwise Φ st (n = 24)
Conflict between morphological and phylogenetic species definitions in Sabatieria and Desmodora
Objective species delimitation is challenging in animal groups where taxonomic information is incomplete and scattered, yet remains fundamental in biodiversity research . For this reason, a combination of several techniques and a conservative method were adopted to delineate species in this study. Congruence in the outcomes of various species delimitation approaches led to the recognition of four species-level lineages for Sabatieria and two for Desmodora. Not all of these coincided with the initial morphologically defined groups, and vice versa (Additional file 1: Table S1.1). In fact, rates of phenotypic and molecular divergence do not always converge , which makes species delimitation all the more tricky. Especially for relatively young species there might be an offset between the process of speciation and the acquisition of secondary properties such as distinct morphology. However, sequence divergence for COI in both genera was substantial (Sabatieria: 25–38%; Desmodora: 23%; Tables 2, 5; Figs. 3b, 4b), making the possibility of dealing with recent speciation less likely in this case. Within the genus Sabatieria, two out of four species differed from the others in morphological appearance (species III and IV), while the other two (species I and II) were not readily distinguishable and might constitute cryptic species. Cryptic speciation is common in marine free-living nematode genera (see  and references therein) and has also been recovered in limno-terrestrial nematodes on the Antarctic continent . Also in other Southern Ocean benthic inhabitants, recent molecular findings have indicated that species which were previously considered eurybathic and/or circum-Antarctic can in fact be partitioned into cryptic species according to depth or geography [22, 25, 68]. Such a pattern has been observed among a variety of benthic invertebrates (e.g., amphipods, isopods, bivalves, crinoids and octopods), indicating that it is a common phenomenon among Southern Ocean species, and may result from isolation of smaller populations undergoing genetic bottlenecks in shelf or slope refugia during glacial times (see  for an overview). Especially for species with low dispersal capacity (such as nematodes), recolonisation of the continental shelf can be a slow process, possibly leading to (cryptic) speciation before secondary contact between previously isolated populations occurs. In contrast to Sabatieria, Desmodora specimens showed no evidence of cryptic speciation. Instead, the opposite phenomenon was observed where (conspicuous) morphological characteristics were not diagnostic in the delimitation of species. This observation of high intraspecific morphological variation for Desmodora casts doubt on previous reports of six different species within the genus based on morphological data for the same locations . Similar high levels of intraspecific variation in morphology have been reported in the deep-sea nematode genus Acantholaimus from the Pacific , and in a Paracanthonchus species rafting on seaweed along the Brazilian coast . In the latter study, these intraspecific morphological differences occurred in the absence of genetic differentiation among its populations. Some nematodes are even capable of resource polyphenism, a situation in which different phenotypes are induced by different thresholds of an environmental cue during their development [70, 71]. As such, relying on morphology alone when discriminating between nematode species may cause substantial bias.
Wide and narrow species ranges in Sabatieria and Desmodora
This study showed that both Sabatieria and Desmodora contained species with wide ranges across the Weddel Sea, as well as species with more limited ranges on either side of the Weddell Sea. This combination of wide and narrow species ranges has been noted in several other Antarctic benthic invertebrates [72, 73], and these contrasting stories have been linked to the species’ dispersive capacities as well as different survival mechanisms during past glacial cycles . Sabatieria species I and II (and also Desmodora species I) were distributed across locations separated by the deep Weddell Sea, indicating a connection at some point in time. Wide and even cosmopolitan species ranges have been reported in marine nematodes (e.g., [2, 9]) and can reflect ongoing dispersal as well as historical connections . Given the fact that nematodes are passive dispersers and that locations in this study are separated by several hundreds of km, historical connectivity might be very important in this case (cf. [74, 75]). High levels of genetic divergence between species (Tables 2, 5) and long branches in tree topologies (Figs 3, 4) seem to support speciation in the distant past. On an evolutionary timescale, the origin of modern Antarctic biota is put shortly after the Gondwana break-up, which marked the onset of vicariance, speciation and diversification [76, 77]. Yet the resulting Antarctic Circumpolar Current (ACC) maintained a certain level of horizontal connectivity between species and populations along the continent, reflected in circum-Antarctic distributions observed in several benthic invertebrate species . The large-scale distribution of both Sabatieria and Desmodora species might have a similar early origin of speciation followed by long-distance dispersal mediated by the presence of large current systems (ACC, ACoC, Weddell gyre) and relatively homogeneous environmental conditions (e.g., seabed temperatures) in the area [20, 78].
High population genetic structure suggests low levels of gene flow in the Southern Ocean
The physical setting of the Southern Ocean – without obvious barriers to gene flow and with the presence of large-scale currents capable of mediating long-distance dispersal – did not change much over the course of history. Combined with the large population sizes of nematodes and the possibility of passive dispersal, this should result in mild genetic differentiation over large distances . Nevertheless, population genetic structuring within Sabatieria and Desmodora species was substantial. Haplotypes were generally confined to a single geographic location or shared between neighbouring sites (only two Sabatieria haplotypes had representatives at both sides of the Weddell Sea; Fig. 2; Additional file 1: Table S1.4), a characteristic of closed populations and not uncommon in taxa that lack pelagic development [25, 74]. Pairwise Φ st values for Sabatieria species I and II were significant in most cases and largest between locations at different sides of the Weddell Sea (Table 4). Similarly large genetic differences between eastern and western Weddell Sea were also revealed by COI and ITS sequences of benthic ostracods in the area . Desmodora species I also showed highly significant pairwise Φ st values (Table 6) but largest differences were situated between populations SG and SO, rather than between eastern and western Weddell Sea locations (Additional file 1: Table S1.5). These high levels of population genetic differentiation can have multiple origins. First, they might reflect poor dispersal capacity  and suggest that contemporary gene flow between populations is strongly limited at the spatial scale considered here. Similar studies for coastal and estuarine nematodes have demonstrated that population genetic structure can be significant at scales of 100 km and less [7, 8, 11], which is well below distances between the different locations for this study. If gene flow is indeed limited due to dispersal limitation, the large observed population genetic differences might point towards a limited efficiency of the ACC and Weddell gyre in homogenising nematode communities over large distances. Second, barriers to gene flow between populations in a marine setting can exist in many forms, such as temperature gradients, depth differences and large areas of unsuitable habitat conditions [11, 79]. The large pairwise differences between populations at both sides of the Weddell Sea and along the Scotia Arc might therefore result from such ‘invisible’ barriers to gene flow rather than true dispersal limitation. Finally, even in the presence of extensive dispersal between habitat patches, populations can show large genetic differences due to differences in the succesful establishment and reproduction of dispersers after settling in a new environment . Local habitat conditions and species-specific niche preferences, followed by rapid adaptation and population growth may result in situations where priority effects, founder effects and genetic bottlenecks result in certain haplotypes being favoured over others . However, such processes are generally assumed to be of less importance at large spatial scales .
Phylogeographic patterns across the Weddell Sea do not support isolation by distance
Strong population genetic structure at large spatial scales (> 300 km) has been observed in many marine species [11, 75, 81], and has often been attributed to an isolation-by-distance mode of genetic differentiation. Yet for all species of Sabatieria and Desmodora with sufficient sample size, no IBD was observed (Additional file 1: Table S1.6). The reason for this is probably related to large variability in genetic divergence between Antarctic Peninsula and Scotia Arc populations. For example, in Sabatieria species I, gene flow was not restricted between populations SG and KG, located approximately 1600 km apart (non-significant small genetic distance; Table 4, Additional file 1: Table S1.5) but was very much so between SG and SO, which are separated by 900 km distance. This pattern was reversed in species III, where pairwise genetic differences between SG and SO were non-significant (Table 4, Additional file 1: Table S1.5). Within Desmodora species I, genetic differences were larger between SG and SO than between either of them and location BX at the other side of the Weddell Sea. Although it has been argued that the tip of the Antarctic Peninsula and Scotia Arc are highly connected due to the ACC system (e.g., ), our population genetic results do not support this. Instead, there seems to be a rather random pattern of genetic structuring between populations at the western Weddell Sea. Genetic structuring in other marine nematodes has shown such a chaotic pattern , which may be linked to oceanographic currents or other environmental variables  posing a certain level of biotic or habitat filtering on dispersing or settling individuals. For example, smaller-scale bottom currents and dynamics might hamper successful settlement, thereby decoupling dispersal from geographic distance (although mainly tested for species with pelagic larvae; ). Such current dynamic data were not assessed at the time of sampling, but could provide an explanation for the random genetic structuring in the stations near the Peninsula. In terms of local environmental conditions for the study locations, these were strongly correlated with geographic distance since largest discrepancies were noted between locations at both sides of the Weddell Sea. AUS and BX clearly had colder bottom temperatures (almost – 2 °C), lower amount of fresh food (assessed as chlorophyll a concentration in the sediment) and coarser sediment than the other three locations in the vicinity of the Antarctic Peninsula (Hauquier, personal observations). Hence, for the study area considered here, measured abiotic variables do not explain the observations of larger population genetic differences between closely located sites. Other variables (e.g., oxygen content) may provide additional explanations, but have not been assessed at the time of sampling. In any case, a better understanding of the complex interactions between species-specific life history traits on one hand, and habitat characteristics and hydrodynamics on the other hand may help to understand the highly variable dispersal patterns through space and time (see review by ), and the substantial patchiness observed in nematode community composition .
Gene flow in the Weddell Sea is strongly reduced in both genera, but more so in the deeper sediment dwelling Sabatieria species
The two genera in this study share a similar endobenthic lifestyle, but population genetic structuring was more pronounced within the Sabatieria species than within Desmodora species I (cf. AMOVA results). This may be the result of their differential vertical distribution and feeding habits. Nematode dispersal is predominantly passive and mediated through hydrodynamic forces, but individuals living in sediment surface layers are more susceptible to resuspension and transportation in the water column, while deeper dwellers are rarely resuspended [36, 38, 84]. Desmodora prefers surface sediments where it can feed on algal particles scraped off the sediment grains, which potentially facilitated contemporary and historical gene flow over larger areas. Dispersal capacity of organisms plays an important role in connectivity between populations, and previous studies have indicated differences in structuring processes between active and passive dispersers (e.g., [75, 85]). Results of this study thus extend this knowledge and support the idea that vertical distribution in the sediment can be an important proxy for dispersal probability in marine nematodes (see also ).
Our results demonstrate the occurrence of cryptic speciation in Antarctic continental shelf nematodes, and provide evidence for different mechanisms underlying spatial genetic structure within surface- and deeper-sediment dwelling nematode taxa. Historically, ocean current systems such as the ACC and Weddell gyre in the area may have served as a transportation route for species across the Weddell Sea, mainly for taxa occurring in surface sediments such as Desmodora, which showed less geographic structure in its distribution than the Sabatieria species. Currently, dispersal limitation in marine nematodes effectively hampers large-scale connectivity between populations across the Weddell Sea. At a smaller spatial scale, population genetic structuring on the western side of the Weddell Sea is rather random.
The authors wish to thank the scientific community, captain and crew of RV Polarstern expedition PS77/ANT-XXVII/3. Andy Vierstraete is thanked for help with primer development. This manuscript has been improved after revision by two anonymous reviewers.
Availability of data and materials
Sequences of the different genetic markers have been deposited to the European Nucleotide Archive (ENA: http://www.ebi.ac.uk/ena) and GenBank and can be found under accession numbers LT577954 – LT578168. All other data are included in this published article and its supplementary information files.
FH acknowledges a doctoral fellowship from the Research Foundation of Flanders (FWO11/ASP/256). This work contributes to project BR/132/A1/vERSO of the Belgian Science Policy (BELSPO/BRAIN).
FH, FL, SD and AV conceived the ideas and concept for this study; FH collected the samples; FH and AR performed the laboratory analyses; FH, FL, AR, SD and AV analysed the data; FH, FL, SD and AV contributed to the writing. All authors read and approved the final version of the manuscript.
The authors declare that they have no competing interests.
Consent for publication
Ethics approval and consent to participate
The manuscript does not report on or involve the use of any vertebrate or human data or tissue, therefore this section is not applicable.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- Appeltans W, Ahyong ST, Anderson G, Angel MV, Artois T, Bailly N, Bamber R, Barber A. The Magnitude of Global Marine Species Diversity. Curr Biol. 2012;22:2189–202.View ArticlePubMedGoogle Scholar
- Bik HM, Thomas WK, Lunt DH, Lambshead PJD. Low endemism, continued deep-shallow interchanges, and evidence for cosmopolitan distributions in free-living marine nematodes (order Enoplida). BMC Evol Biol. 2010;10:389–99.View ArticlePubMedPubMed CentralGoogle Scholar
- De Ley P, De Ley IT, Morris K, Abebe E, Mundo-Ocampo M, Yoder M, Heras J, Waumann D, Rocha-Olivares A, Burr AHJ, Galdwin JG, Thomas WK. An integrated approach to fast and informative morphological vouchering of nematodes for applications in molecular barcoding. Philos Trans R Soc B. 2005;360:1945–58.View ArticleGoogle Scholar
- Bouchet P. The magnitude of marine biodiversity. In: Duarte CM, editor. The Exploration of Marine Biodiversity: Scientific and Technological Challenges. Madrid: Fundación BBVA; 2006. p. 31–64.Google Scholar
- Bucklin A, Steinke D, Blanco-Bercial L. DNA Barcoding of Marine Metazoa. Annu Rev Mar Sci. 2011;3:18.1–18.38.View ArticleGoogle Scholar
- De Oliveira DAS, Decraemer W, Holovachov O, Burr J, De Ley IT, De Ley P, Moens T, Derycke S. An integrative approach to characterize cryptic species in the Thoracostoma trachygaster Hope, 1967 complex (Nematoda: Leptosomatidae). Zool J Linnean Soc. 2012;164:18–35.View ArticleGoogle Scholar
- Derycke S, Remerie T, Vierstraete A, Backeljau T, Vanfleteren J, Vincx M, Moens T. Mitochondrial DNA variation and cryptic speciation within the free-living marine nematode Pellioditis marina. Mar Ecol Prog Ser. 2005;300:91–103.View ArticleGoogle Scholar
- Derycke S, Backeljau T, Vlaeminck C, Vierstraete A, Vanfleteren J, Vincx M, Moens T. Spatiotemporal analysis of population genetic structure in Geomonhystera disjuncta (Nematoda, Monhysteridae) reveals high levels of molecular diversity. Mar Biol. 2007;151:1799–812.View ArticleGoogle Scholar
- Derycke S, Fonseca G, Vierstraete A, Vanfleteren J, Vincx M, Moens T. Disentangling taxonomy within the Rhabditis (Pellioditis) marina (Nematoda, Rhabditidae) species complex using molecular and morphological tools. Zool J Linnean Soc. 2008;152:1–15.View ArticleGoogle Scholar
- Derycke S, De Ley P, De Ley IT, Holovachov O, Rigaux A, Moens T. Linking DNA sequences to morphology: cryptic diversity and population genetic structure in the marine nematode Thoracostoma trachygaster (Nematoda, Leptosomatidae). Zool Scr. 2010;39:276–89.View ArticleGoogle Scholar
- Derycke S, Backeljau T, Moens T. Dispersal and gene flow in free-living marine nematodes. Front Zool. 2013;10:1.View ArticlePubMedPubMed CentralGoogle Scholar
- De Meester N, Derycke S, Bonte D, Moens T. Salinity effects on the coexistence of cryptic species: a case study on marine nematodes. Mar Biol. 2011;158:2717–26.View ArticleGoogle Scholar
- De Meester N, Derycke S, Rigaux A, Moens T. Temperature and salinity induce differential responses in life histories of cryptic nematode species. J Exp Mar Biol Ecol. 2015;472:54–62.View ArticleGoogle Scholar
- Derycke S, De Meester N, Rigaux A, Creer S, Bik H, Thomas WK, et al. Coexisting cryptic species of the Litoditis marina complex (Nematoda) show differential resource use and have distinct microbiomes with high intraspecific variability. Mol Ecol. 2016; doi:10.1111/mec.13597.
- Palumbi SR. Marine speciation on a small planet. Trends Ecol Evol. 1992;7(4):114–8.View ArticlePubMedGoogle Scholar
- Taylor MS, Hellberg ME. Marine radiations at small geographic scales: speciation in neotropical reef gobies (Elacatinus). Evolution. 2005;59(2):374–85.PubMedGoogle Scholar
- Wilke T, Pfenninger M. Separating historic events from recurrent processes in cryptic species: phylogeography of mud snails (Hydrobia spp.). Mol Ecol. 2002;11:1439–51.View ArticlePubMedGoogle Scholar
- Bolnick DI, Fitzpatrick BM. Sympatric speciation: models and empirical evidence. Annu Rev Ecol Evol Syst. 2007;38:459–87.View ArticleGoogle Scholar
- Bowen BW, Rocha LA, Toonen RJ, Karl SA. The origins of tropical marine biodiversity. Trends Ecol Evol. 2013;28:359–66.View ArticlePubMedGoogle Scholar
- Arntz W, Brey T, Gallardo VA. Antarctic zoobenthos. Oceanogr Mar Biol Annu Rev. 1994;32:241–304.Google Scholar
- Deacon GER. The Weddell Gyre. Deep Sea Research Part A. Oceanogr Res Pap. 1979;26(9):981–95.Google Scholar
- Riesgo A, Taboada S, Avila C. Evolutionary patterns in Antarctic marine invertebrates: An update on molecular studies. Mar Genom. 2015;23:1–13. doi:10.1016/j.margen.2015.07.005.View ArticleGoogle Scholar
- Brandt A, Gooday AJ, Brandão SN, Brix S, Brökeland W, Cedhagen T, Choudhury M, Cornelius N, Danis B, De Mesel I, Diaz RJ, Gillan DC, Ebbe B, Howe JA, Janussen D, Kaiser S, Linse K, Malyutina M, Pawlowski J, Raupach M, Vanreusel A. First insights into the biodiversity and biogeography of the Southern Ocean deep sea. Nature. 2007;447:307–11.View ArticlePubMedGoogle Scholar
- Brey T, Dahm C, Gorny M, Klages M, Stiller M, Arntz WE. Do Antarctic benthic invertebrates show an extended level of eurybathy? Antarct Sci. 1996;8(1):3–6.View ArticleGoogle Scholar
- Allcock AL, Strugnell JM. Southern Ocean diversity: new paradigms from molecular ecology. Trends Ecol Evol. 2012;27(9):520–8.View ArticlePubMedGoogle Scholar
- Cowen RK, Gawarkiewicz G, Pineda J, Thorrold SR, Werner FE. Population connectivity in marine systems: An overview. Oceanography. 2007;20(3):14–21.View ArticleGoogle Scholar
- Blouin MS. Molecular prospecting for cryptic species of nematodes: mitochondrial DNA versus internal transcribed spacer. Int J Parasitol. 2002;32:527–31.View ArticlePubMedGoogle Scholar
- Derycke S, Vanaverbeke J, Rigaux A, Backeljau T, Moens T. Exploring the use of cytochrome oxidase c subunit 1 (COI) for DNA barcoding of free-living marine nematodes. PLoS One. 2010;5:e13716. doi:10.1371/journal.pone.0013716.View ArticlePubMedPubMed CentralGoogle Scholar
- Meldal BHM, Debenham NJ, De Ley P, De Ley IT, Vanfleteren JR, Vierstraete AR, Bert W, Borgonie G, Moens T, Tyler PA, Austen MC, Blaxter ML, Rogers AD, Lambshead PJD. An improved molecular phylogeny of the Nematoda with special emphasis on marine taxa. Mol Phylogenet Evol. 2007;42:622–36.View ArticlePubMedGoogle Scholar
- Apolônio Silva de Oliveira D, Decraemer W, Moens T, Amadeu Paiva dos Santos G, Derycke S. Low genetic by high morphological variation over more than 1000 km coastline refutes omnipresence of cryptic diversity in marine nematodes. BMC Evol Biol. 2017;17:71.View ArticlePubMedPubMed CentralGoogle Scholar
- Guilini K, Bezerra TN, Deprez T, Fonseca G, Holovachov O, Leduc D, Miljutin D, Moens T, Sharma J, Smol N, Tchesunov A, Mokievsky V, Vanaverbeke J, Vanreusel A, Vincx M. NeMys: World Database of Free-Living Marine Nematodes. 2016. Accessed at http://nemys.ugent.be on YYYY-MM-DD.
- Ingels J, Vanhove S, De Mesel I, Vanreusel A. The biodiversity and biogeography of the free-living nematode genera Desmodora and Desmodorella (family Desmodoridae) at both sides of the Scotia Arc. Polar Biol. 2006;29:936–49.View ArticleGoogle Scholar
- Ingels J, Hauquier F, Raes M, Vanreusel A. Chapter 5.3. Antarctic free-living marine nematodes. In: De Broyer C, Koubbi P, Griffiths HJ, Raymond B, d’Udekem d’Acoz C, et al., editors. Biogeographic Atlas of the Southern Ocean. Cambridge, Scientific Committee on Antarctic Research; 2014. p. 83–7. isbn:978-0-948277-28-3 http://atlas.biodiversity.aq/index.html.
- Wieser W. Die Beziehung zwischen Mundhöhlenge-stalt, Ernährungsweise und Vorkommen bei freilebenden marinen Nematoden. Arkiv för Zoologi. 1953;4:439–84.Google Scholar
- Hauquier F, Durán Suja L, Gutt J, Vanreusel A. Different oceanographic regimes in the vicinity of the Antarctic Peninsula reflected in benthic nematode communities. PLoS One. 2015;10(9):e0137527. doi:10.1371/journal.pone.0137527.View ArticlePubMedPubMed CentralGoogle Scholar
- Boeckner MJ, Sharma J, Proctor HC. Revisiting the meiofauna paradox: dispersal and colonization of nematodes and other meiofaunal organisms in low- and high-energy environments. Hydrobiologia. 2009;624:91–106.View ArticleGoogle Scholar
- Eskin RA, Palmer MA. Suspension of marine nematodes in a turbulent tidal creek: species patterns. Biol Bull. 1985;169:615–23.View ArticleGoogle Scholar
- Thomas MC, Lana PC. A new look into the small-scale dispersal of free-living marine nematodes. Zoologia. 2011;28:449–56.View ArticleGoogle Scholar
- Slatkin M. Isolation by distance in equilibrium and non-equilibrium populations. Evolution. 1993;47(1):264–79.View ArticleGoogle Scholar
- Wright S. Isolation by distance. Genetics. 1943;28:114–38.PubMedPubMed CentralGoogle Scholar
- Knust R, Gerdes D, Mintenbeck K, editors. The expedition of the research vessel “Polarstern” to the Antarctic in 2011 (ANT-XXVII/3) (CAMBIO). Reports on polar and marine research. Bremerhaven: Alfred Wegener Institute for Polar and Marine Research. 2012;644:202. hdl:10013/epic.39114Google Scholar
- Yoder M, De Ley IT, King I, Mundo-Ocampo M, Mann J, Blaxter M, Poiras L, De Ley P. Dess: a versatile solution for preserving morphology and extractable DNA of nematodes. Nematology. 2006;38:367–76.View ArticleGoogle Scholar
- Heip C, Vincx M, Vranken G. The ecology of marine nematodes. Oceanogr Mar Biol. 1985;23:399–489.Google Scholar
- Vincx M. Meiofauna in marine and freshwater sediments. In: Hall GS, editor. Methods for the examination of organismal diversity in soils and sediments. Wallinford: Cab International; 1996. p. 187–95.Google Scholar
- Williams BD, Schrank B, Huynh C, Shownkeen R, Waterston R. A genetic mapping system in Caenorhabditis elegans based on polymorphic sequence-tagged sites. Genetics. 1992;131:609–24.PubMedPubMed CentralGoogle Scholar
- Larkin MA, Blackshields G, Brown NP, Chenna R, McGettigan PA, McWilliam H, Valentin F, Wallace IM, Wilm A, Lopez R, Thompson JD, Gibson TJ, Higgins DG. Clustal W and Clustal X version 2.0. Bioinformatics. 2007;23:2947–8.View ArticlePubMedGoogle Scholar
- Tamura K, Stecher G, Peterson D, Filipski A, Kumar S. MEGA6: Molecular Evolutionary Genetics Analysis Version 6.0. Mol Biol Evol. 2013;30(12):2725–9.View ArticlePubMedPubMed CentralGoogle Scholar
- Darriba D, Taboada GL, Doallo R, Posada D. jModelTest 2: more models, new heuristics and parallel computing. Nat Methods. 2012;9(8):772.View ArticlePubMedPubMed CentralGoogle Scholar
- Guindon S, Gascuel O. A simple, fast and accurate method to estimate large phylogenies by maximum-likelihood. Syst Biol. 2003;52:696–704.View ArticlePubMedGoogle Scholar
- Stamatakis A. RAxML Version 8: A tool for Phylogenetic Analysis and Post-Analysis of Large Phylogenies. Bioinformatics. 2014;30(9):1312–1313. doi:10.1093/bioinformatics/btu033.View ArticlePubMedPubMed CentralGoogle Scholar
- Drummond AJ, Suchard MA, Xie D, Rambaut A. Bayesian phylogenetics with BEAUti and the BEAST 1.7. Mol Biol Evol. 2012;29:1969–73.View ArticlePubMedPubMed CentralGoogle Scholar
- Rambaut A, Suchard MA, Xie D, Drummond AJ. 2014. Tracer v1.6, available from http://tree.bio.ed.ac.uk/software/Tracer
- Miller MA, Pfeiffer W, Schwartz T. Creating the CIPRES Science Gateway for inference of large phylogenetic trees. New Orleans: Proceedings of the Gateway Computing Environments Workshop (GCE); 2010. p. 1–8.Google Scholar
- Pons J, Barraclough TG, Gomez-Zurita J, Cardoso A, Duran DP, Hazell S, Kamoun S, Sumlin WD, Vogler AP. Sequence-based species delimitation for the DNA taxonomy of undescribed insects. Syst Biol. 2006;55(4):595–609.View ArticlePubMedGoogle Scholar
- R Core Team. R: A language and environment for statistical computing. Vienna: R Foundation for Statistical Computing; 2013. http://www.R-project.org Google Scholar
- Paradis E, Claude J, Strimmer K. APE : analyses of phylogenetics and evolution in R language. Bioinformatics. 2004;20:289–90.View ArticlePubMedGoogle Scholar
- Ezard T, Fujisawa T, Barraclough T. SPLITS: SPecies’ LImits by Threshold Statistics. R package version 1.0–19/r48. 2013. http://R-Forge.R-project.org/projects/splits/.
- Templeton AR, Crandall KA, Sing CF. A cladistic analysis of phenotypic associations with haplotypes inferred from restriction endonuclease mapping and DNA sequence data. III. Cladogram estimation. Genetics. 1992;132:619–33.PubMedPubMed CentralGoogle Scholar
- Clement M, Posada D, Crandall K. TCS: a computer program to estimate gene genealogies. Mol Ecol. 2000;9(10):1657–60.View ArticlePubMedGoogle Scholar
- Clement M, Snell Q, Walke P, Posada D, Crandall K. TCS: estimating gene genealogies. Fort Lauderdale: Proc 16th Int Parallel Distrib Process Symp. 2002;2:184.Google Scholar
- Excoffier L, Lischer HEL. Arlequin suite ver 3.5: A new series of programs to perform population genetics analyses under Linux and Windows. Mol Ecol Resour. 2010;10:564–7.View ArticlePubMedGoogle Scholar
- Holsinger KE, Weir BS. Genetics in geographically structured populations: defining, estimating and interpreting F st. Nat Rev Genet. 2009;10(9):639–50. doi:10.1038/nrg2611.View ArticlePubMedPubMed CentralGoogle Scholar
- Nei M. Molecular evolutionary genetics. New York: Columbia University Press; 1987. p. 512.Google Scholar
- Tajima F. Evolutionary relationship of DNA sequences in finite populations. Genetics. 1983;105(2):437–60.PubMedPubMed CentralGoogle Scholar
- Jensen JL, Bohonak AJ, Kelley ST. Isolation by distance, web service. BMC Genet. 2005;6:13. http://ibdws.sdsu.edu/ View ArticlePubMedPubMed CentralGoogle Scholar
- Fujita MK, Leaché AD, Burbrink FT, McGuire JA, Moritz C. Coalescent-based species delimitation in an integrative taxonomy. Trends Ecol Evol. 2012;27:480–8.View ArticlePubMedGoogle Scholar
- Velasco-Castrillón A, Stevens MI. Morphological and molecular diversity at a regional scale: A step closer to understanding Antarctic nematode biogeography. Soil Biol Biochem. 2014;70:272–84.View ArticleGoogle Scholar
- Hemery LG, Eléaume M, Roussel V, Améziane N, Gallut C, Steinke D, Cruaud C, Couloux A, Wilson NG. Comprehensive sampling reveals circumpolarity and sympatry in seven mitochondrial lineages of the Southern Ocean crinoid species Promachocrinus kerguelensis (Echinodermata). Mol Ecol. 2012;21:2502–18.View ArticlePubMedGoogle Scholar
- Miljutin DM, Miljutina MA. Intraspecific variability of morphological characters in the species-rich deep-sea genus Acantholaimus Allgén, 1933 (Nematoda: Chromadoridae). Nematology. 2016;18(4):455–73. doi:10.1163/15685411-00002970.View ArticleGoogle Scholar
- Fonderie P, Steel H, Moens T, Bert W. Experimental induction of intraspecific morphometric variability in a single population of Halicephalobus cf. gingivalis may surpass total interspecific variability. Nematology. 2013;15(5):529–44.View ArticleGoogle Scholar
- Kiontke K, Fitch DHA. Phenotypic plasticity: different teeth for different feasts. Curr Biol. 2011;20(17):710–2. doi:10.1016/j.cub.2010.07.009.View ArticleGoogle Scholar
- Havermans C, Sonet G, d’Udekem d’Acoz C, Nagy ZT, Martin P, Brix S, Riehl T, Agrawal S, Held C. Genetic and morphological divergences in the cosmopolitan deep-sea amphipod Eurythenes gryllus reveal a diverse abyss and a bipolar species. PLOS ONE. 2013;8(9):e74218. doi:10.1371/journal.pone.0074218.
- Jörger KM, Norenberg JL, Wilson NG, Schrödl M. Barcoding against a paradox? Combined molecular species delineations reveal multiple cryptic lineages in elusive meiofaunal sea slugs. BMC Evol Biol. 2012;12:245.View ArticlePubMedPubMed CentralGoogle Scholar
- Hellberg ME, Burton RS, Neigel JE, Palumbi SR. Genetic assessment of connectivity among marine populations. Bull Mar Sci. 2002;70(1):273–90.Google Scholar
- Pelc RA, Warner RR, Gaines SD. Geographical patterns of genetic structure in marine species with contrasting life histories. J Biogeogr. 2009;36:1881–90.View ArticleGoogle Scholar
- Rogers AD. Evolution and Biodiversity of Antarctic Organisms. In: Rogers AD, Johnston NM, Murphy EJ, Clarke A, editors. Antarctic Ecosystems: An Extreme Environment in a Changing World. Chichester: John Wiley & Sons, Ltd; 2012. doi:10.1002/9781444347241.ch14.View ArticleGoogle Scholar
- Thatje S, Hillenbrand C-L, Larter R. On the origin of Antarctic marine benthic community structure. Trends Ecol Evol. 2005;20(10):534–40.View ArticlePubMedGoogle Scholar
- Griffiths HJ, Barnes DKA, Linse K. Towards a generalized biogeography of the Southern Ocean benthos. J Biogeogr. 2009;36:162–77.View ArticleGoogle Scholar
- Palumbi SR. Genetic divergence, reproductive isolation and marine speciation. Annu Rev Ecol Syst. 1994;25:547–72.View ArticleGoogle Scholar
- Marshall DJ, Monro K, Bode M, Keough MJ, Swearer S. Phenotype-environment mismatches reduce connectivity in the sea. Ecol Lett. 2010;13:128–40.View ArticlePubMedGoogle Scholar
- Selkoe KA, Watson JR, White C, Ben Horin T, Iacchei M, Mitarai S, Siegel DA, Gaines SD, Toonen RJ. Taking the chaos out of genetic patchiness: seascape genetics reveal ecological and oceanographic drivers of genetic patterns in three temperate reef species. Mol Ecol. 2010;19:3708–26.View ArticlePubMedGoogle Scholar
- White C, Selkoe KA, Watson J, Siegel DA, Zacherl DC, Toonen RJ. Ocean currents help explain population genetic structure. Proc R Soc B. 2010; doi:10.1098/rspb.2009.2214.
- Selkoe KA, D’Aloia CC, Crandall ED, Iacchei M, Liggins L, Puritz JB, von der Heyden S, Toonen RJ. A decade of seascape genetics: contributions to basic and applied marine connectivity. Mar Ecol Prog Ser. 2016;554:1–19. doi:10.3354/meps11792.View ArticleGoogle Scholar
- Commito JA, Tita G. Differential dispersal rates in an intertidal meiofauna assemblage. J Exp Mar Biol Ecol. 2002;268:237–56.View ArticleGoogle Scholar
- Bradbury IR, Laurel B, Snelgrove PVR, Bentzen P, Campana SE. Global patterns in marine dispersal estimates: the influence of geography, taxonomic category and life history. Proc R Soc B. 2008;275:1803–9.View ArticlePubMedPubMed CentralGoogle Scholar
- Kimura M. A simple method for estimating evolutionary rate of base substitutions through comparative studies of nucleotide sequences. J Mol Evol. 1980;16:111–20.View ArticlePubMedGoogle Scholar
- Hasegawa M, Kishino K, Yano T. Dating the human-ape splitting by a molecular clock of mitochondrial DNA. J Mol Evol. 1985;22:160–74.View ArticlePubMedGoogle Scholar
- Zharkikh A. Estimation of evolutionary distances between nucleotide sequences. J Mol Evol. 1994;39:315–29.View ArticlePubMedGoogle Scholar
- Tavaré S. Some probabilistic and statistical problems in the analysis of DNA sequences. Some mathematical questions in biology - DNA sequence analysis. Providence, RI. Amer Math Soc. 1986;17:57–86.Google Scholar
- Altschul SF, Madden TL, Schäffer AA, Zhang J, Zhang Z, Miller W, Lipman DJ. Gapped BLAST and PSI-BLAST: a new generation of protein database search programs. Nucleic Acid Res. 1997;25:3389–402.View ArticlePubMedPubMed CentralGoogle Scholar
- Bhadury P, Austen MC, Bilton DT, Lambshead PJD, Rogers AD, Smerdon GR. Evaluation of combined morphological and molecular techniques for marine nematode (Terschellingia spp.) identification. Mar Biol. 2008;154:509–18.View ArticleGoogle Scholar
- Bhadury P, Austen MC, Bilton DT, Lambshead PJD, Rogers AD, Smerdon GR. Development and evaluation of a DNA-barcoding approach for the rapid identification of nematodes. Mar Ecol Prog Ser. 2006;320:1–9.View ArticleGoogle Scholar
- Bik HM, Hawkins LE, Hughes JA, Lambshead PJD. Rapid decline of PCR amplification from genomic extracts of DESS-preserved, slide-mounted nematodes. Nematology. 2009;11(6):827–34.View ArticleGoogle Scholar
- Bhadury P, Bik H, Lambshead PJD, Austen MC, Smerdon GR, Rogers AD. Molecular Diversity of Fungal Phylotypes Co-Amplified Alongside Nematodes from Coastal and Deep-Sea Marine Environments. PLoS One. 2011;6(10):e26445. doi:10.1371/journal.pone.0026445.View ArticlePubMedPubMed CentralGoogle Scholar