Dispersal and gene flow in free-living marine nematodes

Dispersal and gene flow determine connectivity among populations, and can be studied through population genetics and phylogeography. We here review the results of such a framework for free-living marine nematodes. Although field experiments have illustrated substantial dispersal in nematodes at ecological time scales, analysis of the genetic diversity illustrated the importance of priority effects, founder effects and genetic bottlenecks for population structuring between patches <1 km apart. In contrast, only little genetic structuring was observed within an estuary (<50 km), indicating that these small scale fluctuations in genetic differentiation are stabilized over deeper time scales through extensive gene flow. Interestingly, nematode species with contrasting life histories (extreme colonizers vs persisters) or with different habitat preferences (algae vs sediment) show similar, low genetic structuring. Finally, historical events have shaped the genetic pattern of marine nematodes and show that gene flow is restricted at large geographical scales. We also discuss the presence of substantial cryptic diversity in marine nematodes, and end with highlighting future important steps to further unravel nematode evolution and diversity

Salinity effects on competition between cryptic species of the nematode Rhabditis (Pellioditis) marina

Behind the morphological similarity of many species, a substantial hidden genetic diversity can be found. This cryptic diversity has been well documented in the marine nematode, Rhabditis (Pellioditis) marina, where several cryptic species occur sympatrically. Despite the growing knowledge about its cryptic diversity, little is known about potential differences in the ecology and functional roles of the different species. In this experiment four cryptic species of Rhabditis (Pellioditis) marina were reared together, starting from identical initial abundances at two different salinity concentrations to investigate differences in their environmental preferences and how they may affect the outcome of their competitive interactions. Every fifth day of the experiment one third of the adults was removed from the population for analyzing the genetic composition of the population; adult and juvenile population dynamics were also examined. The first results show a difference in population dynamics dependent on the salinity. At low salinity, a population crash occurred after 15 days in all the populations. Only half of the populations recovered from this bottleneck and reached higher abundances than before. Genetic analyses revealed that the composition of these populations consists solely of individuals of one cryptic species. This effect was totally absent in the populations reared at higher salinity. Two explanations are possible: (1) the different cryptic species have different survival rates at different salinities and (2) the competitive interactions between them differ at different salinities. We are currently identifying the sampled adults using a restriction fragment analysis to assess these explanations and to elucidate the effect of salinity on the coexistence of cryptic species of Rhabditis (Pellioditis) marina

Differences in life-histories refute ecological equivalence of cryptic species and provide clues to the origin of bathyal Halomonhystera (Nematoda)

The discovery of morphologically very similar but genetically distinct species complicates a proper understanding of the link between biodiversity and ecosystem functioning. Cryptic species have been frequently observed to co-occur and are thus expected to be ecological equivalent. The marine nematode Halomonhystera disjuncta contains five cryptic species (GD1-5) that co-occur in the Westerschelde estuary. In this study, we investigated the effect of three abiotic factors (salinity, temperature and sulphide) on life-history traits of three cryptic H. disjuncta species (GD1-3). Our results show that temperature had the most profound influence on all life-cycle parameters compared to a smaller effect of salinity. Life-history traits of closely related cryptic species were differentially affected by temperature, salinity and presence of sulphides which shows that cryptic H. disjuncta species are not ecologically equivalent. Our results further revealed that GD1 had the highest tolerance to a combination of sulphides, high salinities and low temperatures. The close phylogenetic position of GD1 to Halomonhystera hermesi, the dominant species in sulphidic sediments of the Hakon Mosby mud volcano (Barent Sea, 1280 m depth), indicates that both species share a recent common ancestor. Differential life-history responses to environmental changes among cryptic species may have crucial consequences for our perception on ecosystem functioning and coexistence of cryptic species

Transcription, signaling receptor activity, oxidative phosphorylation, and fatty acid metabolism mediate the presence of closely related species in distinct intertidal and cold-seep habitats

Bathyal cold seeps are isolated extreme deep-sea environments characterized by low species diversity while biomass can be high. The Hakon Mosby mud volcano (Barents Sea, 1,280 m) is a rather stable chemosynthetic driven habitat characterized by prominent surface bacterial mats with high sulfide concentrations and low oxygen levels. Here, the nematode Halomonhystera hermesithrives in high abundances (11,000 individuals 10 cm(-2)). Halomonhystera hermesi is a member of the intertidal Halomonhystera disjuncta species complex that includes five cryptic species (GD 1-5). GD1-5's common habitat is characterized by strong environmental fluctuations. Here, we compared the transcriptomes of H. hermesi and GD1, H. hermesi's closest relative. Genes encoding proteins involved in oxidative phosphorylation are more strongly expressed in H. hermesi than in GD1, and many genes were only observed in H. hermesi while being completely absent in GD1. Both observations could in part be attributed to high sulfide concentrations and low oxygen levels. Additionally, fatty acid elongation was also prominent in H. hermesi confirming the importance of highly unsaturated fatty acids in this species. Significant higher amounts of transcription factors and genes involved in signaling receptor activity were observed in GD1 (many of which were completely absent in H. hermesi), allowing fast signaling and transcriptional reprogramming which can mediate survival in dynamic intertidal environments. GC content was approximately 8% higher in H. hermesi coding unigenes resulting in differential codon usage between both species and a higher proportion of amino acids with GC-rich codons in H. hermesi. In general our results showed that most pathways were active in both environments and that only three genes are under natural selection. This indicates that also plasticity should be taken in consideration in the evolutionary history of Halomonhystera species. Such plasticity, as well as possible preadaptation to low oxygen and high sulfide levels might have played an important role in the establishment of a cold-seep Halomonhystera population

Phylogenetic clustering and rarity imply risk of local species extinction in prospective deep-sea mining areas of the Clarion-Clipperton Fracture Zone

An understanding of the forces controlling community structure in the deep sea is essential at a time when its pristineness is threatened by polymetallic nodule mining. Because abiotically defined communities are more sensitive to environmental change, we applied occurrence- and phylogeny-based metrics to determine the importance of biotic versus abiotic structuring processes in nematodes, the most abundant invertebrate taxon of the Clarion-Clipperton Fracture Zone (CCFZ), an area targeted for mining. We investigated the prevalence of rarity and the explanatory power of environmental parameters with respect to phylogenetic diversity (PD). We found evidence for aggregation and phylogenetic clustering in nematode amplicon sequence variants (ASVs) and the dominant genus Acantholaimus, indicating the influence of environmental filtering, sympatric speciation, affinity for overlapping habitats and facilitation for community structure. PD was associated with abiotic variables such as total organic carbon, chloroplastic pigments equivalents and/or mud content, explaining up to 57% of the observed variability and providing further support of the prominence of environmental structuring forces. Rarity was high throughout, ranging from 64 to 75% unique ASVs. Communities defined by environmental filtering with a prevalence of rarity in the CCFZ suggest taxa of these nodule-bearing abyssal plains will be especially vulnerable to the risk of extinction brought about by the efforts to extract them

Low genetic but high morphological variation over more than 1000 km coastline refutes omnipresence of cryptic diversity in marine nematodes

Background: The resilience of ecosystems to negative impacts is generally higher when high gene flow, species diversity and genetic diversity are present. Population genetic studies are suitable to investigate genetic diversity and estimate gene flow between populations. Seaweed beds form a dynamic shallow water ecosystem influenced by climate change and human exploitation, as such, seaweed beds are a particularly powerful model to investigate ecosystem resilience in coastal areas. We studied the population genetic structure of the new nematode species Paracanthonchus gynodiporata associated with seaweeds in northeastern Brazil. Nematodes are generally believed to have a limited dispersal capacity because of the lack of planktonic larvae. Yet, they can drift on seaweeds, and water currents might be a natural barrier for their dispersal. Populations of P. gynodiporata were sampled over more than 1000 km coastline in regions across major oceanic currents with and without historical exploitation of seaweed. Results: P. gynodiporata is described in an integrative way using mitochondrial and nuclear sequences and morphological data. The 3D model of the head region shows for the first time a detailed view of the ventrosublateral teeth, a character often overlooked in older taxonomic studies of the genus. A total of 17 mitochondrial COI haplotypes were found with one haplotype representing 63 to 83% of the frequencies in each population. AMOVA showed overall little population genetic structure (F-ST = 0.05204), and no genetic subdivision between the populations under the influence of the two different water currents were found. Effects of historical seaweed exploitation on population genetic diversity were not detected. In contrast, significant differences between populations were found in morphometric characters. This discrepancy in genetic and morphological differentiation between populations across 1000 km of coastline is surprising in view of the frequently observed presence of several cryptic species at small geographical scale in other macroalgal associated nematodes. Conclusions: Our results show that cryptic species are not omnipresent in marine nematode species, suggesting that nematodes associated with seaweeds have been able to disperse over large distances across well-known biogeographic barriers

Population genomic structure of the gelatinous zooplankton species Mnemiopsis leidyi in its nonindigenous range in the North Sea

Nonindigenous species pose a major threat for coastal and estuarine ecosystems. Risk management requires genetic information to establish appropriate management units and infer introduction and dispersal routes. We investigated one of the most successful marine invaders, the ctenophore Mnemiopsis leidyi, and used genotyping-by-sequencing (GBS) to explore the spatial population structure in its nonindigenous range in the North Sea. We analyzed 140 specimens collected in different environments, including coastal and estuarine areas, and ports along the coast. Single nucleotide polymorphisms (SNPs) were called in approximately 40 k GBS loci. Population structure based on the neutral SNP panel was significant (F-ST .02; p < .01), and a distinct genetic cluster was identified in a port along the Belgian coast (Ostend port; pairwise F-ST .02-.04; p < .01). Remarkably, no population structure was detected between geographically distant regions in the North Sea (the Southern part of the North Sea vs. the Kattegat/Skagerrak region), which indicates substantial gene flow at this geographical scale and recent population expansion of nonindigenous M. leidyi. Additionally, seven specimens collected at one location in the indigenous range (Chesapeake Bay, USA) were highly differentiated from the North Sea populations (pairwise F-ST .36-.39; p < .01). This study demonstrates the utility of GBS to investigate fine-scale population structure of gelatinous zooplankton species and shows high population connectivity among nonindigenous populations of this recently introduced species in the North Sea. OPEN RESEARCH BADGES This article has earned an Open Data Badge for making publicly available the digitally-shareable data necessary to reproduce the reported results. The data is available at: The DNA sequences generated for this study are deposited in the NCBI sequence read archive under SRA accession numbers -, and will be made publically available upon publication of this manuscript