Differences in Life-Histories Refute Ecological Equivalence of Cryptic Species and Provide Clues to the Origin of Bathyal Halomonhystera (Nematoda)

Abstract

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 Håkon 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.

About the Authors:

Jelle Van Campenhout

* E-mail: [email protected]

Affiliations Biology Department, Research Group Marine Biology, Ghent University, Ghent, Belgium, Center for Molecular Phylogenetics and Evolution (CeMoFe), Ghent University, Ghent, Belgium

Sofie Derycke

Affiliations Biology Department, Research Group Marine Biology, Ghent University, Ghent, Belgium, Center for Molecular Phylogenetics and Evolution (CeMoFe), Ghent University, Ghent, Belgium

Tom Moens

Affiliations Biology Department, Research Group Marine Biology, Ghent University, Ghent, Belgium, Center for Molecular Phylogenetics and Evolution (CeMoFe), Ghent University, Ghent, Belgium

Ann Vanreusel

Affiliation: Biology Department, Research Group Marine Biology, Ghent University, Ghent, Belgium

Introduction

DNA sequencing of multiple independently evolving gene regions has revealed that distinct taxonomic units, previously classified as single species due to a similar morphology, can be distinguished based on nuclear and mitochondrial sequence markers [1], [2]. Such cryptic diversity is almost evenly distributed among major metazoan taxa and biogeographical regions [2]. However, it has been suggested that cryptic diversity is more common in marine environments because many marine species rely on chemical cues for mate recognition [3]–[5].

The co-occurrence of species with highly similar morphologies has been frequently observed [6]–[8] and indicates that they are adapted to ecologically similar environments. Similarities in body size and life histories of cryptic species may form important equalizing mechanisms which can minimize average fitness differences between species and, consequently, slow down competitive exclusion [9]. However, these equalizing mechanisms are unlikely to lead to stable coexistence [9], [10]. Stable coexistence of species can only exist if both equalizing- and stabilizing mechanisms, such as fitness differences [9], [11], niche differentiation and/or density-dependent life history adjustments [12], [13], are present [14]. Even then, the sufficiency of both mechanisms highly depends on the species composition and the strength of each individual mechanism [9]. To date, ecological characterization of cryptic invertebrate species has focused on nutritional [15] and habitat preferences [16], [17]. However, susceptibility to predation [17], [18], response to abiotic factors [19] and species-specific symbioses [20] are additional examples of stabilizing mechanisms which may also contribute to stable coexistence of cryptic species. Whether and to what extent this implies that cryptic species are ecologically equivalent remains unknown.

In marine sediments, nematodes are usually the most abundant and diverse Metazoa [21], [22]. Cryptic diversity has been reported in marine nematodes belonging to different orders [23]–[25]. These cryptic species can exhibit differential dispersal capacities [26] and differential resource use, and environmental factors such as salinity and temperature may affect the outcome of their competitive interactions when they co-occur [27]. However, at this point, we lack any information on if and how life-history traits of cryptic nematode species differ in response to changes in abiotic factors. Such information is, however, important to understand mechanisms driving the coexistence cryptic species, which in turn is crucial to correctly understand the biodiversity and ecosystem functioning relationship [28].

In this study, we focus on the marine nematode Halomonhystera disjuncta (previously named Geomonhystera disjuncta), a representative of a genus which has a widespread geographical distribution and which has been reported both in shallow-water [23], [29]–[31] and in deep-sea environments [32], [33]. In the Westerschelde estuary (located in the Southwestern part of The Netherlands), five cryptic species (GD1-5) of H. disjuncta have been reported based on nuclear and mitochondrial sequence data; these five species exhibit subtle morphometric differences [23], [34]. Halomonhystera disjuncta has most often been isolated from macroalgal holdfasts and wrack deposits, but has also been frequently observed in the sediment, successfully exploiting organically enriched substrata [21], [35], [36]. In the Westerschelde estuary, the five cryptic species (GD1-5) show sympatric distributions [23]. In view of their mainly intertidal occurrence, these species are subject to strong physical, chemical and biological gradients. Ecological and physiological responses of H. disjuncta from the intertidal zone to changes in salinity, temperature, food quality, food density and heavy metals have been investigated [21], [37]–[42]. However, it is unknown whether the observed broad tolerance to e.g. temperature and heavy metals is the same for different cryptic species.

Halomonhystera is also the dominant nematode genus in the sulphide rich bacterial mats of the Nyegga pockmark (Nordic Norwegian margin) [32] and the Håkon Mosby mud volcano (HMMV, Barent Sea slope) [33], situated at depths of 730 m and 1280 m, respectively. At both locations, fluids escape from the deep-sea floor at lower temperatures and flow rates compared to mid-oceanic ridges, and are, therefore, called cold seeps. The Halomonhystera from both environments were originally identified as H. disjuncta, but the Halomonhystera from the HMMV has morphological and genetic differences compared to the shallow-water H. disjuncta species GD1-5, and has been described as a new species Halomonhystera hermesi [43]. Phylogenetic analysis revealed three intertidal clades (GD1/GD4, GD3, GD2/GD5) and one deep-sea clade Halomonhystera hermesi [44] and showed that H. hermesi is more closely related to GD1 and GD4 than to the other species. A deep-sea invasion from shallow-water regions has hence been hypothesized [44]. Nordic seep colonization from intertidal regions implies that early colonizers had to adapt to low temperatures, higher salinities and the presence of high concentrations of sulphides. If cryptic species of H. disjuncta differ in their tolerance to one or more of these environmental differences, then we could expect that the species with the highest degree of tolerance to bathyal cold seep conditions is phylogenetically closest related to H. hermesi.

Here we investigate the effect of three abiotic factors (salinity, temperature and sulphide) on life-history traits of three cryptic H. disjuncta species (GD1, 2 and 3). Each of the three abiotic factors varied from control conditions (salinity of 25 psu, 16°C and no sulphide) to environmental conditions representative of the HMMV (salinity of 34–35 psu, 2°C and hydrogen sulphide concentrations of ca. 1 mM) [45], [46] in a fully crossed factorial design with three salinity levels (25, 29.5 and 34 psu), three temperatures (16, 10 and 4°C) and two sulphide treatments (no sulphide vs. 1 mM). Because these three cryptic species co-occur in the Westerschelde estuary [23], are phenotypically almost identical [34], and were isolated from the same species of macroalgae (Fucus vesiculosus), they are expected to show a high degree of ecological equivalence with respect to abiotic environmental factors. In contrast, if differences in life history traits would be observed, then we would expect GD1 to better tolerate bathyal cold seep conditions than GD2 and GD3, because it is phylogenetically more closely related to H. hermesi compared to GD2-3. This prediction follows the rationale of the phylogenetic niche conservatism theory (PNC). PNC is the tendency of lineages to retain their niche-related traits through speciation events and over long evolutionary periods [47], [48]. Consequently, phylogenetically closely related species are expected to be ecologically similar [49].

Materials and Methods

Nematode cultures

Halomonhystera disjuncta species GD1 and GD3 were retrieved from decaying Fucus vesiculosus from the Paulina tidal flat in the Westerschelde estuary (The Netherlands, 51° 20′ 56.79′′ N, 3° 43′ 29.56′′ E), while GD2 was isolated from the same decaying macroalgae species from a different location (Kruispolderhaven, 51° 21′ 34.35′′ N, 4° 5′ 52.02′′ E) in the same estuary. All species were collected in April 2012. Permission for the field work was issued by the Provincie Zeeland, the Netherlands (Directie Ruimte, Milieu en Water). Pieces of algae were inoculated on Petri dishes filled with 0.8% nutrient:bacto agar (ratio of 1/7) prepared in artificial seawater [50] with a salinity of 25 psu and placed at a constant temperature of 16°C. These conditions were chosen based on the temperature and salinity values observed in the field locations from where the species were isolated and are perfectly suited for cultivation of H. disjuncta species. Monospecific cultures were established by transferring a single gravid female to a new Petri dish (5.5 cm inner diameter) containing the same agar medium. By doing so, we reduced population variation, which was necessary to ensure monospecific status of the cultures. This, however, implies that differences observed here between cryptic species may partly reflect individual and/or population level differences, for instance resulting from local adaptation, rather than consistent differences between species. The monospecific cultures were maintained at the same temperature (16°C) with frozen-and-thawed Escherichia coli K12 as a food source. Species identity and monospecificity of the cultures was confirmed by PCR and sequencing the ITS gene of ten nematodes per culture. Primers, thermo-cycling conditions and GenBank accession numbers can be found in [23]. Nematodes for the experiment were harvested from the stock cultures.

Experimental design

To test the effect of three abiotic factors (salinity, temperature and sulphide) between cryptic species, a fully crossed design with four factors: salinity (25, 29.5 and 34 psu), temperature (16, 10 and 4°C), species (GD1, GD2 and GD3) and sulphide (present or absent) was set up. For each treatment x species combination, four replicates were prepared, resulting in 216 microcosms. The experiment was conducted on experimental microcosms which consisted of Petri dishes (3.5 cm i.d.) containing 1.5 ml 0.8% nutrient:bacto agar (1/7). Because stock cultures are maintained at a salinity of 25 psu and 16°C, these conditions served as a control. The sulphidic media were prepared from the same agar medium as above, but to which we added sodium thiosulphate (Na₂S₂O₃) and sodium sulphide (Na₂S) in a final concentration of 1 µM for both reagents. Lastly, to each microcosm 350 µl E. coli K12 (ca. 7.6×10⁸ cells/ml) was added as a food source, which is sufficient for the duration of the experiment.

Each replicate microcosm received 10 nematodes. We randomly picked out six female and four male nematodes, in accordance with the sex ratio in the stock cultures. Nematodes were considered adult when all parts of the reproductive system were clearly visible. Only nematodes which had just reached adulthood were randomly selected. Immediately after addition of the nematodes to the Petri dishes, we checked them carefully for specimens which died, were injured or immobilized during manipulation. These were replaced by new specimens of the same sex. The Petri dishes were then closed with parafilm, maintained in temperature-controlled incubators without light, and the experiment was started.

Data collection and analysis

The amounts of eggs, juveniles and adults were counted daily. In addition, we observed mortality of inoculated nematodes. Young F1-adults are easily distinguishable from F0-adults, up to a few days after they reached adulthood, because of their smaller body size. Therefore, the experiment was stopped before F1-adults started to deposit eggs and/or when we were no longer able to distinguish new adults from the inoculated nematodes. This was after ca. 13, 18 and 38 days at 16°C, 10°C and 4°C, respectively. Hence, our results encompassed a single nematode generation.

Four response variables were used for statistical analysis: (1) minimum time until first egg(s) deposition (MEGD) is the minimum time (in days) from the inoculation of the nematodes until we observed the first egg(s); (2) the minimum time for embryonic development (MEMD) was calculated by subtracting MEGD from the time at which we observed the first juvenile; (3) the minimum time for the development of juveniles into adults (MJD) was calculated by subtracting the time at which the first juvenile was observed from the time at which the first F1 adult was detected; (4) the minimum generation time (MGT) was recorded as the time from the start of the experiment (parental generation) until the first new adult (filial generation) was observed.

In addition, fertility was approximated by summing the counts of eggs, juveniles and new adults at the end of the incubation. Fertility was expressed per parental female by dividing total offspring (eggs, juveniles and new adults) of a microcosm, by the number of inoculated females. Furthermore, mortality of the inoculated nematodes was monitored. We, therefore, created an additional variable, the minimum adult lifespan (MALS), which corresponds to the day at which we observed the death of at least one inoculated adult nematode. It must be noted that the experiment for some treatments was stopped before any mortality of inoculated nematodes was observed. Small standard deviations and visual inspection of our data revealed that MALS was not importantly influenced by stochasticity, i.e. the incidental death of an inoculated individual in random replicates and treatments.

Because our data did not meet the assumptions for parametric variance analysis, even after transformation, a Permutational Based Multivariate Analysis of Variance [51], on the basis of Euclidean distances with 9999 permutations, was used. In comparison to ANOVA/MANOVA, both assuming normal distributions and, implicitly, Euclidean distance, PERMANOVA works with any distance measure that is appropriate to the data, and uses permutations to make it distribution free [51]. Each analysis was performed using a fully crossed PERMANOVA design with four fixed factors (salinity, temperature, sulphide presence and species). All analyses were performed within PRIMER v6 with PERMANOVA+ add-on software [51]. Because of the relatively limited level of replication (4 or 3, the latter in few cases where we observed agar dehydration as a result of an imperfect Petri dish closure), a Monte Carlo test was performed. The components of variation, estimated by PERMANOVA, were used to attribute the amount of variation to different factors and interaction terms [52]. A high component value corresponds with a high relative importance of the respective term in the model [51]. Pairwise comparisons were performed on the full dataset and p-values can be found in the supplementary information. If the four-way interaction was significant, we were unable to reliably interpret two-way and three-way interaction terms. In these cases we selected subsets of our data for which only two or three factors varied. However, we tested these two- and three-way interaction effects for all levels and combinations of the factor(s) that were kept constant. PERMDISP [53] was executed to test the homogeneity of multivariate dispersions in order to discriminate between real location or factor effects and effects explained by differences in dispersion for the significant factors.

Figures were made with Graphpad prism 5.

Results

Minimum generation time (MGT)

MGT is the sum of the minimum time until first egg(s) deposition (MEGD), the minimum time for embryonic development (MEMD) and the minimum time for the development of juveniles into adults (MJD). A significant four-way interaction effect on MGT was found (PERMANOVA; p = 0.0001, Table 1) which was not caused by differences in multivariate dispersion (F_{43, 109} = 0.408, p = 0.9994). The estimation of the components of variation showed that temperature was the most important factor influencing MGT (Table 2).

[Figure omitted. See PDF.]

Table 1. P-values of a PERMANOVA test on all possible interaction.

https://doi.org/10.1371/journal.pone.0111889.t001

[Figure omitted. See PDF.]

Table 2. Estimates of components of variation results of the four-way PERMANOVA on single factors and interaction terms for all analyzed variables.

https://doi.org/10.1371/journal.pone.0111889.t002

MGT increased with decreasing temperature (T) for all species (Sp; Table 3). When species are compared in presence and absence of sulphide (S), GD2 generally had the shortest MGT at 16°C and at 4°C. Except at cold seep conditions (at a temperature of 4°C, a salinity of 34 psu and in the presence of sulphide), GD1 had a significantly shorter MGT than GD2 (p = 0.0156, Table S1). A significant shorter MGT for GD1, in comparison to GD2-3, was also observed at 10°C (all p-values <0.0018, Table S1). The difference in MGT between species in relation to temperature was underlined by a significant Sp x T interaction (all p-values = 0.0001, Table 1) independent of salinity and sulphide addition.

[Figure omitted. See PDF.]

Table 3. Minimum generation time of GD1-3 under changing environmental conditions.

https://doi.org/10.1371/journal.pone.0111889.t003

Increasing the salinity (Sal) resulted in an increase of MGT. The response to higher salinities is more pronounced for GD3 (2–10 days) than for GD1-2 (0–6 days; Table 3), in the absence of sulphide. This observation was supported by a significant Sp x Sal interaction at all temperatures (all p-values <0.0468, Table 1).

The effect of adding sulphide to the cultures resulted in the death of all inoculated GD3 nematodes during the first two minutes and no life-history traits could, therefore, be determined. The addition of sulphide resulted in a significant increase in MGT of 1–10 days for GD1 and GD2 at 10°C and 4°C (all p-values <0.0138, Table S2) at all salinity levels (Table 3).

Minimum time until first egg(s) deposition (MEGD)

GD3 was the only species showing an ovoviviparous reproduction strategy throughout the whole experiment. Therefore, we did not include GD3 in the MEGD dataset. For the two other species, abiotic factors significantly affected MEGD (four-way PERMANOVA; p = 0.0065, Table 1). Homogeneity of multivariate dispersions was observed (F_{35, 89} = 0.736, p = 0.8453). The estimation of the components of variation revealed that temperature is the most import factor influencing MEGD (Table 2).

A strong significant increase in MEGD from 10°C to 4°C (all p-values <0.0003, Table S3) was observed for both species (Figure 1). This increase was observed at all salinity levels and in presence/absence of sulphide. Increasing salinity usually resulted in a significant increase in MEGD for both species (p-values in Table S4) and was most pronounced at 4°C for both species, supported by a significant T x Sal interaction term (all p-values <0.0026, Table 1).

[Figure omitted. See PDF.]

Figure 1. Minimum time until first egg(s) deposition of GD1 and GD2 under changing environmental conditions.

Minimum time until first egg(s) deposition in relation to temperature and salinity of cryptic Halomonhystera disjuncta species: GD1 (A, C) and GD2 (B, D), in the absence (No S; A, B) and presence (S; C, D) of sulphide. The minimum egg deposition time corresponds to the day after the start of the experiment on which the first egg was observed within the experimental microcosms. Data shown are mean values ±1stdev of 3 or 4 replicates per treatment.

https://doi.org/10.1371/journal.pone.0111889.g001

The effects of sulphide were heavily dependent on temperature for both species (significant S x T interaction term, all p-values <0.0003 at all salinity levels, Table 1). Addition of sulphide resulted in a significant decrease (all p-values <0.0031, Table S2) in MEGD of GD1 at 16°C, (Figure 1A, 1C) but not of GD2 (Figure 1B, 1C). The observed difference between species was corroborated by significant S x Sp interaction term 16°C at all salinity levels (p-values <0.0383, Table 1). At 10°C, adding sulphide resulted in a significant increase in MEGD of GD1 (all p-values <0.0429, table S2) but not for GD2. At the lowest temperature (4°C), adding sulphide resulted in an increase in MEGD for both species (all p-values <0.0134, Table S2).

Minimum time for embryonic development (MEMD)

Due to GD3’s ovoviviparous reproduction strategy we were unable to estimate MEMD for this species. Similar to the MEGD, a significant four-way interaction (p = 0.0002, Table 1) and homogeneity of multivariate dispersions were found (F_1,123 = 0.042, p = 0.8381). The estimates of components of variation revealed that temperature and the interaction term Sp x T were the most important effects (Table 2).

Comparable with MEGD, lowering the temperature increased MEMD (Figure 2A–D). Both GD1 and GD2 had a stable MEMD between 16°C and 10°C, while a significant increase was observed at 4°C (all p-values <0.0124, Table S3). This increase was much more pronounced for GD1 than for GD2, underlined by a significant Sp x T interaction effect (all p-values <0.001, Table 1) independent of salinity and sulphide.

[Figure omitted. See PDF.]

Figure 2. Minimum time for embryonic development of GD1 and GD2 under changing environmental conditions.

Minimum time for embryonic development in relation to temperature and salinity of cryptic Halomonhystera disjuncta species: GD1 (A, C) and GD2 (B, D), in the absence (No S; A, B) and presence (S; C, D) of sulphide. The minimum time for embryonic development was calculated by subtracting the day we observed the first egg from the day we observed the first juvenile within the experimental microcosms. Data shown are mean values ±1stdev of 3 or 4 replicates per treatment.

https://doi.org/10.1371/journal.pone.0111889.g002

Increasing the salinity had no significant effect on MEMD for GD1, but had an opposite effect on MEMD for GD2 at 4°C between different sulphide treatments (Figure 2B, 2D). This observation is supported by a significant S x Sp x Sal interaction component at 4°C (p-value = 0.001, Table 1). The addition of sulphide had no further consistent or opposing effects on MEMD. Interestingly, at bathyal cold seep conditions, GD2 had a significantly shorter MEMD than GD1 (p<0.0029, Table S1).

Minimum time for the development of juveniles into adults (MJD)

MJD showed a significant four-way interaction effect (p = 0.0433, Table 1) and homogeneity of multivariate dispersions was found (F_43,109 = 0.870, p = 0.6917). According to the estimates of components of variation, temperature affects MJD the most (Table 2).

Lowering temperature resulted in a significant increase in MJD (Figure 3A–E) for both GD2 and GD3 independent of other factors (all p-values <0.0189, Table S3). GD1 also showed a significant increase in MJD from 10°C and 16°C to 4°C (all p-values <0.0024, Table S3; Figure 3A, 3B) but not between 16°C to 10°C (all p-values >0.1092, Table S3). This differential effect of temperature changes on species was confirmed by a significant Sp x T interaction effects at all salinity and sulphide levels (all p-values <0.0003, Table 1). Interestingly, GD3 did not produce offspring at 16°C and a salinity of 34 psu. Even though female nematodes had eggs inside the uterus, which developed into juveniles, they were never deposited.

[Figure omitted. See PDF.]

Figure 3. Minimum time for the development of juveniles into adults of GD1-3 under changing environmental conditions.

Minimum time for the development of juveniles into adults in relation to temperature and salinity of cryptic Halomonhystera disjuncta species: GD1 (A, D), GD2 (B, E) and GD3 (C), in the absence (No S; A–C) and presence (S; D–E) of sulphide. The minimum time for embryonic development was calculated by subtracting the day we observe the first juvenile from the day we observed the first adult within the experimental microcosms. Data shown are mean values ±1stdev of 3 or 4 replicates per treatment. No results were obtained for GD3 in the presence of sulphide because inoculated nematodes died within the first 1–2 minutes after adding Su.

https://doi.org/10.1371/journal.pone.0111889.g003

Our data also shows that at 4°C the effect of salinity interacted with the addition of sulphide and was more pronounced for GD2 than for GD1 (significant S x Sp x Sal interaction term, p = 0.0133, Table 1). In addition, at bathyal cold seep conditions GD1 had a significantly shorter MJD compared to GD2 (p = 0.0009, Table S1).

Minimum adult life span (MALS)

The four-way interaction of MALS was not significant (p = 0.3177, Table 1) and homogeneity of multivariate dispersions (F_{41, 99} = 1.085, p = 0.3638) was observed. Some data points are missing because in some treatments we stopped the experiment before any mortality of inoculated nematodes occurred.

A decrease in temperature generally resulted in a significant increase in MALS for all species (Figure 4). However, no significant increase was observed for GD1 from 16°C to 10°C (all p-values >0.1167, Table S3), in contrast with GD2 and GD3 (all p-values <0.0010, Table S3). At 4°C, in the absence of sulphide, GD2 and GD3 had the shortest and largest MALS, respectively. This difference between species was supported by a significant Sp x T interaction effect at all salinity and sulphide levels (all p-values <0.0002, Table 1).

[Figure omitted. See PDF.]

Figure 4. Minimum adult life span under changing environmental conditions.

Minimum adult life span in relation to temperature and salinity of cryptic Halomonhystera disjuncta species: GD1 (A, D), GD2 (B, E) and GD3 (C), in the absence (No S; A–C) and presence (S; D–E) of sulphide. The minimum adult life span corresponds to the day after the start of the experiment on which the first inoculated nematode died. Data shown are mean values ±1stdev of 3 or 4 replicates per treatment. No results were obtained for GD3 in the presence of sulphide because inoculated nematodes died within the first 1–2 minutes after adding Su. Some data points are missing because the experiment was stopped before the death of any of the inoculated nematodes.

https://doi.org/10.1371/journal.pone.0111889.g004

Addition of sulphide resulted in the death of all inoculated GD3 nematodes within minutes. For both other species the addition of sulphide commonly decreased MALS. At 4°C and at bathyal cold seep conditions, GD1 had a higher MALS than GD2 (all p-values <0.0023, Table S1, Figure 4D–E). Sulphide-induced mortality decreased with decreasing temperature: in the presence of sulphide and at a temperature of 16°C, ca. 99% of both inoculated and new adults of GD1 and GD2 died after one generation. However, this high mortality rate was decreased to ca. 50% at 10°C and was negligible at 4°C.

Offspring per female (OS)

Variances were homogeneous (F_{44, 112} = 0.526, p = 0.9931) and a significant four-way interaction effect was found (p = 0.0001, Table 1). Estimates of components of variation revealed that most factors and interactions had an effect on the OS (Table 2) with temperature being the most important one. Consequently, sub-datasets were analyzed to better understand effects of single factors and their two- and three-way interactions.

All cryptic species had the highest OS at 16°C in the absence of sulphide and at a salinity of 25 psu (except for GD2, Figure 5), and the lowest OS at 4°C in the presence of sulphide and at a salinity of 34 psu. A reduction in OS was observed for all species with decreasing temperature, increasing salinity and addition of sulphide. However, the declination in OS for GD2 was clearly more influenced by temperature than for GD1 and GD3 as was supported by a significant Sp x T interaction component (all p-values = 0.0001, Table 1). Variation in OS, caused by temperature and salinity, was much lower in GD3 in the absence of sulphide (p-value _{Sp x T x Sal} = 0.0001, Table 1), compared to the two other species. However, no offspring was observed at 16°C and a salinity of 34 psu. Generally, adding sulphide resulted in a reduction in OS for both GD1 and GD2. Furthermore, GD1 produced significantly more OS than GD2 at bathyal cold seep conditions (p = 0.006, Table S1).

[Figure omitted. See PDF.]

Figure 5. Offspring per female of GD1-3 under changing environmental conditions.

Offspring per female in relation to temperature and salinity of cryptic Halomonhystera disjuncta species: GD1 (A, D), GD2 (B, E) and GD3 (C), in the absence (No S; A–C) and presence (S; D–E) of sulphide. The offspring corresponds to the amount of eggs, juveniles and new adults at the end of the experiment, averaged across six females. Data shown are mean values ±1stdev of 3 or 4 replicates per treatment. No results were obtained for GD3 in the presence of sulphide because inoculated nematodes died within the first 1–2 minutes after adding sulphide.

https://doi.org/10.1371/journal.pone.0111889.g005

Discussion

The success of nematodes in marine environments relates to their ability to adapt to changing environmental conditions [54], [55]. Because the three cryptic H. disjuncta species commonly occur on decaying macroalgae in the high intertidal, they are regularly exposed to short- and long-term fluctuations in abiotic factors such as salinity and temperature. Moreover, they can be exposed to hydrogen sulphide produced by sulphate reducing bacteria thriving on rotting macroalgae [56]. In summary, the three species were expected to show wide tolerances to the abiotic factors tested in this study. In our experiments, temperature had the strongest impact on life-history traits of all cryptic species. The minimum generation time and adult life span increased with decreasing temperature, which is most likely caused by lower metabolic rates at lower temperatures [57]. Faster generation times at higher temperatures shorten the vulnerable period of the embryo, thus preventing embryonic deformities and/or arrest [54], and allow species to rapidly increase in numbers during favorable episodes, a common feature of opportunistic species [58]. Halomonhystera disjuncta is indeed an opportunistic species capable of quickly colonizing suitable patches such as new algal deposits [35]. Salinity, on the other hand, had a relatively minor impact on the life-history traits studied here. The salinity and temperature ranges tested in our study reflect fluctuations in their natural habitat and indicate a strong capacity for adaptation to natural fluctuations, which is consistent with previous studies on other marine nematodes from very similar intertidal habitats [59], [60].

Cryptic H. disjuncta species are not ecologically equivalent

Co-occurrence of cryptic species raises questions regarding their ecological similarity and the mechanisms influencing their coexistence [15], [19]. The three cryptic species studied here have a very similar morphology, have sympatric occurrences and all thrive on the same habitat, i.e. decaying macroalgae in the Westerschelde estuary. We may, therefore, expect that they have evolved comparable adaptations to similar environmental conditions (equalizing mechanism). The general response to temperature and salinity, i.e. an increase in time for life-history traits and a decrease in number of offspring with lower temperatures and higher salinities, was indeed similar in all cryptic species and points towards ecological equivalence. However, our results also uncovered specific life-history responses to different abiotic factors. Life-history traits of GD1 were more stable in the temperature interval between 10°C and 16°C than those of GD2. Interestingly, GD3 did not survive in the presence of sulphides, whereas GD1 and GD2 not only survived but also reproduced and developed normal offspring for at least one generation. These results reveal that cryptic species are not completely ecologically equivalent. Our results might be partly biased because a single gravid female was used to generate the monospecific cultures. Nevertheless, the observed differential effects of abiotic factors on life-history traits, may provide important stabilizing mechanisms facilitating the coexistence of cryptic species [9].

In addition to its sensitivity to sulphides, the most obvious difference between GD3 and the other two cryptic species was its ovoviviparous reproduction strategy. Intrauterine egg hatching (ovoviviparity) allows a more secure survival and early development of offspring and can be an adaptation strategy in toxic environmental conditions [61]. The low number of GD3 offspring compared to GD1 and GD2 suggests a trade-off between the reproductive strategy and fecundity. This is further supported by the observation that the fecundity of GD3 was less affected by changes in salinity and temperature than that of GD1 and GD2. In addition, the minimum adult life span was always the highest for GD3 at all temperatures. The comparatively smaller impact of temperature on GD3 could also in part explain its occurrence throughout the year and its high relative abundance in winter in the Westerschelde estuary [23]. In contrast to GD3, GD2 was often absent in winter but had the highest relative abundance in spring [23]. Our experimental data revealed that GD2 had the fastest minimum generation time at 16°C and the highest number of offspring at 10°C and 16°C. In spite of this, its number of offspring was more sensitive to temperature changes. These observations may explain the high relative abundance of GD2 in spring and its disappearance in winter.

Differential responses of nematodes to abiotic factors can also affect interspecific interactions [27], [62] and dispersal capacities [26], which could in turn facilitate coexistence. Recently, significant differences in resource use among sympatric cryptic nematode species have been demonstrated, and other potentially discriminating biotic/abiotic factors, such as differential susceptibility to predation and competitive interaction with other taxa, should be taken into account to fully grasp the coexistence and spatial/temporal distribution patterns of these cryptic species because the small morphological differences between GD1-3 [34] may not be sufficient to predict ecological segregation [63].

GD1 and the Nordic seep nematode Halomonhystera hermesi share a common ancestor

It is often suggested that current deep-sea biodiversity has largely resulted from recurrent invasions from bathyal and abyssal depths followed by speciation [64]. The existence of morphologically and/or genetically close relatives from intertidal and deep-sea environments [65] tends to support this contention. In an evolutionary context, understanding the adaptations which have allowed deep-sea colonization is an important prerequisite. GD1 is phylogenetically most closely related to the Nordic seep nematode Halomonhystera hermesi [44] and, therefore, shares a most recent common ancestor with H. hermesi.

In the reduced environment of the Håkon Mosby mud volcano (HMMV), organisms encounter stress levels which selectively favor sulphide tolerant species [66]. In view of the 100% mortality of GD3 in the sulphide treatments, it is very unlikely that this species successfully colonized the HMMV. It must, however, be mentioned that GD3 did survive under lower concentrations of sodium thiosulphate and sodium sulphide (up to 0.25 µM) in a preliminary test. Such a mild exposure could increase resistance to subsequent higher doses (hormesis) of the same or different stressors [67]. Nevertheless, GD3 was unable to cope with high concentrations of sulphides. GD1 and GD2, by contrast, survive and produce offspring at the sulphide concentrations tested here. Such a high sulphide tolerance may be an advantage for colonizing reduced environments [66]. Species living in sulphidic environments often have ecto- or endosymbionts which could either serve as a food source and/or play a role in sulphur detoxification [68], [69]. In addition, sulphur inclusions could temporarily reduce the toxic effect of H₂S [70]. However, no signs of bacterial symbionts in H. hermesi were detected [45], and we did not observe sulphur inclusions in GD1-2. Interestingly, however, low temperature had a positive effect on the survival of GD1 and GD2 (but not GD3) under sulphide exposure. Similarly, low temperatures reduced the mortality rate of H. disjuncta subjected to toxic chromium concentrations [40]. Hence, the low temperatures in the deep sea may have been an important factor enabling Halomonhystera to cope with otherwise toxic sulphide concentrations.

Our results further reveal that GD1 had a higher number of offspring at bathyal cold seep conditions (salinity of 34 psu – temperature of 4°C – presence of sulphides) than GD2. In addition, GD1 had the fastest minimum generation time and a longer minimum life span than GD2 at these conditions. This implies that GD1 is more resistant to a combination of sulphides, low temperature and higher salinities, and is in agreement with the close phylogenetic relationship between GD1 and H. hermesi, supporting the idea that it could have successfully colonized Nordic seeps.

Conclusion

We have shown that life-histories of cryptic were differentially affected by changes in temperature, salinity and presence of sulphides. These observed differences imply that closely related cryptic H. disjuncta species are not necessarily ecologically equivalent and that abiotic factors can strongly impact their life-history traits, which in turn can affect their coexistence. Our results further indicate that Halomonhystera hermesi and GD1 share a common ancestor as the latter appears to be more resistant to bathyal seep conditions compared to other cryptic species. The observed limited ecological equivalence among cryptic species may have important repercussions for our understanding of the link between biodiversity and ecosystem functioning [71].

Supporting Information

[Figure omitted. See PDF.]

Table S1.

P-values of pairwise comparisons between species with the same treatment. PERMANOVA (Euclidean distances, 9999 permutations) was used for the pairwise comparisons. The table part between double lines contains the variable factor species (GD1-3) on which we performed pairwise comparisons for the variables in column 2. Row 1–2 and column 1 contain fixed factor levels for the respective column and row. Significant p-values are indicated with an asterisk.

https://doi.org/10.1371/journal.pone.0111889.s001

(DOCX)

Table S2.

P-values of pairwise comparisons between treatments which differ in presence or absence of sulphide. PERMANOVA (Euclidean distances, 9999 permutations) was used for the pairwise comparisons. The table part between double lines contains the variable factor sulphide (sulphide present = S, sulphide absent = No S) on which we performed pairwise comparisons for the variables in column 2. Significant p-values are indicated with an asterisk.

https://doi.org/10.1371/journal.pone.0111889.s002

(DOCX)

Table S3.

P-values of pairwise comparisons between treatments which differ in one level of the factor temperature. PERMANOVA (Euclidean distances, 9999 permutations) was used for the pairwise comparisons. The table part between double lines contains the variable factor temperature (16, 10 and 4°C) on which we performed pairwise comparisons for the variables in column 2. Significant p-values are indicated with an asterisk.

https://doi.org/10.1371/journal.pone.0111889.s003

(DOCX)

Table S4.

P-values of pairwise comparisons between treatments which differ in one level of the factor salinity. PERMANOVA (Euclidean distances, 9999 permutations) was used for the pairwise comparisons. The table part between double lines contains the variable factor salinity (25, 29.5 and 34 psu) on which we performed pairwise comparisons for the variables in column 2. Significant p-values are indicated with an asterisk.

https://doi.org/10.1371/journal.pone.0111889.s004

(DOCX)

Data S1.

Complete overview of the raw data as was collected during the experiment.

https://doi.org/10.1371/journal.pone.0111889.s005

(XLSX)

Acknowledgments

The data analysis of this research has benefitted from a statistical consult with Ghent University FIRE (Fostering Innovative Research based on Evidence). The authors would like to thank Marti Anderson (Department of Statistics, University of Auckland, New Zealand) and Jan Vanaverbeke (Research Group Marine Biology, Biology Department, Ghent University, Belgium) for statistical guidance in the use of PERMANOVA. S.D. acknowledges a postdoctoral fellowship from the Flemish Fund for Scientific Research (F.W.O).

Author Contributions

Conceived and designed the experiments: JVC SD TM AV. Performed the experiments: JVC. Analyzed the data: JVC SD TM AV. Contributed reagents/materials/analysis tools: JVC SD TM AV. Contributed to the writing of the manuscript: JVC SD TM AV.

Citation: Van Campenhout J, Derycke S, Moens T, Vanreusel A (2014) Differences in Life-Histories Refute Ecological Equivalence of Cryptic Species and Provide Clues to the Origin of Bathyal Halomonhystera (Nematoda). PLoS ONE 9(11): e111889. https://doi.org/10.1371/journal.pone.0111889

References

1. Bickford D, Lohman DJ, Sodhi NS, Ng PKL, Meier R, et al. (2007) Cryptic species as a window on diversity and conservation. Trends Ecol Evol 22: 148–155.

2. Pfenninger M, Schwenk K (2007) Cryptic animal species are homogeneously distributed among taxa and biogeographical regions. Bmc Evol Biol 7: 121.

3. Palumbi SR (1994) Genetic-Divergence, Reproductive Isolation, and Marine Speciation. Annu Rev Ecol Syst 25: 547–572.

4. Lonsdale DJ, Frey MA, Snell TW (1998) The role of chemical signals is copepod reproduction. J Marine Syst 15: 1–12.

5. Stanhope MJ, Connelly MM, Hartwick B (1992) Evolution of a Crustacean Chemical Communication Channel - Behavioral and Ecological Genetic-Evidence for a Habitat-Modified, Race-Specific Pheromone. J Chem Ecol 18: 1871–1887.

6. Knowlton N (2000) Molecular genetic analyses of species boundaries in the sea. Hydrobiologia 420: 73–90.

7. Derycke S, Remerie T, Backeljau T, Vierstraete A, Vanfleteren J, et al. (2008) Phylogeography of the Rhabditis (Pellioditis) marina species complex: evidence for long-distance dispersal, and for range expansions and restricted gene flow in the northeast Atlantic. Mol Ecol 17: 3306–3322.

8. Gomez A, Serra M, Carvalho GR, Lunt DH (2002) Speciation in ancient cryptic species complexes: Evidence from the molecular phylogeny of Brachionus plicatilis (Rotifera). Evolution 56: 1431–1444.

9. Chesson P (2000) Mechanisms of maintenance of species diversity. Annu Rev Ecol Syst 31: 343–366.

10. Hubbell SP (2006) Neutral theory and the evolution of ecological equivalence. Ecology 87: 1387–1398.

11. Adler PB, HilleRisLambers J, Levine JM (2007) A niche for neutrality. Ecol Lett 10: 95–104.

12. Zhang DY, Hanski I (1998) Sexual reproduction and stable coexistence of identical competitors. Journal of Theoretical Biology 193: 465–473.

13. Montero-Pau J, Serra M (2011) Life-Cycle Switching and Coexistence of Species with No Niche Differentiation. PLoS ONE 6: 1399–1410.

14. Leibold MA, McPeek MA (2006) Coexistence of the niche and neutral perspectives in community ecology. Ecology 87: 1399–1410.

15. Gabaldon C, Montero-Pau J, Serra M, Carmona MJ (2013) Morphological Similarity and Ecological Overlap in Two Rotifer Species. PLoS ONE 8: e57087.

16. Ortells R, Gomez A, Serra M (2003) Coexistence of cryptic rotifer species: ecological and genetic characterisation of Brachionus plicatilis. Freshwater Biol 48: 2194–2202.

17. Wellborn GA, Cothran RD (2007) Niche diversity in crustacean cryptic species: complementarity in spatial distribution and predation risk. Oecologia 154: 175–183.

18. Cothran RD, Henderson KA, Schmidenberg D, Relyea RA (2013) Phenotypically similar but ecologically distinct: differences in competitive ability and predation risk among amphipods. Oikos 122: 1429–1440.

19. Montero-Pau J, Ramos-Rodriguez E, Serra M, Gomez A (2011) Long-Term Coexistence of Rotifer Cryptic Species. PLoS ONE 6: e21530.

20. Cunning R, Glynn PW, Baker AC (2013) Flexible associations between Pocillopora corals and Symbiodinium limit utility of symbiosis ecology in defining species. Coral Reefs 32: 795–801.

21. Heip C, Vincx M, Vranken G (1985) The Ecology of Marine Nematodes. Oceanogr Mar Biol 23: 399–489.

22. Lambshead PJD, Boucher G (2003) Marine nematode deep-sea biodiversity - hyperdiverse or hype? J Biogeogr 30: 475–485.

23. Derycke S, Backeljau T, Vlaeminck C, Vierstraete A, Vanfleteren J, et al. (2007) Spatiotemporal analysis of population genetic structure in Geomonhystera disjuncta (Nematoda, Monhysteridae) reveals high levels of molecular diversity. Marine Biology 151: 1799–1812.

24. Derycke S, Remerie T, Vierstraete A, Backeljau T, Vanfleteren J, et al. (2005) Mitochondrial DNA variation and cryptic speciation within the free-living marine nematode Pellioditis marina. Mar Ecol-Prog Ser 300: 91–103.

25. De Oliveira DAS, Decraemer W, Holovachov O, Burr J, De Ley IT, et al. (2012) An integrative approach to characterize cryptic species in the Thoracostoma trachygaster Hope, 1967 complex (Nematoda: Leptosomatidae). Zool J Linn Soc-Lond 164: 18–35.

26. De Meester N, Derycke S, Moens T (2012) Differences in time until dispersal between cryptic species of a marine nematode species complex. Plos ONE 7: e42674.

27. De Meester N, Derycke S, Bonte D, Moens T (2011) Salinity effects on the coexistence of cryptic species: a case study on marine nematodes. Marine Biology 158: 2717–2726.

28. Bik HM, Porazinska DL, Creer S, Caporaso JG, Knight R, et al. (2012) Sequencing our way towards understanding global eukaryotic biodiversity. Trends Ecol Evol 27: 233–243.

29. Mokievsky VO, Filippova KA, Chesunov AV (2005) Nematode fauna associated with detached kelp accumulations in the subtidal zone of the White Sea. Oceanology+ 45: 689–697.

30. Vranken G, Herman PMJ, Heip C (1988) Studies of the Life-History and Energetics of Marine and Brackish-Water Nematodes: 1. Demography of Monhystera disjuncta at Different Temperature and Feeding Conditions. Oecologia 77: 296–301.

31. Trotter D, Webster JM (1983) Distribution and Abundance of Marine Nematodes on the Kelp Macrocystis integrifolia. Marine Biology 78: 39–43.

32. Portnova D, Haflidason H, Todt C (2010) Nematode species distribution patterns at the Nyegga pockmarks. 14th International Meiofauna Conference (FourtIMCo) Ghent, Belgium, pp 174.

33. Van Gaever S, Moodley L, de Beer D, Vanreusel A (2006) Meiobenthos at the Arctic Hakon Mosby Mud Volcano, with a parental-caring nematode thriving in sulphide-rich sediments. Mar Ecol-Prog Ser 321: 143–155.

34. Fonseca G, Derycke S, Moens T (2008) Integrative taxonomy in two free-living nematode species complexes. Biol J Linn Soc 94: 737–753.

35. Derycke S, Van Vynckt R, Vanoverbeke J, Vincx M, Moens T (2007) Colonization patterns of Nematoda on decomposing algae in the estuarine environment: Community assembly and genetic structure of the dominant species Pellioditis marina. Limnol Oceanogr 52: 992–1001.

36. Moens T (1999) Feeding ecology of free-living estuarine nematodes: an experimental approach. Gent: Universiteit Gent. 302 p.

37. Gerlach SA, Schrage M (1971) Life Cycles in Marine Meiobenthos - Experiments at Various Temperatures with Monhystera disjuncta and Theristus pertenuis (Nematoda). Marine Biology 9: 274–280.

38. Vranken G (1987) An autecological study of free-living marine nematodes. Academia Analecta 49: 71–97.

39. Vranken G, Tire C, Heip C (1988) The Toxicity of Paired Metal Mixtures to the Nematode Monhystera disjuncta (Bastian, 1865). Mar Environ Res 26: 161–179.

40. Vranken G, Tire C, Heip C (1989) Effect of Temperature and Food on Hexavalent Chromium Toxicity to the Marine Nematode Monhystera disjuncta. Mar Environ Res 27: 127–136.

41. Vranken G, Vanderhaeghen R, Van Brussel D, Heip C, Hermans D (1984) The toxicity of mercury on the free-living marine nematode Monhystera disjuncta Bastian, 1865. Coomans, A and Heip, CHR (Ed) (1984), Coordinated Research Actions Oceanography: progress report 1983 Ecological and ecotoxicological studies of benthos of the Southern Bight of the North Sea. Brussel: Programmatie van het Wetenschapsbeleid. pp. 186–201.

42. Vranken G, Vanderhaeghen R, Heip C (1985) Toxicity of Cadmium to Free-Living Marine and Brackish Water Nematodes (Monhystera microphthalma, Monhystera disjuncta, Pellioditis marina). Dis Aquat Organ 1: 49–58.

43. Tchesunov AV, Portnova DA, Van Campenhout J Description of two free-living nematode species of Halomonhystera disjuncta complex (Nematoda: Monhysterida) from two peculiar habitats in the sea. Helgoland Mar Res. DOI:10.1007/s10152-014-0416-1.

44. Van Campenhout J, Derycke S, Tchesunov A, Portnova D, Vanreusel A (2013) The Halomonhystera disjuncta population is homogeneous across the Håkon Mosby mud volcano (Barents Sea) but is genetically differentiated from its shallow-water relatives. Journal of Zoological Systematics and Evolutionary Research 52: 203–216.

45. Van Gaever S, Moodley L, Pasotti F, Houtekamer M, Middelburg JJ, et al. (2009) Trophic specialisation of metazoan meiofauna at the Håkon Mosby Mud Volcano: fatty acid biomarker isotope evidence. Marine Biology 156: 1289–1296.

46. Sauter EJ, Muyakshin SI, Charlou JL, Schluter M, Boetius A, et al. (2006) Methane discharge from a deep-sea submarine mud volcano into the upper water column by gas hydrate-coated methane bubbles. Earth Planet Sc Lett 243: 354–365.

47. Ackerly DD (2003) Community assembly, niche conservatism, and adaptive evolution in changing environments. Int J Plant Sci 164: S165–S184.

48. Wiens JJ (2004) Speciation and ecology revisited: Phylogenetic niche conservatism and the origin of species. Evolution 58: 193–197.

49. Harvey PH, Pagel MD (1991) The comparative method in evolutionary biology. Oxford UK: Oxford University Press. 248p.

50. Moens T, Vincx M (1998) On the cultivation of free-living marine and estuarine nematodes. Helgolander Meeresun 52: 115–139.

51. Anderson M, Gorley RN, Clarke KR (2008) PERMANOVA+ for PRIMER: Guide to Software and Statistical Methods: PRIMER-E. Plymouth.

52. Schmidt AL, Wysmyk JKC, Craig SE, Lotze HK (2012) Regional-scale effects of eutrophication on ecosystem structure and services of seagrass beds. Limnol Oceanogr 57: 1389–1402.

53. Anderson M (2004) PERMDISP: a FORTRAN computer program for permutation analysis of multivariate dispersions (for any two-factor ANOVA design) using permutation tests. Department of Statistics, University of Auckland, New Zealand.

54. Tahseen Q (2012) Nematodes in aquatic environments: adaptations and survival strategies. Biodiversity Journal 3: 13–40.

55. Levins R (1968) Evolution in Changing Environments: Some Theoretical Explorations: Princeton University Press. 120p.

56. Bottcher ME, Hespenheide B, Llobet-Brossa E, Beardsley C, Larsen O, et al. (2000) The biogeochemistry, stable isotope geochemistry, and microbial community structure of a temperate intertidal mudflat: an integrated study. Cont Shelf Res 20: 1749–1769.

57. Ferris H, Lau S, Venette R (1995) Population Energetics of Bacterial-Feeding Nematodes - Respiration and Metabolic Rates Based on Co2 Production. Soil Biol Biochem 27: 319–330.

58. Levinton JS (1970) The paleoecological significance of apportunistic species. Lethaia 3: 69–78.

59. Moens T, Vincx M (2000) Temperature and salinity constraints on the life cycle of two brackish-water nematode species. J Exp Mar Biol Ecol 243: 115–135.

60. Moens T, Vincx M (2000) Temperature, salinity and food thresholds in two brackish-water bacterivorous nematode species: assessing niches from food absorption and respiration experiments. J Exp Mar Biol Ecol 243: 137–154.

61. Boffé A (1985) Vergelijkend toxicologisch onderzoek naar sublethale effecten bij de mariene nematode Monhystera disjuncta (Bastian 1865). University Ghent: pp. 70.

62. Lowe CD, Kemp SJ, Diaz-Avalos C, Montagnes DJS (2007) How does salinity tolerance influence the distributions of Brachionus plicatilis sibling species? Marine Biology 150: 377–386.

63. Nicholls B, Racey PA (2006) Habitat selection as a mechanism of resource partitioning in two cryptic bat species Pipistrellus pipistrellus and Pipistrellus pygmaeus. Ecography 29: 697–708.

64. Smith KE, Thatje S (2012) The Secret to Successful Deep-Sea Invasion: Does Low Temperature Hold the Key? PLoS ONE 7: e51219.

65. Bik HM, Thomas WK, Lunt DH, Lambshead PJD (2010) Low endemism, continued deep-shallow interchanges, and evidence for cosmopolitan distributions in free-living marine nematodes (order Enoplida). Bmc Evol Biol 10: 389.

66. Bernardino AF, Levin LA, Thurber AR, Smith CR (2012) Comparative Composition, Diversity and Trophic Ecology of Sediment Macrofauna at Vents, Seeps and Organic Falls. PLoS ONE 7: e33515.

67. Zhao YL, Wang DY (2012) Formation and regulation of adaptive response in nematode Caenorhabditis elegans. Oxid Med Cell Longev 2012: 564093.

68. Ott J, Bright M, Bulgheresi S (2004) Symbioses between marine nematodes and sulfur-oxidizing chemoautotrophic bacteria. Symbiosis 36: 103–126.

69. Levin LA (2005) Ecology of cold seep sediments: Interactions of fauna with flow, chemistry and microbes. Oceanography and Marine Biology - an Annual Review, Vol 43 43: 1–46.

70. Thiermann F, Vismann B, Giere O (2000) Sulphide tolerance of the marine nematode Oncholaimus campylocercoides - a result of internal sulphur formation. Mar Ecol Prog Ser 193: 251–259.

71. Westram AM, Jokela J, Keller I (2013) Hidden Biodiversity in an Ecologically Important Freshwater Amphipod: Differences in Genetic Structure between Two Cryptic Species. PLoS ONE 8: e69576.