Closely related intertidal and deep-sea Halomonhystera species have distinct fatty acid compositions
Helgoland Marine Research volume 70, Article number: 8 (2016)
The deep-sea free-living nematode Halomonhystera hermesi, dominant in the sulphidic sediments of the Håkon Mosby mud volcano (1280 m, Barent sea slope), is part of the mainly estuarine Halomonhystera disjuncta species complex consisting of five cryptic species (GD1–GD5). Cryptic species have a very similar morphology raising questions on their specific environmental differences. This study analyzed total fatty acid (FA) compositions of H. hermesi and GD1, one of H. hermesi’s closest relatives. Additionally, we experimentally investigated the effect of a temperature reduction, salinity increase and sulphide concentrations on GD1’s FA composition. Because nematodes are expected to have low amounts of storage FA, total FA compositions most likely reflect FA contents of cellular membranes. The deep-sea nematode H. hermesi had significantly lower saturation levels and increased highly unsaturated fatty acid (HUFAs) proportions due to the presence of docosahexanoic acid (DHA—22:6ω3) and higher eicosapentaenoic acid (EPA—20:5ω3) proportions. HUFAs were absent in H. hermesi’s food source indicating the ability and need for this nematode to synthesize HUFAs in a deep-sea environment. Our experimental data revealed that only a decrease in temperature resulted in lower saturated fatty acids proportions, indicating that the FA content of H. hermesi is most likely a response to temperature but not to sulphide concentrations or salinity differences. In experimental nematodes, EPA proportions were low and DHA was absent indicating that other factors than temperature, salinity and sulphides mediate the presence of these HUFAs in H. hermesi.
Organisms require adaptive responses when presented to changes within their environment to survive and reproduce successfully. Environmental adaptation can result in morphological, physiological, behavioural or developmental differences. It is known that environmental changes can alter cell membrane fluidity, which could impede membrane function and stability, and performance of membrane-associated enzymes . Cellular membranes function as a physiological barrier with a limited permeability and mainly consist of phospholipids, of which fatty acids are the most variable part . It is, therefore, essential that the liquid-crystalline structure of the membrane is maintained by regulation of the FA composition to acclimatize cells to changing environmental conditions, as is evidenced from temperature effects in poikilotherms . Alterations in the composition of membrane phospholipids in response to environmental changes are of critical importance for the function and integrity of cell membranes, and is a general phenomenon in animals, known as “homeoviscous adaptation” (HVA) or “homeophasic adaptation” [30, 39, 61].
Nematodes, one of the most numerous of all Metazoa in terms of species richness , exist in all habitats that can support life and demonstrate remarkable abilities to withstand stress and adverse conditions [44, 62]. It has been suggested that marine nematodes do not accumulate lipids for energy storage  indicating that nematode fatty acid composition reflects membrane fatty acid composition, however, information remains scant. In this study we will focus on the marine free-living nematode Halomonhystera disjuncta (formerly known as Geomonhystera disjuncta), which is known for its high resistance to environmental stress, especially to heavy metals and temperature changes [25, 75–77]. H. disjuncta and sister species has been reported both in shallow-water [15, 47, 74] and in deep-sea environments [53, 71]. Mitochondrial and nuclear DNA sequences as well as morphometric analysis revealed that H. disjuncta from the Western Scheldt and the North Sea actually consist of five cryptic species (GD1–GD5)  with nearly imperceptible morphometric differences . The same marine nematode has been morphologically identified as the dominant nematode species at the sulphide rich bacterial mats of two methane cold seeps: the Nyegga pockmarks (Nordic Norwegian margin, 730 m water depth)  and Håkon Mosby mud volcano (HMMV, Barents Sea slope, 1280 m water depth) . The HMMV relative is closely related to the intertidal Halomonhystera cryptic species complex but can be distinguished from its shallow-water relatives both morphometrically and with the use of genetic markers , and was, therefore, described as a new species named H. hermesi . Phylogenetic relationships showed that H. hermesi is most closely related to GD1 and GD4. Interestingly, GD1 appeared to be more resistant to simulated cold seep conditions compared to GD2–GD3 . However, lack of any information on the biochemical resources (e.g. FA composition) hampers our understanding of the unique biology of closely related species in different environment.
The deep sea is a relatively harsh environment with constantly low ambient temperatures (~2 °C), high water pressure and absence of light. Furthermore, cold seeps are characterized by methane seepage through sediment and rely on chemosynthetically derived energy. They often have high concentrations of sulphides, as a result of microbial sulphate reduction coupled to anaerobic oxidation of methane . Most studies on adaptive changes of the FA composition to deep-sea conditions have been performed on barotolerant bacteria under different experimental conditions and have revealed increased proportions of mono- and poly-unsaturated FAs in response to increased pressure and/or reduced temperature [13, 67, 80]. Consistent with bacterial responses, a comparison of mitochondrial membranes isolated from shallow and deep-sea fish species has identified similar changes in unsaturation levels [10, 73]. Moreover, highly unsaturated FAs (HUFAs) are believed to be essential in piezophilic adaptation of cellular membranes  and docosahexaenoic acid (DHA—C22:6ω3) was hypothesized to be the key to HVAs to depth and temperature changes in vertically migrating planktonic copepods . As such, we expect higher proportions of unsaturated fatty acids and HUFAs in H. hermesi compared to GD1.
The second part of this study consists of a comparison of the FA composition of GD1 nematodes under simulated estuarine and bathyal seep conditions. We investigated the effect of three abiotic factors salinity, temperature and sulphide on total fatty acid composition. Each of the three abiotic factors varied from control conditions (salinity of 25, 16 °C and no sulphide) to environmental conditions representative of the HMMV (salinity of 34–35, −0.9 °C and hydrogen sulphide concentrations up to 1 mM) [58, 72]. As previously mentioned, temperature decrease has a clear effect on the FA composition in cellular membranes resulting in an increase in the unsaturated fatty acid proportions. Such a response is uniform and we, therefore, anticipate similar outcomes. The effect of increased salinity, however, is an increase of saturated FA levels and decreased proportions of unsaturated FAs as described for the bacteria Shewanella  and Halomonas species , and several plant species from the Ephedra genus , but information remains scant. Moreover, salinity effects were also found to be opposite between different plant varieties of the same species . Lastly, no information on the effects of sulphides on FA composition currently exists. We can, however, consider sulphides to be a toxic compound because H. disjuncta is unable to survive concentrations higher than 1.5 mM (personal observations). In the presence of a toxic compound in sublethal concentrations, cells start to decrease membrane fluidity by increasing the degree of FA saturation [17, 48, 60]. The experimental data, discussed within this paper, allows investigating the individual effects of environmental factors on FA composition. In addition, it can aid in the correct interpretation of observed differences in FA composition between H. hermesi and GD1 sampled from the HMMV and Western Scheldt estuary, respectively.
Intertidal GD1 nematode cultures
Decaying Fucus vesiculosus was collected from the Paulina tidal flat, located along the southern shore of the polyhaline zone of the Western Scheldt estuary (southwestern part of The Netherlands). Macroalgal pieces were transferred to Petri dishes containing 0.8 % nutrient:bacto agar (ratio of 1/7) prepared in artificial seawater  with a salinity of 25 and placed at a constant temperature of 16 °C. Two weeks after inoculation, one single gravid female was transferred to a new Petri dish (5.5 cm inner diameter), containing the same agar medium, to establish monospecific GD1 cultures. Frozen-and-thawed Escherichia coli K12 was added as a food source and cultures were placed at a constant temperature of 16 °C. Monospecificity and species identity were verified, shortly after inoculation, by PCR and sequencing the ITS gene. Primers and thermal-cycling conditions can be found in Derycke et al. .
Intertidal GD1 nematodes were harvested from the monospecific cultures and transferred to rehydrated (salinity of 25), defaunated F. vesiculosus (48 h at 60 °C). They were allowed to grow for ca. 4 weeks in three litterbags (mesh size 200 µm) at 8 °C and a salinity of 25. Litterbags were then placed at the Paulina tidal flat (Western Scheldt estuary, 51°20′56.79″N, 3°43′29.56″E) for 72 h. Following this period, macroalgae were washed over two stacked sieves (top sieve: 1 mm, bottom sieve: 32 µm) and fauna was retained on the bottom sieve. The sieve was rinsed with sieved seawater (32 µm) and replicate samples were frozen in liquid nitrogen and stored at −80 °C awaiting further analysis. Throughout the rest of the study, GD1 samples extracted from the Paulina tidal flat, will be referred to as GD1P.
Deep-sea sediment samples were collected at the HMMV during the Maria S. Merian cruise 2010 using a TV-guided multi-corer. The microbial mats, the preferred habitat of H. hermesi, were selected and sampled from three locations (Station MSM16-2_802-1, 72°0.17′N, 14°43.88′E; Station MSM16-2_829-1, 72°0.16′N, 14°43.94′E; Station MSM16-2_831-1, 72°0.14′N 14°43.94′E). Cores from these locations were immediately sliced on board from 0 to 10 cm with a one cm interval and samples for fatty acid composition were stored at −80 °C.
The top cm of three deep-sea samples and part of the three intertidal samples, were thawed and nematodes were extracted by density gradient centrifugation, using Ludox (a colloidal silica polymer; specific gravity 1.18) as a flotation medium . For each replicate, ca. 333 nematodes of each sex (male and female) and 333 juveniles were manually picked out using a binocular microscope to easily select Halomonhystera species based on their morphology. Per replicate, juveniles and adults were pooled to encompass all life stages and sexes present in the population. Nematodes were then rinsed twice in sterile Mili-Q for at least 1 h to remove most adhering bacteria and transferred to glass 5 ml GC vials in sterile Mili-Q and stored at −80 °C prior to fatty acid extraction.
GD1 experimental set up
The effect of three abiotic factors (salinity, temperature and sulphide presence/absence) on fatty acid composition was tested by a replicated fully crossed design with three factors: a salinity of 25 and 34, temperatures of 16 and 4 °C and presence or absence of sulphide (0 and 1 mM). Because simulated cold deep-sea temperatures resulted in occasional crystallization of the agar medium, a temperature of 4 °C was used. Experimental control conditions were 16 °C and a salinity of 25 because stock cultures are maintained at these conditions. The experiment was conducted on experimental microcosms which consisted of Petri dishes (5.5 cm i.d.) containing 1.5 ml 0.8 % nutrient:bacto agar (1/7). The sulphidic medium consisted of the same agar as described above to which we added sodium thiosulphate (Na2S2O3) and sodium sulphide (Na2S) in a final concentration of 1 mM for both reagents; 350 µl E. coli K12 (ca. 8 × 108 cells/ml) was added as a food source, sufficient for the duration of the experiment. In addition, 350 µl of the E. coli food source was also stored at −80 °C for further fatty acid analysis.
The experiment was started by manually picking out 60 females and 40 males (F0) from the monospecific GD1 stock cultures which were deposited in the microcosms. Petri dishes were closed with parafilm and maintained in temperature-controlled incubators without light for the whole length of the experiment. Two days after the minimum generation times (9–11, 10–15 and 28–39 days at 16, 10 and 4 °C, respectively) , 333 F1 nematodes (111 juveniles, males and females) from each replicate were manually picked out and washed as described above. Because we required 1000 nematodes for further analysis and microcosms did not yield many nematodes, the replicate samples were pooled, frozen in 5 ml GC vials and stored at −80 °C.
Analytic techniques and data treatment
Deep-sea, intertidal and experimental GD1 samples, and the E. coli food source were lyophilized for 24 h before extraction. Extraction and esterification of fatty acids was performed according to a modified one-step method of Abdulkadir and Tsuchiya . Briefly, samples were mixed with 700 µl of a 2.5 % H2SO4–methanol solution in a capped tube and heated at 80 °C for 90 min. After cooling to room temperature, 350 µl hexane and 350 µl of a 0.95 % NaCl solution were added. After stirring, samples were centrifuged for 10 min at 1000 rpm (650 g). Of the two phases which formed, the upper hexane layer, containing the fatty acid methyl esters (FAMEs), was transferred to a capped vial to be injected into the gas chromatographer (GC) for FAME analysis. Five microlitre of the FAMEs solution was injected on a HP88 column (Agilent J&W) at a temperature of 250 °C and separated using gas chromatography (Hewlett Packard 6890N) coupled to a mass spectrometer (HP 5973) to obtain ion profiles for identifications. Oven temperatures were the same as described in De Troch et al. .
FAMEs identification was based on the comparison of the mass spectral fragmentation spectra and retention times to authentic standards and mass spectral libraries (WILEY and in-house libraries) using the MSD ChemStation software (Agilent Technologies). Linear regression of peak areas of external standards (Supelco #47885, Sigma-Aldrich) in known concentrations ranging from 5 to 250 µg ml−1 were used for FAME quantification. The slope of the linear regression of the external standards and amount of hexane is employed to quantify each individual FAME according to De Troch et al. . Ratios between individual FAME and C16:0 peak areas were determined and ratios <0.02 were deemed unreliable and, therefore, removed . In addition, unassigned peaks were removed from the dataset. The nomenclature of the fatty acids follows the universal X:YωZ form, were X represents the number of carbon atoms, Y equals the number of double bonds and Z gives the position of the double bond counting from the terminal methyl group . Because we were unable to reliably distinguish isoforms and anteisoforms, these peaks, of each respective FA, were summed and “(ante)iso” was added to the universal X:YωZ form throughout the rest of this paper. FAME concentrations were standardized to the total amount of nematodes added, resulting in pico grams per nematode. Because organisms have the ability for a genotype to match their phenotype to different environments which can influence the total amount of FAs under different experimental conditions, relative proportions of each quantified FAME to the total amount of pg fatty acids were determined. Similarly, relative FA proportions of the E. coli food source (added to experimental microcosms) were determined as described above.
Statistical analysis on Halomonhystera hermesi and GD1P
To test the difference in relative fatty composition between Halomonhystera hermesi and GD1P, a one factor permutational based multivariate analysis of variance , on the basis of Bray–Curtis similarity with 9999 permutations, was performed. The results were visualized by non-metric multidimensional scaling (MDS). Samples were clustered based on group average and a SIMPER analysis was used to identify which fatty acid accounted for the observed dissimilarity between species. PERMDISP  was executed to test the homogeneity of multivariate dispersions. Because of the limited level of replication, a Monte Carlo test was performed at all times and all p values in this study have been FDR-corrected in R version 3.0.2 .
Comparing relative proportions of FAs based upon the total fatty acid content between species might not truly reflect the composition of SFAs, MUFAs and PUFAs. Therefore, fatty acids were divided into three classes based on the saturation level: saturated fatty acids (SFAs), mono-unsaturated fatty acids (MUFAs) and poly-unsaturated fatty acids (PUFAs). The total values for each class are the sum of values for the following major components: SFAs (C12:0, C13:0, C14:0, (ante)isoC15:0, C15:0, (ante)isoC16:0, C16:0, (ante)isoC17:0, C17:0, C18:0, C20:0 and C22:0), MUFAs (C14:1ω5, C15:1ω5, C16:1ω7, C17:1ω7, C18:1ω12, C18:1ω9, C18:1ω7, C20:1ω9, C22:1ω9 and C24:1ω9) and PUFAs [C16:2ω6, C18:2ω6, C18:3ω3, C18:3ω6, C20:2ω6, C20:3ω6, C20:4ω6, c20:5ω3 (EPA) and C22:6ω3 (DHA)]. Relative values of each FA were determined based upon the total sum of its respective FA class. The difference in SFA-, MUFA- and PUFA-composition, was statistically analysed as described above (Bray–Curtis similarity). The difference between both species of each FA proportion per class was determined using Euclidean distances with 9999 permutations and a PERMANOVA was performed.
In addition, PUFAs with a chain length ≥20 carbons and ≥3 double bonds, known as HUFAs, were summed, SFA/UFA ratio (with UFA defined as the sum of MUFA and PUFA) and the unsaturation index were determined. The unsaturation index is the sum of MUFA and PUFA proportions multiplied by the respective amount of unsaturated bonds of each respective FA. To test the difference between SFA/UFA ratios, unsaturation indices and HUFA proportions, Euclidean distances was used in PERMANOVA. All analyses were performed within PRIMER v6 with PERMANOVA+ add-on software . Figures were made with Graphpad prism 5.
Statistical analysis of the experiment (GD1)
Relative proportions of FAs were determined in the same way as described above and the same parameters were calculated. We have performed a fully crossed PERMANOVA design with three fixed factors (on all parameters separately: total FA-, SFA-, MUFA-, PUFA- and HUFA-composition, SFA/UFA ratios and unsaturation indices). The three factors were temperature (16 and 4 °C), salinity (25 and 34) and sulphide (0 and 1 mM). Because we lack proper replication, we cannot reliably estimate factor interaction effects and, therefore, pooled the two- and three-way interaction terms, thereby overestimating the true residual. This statistical test allowed us to reliably investigate the effect of each individual factor.
FA composition of GD1P and H. hermesi
The total fatty acid composition was significantly different between both species (p = 0.0054), and homogeneity of multivariate dispersions was found (p = 0.1044). Relative values for all FAs, based upon the total FA content can be found in Table 1. Replicate samples from H. hermesi clustered together with a similarity of >90 % and GD1P replicates clustered together with a similarity >80 % (Fig. 1). Total FA composition between both species was 68.65 % similar. SIMPER revealed that MUFA C16:1ω7 was present in significantly higher proportions (p = 0.0032) in H. hermesi (23.56 %) compared to GD1 (9.07 %) and was most important in explaining the dissimilarity (23.11 %) between both species. This analysis further revealed that a total of 17 FAs contributed to explain >90 % of the dissimilarity (Table 2).
Classifying each FA into their respective group (SFAs, MUFAs and PUFAs) showed that MUFAs were the most abundant class for both species followed by PUFAs and SFAs (Fig. 2). SFAs C13:0 and ante(iso)C17:0, and MUFA C15:1ω5 were only observed in GD1P (Table 1; Fig. 3a). However, C13:0 was only identified in one GD1p sample in very low proportions (0.13 %). SFA C22:0, MUFAs C14:1ω5 and C18:1ω12, and PUFAs C16:2ω6, C18:3ω6, C20:3ω6 and C22:6ω3 (DHA) were exclusively observed in H. hermesi (Fig. 3b, c). Of these uniquely observed FAs, DHA had the highest relative proportion (4.23 %) of the total fatty acid content (Fig. 3c). The unsaturation index and SFA/UFA ratio of GD1p was respectively lower (p = 0.0021) and higher (p = 0.0032) in GD1P compared to H. hermesi (Table 3).
The relative amount of SFAs for H. hermesi (8.85 %) was significantly lower than the 28.43 % for GD1P (p = 0. 0023). C16:0 was the dominant SFA (Fig. 3a) in both GD1P (31.0 %) and H. hermesi (42.5 %). In addition, SFA composition was significantly different between both species (p = 0.0017). Except for C13:0 (p = 0.40557), C15:0 (p = 0.26365) and ante(iso)C16:0 (p = 0.05720) all SFAs were significantly different between both species (p values <0.05, Additional file 1: Table S1). C14:0 (p = 0.0102) and ante(iso)C15:0 proportions (p = 0.0023) were significantly higher in GD1P, while C16:0, C17:0, C18:0 and C20:0 were significantly higher in H. hermesi (Fig. 3a).
The low relative SFA proportion in H. hermesi was compensated by significantly higher MUFA—(11.83 %, p = 0.0030) and PUFA (7.75 %, p = 0.0392) proportions in comparison to GD1P (Fig. 2). The MUFA composition in both species was dominated by C16:1ω7 and C18:1ω9 (Fig. 3b) which were responsible for 67.1 and 79.9 % of the total MUFA composition of GD1P and H. hermesi, respectively. The MUFA composition was significantly different between both species (p = 0.0387). MUFAs C18:1ω9, C18:1ω7, C20:1ω9, C24:1ω9 were not significantly different, while all other MUFAs were significantly different between both species (p values <0.05, Additional file 1: Table S1). The MUFA composition of H. hermesi consisted of circa two times higher proportions of C16:1ω7 (p = 0.0102) compared to GD1P (Fig. 3b), and was responsible for 40.2 % of the MUFA dissimilarity between both species.
PUFA compositions were also significantly different between both species (p = 0.0023) and were dominated by C18:2ω6 (Fig. 3c), responsible for 40.14 % of the dissimilarity. The PUFA composition of H. hermesi showed a tendency towards HUFAs (≥20 carbons, ≥3 double bonds) as was indicated by significantly higher HUFA values (p = 0.0032) for H. hermesi (Table 3). This was a result of higher EPA (C20:5ω3—EPA; p = 0.0017) proportions and the presence of DHA (C22:6ω3—DHA), a HUFA that was absent in GD1P (Fig. 3c; Table 1). EPA and DHA were responsible for 11.72 and 9.29 % of the PUFA dissimilarity between both species.
FA composition of GD1 under different experimental conditions
Similar to GD1P, SFAs had the lowest proportions in all experimental microcosms (Fig. 4; Table 1). SFA proportions were significantly affected by temperature (p = 0.0368) resulting in 9.9–14.6 % lower SFA proportions at 4 °C compared to 16 °C (Fig. 4). Similar increases in UFA proportions were observed (Fig. 4) when temperature decreased (p = 0.0069). Thus, the SFA/UFA ratio always significantly decreased (p = 0.0074) with decreasing temperatures while the unsaturation index showed an inverse response (p = 0.0444, Table 3). Increasing the salinity and adding sulphides did not have a significant effect on SFA proportions.
SFAs mainly consisted of (ante)isoform FAs and odd chain length FAs such as C15:0 and C17:0 (Fig. 5a) and the most abundant MUFAs were C16:1ω7 and C18:1ω7 (Fig. 5b). The effect of varying abiotic factors did not reveal any discernible trend in SFA and MUFA composition. PUFA compositions were very similar and highly dominated by C18:2ω6 (71.5–80.7 %) which is similar to GD1P (67.4 %) and H. hermesi (59.3 %). The PUFA composition (Fig. 5c) of experimental nematodes lacked C18:3ω3, compared to GD1P, DHA was not observed in the experimental nematodes and only trace amounts of EPA (<1 %) were present.
FA composition of E. coli
The most abundant fatty acids for the E. coli food source were SFAs (93.5 %) and no PUFAs could be identified (Table 1). SFA composition was similar to experimental nematodes (Fig. 5a) and was dominated by (ante)isoC15:0 (37.3 %) and C16:0 (29.4 %). Meanwhile, the composition of the MUFAs was very different from experimental GD1 nematodes (Fig. 5b) and consisted of C16:1ω7 (32.0 %), C18:1ω9 (51.1 %) and C18:1ω7 (16.8 %).
Fatty acid compositions of H. hermesi and GD1P indicate a homeoviscous adaptation to the deep-sea environment
Survival and successful reproduction of organisms relies on the ability to respond and adapt to the often varying conditions imposed by their environment. Fatty acid composition is known to rapidly (24–72 h) respond to changes in diet  and even daily/hourly fluctuations in FA composition have been observed previously [18, 23]. The barrier function of the cytoplasmic membrane is known to depend critically on the physical state of lipid bilayers , making it susceptible to environmental changes. The deep sea is a constant cold environment. When biological membranes are sufficiently cooled many biochemical reactions may be impaired which results in a decreased fluidity of lipids, a factor of primary importance to cell membrane function . Increased MUFA proportions are most commonly reported as the agents involved in maintaining membrane order, with the addition of further double bonds having little effect on the phase transition properties and hence fluidity of the lipids . Higher 16:1ω7 proportions in H. hermesi could, therefore, be important in HVA to low temperatures. However, the high 16:1ω7 proportions could also have originated from the deep-sea HMMV bacteria Beggiatoa sp., the proposed food source for H. hermesi which is rich in 16:1ω7 . This MUFA has formerly been reported in Beggiatoa species from other cold seeps , suggesting that high 16:1ω7 proportions might be important in HVA to cold-seep environments for Beggiatoa and in extension also for organisms feeding on the bacteria such as H. hermesi. Low temperatures can also result in an increase in PUFA proportions in nematodes , which can explain the increase in PUFA proportion observed in H. hermesi.
High hydrostatic pressure can have similar effect on FA composition as low temperatures, i.e. reduction in SFA proportion and an increase of UFA proportion [6, 13, 14, 82]. Lipids are known to be more compressible than proteins and are, therefore, in comparison more sensitive to pressure effects . Accordingly, many deep-sea organisms modulate their membrane fluidity and composition in response to pressure . Comparison of deep-sea and shallow-water fish revealed increased proportions of membrane fluidizing unsaturated fatty acids . An increase in MUFAs and PUFAs and decrease in SFAs was also observed in vertically migrating marine planktonic copepods with increasing depth . Our results revealed very similar effects in FA proportions, i.e. the proportion of saturated FAs in the deep-sea nematode H. hermesi is nearly 20 % lower compared to its intertidal relative GD1P, which is compensated by higher amounts of MUFAs and PUFAs.
Fatty acid proportions of H. hermesi revealed the presence of the highly unsaturated fatty acid (HUFA) DHA (DHA—C22:6ω3), which was absent in GD1P, and higher concentrations of the HUFA EPA (EPA—C20:5ω3). Marine animals are known to contain high levels of these long-chain ω3 PUFAs [2, 65] and recent data suggest that also marine nematodes are rich in HUFAs  as was observed for H. hermesi. HUFAs are essential for normal functioning of cell membranes and membrane-bound enzyme systems, crucial regulators for animal cell membrane physiology . Moreover, these HUFAs were hypothesized to important in HVA to high pressure and low temperatures . DHA proportions also increased with depth in vertically migrating copepods suggesting its key role in maintaining membrane order in cold, high pressure conditions in the deep sea . The absence of DHA in GDP, contrasting to its presence in H. hermesi, therefore, suggests that DHA might be essential in HVA to a deep sea environment. Moreover, interaction between temperature and pressure on DHA proportions, consistent with homeoviscous adaption, has been recently reported .
Because both species are bacterivorous, one could expect DHA and EPA to originate from their respective food source. Even though DHA is uncommon in most bacteria, it is present in deep-sea bacteria, as is EPA [14, 83]. Studies on deep-sea bacteria [13, 50, 81] and growth experiments of bacterial cultures under high pressure  revealed the importance of increased proportions of the HUFA EPA (EPA—C20:5ω3) under high pressure, which is in agreement with our results. However, a trophoecolocigal study on fatty acid compositions of Beggiatoa sps, the proposed food source of H. hermesi, revealed that HUFAs were absent . Unlike most metazoans, nematodes are, however, able to biosynthesize HUFAs from acetate [7, 56] and from shorter fatty acid precursors . Interestingly, a transcriptome comparison between GD1P and H. hermesi  disclosed an overrepresentation of the KEGG pathway “fatty acid elongation” in H. hermesi. Given the absence of these EPA and DHA in the food source of H. hermesi , the nematode has most likely the ability to synthesize these HUFAs. EPA and DHA in addition can be assimilated into triglyceride forms and stored until they are required for incorporation in the membrane, beta-oxidation, signalizing activity, etc. [55, 79]. It, therefore, remains to be investigated what the primary function of these HUFAs is in H. hermesi.
It must be noted that other environmental differences might influence the FA composition of membranes. High salinities, for example, have an opposite effect on FA composition compared to high pressure and a low temperature [3, 38, 49]. A recent study revealed that increased salinities can also result in increased DHA and EPA proportions . Moreover, synchronous effects of temperature, hydrostatic pressure, and salinity on FA composition have been previously observed . Given the different (a)biotic features between both environments and the similar effect of temperature and pressure on the FA composition it is difficult to assess which abiotic factor is most responsible for the higher UFA proportion in H. hermesi. Therefore, we have experimentally tested the effect of different factors (temperature, salinity and sulphides) on the FA composition of GD1 to elucidate which factors significantly affected UFA proportions.
FA composition of GD1 is affected by temperature, but not by salinity and sulphide changes
Our results revealed that reduced growth temperatures result in greater proportions of unsaturated fatty acids, which is in accordance with previously observed temperature effects in nematodes . However, no effect of salinity and sulphide was observed. This indicates that temperature is most likely the main factor responsible for the difference in UFA proportions between GD1P and H. hermesi, while salinity and sulphides are of lesser importance. Our results further revealed that DHA was absent in all our experimental nematodes and only trace amounts of EPA have been found. This suggests that the high HUFA proportions of H. hermesi, in comparison to GD1P, is most likely a response to other factors such as pressure and less the result of differences in temperature, sulphides and salinities between both environments.
Our experimental data further showed high proportions of C15:0 and C17:0, and especially high (ante)isoform proportions of several FAs. Metazoa almost exclusively synthesize even-chain fatty acids  while bacteria are able to synthesize odd chain length, ante- and iso-forms of FAs, and branched fatty acids . The presence of C17:0 and C15:0 has been previously observed in bacterivorous nematodes, be it in low amounts , in agreement with our study. It is, however, well established that the fatty acid composition of an organism most often reflect their diet [45, 57, 72]. Meanwhile, transfer up the food chain of C15:0 and C17:0 have been seen in copepods feeding on bacterivorous ciliates  and an unusually high amount (up to 50 %) of odd chain length FAs have been observed in the amphipod Pontoporeia femorata . Both observations suggest that an organism is able to incorporate bacterial FAs. Odd chain length FAs and ante(iso) FAs are indeed present in the E. coli food source for GD1 and reflect the SFA composition of GD1 suggesting that the latter incorporates fatty acids from its food source to use in its membranes. In addition, absence of PUFAs in E. coli does clearly indicate and confirm that nematodes are able to biosynthesize long unsaturated fatty acids from shorter chain fatty acids as has previously been observed .
In this study we have shown that the fatty acid composition of closely related deep-sea and intertidal nematodes (GD1P) clearly differs. The deep-sea nematode H. hermesi had higher unsaturated fatty acid proportions which could be attributed to the cold deep-sea environment. In addition, the saturated-, monounsaturated-, and poly-unsaturated fatty acid composition of both species is different. H. hermesi showed significantly higher 16:1ω7 proportions, which might be important in HVA to cold-seep environments. The deep-sea nematode had elevated HUFAs such as EPA and DHA, suggesting that these HUFAs are important in HVA to the cold deep-sea environment as was previously reported for other deep-sea organisms. Moreover, the lack of DHA and trace amounts of EPA in our experimental microcosms suggest that these HUFAs are a response to factors such as pressure and not a consequence of differences in temperature, salinity and sulphide concentrations. Subjecting GD1 to changing abiotic factors such as temperature, salinity and sulphides, revealed that only temperature significantly increased the unsaturated fatty acid proportions. This indicates that temperature is important in explaining differences in UFA proportions between GD1P and H. hermesi. Finally, we were able to show that nematodes have the ability to synthesize polyunsaturated fatty acids which are not present in their food source.
Abdulkadir S, Tsuchiya M. One-step method for quantitative and qualitative analysis of fatty acids in marine animal samples. J Exp Mar Biol Ecol. 2008;354(1):1–8.
Ackman RG. Fatty acids. In: Ackman RG, editor. Marine biogenic lipids, fats, and oils. Boca Raton: CRC Press; 1989. p. 103–37.
Alqarawi AA, Hashem A, Abd Allah EF, Alshahrani TS, Huqail AA. Effect of salinity on moisture content, pigment system, and lipid composition in Ephedra alata Decne. Acta Biol Hung. 2014;65(1):61–71.
Anderson MJ. PERMDISP: a FORTRAN computer program for permutation analysis of multivariate dispersions (for any two-factor ANOVA design) using permutation tests. Auckland: Department of Statistics, University of Auckland; 2004.
Anderson MJ, Gorley RN, Clarke KR. PERMANOVA+ for PRIMER: guide to software and statistical methods. Plymouth: PRIMER-E Ltd.; 2008.
Bartlett DH. Pressure effects on in vivo microbial processes. BBA Protein Struct Mol Enzymol. 2002;1595(1–2):367–81.
Bolla R. Nematode energy metabolism. In: Zuckermann BM, editor. Nematodes as biological models, vol. 2. New York: Academic Press; 1980. p. 165–92.
Brett M, Müller-Navarra D. The role of highly unsaturated fatty acids in aquatic foodweb processes. Freshw Biol. 1997;38(3):483–99.
Cossins AR, Macdonald AG. The adaptation of biological membranes to temperature and pressure: fish from the deep and cold. J Bioenerg Biomembr. 1989;21(1):115–35.
Cossins AR, Macdonald AG. Homeoviscous adaptation under pressure 3. The fatty-acid composition of liver mitochondrial phospholipids of deep-sea fish. Biochim Biophys Acta. 1986;860(2):325–35.
Cronan JE, Gennis RB, Maloy SR. Cytoplasmic membrane. In: Neidhardt FC, Ingraham JL, Low KB, Magasanik B, Schaechter M, Umbarger HE, editors. Escherichia coli and Salmonella: cellular and molecular biology. Washington, D.C: American Society for Microbiology; 1987. p. 31–55.
De Troch M, Boeckx P, Cnudde C, Van Gansbeke D, Vanreusel A, Vincx M, Caramujo MJ. Bioconversion of fatty acids at the basis of marine food webs: insights from a compound-specific stable isotope analysis. Mar Ecol Prog Ser. 2012;465:53–67.
DeLong EF, Yayanos AA. Adaptation of the membrane lipids of a deep-sea bacterium to changes in hydrostatic pressure. Science. 1985;228(4703):1101–3.
DeLong EF, Yayanos AA. Biochemical function and ecological significance of novel bacterial lipids in deep-sea procaryotes. Appl Environ Microbiol. 1986;51(4):730–7.
Derycke S, Backeljau T, Vlaeminck C, Vierstraete A, Vanfleteren J, Vincx M, Moens T. Spatiotemporal analysis of population genetic structure in Geomonhystera disjuncta (Nematoda, Monhysteridae) reveals high levels of molecular diversity. Mar Biol. 2007;151(5):1799–812.
Drijber RA, Doran JW, Parkhurst AM, Lyon DJ. Changes in soil microbial community structure with tillage under long-term wheat-fallow management. Soil Biol Biochem. 2000;32(10):1419–30.
Duldhardt I, Gaebel J, Chrzanowski L, Nijenhuis I, Hartig C, Schauer F, Heipieper HJ. Adaptation of anaerobically grown Thauera aromatica, Geobacter sulfurreducens and Desulfococcus multivorans to organic solvents on the level of membrane fatty acid composition. Microb Biotechnol. 2010;3(2):201–9.
Dumm INTD, Dealaniz MJT, Brenner RR. Daily variations of the biosynthesis and composition of fatty-acids in rats fed on complete and fat-free diets. Lipids. 1984;19(2):91–5.
Ederington MC, Mcmanus GB, Harvey HR. Trophic transfer of fatty-acids, sterols, and a triterpenoid alcohol between bacteria, a ciliate, and the copepod Acartia tonsa. Limnol Oceanogr. 1995;40(5):860–7.
Elvert M, Boetius A, Knittel K, Jorgensen BB. Characterization of specific membrane fatty acids as chemotaxonomic markers for sulfate-reducing bacteria involved in anaerobic oxidation of methane. Geomicrobiol J. 2003;20(4):403–19.
Fonseca-Madrigal J, Pineda-Delgado D, Martinez-Palacios C, Rodriguez C, Tocher DR. Effect of salinity on the biosynthesis of n-3 long-chain polyunsaturated fatty acids in silverside Chirostoma estor. Fish Physiol Biochem. 2012;38(4):1047–57.
Fonseca G, Derycke S, Moens T. Integrative taxonomy in two free-living nematode species complexes. Biol J Linn Soc. 2008;94(4):737–53.
Fozo EA, Quivey RG. Shifts in the membrane fatty acid profile of Streptococcus mutans enhance survival in acidic environments. Appl Environ Microbiol. 2004;70(2):929–36.
Fulco AJ. Fatty-acid metabolism in bacteria. Prog Lipid Res. 1983;22(2):133–60.
Gerlach SA, Schrage M. Life cycles in marine meiobenthos—experiments at various temperatures with Monhystera disjuncta and Theristus pertenuis (Nematoda). Mar Biol. 1971;9(3):274–80.
Guckert JB, Antworth CP, Nichols PD, White DC. Phospholipid, ester-linked fatty-acid profiles as reproducible assays for changes in prokaryotic community structure of estuarine sediments. FEMS Microbiol Ecol. 1985;31(3):147–58.
Guilini K, Veit-Kohler G, De Troch M, Van Gansbeke D, Vanreusel A. Latitudinal and temporal variability in the community structure and fatty acid composition of deep-sea nematodes in the Southern Ocean. Prog Oceanogr. 2013;110:80–92.
Guschina IA, Harwood JL. Mechanisms of temperature adaptation in poikilotherms. FEBS Lett. 2006;580(23):5477–83.
Hazel JR. Thermal adaptation in biological membranes: is homeoviscous adaptation the explanation? Annu Rev Physiol. 1995;57(1):19–42.
Hazel JR, Williams EE. The Role of alterations in membrane lipid-composition in enabling physiological adaptation of organisms to their physical-environment. Prog Lipid Res. 1990;29(3):167–227.
Heip C, Vincx M, Vranken G. The ecology of marine nematodes. Oceanogr Mar Biol. 1985;23:399–489.
Hutzell PA, Krusberg LR. Fatty acid compositions of Caenorhabditis elegans and C. briggsae. Comp Biochem Physiol Part B Comp Biochem. 1982;73(3):517–20.
Jagdale GB, Gordon R. Effect of temperature on the composition of fatty acids in total lipids and phospholipids of entomopathogenic nematodes. J Therm Biol. 1997;22(4–5):245–51.
Javed S, Bukhari SA, Ashraf MY, Mahmood S, Iftikhar T. Effect of salinity on growth, biochemical parameters and fatty acid composition in safflower (Carthamus tinctorius L.). Pak J Bot. 2014;46(4):1153–8.
Kamimura K, Fuse H, Takimura O, Yamaoka Y. Effects of growth pressure and temperature on fatty acid composition of a barotolerant deep-sea bacterium. Appl Environ Microbiol. 1993;59(3):924–6.
Kaneda T. Iso-fatty and anteiso-fatty acids in bacteria—biosynthesis, function, and taxonomic significance. Microbiol Rev. 1991;55(2):288–302.
Kaye JZ, Baross JA. Synchronous effects of temperature, hydrostatic pressure, and salinity on growth, phospholipid profiles, and protein patterns of four Halomonas species isolated from deep-sea hydrothermal-vent and sea surface environments. Appl Environ Microbiol. 2004;70(10):6220–9.
Kaye JZ, Marquez MC, Ventosa A, Baross JA. Halomonas neptunia sp. nov., Halomonas sulfidaeris sp. nov., Halomonas axialensis sp. nov. and Halomonas hydrothermalis sp. nov.: halophilic bacteria isolated from deep-sea hydrothermal-vent environments. Int J Syst Evol Microbiol. 2004;54(2):499–511.
Koštál V. Cell structural modifications in insects at low temperatures. In: Lee J, Lee Jr RE, editors. Low temperature biology of insects. New York: Cambridge University Press; 2010. p. 116–40.
Kreps EM. Lipids of cell membranes. Leningrad: Nauka; 1981.
Lambshead PJD. Marine nematode biodiversity. In: Chen ZX, Chen SY, Dickson DW, editors. Nematology: advances and perspectives, vol 1: Nematode morphology, physiology and ecology. Wallingford: International CABI; 2004. p. 436–67.
Leduc D, Probert PK. The effect of bacterivorous nematodes on detritus incorporation by macrofaunal detritivores: a study using stable isotope and fatty acid analyses. J Exp Mar Biol Ecol. 2009;371(2):130–9.
Levin LA. Ecology of cold seep sediments: interactions of fauna with flow, chemistry and microbes. Oceanogr Mar Biol Annu Rev. 2005;43:1–46.
Levins R. Evolution in changing environments: some theoretical explorations. Princeton: Princeton University Press; 1968.
MacAvoy SE, Macko SA, Joye SB. Fatty acid carbon isotope signatures in chemosynthetic mussels and tube worms from gulf of Mexico hydrocarbon seep communities. Chem Geol. 2002;185(1–2):1–8.
Moens T, Vincx M. On the cultivation of free-living marine and estuarine nematodes. Helgol Meeresunters. 1998;52(2):115–39.
Mokievsky VO, Filippova KA, Chesunov AV. Nematode fauna associated with detached kelp accumulations in the subtidal zone of the White Sea. Oceanology+. 2005;45(5):689–97.
Mrozik A, Piotrowska-Seget Z, Łabużek S. Changes in whole cell-derived fatty acids induced by naphthalene in bacteria from genus Pseudomonas. Microbiol Res. 2004;159(1):87–95.
Nichols DS, Olley J, Garda H, Brenner RR, McMeekin TA. Effect of temperature and salinity stress on growth and lipid composition of Shewanella gelidimarina. Appl Environ Microbiol. 2000;66(6):2422–9.
Nogi Y, Kato C, Horikoshi K. Taxonomic studies of deep-sea barophilic Shewanella strains and description of Shewanella violacea sp. nov. Arch Microbiol. 1998;170(5):331–8.
Paradis M, Ackman RG. Localization of a marine source of odd chain-length fatty-acids. 1. Amphipod Pontoporeia femorata (Kroyer). Lipids. 1976;11(12):863–70.
Pond DW, Tarling GA, Mayor DJ. Hydrostatic pressure and temperature effects on the membranes of a seasonally migrating marine copepod. PLoS One. 2014;9(10):e111043.
Portnova D, Haflidason H, Todt C. Nematode species distribution patterns at the Nyegga pockmarks. In: 14th international meiofauna conference (FourtIMCo), Ghent, Belgium; 2010. p. 174.
R Core Team. R: a language and environment for statistical computing. Vienna: R Foundation for Statistical Computing; 2013. http://www.R-project.org/.
Reisner K. Physiological roles of fatty acids and endocannabinoids in the nematode Caenorhabditis elegans. Ph.D. dissertation, University of Easten Finland; 2012.
Rothstein M, Gotz P. Biosynthesis of fatty acids in the free-living nematode, Turbatrix aceti. Arch Biochem Biophys. 1968;126(1):131–40.
Ruess L, Haggblom MM, Zapata EJG, Dighton J. Fatty acids of fungi and nematodes—possible biomarkers in the soil food chain? Soil Biol Biochem. 2002;34(6):745–56.
Sauter EJ, Muyakshin SI, Charlou J-L, Schlüter M, Boetius A, Jerosch K, Damm E, Foucher J-P, Klages M. Methane discharge from a deep-sea submarine mud volcano into the upper water column by gas hydrate-coated methane bubbles. Earth Planet Sci Lett. 2006;243(3–4):354–65.
Schlechtriem C, Tocher DR, Dick JR, Becker K. Incorporation and metabolism of fatty acids by desaturation and elongation in the nematode, Panagrellus redivivus. Nematology. 2004;6:783–95.
Sikkema J, Debont JAM, Poolman B. Mechanisms of membrane toxicity of hydrocarbons. Microbiol Rev. 1995;59(2):201–22.
Sinensky M. Homeoviscous adaptation—a homeostatic process that regulates the viscosity of membrane lipids in Escherichia coli. Proc Natl Acad Sci USA. 1974;71(2):522–5.
Tahseen Q. Nematodes in aquatic environments: adaptations and survival strategies. Biodivers J. 2012;3(1):13–40.
Tanaka T, Ikita K, Ashida T, Motoyama Y, Yamaguchi Y, Satouchi K. Effects of growth temperature on the fatty acid composition of the free-living nematode Caenorhabditis elegans. Lipids. 1996;31(11):1173–8.
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. Helgol Mar Res. 2015;69:57–85.
Tocher DR. Metabolism and functions of lipids and fatty acids in teleost fish. Rev Fish Sci. 2003;11(2):107–84.
Turner JP, Rooker JR. Effect of dietary fatty acids on the body tissues of larval and juvenile cobia and their prey. J Exp Mar Biol Ecol. 2005;322(1):13–27.
Valentine RC, Valentine DL. Omega-3 fatty acids in cellular membranes: a unified concept. Prog Lipid Res. 2004;43(5):383–402.
Van Campenhout J, Derycke S, Moens T, Vanreusel A. Differences in life-histories refute ecological equivalence of cryptic species and provide clues to the origin of bathyal Halomonhystera (Nematoda). PLoS One. 2014;9(11):e111889.
Van Campenhout J, Derycke S, Tchesunov A, Portnova D, Vanreusel A. The Halomonhystera disjuncta population is homogeneous across the Håkon Mosby mud volcano (Barents Sea) but is genetically differentiated from its shallow-water relatives. J Zool Syst Evol Res. 2013;52(3):203–16.
Van Campenhout J, Vanreusel A, Van Belleghem S, Derycke S. Transcription, signaling receptor activity, oxidative phosphorylation, and fatty acid metabolism mediate the presence of closely related species in distinct intertidal and cold-seep habitats. Genome Biol Evol. 2016;8(1):51–69.
Van Gaever S, Moodley L, de Beer D, Vanreusel A. Meiobenthos at the Arctic Hakon Mosby Mud Volcano, with a parental-caring nematode thriving in sulphide-rich sediments. Mar Ecol Prog Ser. 2006;321:143–55.
Van Gaever S, Moodley L, Pasotti F, Houtekamer M, Middelburg JJ, Danovaro R, Vanreusel A. Trophic specialisation of metazoan meiofauna at the Håkon Mosby Mud Volcano: fatty acid biomarker isotope evidence. Mar Biol. 2009;156(6):1289–96.
Vettier A, Labbe C, Amerand A, Da Costa G, Le Rumeur E, Moisan C, Sebert P. Hydrostatic pressure effects on eel mitochondrial functioning and membrane fluidity. Undersea Hyperb Med. 2006;33(3):149–56.
Vranken G, Herman PMJ, Heip C. 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. 1988;77(3):296–301.
Vranken G, Tire C, Heip C. Effect of temperature and food on hexavalent chromium toxicity to the marine nematode Monhystera disjuncta. Mar Environ Res. 1989;27(2):127–36.
Vranken G, Tire C, Heip C. The toxicity of paired metal mixtures to the nematode Monhystera disjuncta (Bastian, 1865). Mar Environ Res. 1988;26(3):161–79.
Vranken G, Vanderhaeghen R, Heip C. Toxicity of cadmium to free-living marine and brackish water nematodes (Monhystera microphthalma, Monhystera disjuncta, Pellioditis marina). Dis Aquat Organ. 1985;1(1):49–58.
Weber G, Drickamer HG. The effect of high-pressure upon proteins and other biomolecules. Q Rev Biophys. 1983;16(1):89–112.
Williams CM, Burdge G. Long-chain n-3 PUFA: plant v. marine sources. Proc Nutr Soc. 2006;65(1):42–50.
Yano Y, Nakayama A, Ishihara K, Saito H. Adaptive changes in membrane lipids of barophilic bacteria in response to changes in growth pressure. Appl Environ Microbiol. 1998;64(2):479–85.
Yano Y, Nakayama A, Yoshida K. Distribution of polyunsaturated fatty acids in bacteria present in intestines of deep-sea fish and shallow-sea poikilothermic animals. Appl Environ Microbiol. 1997;63(7):2572–7.
Yayanos AA, DeLong EF. Deep-sea bacterial fitness to environmental temperatures and pressures. In: Jannasch HW, Marquis RE, Zimmerman AM, editors. Current perspectives in high pressure biology. New York: Academic Press; 1987. p. 17–32.
Yazawa K, Araki K, Okazaki N, Watanabe K, Ishikawa C, Inoue A, Numao N, Kondo K. Production of eicosapentaenoic acid by marine bacteria. J Biochem. 1988;103(1):5–7.
JVC and AV conceived the study; JVC performed the experiment, sampling and all the analyses and data handling. All authors contributed to data interpretation. JVC wrote the manuscript with significant input of AV. All authors read and approved the final manuscript.
We acknowledge chief scientist A. Boetius and project LOOME coordinator D. De Beer for the sampling opportunity during the RV Maria S. Merian campaign MSM 16/2. Special thanks to the captain and the crew of the RV ‘Maria S. Merian’, IFREMER (France). We thank Dirk Van Gansbeke for the fatty acid extractions and identifications. The data analysis of this research has benefitted from a statistical consult with Ghent University FIRE (Fostering Innovative Research based on Evidence). Additionally, we thank Bob Clarke (Plymouth marine laboratory) for statistical guidance in the use of PERMANOVA. This project was financed by the EU project HERMIONE (No. 226354), the Flemish Fund for Scientific Research (F.W.O) deep-sea project (No. 3G083512) and by the DFG for the Maria S. Merian cruise MSM 16/2.
The authors declare that they have no competing interests.
An erratum to this article is available at http://dx.doi.org/10.1186/s10152-016-0473-8.
About this article
Cite this article
Van Campenhout, J., Vanreusel, A. Closely related intertidal and deep-sea Halomonhystera species have distinct fatty acid compositions. Helgol Mar Res 70, 8 (2016). https://doi.org/10.1186/s10152-016-0467-6
- Deep sea
- Fatty acids
- Hakon mosby mud volcano