Nematodes alter the taxonomic and functional profiles of benthic bacterial communities: A metatranscriptomic approach.

Marine sediments cover 70% of the Earth's surface, and harbour diverse bacterial communities critical for marine biogeochemical processes, which affect climate change, biodiversity and ecosystem functioning. Nematodes, the most abundant and species-rich metazoan organisms in marine sediments, in turn, affect benthic bacterial communities and bacterial-mediated ecological processes, but the underlying mechanisms by which they affect biogeochemical cycles remain poorly understood. Here, we demonstrate using a metatranscriptomic approach that nematodes alter the taxonomic and functional profiles of benthic bacterial communities. We found particularly strong stimulation of nitrogen-fixing and methane-oxidizing bacteria in the presence of nematodes, as well as increased functional activity associated with methane metabolism and degradation of various carbon compounds. This study provides empirical evidence that the presence of nematodes results in taxonomic and functional shifts in active bacterial communities, indicating that nematodes may play an important role in benthic ecosystem processes.

Environmental filtering along a bathymetric gradient: A metabarcoding meta-analysis of free-living nematodes.

Identifying and understanding patterns of biological diversity is crucial at a time when even the most remote and pristine marine ecosystems are threatened by resource exploitation such as deep-seabed mining. Metabarcoding provides the means through which one can perform comprehensive investigations of diversity by examining entire assemblages simultaneously. Nematodes commonly represent the most abundant infaunal metazoan group in marine soft sediments. In this meta-analysis, we compiled all publicly available metabarcoding datasets targeting the 18S rRNA v1-v2 region from sediment samples to conduct a global-scale examination of nematode amplicon sequence variant (ASV) alpha diversity patterns and phylogenetic community structure at different depths and habitats. We found that nematode ASV richness followed a parabolic trend, increasing from the intertidal to the shelf, reaching a maximum in the bathyal and decreasing in the abyssal zone. No depth- or habitat-specific assemblages were identified as a large fraction of genera were shared. Contrastingly, the vast majority of ASVs were unique to each habitat and/or depth zone; genetic diversity was thus highly localized. Overwhelmingly, nematode ASVs in all habitats exhibited phylogenetic clustering, pointing to environmental filtering as the primary force defining community assembly rather than competitive interactions. This finding stresses the importance of habitat preservation for the maintenance of marine nematode diversity.

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.

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

BACKGROUND: Dispersal ability, population genetic structure and species divergence in marine nematodes are still poorly understood, especially in remote areas such as the Southern Ocean. We investigated genetic differentiation of species and populations of the free-living endobenthic nematode genera Sabatieria and Desmodora using nuclear 18S rDNA, internal transcribed spacer (ITS) rDNA, and mitochondrial cytochrome oxidase I (COI) gene sequences. Specimens were collected at continental shelf depths (200-500 m) near the Antarctic Peninsula, Scotia Arc and eastern side of the Weddell Sea. The two nematode genera co-occurred at all sampled locations, but with different vertical distribution in the sediment. A combination of phylogenetic (GMYC, Bayesian Inference, Maximum Likelihood) and population genetic (AMOVA) analyses were used for species delimitation and assessment of gene flow between sampling locations. RESULTS: Sequence analyses resulted in the delimitation of four divergent species lineages in Sabatieria, two of which could not be discriminated morphologically and most likely constitute cryptic species. Two species were recognised in Desmodora, one of which showed large intraspecific morphological variation. Both genera comprised species that were restricted to one side of the Weddell Sea and species that were widely spread across it. Population genetic structuring was highly significant and more pronounced in the deeper sediment-dwelling Sabatieria species, which are generally less prone to resuspension and passive dispersal in the water column than surface Desmodora species. CONCLUSIONS: Our results indicate that gene flow is restricted at large geographic distance in the Southern Ocean, which casts doubt on the efficiency of the Weddell gyre and Antarctic Circumpolar Current in facilitating circum-Antarctic nematode species distributions. We also show that genetic structuring and cryptic speciation can be very different in nematode species isolated from the same geographic area, but with different habitat preferences (surface versus deeper sediment layers).

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.

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.

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.

Genetic structure of Xiphinema pachtaicum and X. index populations based on mitochondrial DNA variation.

The dagger nematodes Xiphinema pachtaicum and X. index are two of the most widespread and frequently occurring Xiphinema spp. co-infesting vineyards and other crops and natural habitats worldwide. Sexual reproduction is rare in these species. The primary objective of this study was to determine the genetic structure of X. pachtaicum and X. index populations using eight and seven populations, respectively, from different "wine of denomination of origin (D.O.) zones" in Spain and Sardinia (Italy), by studying mitochondrial (cytochrome oxidase c subunit 1 or COI) and nuclear (D2-D3 expansion segments of 28S rDNA) markers. Both Xiphinema spp. showed low intraspecific divergence among COI sequences, ranging from 0.2% (1 base substitution) to 2.3% (10 substitutions) in X. pachtaicum and from 0.2% (1 base substitution) to 0.4% (2 substitutions) in X. index. Population genetic structure was strong for both species. Nevertheless, molecular differences among grapevine-growing areas were not significant, and intrapopulation diversity was very low. It is hypothesized that this genetic homogeneity in the nematode populations reflects their predominant parthenogenetic reproduction mode and low dispersal abilities. Our results also show that X. pachtaicum populations in Spain have possibly been established from two different populations of origin. Results also demonstrated that the two DNA regions studied are suitable diagnostic markers for X. index and X. pachtaicum.

Colonization of macroalgal deposits by estuarine nematodes through air and potential for rafting inside algal structures.

Dispersal is an important life-history trait. In marine meiofauna, and particularly in nematodes, dispersal is generally considered to be mainly passive, i.e. through transport with water currents and bedload transport. Because nematodes have no larval dispersal stage and have a poor swimming ability, their per capita dispersal capacity is expected to be limited. Nevertheless, many marine nematode genera and even species have near-cosmopolitan distributions, and at much smaller spatial scales, can rapidly colonise new habitat patches. Here we demonstrate that certain marine nematodes, like the morphospecies Litoditis marina, can live inside macroalgal structures such as receptacula and-to a lesser extent-floating bladders, which may allow them to raft over large distances with drifting macroalgae. We also demonstrate for the first time that these nematodes can colonize new habitat patches, such as newly deposited macroalgal wrack in the intertidal, not only through seawater but also through air. Our experimental set-up demonstrates that this aerial transport is probably the result of hitchhiking on vectors such as insects, which visit, and move between, the patches of deposited algae. Transport by wind, which has been observed for terrestrial nematodes and freshwater zooplankton, could not be demonstrated. These results can be important for our understanding of both large-scale geographic distribution patterns and of the small-scale colonization dynamics of habitat patches by marine nematodes.

Identification and semi-quantification of Atlantic salmon in processed and mixed seafood products using Droplet Digital PCR (ddPCR).

Fishery products are often subject to substitution fraud, which is hard to trace due to a lack of morphologic traits when processed, gutted, or decapitated. Traditional molecular methods (DNA barcoding) fail to identify products containing multiple species and cannot estimate original weight percentages. As a proof of concept, an Atlantic salmon (Salmo salar) specific ddPCR assay was designed to authenticate mixed food products. The method proved to be specific and able to accurately quantify S. salar when using DNA extracts, even in the presence of DNA from closely related salmon species. The ddPCR estimates correlated well with the percentage of S. salar in artificially assembled tissue mixtures. The effect of common salmon processing techniques (freezing, smoking, poaching with a "Bellevue" recipe and marinating with a 'Gravad lax' recipe) on the ddPCR output was investigated and freezing and marinating appeared to lower the copies detected by the ddPCR. Finally, the assay was applied to 46 retail products containing Atlantic or Pacific salmon, and no indications of substitution fraud were detected. The method allows for a semi-quantitative evaluation of the S. salar content in processed food products and can rapidly screen Atlantic salmon products and flag potentially tampered samples for further investigation.

A fast HRMA tool to authenticate eight salmonid species in commercial food products.

Atlantic and Pacific salmon are frequently consumed species with very different economic values: farmed Atlantic salmon is cheaper than wild-caught Pacific salmons. Species replacements occur with the high valued Pacific species (Oncorhynchus keta, O. gorbuscha, O. kisutch, O. nerka and O. tshawytscha) substituted by cheaper farmed Atlantic salmon (Salmo salar) and Atlantic salmon by rainbow trout (Oncorhynchus mykiss) and brown trout (Salmo trutta). Here we use High-Resolution Melting Analysis (HRMA) to identify eight salmonid species. We designed primers to generate short amplicons of 72 and 116 bp from the fish barcode genes CO1 and CYTB. The time of analysis was under 70 min, after DNA extraction. Food processing of Atlantic salmon (fresh, "Bellevue", "gravadlax", frozen and smoked) did not impact the HRMA profiles allowing reliable identification. A blind test was conducted by three different institutes, showing correct species identifications irrespective of the laboratory conducting the analysis. Finally, a total of 82 retail samples from three European countries were analyzed and a low substitution rate of 1.2% was found. The developed tool provides a quick way to investigate salmon fraud and contributes to safeguard consumers.

Species boundaries and phylogenetic relationships between Atlanto-Mediterranean shallow-water and deep-sea coral associated Hexadella species (Porifera, Ianthellidae).

Coral reefs constitute the most diverse ecosystem of the marine realm and an increasing number of studies are focusing on coral species boundaries, distribution, and on processes that control species ranges. However, less attention has been paid to coral associated species. Deep-sea sponges dominate cold-water coral ecosystems, but virtually nothing is known about their molecular diversity. Moreover, species boundaries based on morphology may sometimes be inadequate, since sponges have few diagnostic characters. In this study, we investigated the molecular diversity within the genus Hexadella (Porifera, Demospongiae, Verongida, Ianthellidae) from the European shallow-water environment to the deep-sea coral ecosystems. Three molecular markers were used: one mitochondrial (COI) and two nuclear gene fragments (28S rDNA and the ATPS intron). Phylogenetic analyses revealed deeply divergent deep-sea clades congruent across the mitochondrial and nuclear markers. One clade contained specimens from the Irish, the Scottish, and the Norwegian margins and the Greenland Sea (Hexadella dedritifera) while another clade contained specimens from the Ionian Sea, the Bay of Biscay, and the Irish margin (H. cf. dedritifera). Moreover, these deeply divergent deep-sea clades showed a wide distribution suggesting a connection between the reefs. The results also point to the existence of a new deep-sea species (Hexadella sp.) in the Mediterranean Sea and of a cryptic shallow-water species (Hexadella cf. pruvoti) in the Gorringe Bank. In contrast, low genetic differentiation between H. cf. dedritifera and H. pruvoti from the Mediterranean Sea was observed. All Hexadella racovitzai specimens from the Mediterranean Sea (shallow and deep) to the Atlantic formed a monophyletic group.

A high-quality genetic reference database for European commercial fishes reveals substitution fraud of processed Atlantic cod (Gadus morhua) and common sole (Solea solea) at different steps in the Belgian supply chain.

Seafood is an important component of the human diet. With depleting fish stocks and increasing prices, seafood is prone to fraudulent substitution. DNA barcoding has illustrated fraudulent substitution of fishes in retail and restaurants. Whether substitution also occurs in other steps of the supply chain remains largely unknown. DNA barcoding relies on public reference databases for species identification, but these can contain incorrect identifications. The creation of a high quality genetic reference database for 42 European commercially important fishes was initiated containing 145 Cytochrome c oxidase subunit I (COI) and 152 Cytochrome b (cytB) sequences. This database was used to identify substitution rates of Atlantic cod (Gadus morhua) and common sole (Solea solea) along the fish supply chain in Belgium using DNA barcoding. Three out of 132 cod samples were substituted, in catering (6%), import (5%) and fishmongers (3%). Seven out of the 41 processed sole samples were substituted, in wholesale (100%), food services (50%), retailers (20%) and catering (8%). Results show that substitution of G. morhua and S. solea is not restricted to restaurants, but occurs in other parts of the supply chain, warranting for more stringent controls along the supply chain to increase transparency and trust among consumers.

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 SRR6950721-SRR6950884, and will be made publically available upon publication of this manuscript.

Metabarcoding free-living marine nematodes using curated 18S and CO1 reference sequence databases for species-level taxonomic assignments.

High-throughput sequencing has the potential to describe biological communities with high efficiency yet comprehensive assessment of diversity with species-level resolution remains one of the most challenging aspects of metabarcoding studies. We investigated the utility of curated ribosomal and mitochondrial nematode reference sequence databases for determining phylum-specific species-level clustering thresholds. We compiled 438 ribosomal and 290 mitochondrial sequences which identified 99% and 94% as the species delineation clustering threshold, respectively. These thresholds were evaluated in HTS data from mock communities containing 39 nematode species as well as environmental samples from Vietnam. We compared the taxonomic description of the mocks generated by two read-merging and two clustering algorithms and the cluster-free Dada2 pipeline. Taxonomic assignment with the RDP classifier was assessed under different training sets. Our results showed that 36/39 mock nematode species were identified across the molecular markers (18S: 32, JB2: 19, JB3: 21) in UClust_ref OTUs at their respective clustering thresholds, outperforming UParse_denovo and the commonly used 97% similarity. Dada2 generated the most realistic number of ASVs (18S: 83, JB2: 75, JB3: 82), collectively identifying 30/39 mock species. The ribosomal marker outperformed the mitochondrial markers in terms of species and genus-level detections for both OTUs and ASVs. The number of taxonomic assignments of OTUs/ASVs was highest when the smallest reference database containing only nematode sequences was used and when sequences were truncated to the respective amplicon length. Overall, OTUs generated more species-level detections, which were, however, associated with higher error rates compared to ASVs. Genus-level assignments using ASVs exhibited higher accuracy and lower error rates compared to species-level assignments, suggesting that this is the most reliable pipeline for rapid assessment of alpha diversity from environmental samples.

Living apart-together: Microhabitat differentiation of cryptic nematode species in a saltmarsh habitat.

Coexistence of highly similar species is at odds with ecological theory of competition; coexistence, then, requires stabilizing mechanisms such as differences in ecological niche. In the bacterivore nematode Litoditis marina species complex, which occurs associated with macro-algae, four cryptic lineages (Pm I-IV) co-occur in the field along the south-western coast and estuaries of The Netherlands. Here we investigate the temporal and/or spatial niche differentiation in their natural environment using a qPCR-based detection and relative quantification method. We collected different algal species (i.e. two Fucus species and Ulva sp.) and separated algal structures (i.e. receptacula, thalli, non-fertile tips and bladders) at different sampling months and times (i.e. twice per sampling month), to examine differences in microhabitat use between coexisting L. marina species. Results demonstrate that the cryptic species composition varied among different algal species and algal structures, which was also subject to temporal shifts. Pm I dominated on Fucus spp., Pm II showed dominance on Ulva sp., while Pm III overall had the lowest frequencies. Microhabitat partitioning was most pronounced between the two cryptic species which had similar microbiomes (Pm I and Pm II), and less so between the two species which had significantly different microbiomes (Pm I and Pm III), suggesting that species which share the same microhabitats may avoid competition through resource partitioning. The interplay of microhabitat differentiation and temporal dynamics among the cryptic species of L. marina implies that there is a complex interaction between biotic components and abiotic factors which contributes to their coexistence in the field.

The effect of bio-irrigation by the polychaete Lanice conchilega on active denitrifiers: Distribution, diversity and composition of nosZ gene.

The presence of large densities of the piston-pumping polychaete Lanice conchilega can have important consequences for the functioning of marine sediments. It is considered both an allogenic and an autogenic ecosystem engineer, affecting spatial and temporal biogeochemical gradients (oxygen concentrations, oxygen penetration depth and nutrient concentrations) and physical properties (grain size) of marine sediments, which could affect functional properties of sediment-inhabiting microbial communities. Here we investigated whether density-dependent effects of L. conchilega affected horizontal (m-scale) and vertical (cm-scale) patterns in the distribution, diversity and composition of the typical nosZ gene in the active denitrifying organisms. This gene plays a major role in N2O reduction in coastal ecosystems as the last step completing the denitrification pathway. We showed that both vertical and horizontal composition and richness of nosZ gene were indeed significantly affected when large densities of the bio-irrigator were present. This could be directly related to allogenic ecosystem engineering effects on the environment, reflected in increased oxygen penetration depth and oxygen concentrations in the upper cm of the sediment in high densities of L. conchilega. A higher diversity (Shannon diversity and inverse Simpson) of nosZ observed in patches with high L. conchilega densities (3,185-3,440 ind. m-2) at deeper sediment layers could suggest a downward transport of NO3- to deeper layers resulting from bio-irrigation as well. Hence, our results show the effect of L. conchilega bio-irrigation activity on denitrifying organisms in L. conchilega reefs.

Limited possibilities for prezygotic barriers in the reproductive behaviour of sympatric Ophthalmotilapia species (Teleostei, Cichlidae).

Since prezygotic rather than postzygotic barriers are believed to maintain the diversity of closely related sympatric cichlids, differences in phenotypic traits and reproductive behaviours are likely involved in maintaining species boundaries. Here, we focused on the reproductive behaviour of three Ophthalmotilapia species with distributions that only overlap on a small stretch of the shore line of Lake Tanganyika. Repeated introgression of mitochondrial DNA between these species was previously reported, which suggested they can hybridise. Our aim is to test the hypothesis that reproductive behaviour acts as a prezygotic barrier that prevents frequent hybridisation in sympatric Ophthalmotilapia species. We performed a quantitative analysis of twelve reproductions (four for O. ventralis, six for O. nasuta, one for O. boops, and one between a female O. ventralis and a male O. nasuta). Although similar ethograms were obtained for these reproductions, the O. ventralis and O. boops males displayed a behaviour that was never performed by O. nasuta males. This behaviour was displayed during courtship and we called it 'invite'. In O. ventralis, we could show that it was associated with the emission of a single pulse sound. The comparison of O. nasuta and O. ventralis reproductive behaviours also revealed some quantitative differences: O. ventralis males showed the location of the bower more often to the female, whereas O. ventralis females followed the male more often. The similarity between the reproductive behaviours in O. ventralis and O. nasuta could explain the occurrence of the heterospecific spawning event recorded between an O. nasuta male and an O. ventralis female. Importantly, few eggs were laid and the maternal mouthbrooding that resulted from this heterospecific reproduction only lasted for two days, which suggested the abortion of egg development. Hence, in the absence of conspecifics, courtship and mating behaviours alone do not constitute perfect prezygotic barriers between these two species.

Neurogenomic Profiling Reveals Distinct Gene Expression Profiles Between Brain Parts That Are Consistent in Ophthalmotilapia Cichlids.

The detection of external and internal cues alters gene expression in the brain which in turn may affect neural networks that underly behavioral responses. Previous studies have shown that gene expression profiles differ between major brain regions within individuals and between species with different morphologies, cognitive abilities and/or behaviors. A detailed description of gene expression in all macroanatomical brain regions and in species with similar morphologies and behaviors is however lacking. Here, we dissected the brain of two cichlid species into six macroanatomical regions. Ophthalmotilapia nasuta and O. ventralis have similar morphology and behavior and occasionally hybridize in the wild. We use 3' mRNA sequencing and a stage-wise statistical testing procedure to identify differential gene expression between females that were kept in a social setting with other females. Our results show that gene expression differs substantially between all six brain parts within species: out of 11,577 assessed genes, 8,748 are differentially expressed (DE) in at least one brain part compared to the average expression of the other brain parts. At most 16% of these DE genes have |log2FC| significantly higher than two. Functional differences between brain parts were consistent between species. The majority (61-79%) of genes that are DE in a particular brain part were shared between both species. Only 32 genes show significant differences in fold change across brain parts between species. These genes are mainly linked to transport, transmembrane transport, transcription (and its regulation) and signal transduction. Moreover, statistical equivalence testing reveals that within each comparison, on average 89% of the genes show an equivalent fold change between both species. The pronounced differences in gene expression between brain parts and the conserved patterns between closely related species with similar morphologies and behavior suggest that unraveling the interactions between genes and behavior will benefit from neurogenomic profiling of distinct brain regions.

Effects of sublethal abiotic stressors on population growth and genetic diversity of Pellioditis marina (Nematoda) from the Westerschelde estuary.

Understanding the effects of anthropogenic pollutants at the ecosystem level requires a proper understanding of the toxicological effects at the population level. Species living in estuaries resist highly fluctuating conditions, and are often exposed to sublethal concentrations of pollutants coming from industrial and domestic wastes. In the Westerschelde estuary, the most upstream sampled population of the nematode Pellioditis marina is genetically less diverse than elsewhere. It experiences lower salinities and higher Cd concentrations than more downstream populations in the estuary. In the present study, we investigate whether these environmental conditions may explain the lower genetic diversity in the most upstream location. To this end we followed the development of genetically diverse P. marina populations under experimental conditions during 14 days. Genetic diversity was assessed in the F1, F2 and F5 generation by screening mitochondrial cytochrome oxidase c subunit 1 variation with the single-strand conformation polymorphism method (SSCP) and nucleotide sequencing. Our results show that sublethal Cd concentrations reduce population development of P. marina at suboptimal salinities, and that low salinity conditions induce responses at the genetic level. Nevertheless, the genetic effects were not persistent over generations, which emphasize the need for longer multigenerational experiments.