- Original Article
Trophic niche partitioning of littoral fish species from the rocky intertidal of Helgoland, Germany
Helgoland Marine Research volume 69, pages 385–399 (2015)
During a 3-year field study, interspecific and interannual differences in the trophic ecology of littoral fish species were investigated in the rocky intertidal of Helgoland island (North Sea). We investigated trophic niche partitioning of common coexisting littoral fish species based on a multi-tracer approach using stable isotope and fatty acids in order to show differences and similarities in resource use and feeding modes. The results of the dual-tracer approach showed clear trophic niche partitioning of the five target fish species, the goldsinny wrasse Ctenolabrus rupestris, the sand goby Pomatoschistus minutus, the painted goby Pomatoschistus pictus, the short-spined sea scorpion Myoxocephalus scorpius and the long-spined sea scorpion Taurulus bubalis. Both stable isotopes and fatty acids showed distinct differences in the trophic ecology of the studied fish species. However, the combined use of the two techniques added an additional resolution on the interannual scale. The sand goby P. minutus showed the largest trophic plasticity with a pronounced variability between years. The present data analysis provides valuable information on trophic niche partitioning of fish species in the littoral zones of Helgoland and on complex benthic food webs in general.
Understanding the trophic interactions in complex and diverse marine ecosystems is a significant challenge. In order to obtain a profound knowledge on complex ecosystems, it is important to address trophodynamic interactions and to account for trophic niche partitioning of coexisting species in individual systems.
While stable carbon isotopes (13C/12C) provide a measure of the carbon sources used by consumers, nitrogen isotopes (15N/14N) serve as an indicator of the relative trophic position of the organism within the food web (DeNiro and Epstein 1978; Minagawa and Wada 1984; Peterson and Fry 1987). In contrast to e.g. quantitative gut content analyses which provide only ‘snapshots’ of ingested prey items used by a specific consumer, stable isotopes (SI) provide long-term, time-integrated information on the food sources, thus allowing estimates on the trophic ecology of individual species (Polacik et al. 2014).
The use of δ 13C and δ 15N in food web studies is considered as a reliable tool for estimates on community structure and niche width (Layman et al. 2007; Jackson et al. 2011), as well as for intra- and interspecific variability and feeding modes (Bearhop et al. 2004; Aberle et al. 2005; Newsome et al. 2007; Polacik et al. 2014) in complex ecosystems. To avoid misinterpretations of SI data, it is beneficial to combine it with additional approaches, e.g. fatty acid (FA) analysis—providing supplementary dietary information such as a clear indication of the food sources used (Herman et al. 2005; Perga et al. 2006; Inger and Bearhop 2008). Herman et al. (2005) showed that, based on SI data, killer whales Orcinus orca predominantly prey marine mammals, while FA profiles pointed at a more fish-based diet. Here, the combination of different tracer methods widened the range of dietary information in this specific field study. Consequently, the combined use of SI and FA analyses is a promising approach to characterize complex food web interactions and to define specific trophic niches in a given ecosystem (Peterson and Fry 1987; Fry et al. 1999).
In the present study, we analysed the trophic ecology of the common littoral marine fish species goldsinny wrasse, sand goby, painted goby, long-spined and short-spined scorpion due to their wide distribution ranges in a variety of habitats along the European coastline (i.e. Wheeler and Du Heaume 1969; Hureau et al. 1986). The study focussed on the rocky intertidal fish community of Helgoland, a small and isolated island in the German Bight, North Sea. The aim was to disentangle food web interactions of our target fish species and to define their trophic niches and position of coexisting species over a period of three consecutive years. Due to the ecological relevance of Gobiidae in rocky and soft-bottom habitats (Wennhage and Pihl 2002), the specific role of the sand goby P. minutus was considered in detail. The crucial role of P. minutus in marine and estuarine ecosystems has been stressed in the literature since this species acts both as predator as well as prey (e.g. Leitao et al. 2006 and references therein), serving as a preferred prey item of coexisting fish species e.g. T. bubalis and M. scorpius (Ehrenbaum 1936; Wheeler and Du Heaume 1969). Therefore, we considered this goby species more specifically and tried to define its trophic role among other littoral fish species within the same littoral community.
We hypothesized that we would find (1) interspecific variations in trophic niches of coexisting littoral fish species and (2) annual patterns in trophic niche partitioning (e.g. changes in prey availability).
Materials and methods
Study area and sample collection
A 3-year sampling campaign of littoral fish species was conducted (2009–2011) around the island of Helgoland (54°11′N, 7°53′O), situated in the North Sea 50 km off the German west coast (Fig. 1) to investigate the trophic ecology and trophic niche partitioning of littoral marine fish species.
Target fish species were: goldsinny wrasse Ctenolabrus rupestris (L.), sand goby Pomatoschistus minutus (Pallas), painted goby Pomatoschistus pictus (Malm), short-spined sea scorpion Myoxocephalus scorpius (L.) and long-spined sea scorpion Taurulus bubalis (Euphrasen) (Table 1). Specimen were collected from late spring to autumn (June to November) over a period of three consecutive years (2009–2011) using beach seining (20 m × 1.20 m, 6-mm codend), fyke nets and hand netting with suitable mesh sizes as fishing techniques. Furthermore, bivalve (Mytilus edulis) and gastropod (Littorina littorea) samples were collected each year at monthly intervals from June to November to provide a seasonal baseline for SI values of fish species. In general, 10–20 individuals of each species (incl. M. edulis and L. littorea) were collected based on the availability at each sampling season. Subsamples were frozen immediately after processing and stored at −80 °C until analyses.
Fatty acid (FA) analysis
The extraction of FAs was carried out on each individual fish applying a modification of the methods following the analytical procedure by Malzahn et al. (2007). For all fish species, we used dorsal muscle samples dissected directly from posterior to the head. Fatty acids were analysed by gas chromatography (Varian CP 8400 gas chromatograph equipped with a DB-225 column [J&W Scientific, 30-m length, 0.25-mm inner-diameter (ID), 0.25-μm film]. Quantification of the fatty acid methyl esters (FAME, i.e. a type of fatty acid ester that are derived by transesterification of fats with methanol) was done by using GALAXIE Workstation software (Galaxie 1.9 SP2b, Varian, Inc., USA), and each chromatogram was manually assessed for correct peak identification (and also reintegrated when it was necessary). A known amount of C23:0 as internal standard was added additionally before starting the FAME extraction.
All FA data were reported in μg FA/µg C referring to the carbon (C) amount of defatted muscle tissue from the SI measurements (Table 2). The trophic tracers of dietary fatty acids used in this study are DHA/EPA, 18:1n9/18:1n7 and PUFA/SFA (El-Sabaawi et al. 2009).
Samples showing a total FA content higher than 0.15 µg FA/µg C were not included into the analysis since analytical errors were expected.
Stable isotope (SI) analysis
SI analysis was performed using solely muscle tissue of fish due to discrepancies in synthesis and replacement processes in different animal tissues (Pinnegar and Polunin 1999; Bearhop et al. 2002). Each tissue sample was freeze-dried and defatted with a dichloromethane–methanol (2:1 vol/vol) mixture, treated in an ultrasound bath for 30 min and dried again before ±1 mg of sample was weighed into a tin capsule (HEKAtech GmbH) for further SI analysis (δ 15N and δ 13C). While lipid extraction can alter the δ 13C of consumers (Kiljunen et al. 2006), δ 15N usually remains unaffected by the defatting procedure (Ingram et al. 2007).
Freeze-dried bivalve and gastropod tissue was scraped from the shells prior to grinding. An adequate amount of ground tissue was defatted for 24 h before preparing subsamples for SI analysis.
13C/12C and 15N/14N ratios were determined using continuous flow isotope ratio mass spectrometry (PDZ Europa ANCA-GSL elemental analyser interfaced to a PDZ Europa 20–20 isotope ratio mass spectrometer, Sercon Ltd., Cheshire, UK) with an analytical precision of 0.2 ‰ for δ 13C and 0.3 ‰ for δ 15N carried out by UC Davis Stable Isotope Facility, University of California, USA. The ratios (R) of 13C/12C and 15N/14N are expressed in conventional delta notation (δ) as parts per thousand relative to the international carbon standard Pee Dee Belemnite and atmospheric N2 for nitrogen:
where X: equal to 13C or 15N and R: corresponding ratio 13C/12C or 15N/14N.
Calculations on the trophic position (TP) of our fish species relied in general on the baseline SI of our benthic grazer L. littorea. Only P. minutus from 2010 showed a stronger dependency on 13C derived from pelagic production, and therefore, M. edulis served as baseline organism in this specific case (Table 3).
In order to estimate the trophic position (TP) of each individual, the following equation was used (Post 2002):
assuming that 2.9 ‰ describes the averaged 15N trophic enrichment factor after McCutchan et al. (2003). 2 represents the trophic position of the baseline organisms and enables a legitimate tool for comparability with other trophic level calculations from literature.
Length–weight relationship of P. minutus
One of the most frequently used morphometric indices is Fulton’s condition factor K that can be used to describe interannual variations in P. minutus condition (Nash et al. 2006). Due to its conflicting results in preceding studies (i.e. Trippel and Hubert 1989; Davidson and Marshall 2010; Mozsar et al. 2015), we used the relative condition factor K n , with:
where W: weight (g), SL: standard length (mm), a and n are the exponential form of the intercept and slope, respectively, of the logarithmic (log) length–weight equation.
To assess the trophic ecology and interactions between fish species, we analysed the degree of similarity between FAs and SI compositions.
Patterns in FA composition of individual fish were analysed using principal component analysis (PCA) with PRIMER 6. We selected 20 FAs that were the most abundant across all species, including 16:0, 16:1n7, 18:0, 18:1n7, 18:1n9, 18:2n6, 18:3n6, 18:3n4, 20:0, 20:1n9, 20:4n6, 20:4n3, 20:5n3 (EPA), 22:0, 22:1n9, 22:2n6, 22:5n3 (DPA), 22:6n3 (DHA), 24:0 and 24:1n9. Prior to analysis, each FA content was ln-transformed in order to reduce the relationship between the mean and the variance and to linearize the relationship between variables. A standardized PCA based on the correlation matrix was performed on the transformed values, and the sample scores were shown on the first factorial plane (PC1 vs. PC2). In order to further test for interspecific and interannual FA patterns, a redundancy analysis (RDA; Legendre and Legendre 1998) was performed. The significance of the relationship between FA and the two factors as well as their interactive effects was investigated by Monte Carlo permutation tests. First, a general test was applied to test the null hypothesis of independence between datasets, the response and the explanatory variables. Second, a forward selection process was performed in order to identify the most important variables. Here, the adjusted R 2-function is maximized, similar to the selection process in stepwise multiple regression. RDA analyses were performed with the software R 3.0.2 (R Core Team 2013) for statistical computing and the related package vegan (Oksanen et al. 2013).
To test for the suitability of carnivory markers as well as the capacity of FA indices as indicators for trophic position, we performed linear regressions between the baseline-corrected δ 15N of each fish with the three FA markers for carnivory (PUFA/SFA, DHA/EPA and 18:1n9/18:1n7).
To account for interspecific/interannual differences and potential interactive effects in SI, analysis of variance (ANOVA) were performed with Statistica 9.1, followed by Tukey’s honest significant difference test (HSD) for unequal sample size (N). Interannual changes in the condition of P. minutus were investigated with analysis of covariance (ANCOVA) by comparing the logarithmic length–weight relationships between sampling years. Prior to these analyses, data were tested for the analytical assumptions of variance homogeneity and normality.
Fatty acids (FA)
Applying a principal component analysis (PCA), the most important 20 FAs of all samples resulted in 69.1 % explained variability on the first two axes (PC1 = 47 %, PC2 = 22.1 %), and only little further information was added by subsequent PCs (e.g. PC3 = 8.6 %, Fig. 2). Most of the variability in the FA composition of the five littoral fish species was caused by FAs listed in Table 4. Among these FAs, the trophic markers 18:1n7, 18:1n9 and 22:6n3 (DHA) displayed high loadings on the first two axes contributing to the clustering effects in Fig. 2a, b. We found distinct annual FA compositions with a slight overlap between the samples collected in 2010 and 2011 (Fig. 2a). The interannual differences were mainly described by the second axis (PC2) with negative scores for 2010 and 2011 and positive scores for 2009. High positive correlations with the second PC were found e.g. for the FAs 18:3n6 and 18:3n4 and a negative correlation with 16:0 (see Table 4), meaning that the first two FAs were found in above-average concentrations in 2009 and in below-average concentrations in 2010 and 2011. In contrast to this, interspecific differences were less pronounced and species-specific samples overlapped to a high degree on the first factorial plane (Fig. 2b). However, the highest variability between species was found on the first PC, to which all FAs displayed a negative correlation (Table 4). This means that samples with positive PC scores had below-average concentrations for most of these FAs, whereas samples with negative PC scores showed above-average FA concentrations. Highest scores were found for P. minutus, lowest for P. pictus. The other three species (C. rupestris, M. scorpius and T. bubalis) were not clearly separated on the first two axes.
To elucidate the differences among individual FA patterns, redundancy analysis (RDA) was used in order to test for species- and annual-specific effects. Both factors explained a significant amount of the variability in the FA contents (p < 0.005) with a total of 58.3 % variability explained on all axes and 42.3 % on the first factorial plane (Table 5). The interaction between both factors, species and sampling year, was significant on the p < 0.01 level. When accounting for factors separately as well as sequentially, all effects turned out to be significant with a slightly more variation explained by sampling year than by species (Table 5). In order to assess the importance of each factor level, we additionally repeated the analysis with binary-coded variables for each species and sampling year. In the selection process, the year 2010 and the species P. pictus were revealed to be the variables explaining most of the variation in the FA patterns. This suggests that the gobies caught in 2010 showed strong deviations from the general FA pattern in comparison with all other species.
Stable isotopes (SI)
Compared to FA profiles, SI showed less pronounced interannual differences in fish species (Fig. 3a; Table 2). Significant interannual differences were observed in δ 13C values with a smaller difference observed between 2009 and 2010 (ANOVA, p < 0.05) and a bigger difference between 2010 and 2011 (ANOVA, p < 0.01). In contrast to this, annual mean δ 15N values showed no clear difference between years, but δ 15N were extremely variable in 2010, ranging from 14.5 ‰ for P. minutus to 19.9 ‰ for M. scorpius. This variability was mainly caused by P. minutus that moreover showed significant differences in 2010 compared to P. minutus sampled in 2009 and 2011 with 17.8 ‰ and 17.3 ‰, respectively, (2009–2010, 2010–2011 both with an ANOVA p value of <0.01, Fig. 3). Hence, some species also experienced differences in SI throughout the sampling years (i.e. T. bubalis showed a different SI in 2010), but P. minutus displayed the strongest variability in calculated SI ± standard deviation (δ 15N: 16 ± 2 ‰ and δ 13C: −16 ± 2 ‰).
With regard to δ 15N, the two gobiid species P. minutus and P. pictus were clearly separated from the remaining species (ANOVA, p < 0.01), with distinct δ 15N values of ~16 ‰ for P. minutus and ~18 ‰ for P. pictus. Differences in δ 13C were less pronounced and only the δ 13C of C. rupestris and M. scorpius differed significantly from each other (ANOVA, p < 0.01). A notable distinction in the 13C-enrichment was found between M. scorpius and the remaining species C. rupestris, P. pictus and T. bubalis. (−15.6 ‰ for M. scorpius vs. ~−17 ‰, respectively). The δ 13C values of P. minutus varied strongly between years. The 2010 samples showed significantly depleted δ 13C values compared to samples from 2009 and 2011 (−18, −15 and −15 ‰, respectively, ANOVA, p < 0.05). The mean trophic position (TP) among the two gobiid species did not differ significantly (TP of 3.0–3.5), but the gobies differed from all other species (TP 4) resulting in a lower trophic position than C. rupestris, M. scorpius and T. bubalis (Table 2). However, TP of P. minutus did not differ significantly among sampling years when baseline corrections were applied for specimen sampled in 2010.
We used a PCA to combine the results of SI and FA analyses (Fig. 4; Table 6). Here, 63.4 % of the total variation was explained by the first two principal components (PC1: 40.5 %, PC2: 22.9 %). Moreover, PC3 explained 17.2 % of the variance and δ 15N showed the highest loading therein. For clarity, only the species-specific annual mean PC scores were visualized in Fig. 4. PC2 mainly represents the FA trophic marker PUFA/SFA ratio and accounts for the differences that were found between 2009 and the other two sampling years (2010 and 2011). These were, with the exception of samples from P. minutus in 2010, similar according to their PC2 score. Species scores are separated on the first axis, with FAs (i.e. DHA/EPA and 18:1n9/18:1n7) loading positively, SIs loading negatively on PC1. C. rupestris, P. pictus and the 2010 population of P. minutus have above-average values for FA biomarkers and below-average values for SIs. In contrast, samples of M. scorpius, T. bubalis and P. minutus (2009 + 2010) are characterized by negative PC1 values. Thus, lower levels of 18:1n9/18:1n7 and DHA/EPA and higher levels of SI biomarkers were observed. Furthermore, species-specific scores of T. bubalis resembled those of M. scorpius, over all sampling years. P. minutus is the only species where interannual differences were pronounced not only on the second, but also on the first PC forming a separate cluster of the 2010 population in contrast to 2009 and 2011. The latter are presented by an opposing trend caused by a significant lower DHA/EPA ratios and remarkably higher δ 13C as well as δ 15N values (Fig. 4). None of the three FA markers—suggested being representative for carnivory (El-Sabaawi et al. 2009)—correlated significantly with the baseline-corrected δ 15N of fish. However, we still decided to use these markers in our analysis in order to separate fish species from each other.
Length–weight relationship of P. minutus
Due to its distinct role, we compared variations in P. minutus condition by logarithmic length-weight relationships among sampling years (Fig. 5). No significant interactions between log length and sampling year were found, meaning that the slope of the regression line remained the same over the entire investigation period. However, when comparing the intercepts and calculating marginal means, the P. minutus population of 2010 showed considerably lower weights at the same length compared to the populations from 2009 and 2011 (p < 0.001, ANCOVA), indicating overall a lower condition.
The present study aimed at analysing trophic variation and trophic niche partitioning among littoral fish species within the rocky intertidal community off Helgoland (German Bight).
The fauna around Helgoland is characterized by a high diversity and complex trophic relations (Janke 1990). The habitat structure is of considerable relevance when the ranges of trophic positions within a given food web are addressed (Bearhop et al. 2004; Layman et al. 2007). In general, trophic niche partitioning is considered as a major strategy for fish assemblages which might be of even greater relevance than temporal and spatial separation (Ross 1986). In this study, we used FA and SI analyses since these tracers allow an integration of diet preferences over a longer period of time prior to sampling. SIs and FAs analyses allowed us to characterize the trophic ecology of the five fish species and provided valuable information on the trophic positions as well as interannual and interspecific variations in our local fish community.
Fatty acid composition of target fish species
Fatty acid (FA) analyses are considered as a reliable tool to trace food sources within a given food web by identifying characteristic, taxa-specific FA pattern at the base of the food web, e.g. aquatic primary producers which are in turn transferred to higher trophic levels (Dalsgaard et al. 2003). However, this approach has recently been criticized pointing at specific uncertainties in using FA as trophic markers when studying the trophic ecology of consumers especially in benthic communities (Kelly and Scheibling 2012; McLeod et al. 2013). Despite potential flaws, FAs are still considered to provide useful dietary information especially with regard to higher-order consumers, e.g. vertebrates, mammals (Kelly and Scheibling 2012).
FA profiles of the examined fish species provided insights into the trophic ecology within the Helgoland rocky intertidal food web, suggesting that FA patterns are driven by differences in dietary preferences, associated with interannual variability. The unsaturated FAs 18:1n7, 18:1n9 and 22:6n3 (DHA) were identified as the most important drivers characterizing trophic relations. These FAs are considered as essential components in the diets of fish (Sargent et al. 1999b) acting as reliable indicators, i.e. functioning as dietary markers by providing valuable information on trophic niche partitioning of fish species. The conservative character of these FAs leads to a valid detection of the transfer from ingested dietary sources to higher trophic levels (Bell et al. 1986; Dalsgaard and St John 2004). The occurrence of FAs with other functions (i.e. metabolic processes and/or storage function in organisms) like 16:0 further provide information on species-specific markers that are of considerable relevance when describing complex food webs. Overall, the present results support the assumption that diet variation results in broad FA spectra of consumers (Budge et al. 2002).
In this study, interannual differences in FA signatures were slightly more pronounced than differences between species, but both factors were highly significant. In addition, we found some evidence for intraspecific differences in FA signatures, but interspecific variations played a more pronounced role. This is in line with observations from other fish assemblages at e.g. the Scotian Shelf, Georges Bank, Gulf of Lawrence or the North-western Hawaiian Islands where strong interspecific variations in FA were reported (Budge et al. 2002; Piche et al. 2010). Overall, the factors sampling year and species explained most of the FA variation of our five fish species with slightly more variation explained by interannual changes in FA profiles.
The observed interannual variations might be explained by changes in the foraging strategy of our fish species. We sampled all specimens with regard to availability, and therefore, spatial and temporal differences between sampling years and seasonal changes in prey availability and/or consumers’ preferences have to be taken into consideration (e.g. El-Sabaawi et al. 2009). Further, intrapopulation variations in terms of resource use and habitat need to be addressed when evaluating food web relations (Bolnick et al. 2003). The observed intraspecific and interannual variability in this study can most likely be explained by a strong flexibility in the foraging behaviour, since changes in feeding modes are a common pattern in fish and diet switches resulting from changes in prey availability and quality commonly occur (Hart 1993).
Stable isotope composition of target fish species
The δ 15N and δ 13C signatures of our fish species showed distinct niche separation patterns. While δ 15N reflects the ‘stepwise enrichment with trophic transfers` (Fry 2006; Inger and Bearhop 2008), δ 13C is an indicator for the diet preferences of consumers (Peterson and Fry 1987; Post 2002). Strong overlaps in δ 13C indicate a similar pool of carbon sources used, while species with a similar δ 15N are considered to occupy similar trophic positions. In our study, we found a close relationship between the two gobiid species P. pictus and P. minutus resulting in similar trophic positions (TP) and diet preferences (δ 13C) with a tendency towards a more generalistic feeding mode of P. minutus. Leitao et al. (2006) characterized P. minutus as an opportunistic carnivorous species, with a feeding mode which is confirmed by the wide δ 13C ranges found in our field study. This is in line with studies reporting on similar feeding preferences of both species with only slightly different dietary sources (Wheeler and Du Heaume 1969; Hureau et al. 1986) and minor differences in habitat preference in terms of microhabitat utilization: semi-hyperbenthic (P. pictus) versus epibenthic (P. minutus) guilds identified by Wilkins and Myers (1992). We are aware that the different sampling sites in this study might have been influenced by differences in microhabitat utilization. Due to the close vicinity of the sampling sites, we are convinced that our study still provides representative data of the entire rocky intertidal fish community. Moreover, differences in mouth gape opening between similar-sized specimens of both species could furthermore explain the observed moderate differences in prey preference/consumption between P. pictus and P. minutus.
The wide range for δ 13C in M. scorpius implies a generalistic feeding mode which is in good agreement with other field studies (Norderhaug et al. 2005), while C. rupestris, commonly considered as a generalist (Fjosne and Gjosaeter 1996), showed a more specialized feeding mode in our study. The closely related Cottidae T. bubalis and M. scorpius are not only characterized by a strong morphological similarity, but also by a specific overlap in their preferred diets with a varying degree of similarity across rocky- and soft-bottom habitats (Wennhage and Pihl 2002). In our field study, distinct prey-derived carbon sources were observed assuming different preferences in their prevalent prey taxa. With regard to carbon isotopes, planktonic consumers in marine coastal environments such as M. edulis are in general considered to show a lower 13C/12C ratio (i.e. reflected by more negative δ 13C values) relative to benthic consumers (France 1995) such as L. littorea.
The five fish species analysed in this study were 13C enriched relative to the mussels assuming the use of carbon sources derived from benthic production (with the exception of P. minutus 2010). This observation is in line with trophic information given by, e.g. Wennhage and Pihl (2002) for benthic fish communities in the rocky intertidal along the Swedish West coast as well as Norderhaug et al. (2005) for a Norwegian kelp forest habitat. However, M. scorpius showed a slightly stronger enrichment in 13C relative to C. rupestris, P. pictus and T. bubalis pointing at strictly benthic-derived carbon sources, whereas P. minutus populations from 2010 as appeared to be the only component which was influenced by pelagic-derived carbon sources.
Moreover, the SI of C. rupestris, M. scorpius, T. bubalis point at a higher trophic position when compared to the gobiid species. Generally, the TP estimates of our fish species are in a good agreement with the predicted ones published in the literature (see Froese and Pauly 2000). In this context, especially the role of P. minutus needs further consideration due to its role as a key component of littoral fish assemblages and its relevance in marine and estuarine food webs. It has been stressed that P. minutus acts both as an intermediate predator and a preferred prey organism for other fish species (Ehrenbaum 1936; Wheeler and Du Heaume 1969; Leitao et al. 2006). Throughout our study, P. minutus and P. pictus showed lower TPs compared with C. rupestris, M. scorpius and T. bubalis. The lower TP in our gobiids might support previous observations pointing at a strong predation on gobiids by C. rupestris, M. scorpius and T. bubalis as reported by previous studies (Hilldén 1978; Wennhage and Pihl 2002; Norderhaug et al. 2005). With regard to the trophic ecology of P. minutus, Dolbeth et al. (2008) reported a more specialized feeding strategy which indicates that the diet ranges in P. minutus can also be restricted depending on the life history variations and habitat-specific differences of a given population. The above-mentioned conclusion could explain the wide range of δ 13C and δ 15N values measured from P. minutus samples between the 3 years observed in this study. Furthermore, the observed distinction among years could point to a more plankton-derived C influence in 2010 (i.e. caused by a massive decay of jellyfish blooms, pers. observation) which could have affected the local food web.
Various publications have suggested the combined use of different biomarkers in order to gain more reliable results on trophodynamic interactions in food web studies (Alfaro et al. 2006; Fry 2006; El-Sabaawi et al. 2009; Kelly and Scheibling 2012). In this study, we used a tracer approach where we combined SI and FA biomarkers. On the one hand, we reduced our dataset of all measured FA up to a specific group of relevant diet-specific FA, the so-called fatty acid trophic markers (FATM after El-Sabaawi et al. (2009), assuming that these dietary components act as reliable indicators for a trophic niche partitioning of our fish species. We are aware that our FATM selection might have caused a simplification of the actual trophic interaction patterns. On the other hand, differences in FA might be obscured either by FATM composition of prey and/or by consumer’s metabolic activities, and therefore, food sources and their differences in quality could modify the observed trophic interactions within the investigated community. Furthermore, we added the SI markers δ 15N and δ 13C to our PCA as response variables. It is known that not only fatty acid but also isotopic variation is not simply related to dietary differences among individuals in a consumer population, but likely to be a function of a number of processes such as nutrient availability, starvation or nutrient-specific assimilation efficiency (Aberle and Malzahn 2007; Inger and Bearhop 2008). Given that SI analysis has limitations too, the use of a combined tracer approach minimizes susceptibility that emerges from single tracer techniques (Kelly and Scheibling 2012). Thus, the proximity in PCA of T. bubalis and M. scorpius uncovered a distinct mutuality in feeding guilds (Wennhage and Pihl 2002) showing similar patterns in feeding strategies and life histories of these two, closely related fish species (e.g. Ehrenbaum 1936; Wheeler and Du Heaume 1969).
In general, the DHA/EPA ratio is used as an index for carnivory (Dalsgaard et al. 2003; Piche et al. 2010). With regard to C. rupestris, a preference for crustaceans and mussels described by Hilldén (1978) and Wennhage and Pihl (2002) is supported by the high values in both FA 18:1n9/18:1n7 and DHA/EPA, indicating carnivorous feeding modes. In contrast, M. scorpius and T. bubalis commonly considered as carnivores in coastal food webs (Wheeler and Du Heaume 1969), showed comparatively low DHA/EPA values in our study. Using SI as an additional trophic marker allowed us to confirm carnivory in M. scorpius and T. bubalis by displaying the highest δ 15N in the present study.
Although FAs were reduced to only three dietary biomarkers, we were still able to confirm interannual variability patterns which were mainly driven by consistently higher PUFA/SFA ratios in 2009. Moreover, the interannual variability of P. minutus, characterized by a high DHA/EPA and low δ 15N as well as δ 13C in 2010, was even more pronounced than in previous analyses. The combination of FA and SI markers allowed us to describe clear trophic patterns: one cluster of P. minutus (except 2010), M. scorpius and T. bubalis feeding on prey relying exclusively on benthic-derived carbon sources and another cluster of C. rupestris and P. pictus, which are characterized by the consumption of prey items using carbon sources derived from planktonic production. Furthermore, SI biomarkers rather than FA biomarkers enabled a clear segregation of P. minutus, M. scorpius and T. bubalis from the remaining species thus allowing a detailed view on trophic relations within our fish community. While FA signatures of consumers are known to become harder to differentiate with increasing trophic levels due to a wide variety in their dietary sources (Dalsgaard et al. 2003), the additional use of SI filled this gap and provided valuable insights into our local food web. It has to be noted that FA measures like DHA/EPA, PUFA/SFA and 18:1n9/18:1n7, which were suggested being representative for carnivory in the literature, did not correlate with the trophic position of the fish in our study. Initially, these markers were reported as a measure for trophic relations among copepods (El-Sabaawi et al. 2009) as well as other crustaceans, e.g. amphipods and krill (Stevens et al. 2004; Nyssen et al. 2005; Schmidt et al. 2006; Hagen et al. 2007). It can be speculated that taxonomic differences and differences in fatty acid metabolism between benthic versus pelagic consumers as well as vertebrate versus invertebrates might be responsible for the contrasting results.
Identifying substructures in apparently homogeneous populations is another advantage of multi-tracer approaches. In particular, the trophodynamic patterns in P. minutus were significantly different in 2010—suggesting that FA and SI reflect differences in feeding modes and/or prey availability/quality. In general, a high interannual heterogeneity can be expected in gobiids, as these small fish are well known for their “great ecological diversity” (Bouchereau and Guelorget 1998). Furthermore, habitat selection in gobiids implies that their distribution ranging from estuaries, coastal areas to lagoons (Wheeler and Du Heaume 1969; Jaquet and Raffaelli 1989; Leitao et al. 2006) is wider than the one from other more stationary fishes. In this context, the low δ 13C and δ 15N of P. minutus in 2010 could result from differences in feeding habits or shifts in prey type and abundance. The pronounced interannual variability is supported by FA biomarkers which showed a consistently high PUFA/SFA ratio in 2009 and a distinct DHA/EPA ratio in 2010. Since we did not find any significant shifts in the baseline, such interannual variations in trophic markers can most likely be related to dietary shifts or changes in metabolic processes/activities (Sargent et al. 1999a).
In order to elucidate the differences among years for P. minutus, the body condition in terms of the interannual length–weight relationship as morphometric condition index for estimating the magnitude of energy reserve (Davidson and Marshall 2010) was further examined. The 2010 population turned out to have a significantly lower condition compared with the other two sampling years, and hence ranked the 2010 population to a lower nutritional condition. It can be assumed that the relative contribution of benthic and pelagic carbon sources caused the observed alterations in length–weight relationship in P. minutus. With regard to the gobies, a carbon source derived from pelagic production seems to be less beneficial for the condition of P. minutus.
This aspect would need further consideration in the future when addressing the pros and cons of differently derived carbon sources on the growth and condition of littoral fish species.
Trophic markers enable the characterization of diets, predator–prey relationships and trophic relations in food webs, especially when higher trophic levels are involved.
The observed differences in our target fish species were based on differences in SI as well as FA profiles, and hence the trophic position of each species within the local food web was highlighted. Furthermore, the present data enabled insights into the functional diversity of the Helgoland intertidal food web in the North Sea, thus allowing a characterization of this specific community. Furthermore, we were able to highlight the advantage of using a multi-dimensional trophic marker approach, thus strengthening the evaluation in the light of trophic niche partitioning. This permitted a detailed description of trophic structures in our local food web. Further, among- and within-species interactions and the advantages of single trophic markers were investigated in detail. Here, especially the focus on intraspecific interactions within P. minutus elucidated a high ecological diversity in this community thus pointing at the great advantage of multi-tracer approaches in food web studies.
Aberle N, Malzahn AM (2007) Interspecific and nutrient-dependent variations in stable isotope fractionation: experimental studies simulating pelagic multitrophic systems. Oecologia 154:291–303
Aberle N, Hillebrand H, Grey J, Wiltshire KH (2005) Selectivity and competitive interactions between two benthic invertebrate grazers (Asellus aquaticus and Potamopyrgus antipodarum): an experimental study using 13C- and 15N-labelled diatoms. Freshw Biol 50:369–379
Alfaro AC, Thomas F, Sergent L, Duxbury M (2006) Identification of trophic interactions within an estuarine food web (northern New Zealand) using fatty acid biomarkers and stable isotopes. Estuar Coast Shelf Sci 70:271–286
Bearhop S, Waldron S, Votier SC, Furness RW (2002) Factors that influence assimilation rates and fractionation of nitrogen and carbon stable isotopes in avian blood and feathers. Physiol Biochem Zool 75:451–458
Bearhop S, Adams CE, Waldron S, Fuller RA, Macleod H (2004) Determining trophic niche width: a novel approach using stable isotope analysis. J Anim Ecol 73:1007–1012
Bell M, Henderson R, Sargent J (1986) The role of polyunsaturated fatty acids in fish. Comp Biochem Physiol A: Mol Integr Physiol 83:711–719
Bolnick DI, Svanbäck R, Fordyce JA, Yang LH, Davis JM, Hulsey CD, Forister ML (2003) The ecology of individuals: incidence and implications of individual specialization. Am Nat 161:1–28
Bouchereau JL, Guelorget O (1998) Comparison of three Gobiid ae (Teleostei) life history strategies over their geographical range. Oceanol Acta 21:503–517
Budge SM, Iverson SJ, Bowen WD, Ackman RG (2002) Among- and within-species variability in fatty acid signatures of marine fish and invertebrates on the Scotian Shelf, Georges Bank, and southern Gulf of St. Lawrence. Can J Fish Aquat Sci 59:886–898
Dalsgaard J, St John M (2004) Fatty acid biomarkers: validation of food web and trophic markers using 13C-labelled fatty acids in juvenile sandeel (Ammodytes tobianus). Can J Fish Aquat Sci 61:1671–1680
Dalsgaard J, St John M, Kattner G, Muller-Navarra D, Hagen W (2003) Fatty acid trophic markers in the pelagic marine environment. Adv Mar Biol 46(46):225–340
Davidson D, Marshall CT (2010) Are morphometric indices accurate indicators of stored energy in herring Clupea harengus? J Fish Biol 76:913–929
DeNiro M, Epstein S (1978) Influence of diet on the distribution of carbon isotopes in animals. Geochim Cosmochim Acta 42:495–506
Dolbeth M, Martinho F, Leitao R, Cabral H, Pardal MA (2008) Feeding patterns of the dominant benthic and demersal fish community in a temperate estuary. J Fish Biol 72:2500–2517
Ehrenbaum E (1936) Naturgeschichte und wissenschaftliche Bedeutung der Seefische Nordeuropas. E. Schweitzerbart’sche Verlagsbuchhandlung (Erwin Naegele), Stuttgart:337
El-Sabaawi R, Dower JF, Kainz M, Mazumder A (2009) Characterizing dietary variability and trophic positions of coastal calanoid copepods: insight from stable isotopes and fatty acids. Mar Biol 156:225–237
Fjosne K, Gjosaeter J (1996) Dietary composition and the potential of food competition between 0-group cod (Gadus morhua L.) and some other fish species in the littoral zone. ICES J Mar Sci 53:757–770
France RL (1995) Carbon-13 enrichment in benthic compared to planktonic algae: foodweb implications. Mar Ecol Prog Ser 124:307–312
Froese R, Pauly D (2000) FishBase 2000: concepts, design and data sources. WorldFish, Penang
Fry B (2006) Stable isotope ecology. Springer, Berlin
Fry B, Mumford PL, Tam F, Fox DD, Warren GL, Havens KE, Steinman AD (1999) Trophic position and individual feeding histories of fish from Lake Okeechobee, Florida. Can J Fish Aquat Sci 56:590–600
Hagen W, Yoshida T, Virtue P, Kawaguchi S, Swadling KM, Nicol S, Nichols PD (2007) Effect of a carnivorous diet on the lipids, fatty acids and condition of Antarctic krill, Euphausia superba. Antarct Sci 19:183–188
Hart PJB (1993) Teleost foraging: facts and theories
Herman DP, Burrows DG, Wade PR, Durban JW, Matkin CO, LeDuc RG, Barrett-Lennard LG, Krahn MM (2005) Feeding ecology of eastern North Pacific killer whales Orcinus orca from fatty acid, stable isotope, and organochlorine analyses of blubber biopsies. Mar Ecol Prog Ser 302:275–291
Hilldén N-O (1978) On the feeding of the goldsinny, Ctenolabrus rupestris L. (Pisces, Labridae). Ophelia 17:195–198
Hureau J, Bauchot M, Nielsen J, Tortonese E (1986) Fishes of the North-eastern Atlantic and the Mediterranean, vol 3. Unesco, Paris
Inger R, Bearhop S (2008) Applications of stable isotope analyses to avian ecology. Ibis 150:447–461
Ingram T, Matthews B, Harrod C, Stephens T, Grey J, Markel R, Mazumder A (2007) Lipid extraction has little effect on the delta N-15 of aquatic consumers. Limnol Oceanogr Methods 5:338–343
Jackson AL, Inger R, Parnell AC, Bearhop S (2011) Comparing isotopic niche widths among and within communities: SIBER—stable isotope bayesian ellipses in R. J Anim Ecol 80:595–602
Janke K (1990) Biological interactions and their role in community structure in the rocky intertidal of Helgoland (German Bight, North Sea). Helgolaender Meeresunters 44:219–263
Jaquet N, Raffaelli D (1989) The ecological importance of the sand goby Pomatoschistus minutus (Pallas). J Exp Mar Biol Ecol 128:147–156
Kelly JR, Scheibling RE (2012) Fatty acids as dietary tracers in benthic food webs. Mar Ecol Prog Ser 446:1–22
Kiljunen M, Grey J, Sinisalo T, Harrod C, Immonen H, Jones RI (2006) A revised model for lipid-normalizing delta C-13 values from aquatic organisms, with implications for isotope mixing models. J Appl Ecol 43:1213–1222
Layman CA, Arrington DA, Montana CG, Post DM (2007) Can stable isotope ratios provide for community-wide measures of trophic structure? Ecology 88:42–48
Legendre P, Legendre L (1998) Numerical ecology (2nd English ed.). Elsevier Science BV, Amsterdam
Leitao R, Martinho E, Neto JM, Cabral H, Marques JC, Pardal MA (2006) Feeding ecology, population structure and distribution of Pomatoschistus microps (Kroyer, 1838) and Pomatoschistus minutus (Pallas, 1770) in a temperate estuary, Portugal. Estuar Coast Shelf Sci 66:231–239
Malzahn AM, Aberle N, Clemmesen C, Boersma M (2007) Nutrient limitation of primary producers affects planktivorous fish condition. Limnol Oceanogr 52:2062–2071
McCutchan JH, Lewis WM, Kendall C, McGrath CC (2003) Variation in trophic shift for stable isotope ratios of carbon, nitrogen, and sulfur. Oikos 102:378–390
McLeod RJ, Hyndes GA, Hurd CL, Frew RD (2013) Unexpected shifts in fatty acid composition in response to diet in a common littoral amphipod. Mar Ecol Prog Ser 479:1–12
Minagawa M, Wada E (1984) Stepwise enrichment of 15N along food-chains—further evidence and the relation between δ 15N and animal age. Geochim Cosmochim Acta 48:1135–1140
Mozsar A, Boros G, Saly P, Antal L, Nagy SA (2015) Relationship between Fulton’s condition factor and proximate body composition in three freshwater fish species. J Appl Ichthyol 31:315–320
Nash RDM, Valencia AH, Geffen AJ (2006) The origin of Fulton’s condition factor—setting the record straight. Fisheries 31:236–238
Newsome SD, del Rio CM, Bearhop S, Phillips DL (2007) A niche for isotopic ecology. Front Ecol Environ 5:429–436
Norderhaug K, Christie H, Fossa J, Fredriksen S (2005) Fish-macrofauna interactions in a kelp (Laminaria hyperborea) forest. J Mar Biol Assoc UK (JMBA) 85:1279–1286
Nyssen F, Brey T, Dauby P, Graeve M (2005) Trophic position of Antarctic amphipods—enhanced analysis by a 2-dimensional biomarker assay. Mar Ecol Prog Ser 300:135–145
Oksanen J, Blanchet F, Kindt R, Legendre P, O’Hara R, Simpson G, Solymos P, Stevens M, Wagner H (2013) Vegan: community ecology package. 2013. R package version 2.0-7
Perga ME, Kainz M, Matthews B, Mazumder A (2006) Carbon pathways to zooplankton: insights from the combined use of stable isotope and fatty acid biomarkers. Freshw Biol 51:2041–2051
Peterson BJ, Fry B (1987) Stable isotopes in ecosystem studies. Annu Rev Ecol Syst 18:293–320
Piche J, Iverson SJ, Parrish FA, Dollar R (2010) Characterization of forage fish and invertebrates in the Northwestern Hawaiian Islands using fatty acid signatures: species and ecological groups. Mar Ecol Prog Ser 418:1-U410
Pinnegar JK, Polunin NVC (1999) Differential fractionation of δ 13C and δ 15N among fish tissues: implications for the study of trophic interactions. Funct Ecol 13:225–231
Polacik M, Harrod C, Blazek R, Reichard M (2014) Trophic niche partitioning in communities of African annual fish: evidence from stable isotopes. Hydrobiologia 721:99–106
Post DM (2002) Using stable isotopes to estimate trophic position: models, methods, and assumptions. Ecology 83:703–718
R Core Team (2013) A language and environment for statistical computing. R foundation for statistical computing, Vienna
Ross ST (1986) Resource partitioning in fish assemblages—a review of field studies. Copeia 352–388
Sargent J, Bell G, McEvoy L, Tocher D, Estevez A (1999a) Recent developments in the essential fatty acid nutrition of fish. Aquaculture 177:191–199
Sargent J, McEvoy L, Estevez A, Bell G, Bell M, Henderson J, Tocher D (1999b) Lipid nutrition of marine fish during early development: current status and future directions. Aquaculture 179:217–229
Schmidt K, Atkinson A, Petzke K-J, Voss M, Pond DW (2006) Protozoans as a food source for Antarctic krill, Euphausia superba: complementary insights from stomach content, fatty acids, and stable isotopes. Limnol Oceanogr 51:2409–2427
Stevens CJ, Deibel D, Parrish CC (2004) Incorporation of bacterial fatty acids and changes in a wax ester-based omnivory index during a long-term incubation experiment with Calanus glacialis Jaschnov. J Exp Mar Biol Ecol 303:135–156
Trippel EA, Hubert JJ (1989) Common statistical errors in fishery research. Department of Mathematics and Statistics, University of Guelph, Guelph
Wennhage H, Pihl L (2002) Fish feeding guilds in shallow rocky and soft bottom areas on the Swedish west coast. J Fish Biol 61:207–228
Wheeler A, Du Heaume V (1969) The fishes of the British Isles and north-west Europe. Macmillan, London
Wilkins HKA, Myers AA (1992) Microhabitat utilization by an assemblage of temperate Gobiidae (Pisces: Teleostei). Mar Ecol Prog Ser 90:103–112
Special thanks to the technicians of the Alfred Wegener Institute, Helmholtz Centre for Polar and Marine Research, Biological Station Helgoland for their support in sample processing and lab analyses. We are grateful to the captains and crew members of RV Uthörn and RV Aade for providing fish samples. Further, our colleagues and volunteers within the AWI food web group are thanked for their kind support during beach seining campaigns regardless of wind and weather conditions. This project was funded by ERAnet (MarinERA project reference no. 189570, partner 4: AWI, PTJ Förderzeichen 03F0600A). In memory of our highly esteemed project leader, colleague and friend Vitor Almada.
Communicated by H.-D. Franke.
About this article
Cite this article
Hielscher, N.N., Malzahn, A.M., Diekmann, R. et al. Trophic niche partitioning of littoral fish species from the rocky intertidal of Helgoland, Germany. Helgol Mar Res 69, 385–399 (2015). https://doi.org/10.1007/s10152-015-0444-5
- Fatty acids
- Stable isotopes
- Ctenolabrus rupestris
- Pomatoschistus minutus
- Pomatoschistus pictus
- Myoxocephalus scorpius
- Taurulus bubalis
- Trophic ecology
- Feeding ecology