Investigating the ecology and evolution of cryptic marine nematode species through quantitative real-time PCR of the ribosomal ITS region

Temperature-driven dynamics: unraveling the impact of climate change on cryptic species interactions within the Litoditis marina complex

Anthropogenic climate change and the associated increase in sea temperatures are projected to greatly impact marine ecosystems. Temperature variation can influence the interactions between species, leading to cascading effects on the abundance, diversity and composition of communities. Such changes in community structure can have consequences on ecosystem stability, processes and the services it provides. Therefore, it is important to better understand the role of species interactions in the development of communities and how they are influenced by environmental factors like temperature. The coexistence of closely related cryptic species, with significant biological and ecological differences, makes this even more complex. This study investigated the effect of temperature on species growth and both intra- and interspecific interactions of three species within the free-living nematode Litoditis marina complex. To achieve this, closed microcosm experiments were conducted on the L. marina species Pm I, Pm III and Pm IV in monoculture and combined cultures at two temperature treatments of 15 °C and 20 °C. A population model was constructed to elucidate and quantify the effects of intra- and interspecific interactions on nematode populations. The relative competitive abilities of the investigated species were quantified using the Modern Coexistence Theory (MCT) framework. Temperature had strong and disparate effects on the population growth of the distinct L. marina species. This indicates temperature could play an important role in the distribution of these cryptic species. Both competitive and facilitative interactions were observed in the experiments. Temperature affected both the type and the strength of the species interactions, suggesting a change in temperature could impact the coexistence of these closely related species, alter community dynamics and consequently affect ecosystem processes and services.

Selection of a reference gene for studies on lipid-related aquatic adaptations of toothed whales (Grampus griseus)

Toothed whales are one group of marine mammals that has developed special adaptations, such as echolocation for predation, to successfully live in a dynamic aquatic environment. Their fat metabolism may differ from that of other mammals because toothed whales have acoustic fats. Gene expression in the metabolic pathways of animals can change with respect to their evolution and environment. A real-time quantitative polymerase chain reaction (RT-qPCR) is a reliable technique for studying the relative expressions of genes. However, since the accuracy of RT-qPCR data is totally dependent on the reference gene, the selection of the reference gene is an essential step. In this study, 10 candidate reference genes (ZC3H10, FTL, LGALS1, RPL27, GAPDH, FTH1, DCN, TCTP, NDUS5, and UBIM) were initially tested for amplification efficiency using RT-qPCR. After excluding DCN, the remaining nine genes, which are nearly 100% efficient, were selected for the gene stability analysis. Stable reference genes across eight different fat tissue, liver, and muscle samples from Grampus griseus were identified by four algorithms, which were provided in Genorm, NormFinder, BestKeeper, and Delta CT. Finally, a RefFinder comprehensive ranking was performed based on the stability values, and the nine genes were ranked as follows: LGALS1 > FTL > GAPDH > ZC3H10 > FTH1 > NDUS5 > TCTP > RPL27 > UBIM. The LGALS1 and FTL genes were identified as the most stable novel reference genes. The third-ranked gene, GAPDH, is a well-known housekeeping gene for mammals. Ultimately, we suggest the use of LGALS1 as a reliable novel reference gene for genomics studies on the lipid-related aquatic adaptations of toothed whales.

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.

Insights on foodborne zoonotic trematodes in freshwater snails in North and Central Vietnam

Foodborne zoonotic trematode (FZT) infections are common neglected tropical diseases in Southeast Asia. Their complicated life cycles involve freshwater snails as intermediate hosts. A cross-sectional study was conducted in Yen Bai and Thanh Hoa provinces in North and Central Vietnam, to investigate the diversity of cercariae of potential FZT and to construct the phylogenetic relationship of trematode cercariae based on the Internal Transcribed Spacer 2 (ITS2) region. Among 17 snail species collected from various habitats, 13 were infected by 10 cercarial groups among which parapleurolophocercous, pleurolophocercous, and echinostome cercariae were of zoonotic importance. The monophyletic tree separated cercarial sequences into different groups following the description of the cercariae families in which Haplorchidae, Opisthorchiidae, Echinochasmidae, and Echinostomatidae are important families of FZT. The overall prevalence was different among snail species and habitats and showed a seasonal trend. Parapleurolophocercous and echinostome cercariae emerged as the most common cercariae in snails in Yen Bai, while monostome, echinostome, and megalura cercariae were most common in Thanh Hoa. Using a molecular approach, we identified Parafossarulus striatulus as the first intermediate snail host of Clonorchis sinensis in Thac Ba Lake. Melanoides tuberculata and Bithynia fuchsiana were we identified preferred intermediate snail hosts of a diverse range of trematode species including intestinal flukes (i.e., Haplorchis pumilio and Echinochasmus japonicus) in Yen Bai and Thanh Hoa, respectively.

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.

Imagining Sisyphus happy: DNA barcoding and the unnamed majority

The vast majority of life on the Earth is physically small, and is classifiable as micro- or meiobiota. These organisms are numerically dominant and it is likely that they are also abundantly speciose. By contrast, the vast majority of taxonomic effort has been expended on ‘charismatic megabionts’: larger organisms where a wealth of morphology has facilitated Linnaean species definition. The hugely successful Linnaean project is unlikely to be extensible to the totality of approximately 10 million species in a reasonable time frame and thus alternative toolkits and methodologies need to be developed. One such toolkit is DNA barcoding, particularly in its metabarcoding or metagenetics mode, where organisms are identified purely by the presence of a diagnostic DNA sequence in samples that are not processed for morphological identification. Building on secure Linnaean foundations, classification of unknown (and unseen) organisms to molecular operational taxonomic units (MOTUs) and deployment of these MOTUs in biodiversity science promises a rewarding resolution to the Sisyphean task of naming all the world's species.

This article is part of the themed issue ‘From DNA barcodes to biomes’.

Mitogenomics reveals high synteny and long evolutionary histories of sympatric cryptic nematode species

Species with seemingly identical morphology but with distinct genetic differences are abundant in the marine environment and frequently co-occur in the same habitat. Such cryptic species are typically delineated using a limited number of mitochondrial and/or nuclear marker genes, which do not yield information on gene order and gene content of the genomes under consideration. We used next-generation sequencing to study the composition of the mitochondrial genomes of four sympatrically distributed cryptic species of the Litoditis marina species complex (PmI, PmII, PmIII, and PmIV). The ecology, biology, and natural occurrence of these four species are well known, but the evolutionary processes behind this cryptic speciation remain largely unknown. The gene order of the mitochondrial genomes of the four species was conserved, but differences in genome length, gene length, and codon usage were observed. The atp8 gene was lacking in all four species. Phylogenetic analyses confirm that PmI and PmIV are sister species and that PmIII diverged earliest. The most recent common ancestor of the four cryptic species was estimated to have diverged 16 MYA. Synonymous mutations outnumbered nonsynonymous changes in all protein-encoding genes, with the Complex IV genes (coxI-III) experiencing the strongest purifying selection. Our mitogenomic results show that morphologically similar species can have long evolutionary histories and that PmIII has several differences in genetic makeup compared to the three other species, which may explain why it is better adapted to higher temperatures than the other species.

Unraveling the intraguild competition between Oscheius spp. nematodes and entomopathogenic nematodes: Implications for their natural distribution in Swiss agricultural soils

Entomopathogenic nematodes (EPN) are excellent biological control agents to fight soil-dwelling insect pests. In a previous survey of agricultural soils of Switzerland, we found mixtures of free-living nematodes (FLN) in the genus Oscheius, which appeared to be in intense competition with EPN. As this may have important implications for the long-term persistence of EPN, we studied this intraguild competition in detail. We hypothesized that (i) Oscheius spp. isolates act as scavengers rather than entomopathogens, and (ii) cadavers with relatively small numbers of EPN are highly suitable resources for Oscheius spp. reproduction. To study this, we identified Oscheius spp. isolated from Swiss soils, quantified the outcome of EPN/Oscheius competition in laboratory experiments, developed species-specific primers and probe for quantitative real-time PCR, and evaluated their relative occurrence in the field in the context of the soil food web. Molecular analysis (ITS/D2D3) identified MG-67/MG-69 as Oscheius onirici and MG-68 as O. tipulae (Dolichura-group). Oscheius spp. indeed behaved as scavengers, reproducing in ∼64% of frozen-killed cadavers from controlled experiments. Mixed infection in the laboratory by Oscheius spp. with low (3 IJs) or high (20 IJs) initial EPN numbers revealed simultaneous reproduction in double-exposed cadavers which resulted in a substantial reduction in the number of EPN progeny from the cadaver. This effect depended on the number of EPN in the initial inoculum and differed by EPN species; Heterorhabditis megidis was better at overcoming competition. This study reveals Oscheius spp. as facultative kleptoparasites that compete with EPN for insect cadavers. Using real-time qPCR, we were able to accurately quantify this strong competition between FLN and EPN in cadavers that were recovered after soil baiting (∼86% cadavers with >50% FLN production). The severe competition within the host cadavers and the intense management of the soils in annual crops readily explain the low EPN numbers in Swiss field samples. The developed molecular tools can be used to elucidate the extent to which the competitive interactions affect EPN populations. This can help to develop strategies to achieve good persistence and natural EPN recycling, in particular in systems where native EPN levels are low, such as annual crops.

Daily Temperature Fluctuations Alter Interactions between Closely Related Species of Marine Nematodes

In addition to an increase in mean temperature, climate change models predict decreasing amplitudes of daily temperature fluctuations. In temperate regions, where daily and seasonal fluctuations are prominent, such decreases in daily temperature fluctuations can have a pronounced effect on the fitness of species and on the outcome of species interactions. In this study, the effect of a temperature regime with daily fluctuations versus a constant temperature on the fitness and interspecific interactions of three cryptic species of the marine nematode species complex of Litoditis marina (Pm I, Pm III and Pm IV) were investigated. In a lab experiment, different combinations of species (monospecific treatment: Pm I and Pm IV and Pm III alone; two-species treatment: Pm I + Pm IV; three-species treatment: Pm I + Pm IV + Pm III) were subjected to two different temperature regimes: one constant and one fluctuating temperature. Our results showed that fluctuating temperature had minor or no effects on the population fitness of the three species in monocultures. In contrast, interspecific interactions clearly influenced the fitness of all three species, both positively and negatively. Temperature regime did have a substantial effect on the interactions between the species. In the two-species treatment, temperature regime altered the interaction from a sort of mutualism to commensalism. In addition, the strength of the interspecific interactions changed depending on the temperature regime in the three-species treatment. This experiment confirms that interactions between the species can change depending on the abiotic environment; these results show that it is important to incorporate the effect of fluctuations on interspecific interactions to predict the effect of climate change on biodiversity.

First molecular phylogeny of the circumtropical bivalve family Pinnidae (Mollusca, Bivalvia): evidence for high levels of cryptic species diversity

The family Pinnidae Leach, 1819, includes approximately 50 species of large subtidal and coastal marine bivalves. These commercially important species occur in tropical and temperate waters around the world and are most frequently found in seagrass meadows. The taxonomy of the family has been revised a number of times since the early 20th Century, the most recent revision recognizing 55 species distributed in three genera: Pinna, Atrina and Streptopinna, the latter being monotypic. However, to date no phylogenetic analysis of the family has been conducted using morphological or molecular data. The present study analyzed 306 pinnid specimens from around the world, comprising the three described genera and ca. 25 morphospecies. We sequenced the mitochondrial genes 16S rRNA and cytochrome c oxidase subunit I, and the nuclear ribosomal genes 18S rRNA and 28S rRNA. Phylogenetic analysis of the data revealed monophyly of the genus Atrina but also that the genus Streptopinna is nested within Pinna. Based on the strong support for this relationship we propose a new status for Streptopinna Martens, 1880 and treat it as a subgenus (status nov.) of Pinna Linnaeus, 1758. The phylogeny and the species delimitation analyses suggest the presence of cryptic species in many morphospecies displaying a wide Indo-Pacific distribution, including Pinna muricata, Atrina assimilis, A. exusta and P. (Streptopinna) saccata but also in the Atlantic species A. rigida. Altogether our results highlight the challenges associated with morphological identifications in Pinnidae due to the presence of both phenotypic plasticity and morphological stasis and reveal that many pinnid species are not as widely distributed as previously thought.

SSU ribosomal DNA-based monitoring of nematode assemblages reveals distinct seasonal fluctuations within evolutionary heterogeneous feeding guilds

Soils are among the most complex, diverse and competitive habitats on Earth and soil biota are responsible for ecosystem services such as nutrient cycling, carbon sequestration and remediation of freshwater. The extreme biodiversity prohibits the making of a full inventory of soil life. Hence, an appropriate indicator group should be selected to determine the biological condition of soil systems. Due to their ubiquity and the diverse responses to abiotic and biotic changes, nematodes are suitable indicators for environmental monitoring. However, the time-consuming microscopic analysis of nematode communities has limited the scale at which this indicator group is used. In an attempt to circumvent this problem, a quantitative PCR-based tool for the detection of a consistent part of the soil nematofauna was developed based on a phylum-wide molecular framework consisting of 2,400 full-length SSU rDNA sequences. Taxon-specific primers were designed and tested for specificity. Furthermore, relationships were determined between the quantitative PCR output and numbers of target nematodes. As a first field test for this DNA sequence signature-based approach, seasonal fluctuations of nematode assemblages under open canopy (one field) and closed canopy (one forest) were monitored. Fifteen taxa from four feeding guilds (covering ∼ 65% of the free-living nematode biodiversity at higher taxonomical level) were detected at two trophic levels. These four feeding guilds are composed of taxa that developed independently by parallel evolution and we detected ecologically interpretable patterns for free-living nematodes belonging to the lower trophic level of soil food webs. Our results show temporal fluctuations, which can be even opposite within taxa belonging to the same guild. This research on nematode assemblages revealed ecological information about the soil food web that had been partly overlooked.

Differences in time until dispersal between cryptic species of a marine nematode species complex

Co-occurrence of closely related species may be achieved in environments with fluctuating dynamics, where competitively inferior species can avoid competition through dispersal. Here we present an experiment in which we compared active dispersal abilities (time until first dispersal, number and gender of dispersive adults, and nematode densities at time of dispersal) in Litoditis marina, a common bacterivorous nematode species complex comprising four often co-occurring cryptic species, Pm I, II, III, and IV, as a function of salinity and food distribution. The experiment was conducted in microcosms consisting of an inoculation plate, connection tube, and dispersal plate. Results show species-specific dispersal abilities with Pm I dispersing almost one week later than Pm III. The number of dispersive adults at time of first dispersal was species-specific, with one dispersive female in Pm I and Pm III and a higher, gender-balanced, number in Pm II and Pm IV. Food distribution affected dispersal: in absence of food in the inoculation plate, all species dispersed after ca four days. When food was available Pm I dispersed later, and at the same time and densities irrespective of food conditions in the dispersal plate (food vs no food), suggesting density-dependent dispersal. Pm III dispersed faster and at a lower population density. Salinity affected dispersal, with slower dispersal at higher salinity. These results suggest that active dispersal in Litoditis marina is common, density-dependent, and with species, gender- and environment-specific dispersal abilities. These differences can lead to differential responses under suboptimal conditions and may help to explain temporary coexistence at local scales.

Wide interguild relationships among entomopathogenic and free-living nematodes in soil as measured by real time qPCR

Entomopathogenic nematodes (EPNs) are promising biological control agents of soil-dwelling insect pests of many crops. These nematodes are ubiquitous in both natural and agricultural areas. Their efficacy against arthropods is affected directly and indirectly by food webs and edaphic conditions. It has long been suggested that a greater understanding of EPN ecology is needed to achieve consistent biological control by these nematodes and the development of molecular tools is helping to overcome obstacles to the study of cryptic organisms and complex interactions. Here we extend the repertoire of molecular tools to characterize soil food webs by describing primers/probe set to quantify certain free-living, bactivorous nematodes (FLBNs) that interact with EPNs in soil. Three FLBN isolates were recovered from soil baited with insect larvae. Morphological and molecular characterization confirmed their identities as Acrobeloides maximum (RT-1-R15C and RT-2-R25A) and Rhabditis rainai (PT-R14B). Laboratory experiments demonstrated the ability of these FLBNs to interfere with the development of Steinernema diaprepesi, Steinernema riobrave and Heterorhabditis indica parasitizing the weevil Diaprepes abbreviatus (P<0.001), perhaps due to resource competition. A molecular probe was developed for the strongest competitor, A. maximum. We selected the highly conserved SSU rDNA sequence to design the primers/probe, because these sequences are more abundantly available for free-living nematodes than ITS sequences that can likely provide better taxonomic resolution. Our molecular probe can identify organisms that share ⩾98% similarity at this locus. The use of this molecular probe to characterize soil communities from samples of nematode DNA collected within a citrus orchard revealed positive correlations (P<0.01) between Acrobeloides-group nematodes and total numbers of EPNs (S. diaprepesi, H. indica and Heterorhabditis zealandica) as well as a complex of nematophagous fungi comprising Catenaria sp. and Monachrosporium gephyropagum that are natural enemies of EPNs. These relationships can be broadly interpreted as supporting Linford's hypothesis, i.e., decomposition of organic matter (here, insect cadavers) greatly increases bactivorous nematodes and their natural enemies.