Skip to main content
Advertisement
Browse Subject Areas
?

Click through the PLOS taxonomy to find articles in your field.

For more information about PLOS Subject Areas, click here.

  • Loading metrics

A non-lethal stable isotope analysis of valued freshwater predatory fish using blood and fin tissues as alternatives to muscle tissue

  • Lukáš Vejřík ,

    Contributed equally to this work with: Lukáš Vejřík, Ivana Vejříková, Mojmír Vašek

    Roles Conceptualization, Investigation, Methodology, Project administration, Visualization, Writing – original draft, Writing – review & editing

    Affiliation Biology Centre of the Czech Academy of Sciences, Institute of Hydrobiology, České Budějovice, Czech Republic

  • Ivana Vejříková ,

    Contributed equally to this work with: Lukáš Vejřík, Ivana Vejříková, Mojmír Vašek

    Roles Conceptualization, Writing – original draft, Writing – review & editing

    ivana.vejrikova@hbu.cas.cz

    Affiliation Biology Centre of the Czech Academy of Sciences, Institute of Hydrobiology, České Budějovice, Czech Republic

  • Zuzana Sajdlová ,

    Roles Investigation, Visualization, Writing – review & editing

    ‡ ZS, LK, TK, DB, TJ, PB, JP, and MC also contributed equally to this work.

    Affiliation Biology Centre of the Czech Academy of Sciences, Institute of Hydrobiology, České Budějovice, Czech Republic

  • Luboš Kočvara ,

    Roles Investigation, Writing – review & editing

    ‡ ZS, LK, TK, DB, TJ, PB, JP, and MC also contributed equally to this work.

    Affiliation Biology Centre of the Czech Academy of Sciences, Institute of Hydrobiology, České Budějovice, Czech Republic

  • Tomáš Kolařík ,

    Roles Investigation, Writing – review & editing

    ‡ ZS, LK, TK, DB, TJ, PB, JP, and MC also contributed equally to this work.

    Affiliation Biology Centre of the Czech Academy of Sciences, Institute of Hydrobiology, České Budějovice, Czech Republic

  • Daniel Bartoň ,

    Roles Investigation, Writing – review & editing

    ‡ ZS, LK, TK, DB, TJ, PB, JP, and MC also contributed equally to this work.

    Affiliation Biology Centre of the Czech Academy of Sciences, Institute of Hydrobiology, České Budějovice, Czech Republic

  • Tomáš Jůza ,

    Roles Investigation, Writing – review & editing

    ‡ ZS, LK, TK, DB, TJ, PB, JP, and MC also contributed equally to this work.

    Affiliation Biology Centre of the Czech Academy of Sciences, Institute of Hydrobiology, České Budějovice, Czech Republic

  • Petr Blabolil ,

    Roles Investigation, Visualization, Writing – review & editing

    ‡ ZS, LK, TK, DB, TJ, PB, JP, and MC also contributed equally to this work.

    Affiliations Biology Centre of the Czech Academy of Sciences, Institute of Hydrobiology, České Budějovice, Czech Republic, Faculty of Science, University of South Bohemia in České Budějovice, České Budějovice, Czech Republic

  • Jiří Peterka ,

    Roles Funding acquisition, Investigation, Writing – review & editing

    ‡ ZS, LK, TK, DB, TJ, PB, JP, and MC also contributed equally to this work.

    Affiliation Biology Centre of the Czech Academy of Sciences, Institute of Hydrobiology, České Budějovice, Czech Republic

  • Martin Čech ,

    Roles Funding acquisition, Investigation, Writing – review & editing

    ‡ ZS, LK, TK, DB, TJ, PB, JP, and MC also contributed equally to this work.

    Affiliation Biology Centre of the Czech Academy of Sciences, Institute of Hydrobiology, České Budějovice, Czech Republic

  • Mojmír Vašek

    Contributed equally to this work with: Lukáš Vejřík, Ivana Vejříková, Mojmír Vašek

    Roles Investigation, Writing – review & editing

    Affiliation Biology Centre of the Czech Academy of Sciences, Institute of Hydrobiology, České Budějovice, Czech Republic

Abstract

Stable isotope analysis (SIA) is widely used to study trophic ecology and food webs in aquatic ecosystems. In the case of fish, muscle tissue is generally preferred for SIA, and the method is lethal in most cases. We tested whether blood and fin clips can be used as non-lethal alternatives to muscle tissue for examining the isotopic composition of two freshwater predatory fish, European catfish (Silurus glanis) and Northern pike (Esox lucius), species of high value for many freshwater systems as well as invasive species in many others. Blood samples from the caudal vein, anal fin clips, and dorsal muscle obtained by biopsy punch were collected from four catfish and pike populations (14–18 individuals per population). Subsequently, these samples were analyzed for δ13C and δ15N. The effects of alternative tissues, study site, and fish body mass on the isotopic offset were investigated. Both species showed a correlation between the isotopic offset and the tissue type, as well as the study site, but no significant relationship with the body mass. The isotopic offsets between tissues were used to calculate the conversion equations. The results demonstrated that both blood and fin clips are suitable and less invasive alternative to muscle in SIA studies focused on European catfish and Northern pike. Blood provided better correspondence to muscle isotope values. However, our results clearly demonstrated that isotopic offsets between tissues vary significantly among populations of the same species. Therefore, obtaining a muscle biopsy from several individuals in any population is advisable to gain initial insights and establish a possible population-specific inter-tissue conversion.

Introduction

Stable isotopes are widely used in aquatic science [13], primarily δ13C and δ15N [3, 4]. They are used in a broad spectrum of studies such as indicating nutrient pollution [5, 6], determining contaminant bioaccumulation [7, 8], tracking changes in carbon cycle over time [9], investigating aquatic food webs [1012], and even assessing individual specialization [13]. Muscle tissue is commonly used for stable isotope analysis (SIA) of vertebrates [14], specifically the dorsal white muscle in fish studies [15]. Approximately 0.5 mg of dry mass is necessary for one analysis which constitutes a relatively substantial and often even lethal intrusion into the fish body [16]. Fish typically occupy higher positions within food webs underscoring the significance of their inclusion in the analysis of trophic structure and energy flow within aquatic ecosystems [17]. Nevertheless, the considerable size and high trophic position of predatory fish species impose constraints on sample sizes [12, 18]. Predatory fish exhibit notably lower population densities, and capturing large individuals through conventional sampling methods is usually challenging [12]. The substantial reduction of key species populations can also have a significant impact on ecosystem functioning [19]. Additionally, predatory fish are intentionally introduced into various systems (e.g., drinking water reservoirs) to exert a top-down effect on lower trophic levels. Therefore, it is preferable not to reduce their population density by lethal sampling [20]. In certain locations, such actions are even prohibited by the regulations of respective authorities.

The use of non-lethal sampling in stable isotope studies of fish, aimed at preventing unnecessary mortality of sampled animals, has recently received significant attention [16, 17, 21, 22]. Recent studies have explored non-lethal and minimally invasive alternatives, including scales, mucus and fins [17, 2225]. The use of fin tissue has been the most prevalent [12], and such tissue is capable of regeneration [26, 27]. In contrast, collecting and employing blood for fish studies is relatively infrequent [28], despite its common application in various analyses, including SIA of higher vertebrates [29, 30]. The sampling is straightforward, particularly for large individuals, and results in minimal injury [13] compared to the muscle tissue sampling through biopsy punch [17, 31]. Although a biopsy itself is considered a non-lethal method for obtaining tissue, and according to Henderson et al. [31], it can even be performed in fish under 30 cm in size, the primary advantage of blood tissue collection over biopsy lies in the easier and gentler collection process. Additionally, there is a lower risk of introducing infections due to minimal external injury [32, 33]. This method is more comfortable for both the sampled individual and the person conducting the sampling. Despite the fact that biopsies often disrupt vessels in the muscle, causing bleeding in volumes greater than required for SIA blood analysis (personal observation), blood collection from fish results in complete healing within 2–3 weeks, demonstrating a 100% survival rate in teleost blood draws of 1 μL g-1 with no post-treatment [34]. Approximately 200 μL of blood is sufficient for SIA, making this collection method safe for fish as small as 200 g of body mass [32, 34]. While it is possible to reduce the risk of introducing infections into biopsy wounds by using techniques such as applying Fish BandageTM and a non-toxic, non-allergenic cellulose-based powder that forms a clear viscous gel upon contact with water and can treat skin ulcers and bind open wounds in fish when combined with an antiseptic [31], this represents an additional step in the sampling process, consequently increasing the sampling time.

Several studies have demonstrated correlations between isotopic signals from ’non-lethal tissues’ and those from muscle tissue [16, 17, 22, 35, 36]. Nevertheless, it is essential to consider that the isotopic offsets, defined as the differences in isotope value between individual tissues, are not solely attributed to variations in tissue protein composition [37]. Stable isotope biokinetics also play a role in this regard [13, 38], and although often overlooked, they likely have a significant impact on differences among populations and various life stages within the same species [21, 39]. The biokinetics of stable isotopes in various tissue types, species, and environments is a complex issue. Unfortunately, current findings do not provide clear conclusions about the isotopic turnover rate of different tissues [13, 38, 40, 41]. Turnover depends not only on tissue type but also on individual size and environmental temperature [38]. For an average organism weighing 1 mg, the isotopic half-life at 10°C is nine days, and at 40°C, it is only three days for both 13C and 15N (full turnover is about four to five times longer). However, for a 100 g organism, the half-life at 10°C is 80 days for 13C and 84 days for 15N. At 40°C, these values decrease to 28 days for 13C and 30 days for 15N. In the case of a 100 kg organism, the half-life at 10°C is 303 days for 13C and 321 days for 15N. At 40°C, these values decrease significantly to 105 days for 13C and 115 days for 15N [38]. As shown above, isotopic turnover rates vary among individual elements. It appears that 15N turnover rates are slightly longer than those of 13C [38, 40]. However, in the case of summer flounder (Paralichthys dentatus), no difference in the half-life between 15N and 13C was observed [41].

Based on several studies, the order of tissues in terms of isotopic half-life, from fastest to slowest, is approximately as follows (with slight variations in individual studies): Plasma → Liver → Fin → Heart → Mucus → Blood → Red blood cells → Bone collagen → Scale → Muscle [13, 38, 40, 41]. Another important factor that can influence isotopic turnover is the type of diet. It appears that different body tissues respond to dietary changes by altering their isotopic turnover rates. Tissues with longer isotopic turnover tend to exhibit greater differences depending on the type of diet [40].

Although the notable increase and growing trend towards non-invasive sampling are evident, there remains a lack of information concerning the relationship between isotope signatures of muscle and less harmful tissues for numerous species and populations [17, 22]. This study aims to contribute to the existing knowledge by examining the use of fin clips and blood in comparison to muscle tissue for two key apex predator fish species: European catfish (Silurus glanis) and Northern pike (Esox lucius) [12, 42]. The studied species are among the most widespread predatory fish in Europe and thus play a crucial role in maintaining ecosystem stability through top-down control [43]. To mitigate the influence of any single location, samples were collected from four distinct study sites. The primary aim was to investigate the feasibility of using fin and blood tissues as less harmful alternatives for δ15N and δ13C analyses in place of muscle tissue. The specific aims were as follows: (i) to examine the correlation between blood and fin tissues isotope signatures and those of muscle tissue, (ii) to assess the suitability of blood or fin tissue as a viable alternative, (iii) to develop appropriate conversion equations to enable the interchangeable use of different tissues in future stable isotope studies, and (iv) to explore whether variations in isotope signatures among tissues correlate with fish body mass differences.

Methods

Study design

Fishes were treated in accordance with the Experimental Animal Welfare Commission under the Ministry of Agriculture of the Czech Republic guidelines (Ref. No. CZ 01679). The work was approved by the Ethics Committee of the Czech Academy of Sciences.

The study was conducted in 2017 at four water bodies of similar size and fish species composition, but varying in their trophic states: two oligotrophic post-mining lakes, Most and Milada, and two meso-eutrophic reservoirs, Žlutice and Římov, Czech Republic (Fig 1). For basic parameters of the study sites see Table 1, and for further details see [20, 44].

thumbnail
Fig 1. Map of the Czech Republic with four study sites, Most and Milada artificial lakes, and Žlutice and Římov reservoirs.

Lakes and reservoirs share the same scale.

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

thumbnail
Table 1. Basic hydrological and geographical parameters of the four study sites.

Total phosphorus (TP) presents an average from longitudinal profiles of three areas sampled four times per year in 2013–2017.

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

Electrofishing (method description by [45]) and long-lines (method description by [20, 46, 47]) were used in 2017 to capture adult individuals of European catfish and Northern pike on 4–7 September in Most Lake, 11–14 September in Milada Lake, 25–28 July in Žlutice Reservoir, and 18–21 July in Římov Reservoir.

Body size and mass of European catfish ranged from 630 to 1,500 mm SL (mean±SD: 1046±214 mm) and 1.6 to 21.8 kg (mean±SD: 8.35±5.18 kg), and Northern pike body size ranged from 450 to 1,200 mm SL (729±202 mm), and mass ranged from 0.6 to 13.8 kg (3.75±3.06 kg). All individuals were weighed to the nearest gram, and samples of three body tissues were collected for SIA under anesthesia, a bath containing clove oil. Specifically, a small sample of the anal fin (1 cm2, resection), 0.2 cm3 of dorsal muscle (biopsy punch Miltex, skin removed), and 1 mL of blood from the caudal vein using an 18-gauge needle Sterican were taken. The three tissue samples from each individual were placed on ice and transferred to a laboratory freezer for SIA. Subsequently, all European catfish and Northern pike were released back into the water bodies.

Stable isotope analysis (SIA)

All frozen samples for SIA were subsequently dried at 60°C for 48h and ground into a uniform powder using a Retsch MM 200 ball-mill (Retsch GmbH, Haan, Germany). Minor subsamples (0.52–0.77 mg) were carefully placed into tin cups for δ13C and δ15N analysis. All SIA were performed using a FlashEA 1112 elemental analyzer coupled to a Thermo Finnigan DELTAplus Advantage mass spectrometer (Thermo Fisher Scientific Corporation, USA) at the University of Jyväskylä, Finland. Carbon and nitrogen isotope ratios are presented as δ13C and δ15N values relative to the international standards for carbon (Vienna PeeDeeBelemnite, Austria) and nitrogen (atmospheric nitrogen). The precision of the analytical results was determined to be ±0.20‰ through repeated analyses of a consistent standard (Northern pike white muscle tissue), which was inserted in each analytical run following every five samples. The data were not adjusted for lipid content due to elemental carbon–nitrogen (C:N) ratios observed in all tissues being less than 4 (except for a single catfish blood sample at Římov Reservoir) indicating a low lipid content [48].

Data analysis

Prior the analyses, the data distribution was visually checked by histogram plotting. No obvious skewness distribution was observed. Paired t-tests were conducted to assess whether blood and fin isotope values differ from those of muscle, and also to determine whether C:N ratios in blood and fin of individuals at each study site differ compared to those muscle tissue. The isotopic offsets between tissues, defined as the differences between blood and muscle or fin and muscle stable isotope values, were compared among study sites (Milada, Most, Žlutice, Římov) for European catfish and Northern pike using a one-way repeated measures analysis of variance (ANOVA) with multiple comparisons (Tukey HSD test). Additionally, the one-way ANOVA was used to assess the differences in C:N ratios among all tissue types across all study sites. The statistical analyses were performed in STATISTICA 9.1 [49] and used the standard level of significance of α = 0.05. Linear model with the random effect of site was applied to model the relationships between mussel isotopes with the isotopes in fin and blood and between the isotope offsets and fish mass. The models were prepared for each species and C and N isotope separately. After fitting the models, regression diagnostics was applied by plotting the residual vs. fitted values and Normal QQ Plot to check the residuals variance. The models were developed using R packages “lme4” [50] and “nlme” [51] in R software [52].

Results

The δ13C and δ15N values in blood and fins of European catfish and Northern pike were generally found to be significantly different from those in their muscle tissue. However, the difference in δ13C between the blood and muscle was not significant for European catfish at one study site and for Northern pike at two study sites. Regarding δ15N, there was no significant difference between the blood and muscle of Northern pike in two study sites, and between the fin and muscle of Northern pike in one study site (Table 2).

thumbnail
Table 2. Mean ± SE values of δ13C and δ15N, in muscle (M), blood (B) and fin (F) tissues of European catfish (top) and Northern pike (bottom) at four study sites (Most, Milada, Žlutice and Římov).

Significant differences between offset (means how given tissue is enriched or depleted compared to muscle tissue) of blood and muscle (B) and fin and muscle (F) (paired t-test) were examined. The C:N ratios were measured in three tissues of European catfish (top) and Northern pike (bottom). Asterisks indicate a significance (paired t-test, p > 0.05 = -, p < 0.05 = *, p < 0.01 = **, p < 0.001 = ***).

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

Blood generally showed neither clear enrichment nor depletion of δ13C compared to muscle tissue. In European catfish, the blood exhibited a slight δ13C enrichment at three study sites and a slight depletion at one study site. For Northern pike, the blood displayed slight δ13C enrichment at only one study site and slight depletion at three study sites. The fin was consistently and significantly enriched in δ13C in comparison to the muscle tissue in both species across all study sites (Table 2). In both European catfish and Northern pike, both the blood and fins were consistently and significantly enriched in δ15N when compared to the muscle tissue, observed across most study sites. The degree of these differences was more pronounced in the blood compared with fins of European catfish. In the case of Northern pike, on average, the degree of these differences between blood and fin was similar. The difference in both δ13C and δ15N between less harmful tissues and muscle tissue tended to be greater in European catfish when compared to Northern pike (Table 2).

For European catfish, both δ13C and δ15N offsets between both blood and muscle, and fin and muscle significantly differed among study sites. The statistical values for δ13C offset were: ANOVA: F3,62 = 6.87, p < 0.001 for blood vs. muscle, and ANOVA: F3,62 = 4.26, p < 0.01 for fin vs. muscle. The statistical values for δ15N offset were: ANOVA: F3,62 = 12.31, p < 0.001 for blood vs. muscle, and ANOVA: F3,62 = 4.30, p < 0.01 for fin vs. muscle.

Also for Northern pike, the differences in both δ13C and δ15N offsets between both blood and muscle, and fin and muscle were significant among study sites. The statistical values for δ13C offset were: ANOVA: F3,56 = 5.67, p < 0.01 for blood vs. muscle, and ANOVA: F3,56 = 9.77, p < 0.001 for fin vs. muscle. The statistical values for δ15N offset were: ANOVA: F3,56 = 3.81, p < 0.05 for blood vs. muscle, and ANOVA: F3,56 = 7.44, p < 0.001 for fin vs. muscle.

Tukey HSD comparisons of pairwise offset differences between study sites revealed significant distinction. Specifically, there were significant differences in δ13C of European catfish between blood and muscle in two cases, between fin and muscle in one case. Additionally, for δ15N, significant differences were between blood and muscle in three cases and between fin and muscle in one case (Table 3). For Northern pike, Tukey HSD comparisons of pairwise offset differences between study sites revealed significant differences in δ13C between blood and muscle in one case, between fin and muscle in three cases, and for δ15N between blood and muscle in one case and between fin and muscle in two cases (Table 3).

thumbnail
Table 3. Tukey HSD statistical values (p) of the difference in δ13C and δ15N offsets between blood and muscle (B—M), and fin and muscle (F—M) among study sites (Most, Milada, Žlutice, Římov).

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

The C:N ratios were, on average, lowest in the muscle tissue of both European catfish and Northern pike. In contrast, blood reached the highest average values (Fig 2). In the case of European catfish, both blood and fin reached significantly higher C:N values than muscle in Milada and Most lakes (paired t-test, p < 0.001). Also in Žlutice Reservoir, C:N ratio was significantly higher in blood than in muscle of European catfish (paired t-test, p < 0.05), however no difference was found between fin and muscle (paired t-test t17 = -0.75, p = 0.46). In Římov Reservoir, there was no difference in C:N ratio between both blood and muscle (paired t-test t15 = -0.80, P = 0.44) and fin and muscle (paired t-test t15 = -0.04, p = 0.97) of European catfish. For Northern pike, the C:N ratio was consistently and significantly higher in both blood and fin compared to muscle, across all study sites (paired t-test, p < 0.001) (Fig 2). Among study sites, C:N ratios were significantly different in muscle and fin of European catfish (muscle: ANOVA: F3,62 = 7.98, p < 0.001; fin: ANOVA: F3,62 = 5.69, p < 0.001) and Northern pike (muscle: ANOVA: F3,56 = 6.30, p < 0.001; fin: ANOVA: F3,56 = 6.60, p < 0.001). In contrast, there was no significant difference in C:N ratio in blood among the study sites either for European catfish (ANOVA: F3,62 = 0.16, p = 0.92) or Northern pike (ANOVA: F3,56 = 1.03, p = 0.39).

thumbnail
Fig 2.

C:N ratios of European catfish (a) and Northern pike (b) in tissues (muscle: blue, fin: gray, blood: red) at four study sites (Most, Milada, Žlutice, Římov). Box and whiskers plots: upper and lower quartiles (boxes), median values (line inside the boxes), mean value (crosses), maximum and minimum values (whiskers), and outliers (circles) are shown.

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

The linear regression revealed significant relationships between the δ15N and δ13C values in muscle, fin and blood of both European catfish and Northern pike (Figs 3 and 4). Thus, both less harmful tissues served as reliable predictors for muscle. Specifically, blood demonstrated strong predictive capabilities for both δ13C and δ15N values in the muscle of both European catfish and Northern pike. For European catfish, the coefficients of determination (R2) for δ13C were, on average, 0.83 and 0.79 for blood and fin, respectively (Fig 3). Similarly, for δ15N, the R2 values were 0.91 and 0.84 for blood and fin, respectively (Fig 4). Regarding Northern pike, the R2 values for δ13C averaged 0.94 and 0.86 for blood and fin, respectively (Fig 3), and for δ15N the R2 values were 0.92 and 0.88 in blood and fin, respectively (Fig 4). The regression diagnostics of isotope linear model with the random effect of site demonstrate uniform residuals variance and lack of high leverage effects. Models intercepts and standard errors in fixed effects were higher for European catfish compared to Northern pike and slopes were comparable with exception of δ13C for European catfish with lower slope (Table 4). Similarly, random effects were higher for European catfish compared to Northern pike (Table 4).

thumbnail
Fig 3.

Linear regressions between δ13C of less harmful tissues (blood: black circles and line, fin: gray circles and line) and muscle tissue for European catfish (left column) and Northern pike (right column) in Most lake (a, e), Milada lake (b, f), Žlutice reservoir (c, g) and Římov reservoir (d, h). The same scale is maintained on both axes in all images. The linear regression equations for European catfish and Northern pike to convert δ13C values in blood (black equations) and fin (gray equations) to those in muscle are presented with added statistical significances of the linear regression and the 95% confidence intervals of the slopes.

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

thumbnail
Fig 4.

Linear regressions between δ15N of less harmful tissues (blood: black circles and line, fin: gray circles and line) and muscle tissue for European catfish (left column) and Northern pike (right column) in Most lake (a, e), Milada lake (b, f), Žlutice reservoir (c, g) and Římov reservoir (d, h). The same scale is maintained on both axes in all images. The linear regression equations for European catfish and Northern pike to convert δ13C values in blood (black equations) and fin (gray equations) to those in muscle are presented with added statistical significances of the linear regression and the 95% confidence intervals of the slopes.

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

thumbnail
Table 4. Parameters and coefficients of isotope linear model with the random effect of site.

https://doi.org/10.1371/journal.pone.0297070.t004

No distinct trends were found through linear regression analysis between isotopic tissue offsets and fish body mass in either species or in any study site. For European catfish, the R2 reached mean values of 0.11 ± 0.17 SD for the blood and muscle δ13C offset, and 0.17 ± 0.17 SD for the fin and muscle δ13C offset in relation to the body mass. Regarding δ15N, the R2 reached mean values of 0.24 ± 0.22 SD for the blood and muscle offset, and 0.05 ± 0.04 SD for the fin and muscle offset in relation to the body mass. For Northern pike, the R2 reached mean values of 0.10 ± 0.11 SD for the blood and muscle δ13C offset, and 0.01 ± 0.02 SD for the fin and muscle δ13C offset. In the case of δ15N, the R2 reached mean values of 0.07 ± 0.07 SD for the blood and muscle offset, and 0.04 ± 0.05 SD for the fin and muscle offset in relation to the body mass. The regression diagnostics of offset isotope linear model with the random effect of site demonstrate uniform residuals variance and lack of high leverage effects. Models’ intercepts were the highest for European catfish offset δ15N for all tissues combination, offset δ13C fin to body mass and Northern pike offset δ13C fin to body mass comparison (Table 5). The slopes values were generally low with the highest values for European catfish offset δ13C. The random effect intercept and residual variances were the highest for European catfish offset δ13C (Table 5).

thumbnail
Table 5. Parameters and coefficients of isotope offsets linear model with the random effect of site.

https://doi.org/10.1371/journal.pone.0297070.t005

Discussion

Sampling white muscle typically involves sacrificing the fish [16]. However, in our case, none of the predatory fish were directly sacrificed. We employed a biopsy punch to extract the muscle, a method feasible for larger fish species [33]. Nevertheless, this procedure still constitutes a relatively invasive intervention into the fish body, and the potential consequences upon its release into the water are generally unknown [31]. In contrast, methods such as fin clipping [22, 53] and blood collection [54] are significantly less disruptive. Blood sampling through the caudal vasculature is widely employed in fish biology for investigating health and physiology [54]. In live fish, it offers a rapid, uncomplicated, and comparatively less harmful approach for obtaining tissue, in contrast to the invasive and even more intricate nature of biopsies [31, 54]. The fin clip collection is a routine and widely used technique, employed either to acquire tissue samples for various analyses [35] or to mark fish [55].

In general, the differences in δ13C and δ15N isotope values between muscle and less harmful tissues in both predatory fish were smaller in the case of blood than in the fin. This is likely to be due to a more comparable protein composition between muscle and blood than between muscle and fin [37]. However, it appears that this effect is not related to the biokinetics of stable isotopes.

In our previous study [13], we calculated isotopic half-life [38] for European catfish and Northern pike in the presented study sites. We considered both the mean mass of the individuals and the temperature at which they were found, as determined by telemetry measurements [42]. Based on these calculations, the average isotopic half-life for blood and fin were quite similar in both species. Conversely, the half-life of these tissues in both species significantly differed from the notably longer half-life of muscle. For European catfish, the average isotopic half-life values were 39, 44 and 153 days for blood, fin, and muscle, respectively. In the case of Northern pike, the values were 39, 41 and 139 days for blood, fin and muscle [13]. However, as mentioned in the introduction, the issue of biokinetics is highly complex and can be influenced by various other factors, such as the diet composition, which is unknown in this case [40]. Nevertheless, significant variations in both δ13C and δ15N were observed in most instances between muscle and blood, or between muscle and fin. The fin was enriched in δ13C compared to muscle in both European catfish and Northern pike across all study sites. This pattern is consistent with findings from several other studies focusing on diverse fish species [22, 53, 56] as well as various Northern pike populations [57, 58]. The fin tissue was enriched also in δ15N in both species at nearly all study sites. This is noteworthy since many previous studies have reported δ15N depletion in fin tissue compared to muscle in other species [22, 53] or have identified an inconclusive trend between fin and muscle [17, 33, 56]. However, for example, in the case of callop (Macquaria ambigua) in the Australian Darling River, the fins were also enriched in δ15N compared to the muscle [17]. A similar pattern was observed with the bleak (Alburnus alburnus) in Spain [53]. Therefore, it appears that, unlike δ13C, the trend in δ15N is not as clear and may vary among species. Notably, Syväranta et al. [59] determined that the δ15N of European catfish fin tissue exhibited a slight but insignificant depletion in relation to muscle tissue. However, in order to compare isotopic signatures between tissues, these authors aggregated data from catfish individuals collected across four distinct locations. Consequently, it became impractical to draw conclusions about population-specific disparities in δ13C and δ15N between muscle and fin tissues. The two prior studies on Northern pike that investigated the isotopic relationship between fin and muscle [58, 60] revealed that fin δ15N was significantly depleted compared to muscle, whereas our observations demonstrated an opposing fin-muscle δ15N difference in three out of four pike populations (with no significant difference between the two tissues in δ15N within one studied population). Our findings undeniably illustrate that isotopic offsets between tissues can considerably differ among populations of the same species. Thus, it is unrealistic to assume the feasibility of devising a singular species-specific conversion equation that would universally apply to all populations of the given species. Rather, aligning with the perspectives of other researchers [39, 61, 62], we propose that a tissue-to-tissue conversion equation should be developed distinctly for each specific population of interest. The isotopic signal of individuals and their various tissues partly reflects the distinct protein composition of the tissues [37]. However, it is primarily influenced by the dietary composition of individuals or the entire populations. This diet composition can vary considerably within a species, both seasonally and across different localities [1113]. The diet of apex predators in our study sites was extensively examined, based on both SIA and the stomach content analysis [12, 63]. The study sites differed significantly from each other in terms of the diet composition of apex predators and their growth rates [12, 63]. In the oligotrophic sites of Milada and Most, apex predators exhibited significantly slower growth compared to the meso-eutrophic reservoirs Žlutice and Římov [63]. Moreover, at Milada and Most lakes, the diet of apex predators displayed high variety, including invertebrates, fish, and semi-aquatic prey such as waterfowl, mammals, and amphibians. In contrast, at the Římov and Žlutice, their diet consisted almost exclusively of fish [12, 63]. Particularly, semi-aquatic prey significantly differed in the SIA signal from other food sources [12, 59]. This distinct diet composition among study sites, along with varying growth rates, which likely affect the isotopic turnover rate, could be the main reason for the offset differences between tissues at different locations. This is supported by the observation that differences were evident in nine cases between oligotrophic and meso-eutrophic sites, but only in three cases between meso-eutrophic sites and once between the oligotrophic sites. In addition to the different protein composition of tissues [37], the dietary changes during the season are probably the main factor contributing to the offset differences between tissues of individual species within one population at one study site. This is further supported by the observation that the European catfish, which exhibits higher seasonal diet plasticity than the Northern pike [12, 13], shows a more pronounced difference in offset between tissues. The better body condition of catfish in the more productive Římov and Žlutice may also have caused the higher C:N ratio variability in muscle than in blood and fin due to a potentially higher level of lipid reserves stored in muscle tissue of some predators. The muscle tissue, being metabolically active, undergoes continuous metabolic changes such as growth, repair, and energy storage [64] that can result in a more variable range of C:N ratios compared to more stable tissues like blood and fin.

To the best of our knowledge, blood has not yet been proposed as an alternative tissue to muscle, and our study is also the first ever to compare the isotope values between blood and muscle in European catfish and Northern pike. Within the studied populations of the two predatory species, the isotopic offsets between blood and muscle were usually smaller compared to the isotopic offsets observed between fin and muscle. Particularly in the case of Northern pike, isotope values between blood and muscle did not significantly differ in half of the studied populations. Even in those populations where these values differed significantly, the offset between tissues was not substantial (ranging from 0.10 to 0.43‰; Table 2). These results indicate that in certain species and populations the isotopic offsets between blood and muscle of δ13C and δ15N may be negligible. However, predicting routinely in which populations such a situation will occur is challenging. Thus, it is essential to asses initially the isotopic offsets between tissues for each population of interest. If necessary, a population-specific inter-tissue conversion can then be calculated.

For both species, the offsets between tissues differed significantly across some of the study sites. We calculated the conversion equations for each study site individually, both for δ13C and δ15N in both blood and fin of both studied species. The conversion equations from blood to muscle and from fin to muscle were all significant. Their coefficients of determination were notably high: R2 for δ13C exceeded 0.77 and 0.46 for blood and fin, respectively, and for δ15N it exceeded 0.87 and 0.74 for blood and fin, respectively. The isotopic relationships between tissues in individual populations of European catfish and Northern pike exhibited better results on average than those observed in other fish species [17, 22, 39, 33, 58]. This could possibly be attributed to our study’s focus on adult individuals of two large fish species, as isotopic turnover rates generally decrease with body mass [30]. Thus, the δ13C and δ15N values between tissues may have been in better agreement in these larger adult fish. As such, the isotopic offsets between blood and muscle, as well as fin and muscle, can be used to formulate a tissue conversion model, allowing collected blood and fin tissue values to be converted to equivalent muscle tissue values for integration into food web analyses. However, establishing a species- or population-specific conversion model might require a muscle biopsy from a small number of fish, as the models for both species and both tissues varied to some extent among the study sites. A number of prior studies focusing on inter-tissue isotopic differences in other fish species also suggest that fractionation between tissues could differ among populations of the same species, potentially necessitating population-specific conversion equations [22, 39, 61, 62]. We did not observe any evident trend in the effect of body mass on the isotopic differences between tissues for either European catfish or Northern pike at any study sites. This lack of trend might indicate that inter-tissue isotope fractionations remain relatively consistent during the adult life-history stage. Regarding lipid content (approximated by C:N ratio), both species exhibited a clear pattern. The C:N ratio was consistently lowest in muscle, intermediate in fin, and highest in blood across all cases. Although significant differences in C:N ratio were evident in muscle and fin across study sites, the blood consistently displayed an unexpectedly stable C:N ratio, with no noticeable variation among study sites.

In conclusion, both δ13C and δ15N values for blood and fin serve as highly accurate predictors of respective isotope values for muscle tissue in adult apex predators like European catfish and Northern pike. Generally, blood proves to be a more suitable tissue than fin as it corresponds more closely to muscle tissue. Although blood sampling is time-consuming and requires technical proficiency, with practice, it is feasible to routinely draw blood from larger fish using minimally invasive methods [54]. However, the routine collection of fin clips stands out as a simpler and more commonly employed approach [17, 53, 55]. The variations in offsets among populations were notably different, not only for the studied European catfish and Northern pike but also for other species examined previously [22, 61, 62]. These differences likely stem from site-specific environmental conditions, population-related physiological factors, and mainly due to the varying diet composition. Thus, it is advisable to determine isotopic offsets between tissues from a small sample e.g. 10 individuals within each population of interest, in order to ascertain the necessary conversion correction.

Finally, our optimal guide for SIA, not only for valuable fish but ideally for any other species, suggests avoiding the need to kill individuals. Instead, if their size allows, it is preferable to conduct a muscle tissue biopsy in a necessary subsample of approximately 10 individuals from each studied locality. Concurrently, for these and all other individuals, opt for less harmful tissue sampling methods, such as blood or fin. Based on the 10 individuals, perform a regression analysis between muscle tissue and the less harmful tissues. Subsequently, these conversion-corrected data could be used. This process eliminates the need to harm or kill more individuals.

Supporting information

S1 File. Dataset.

The dataset analyzed during the present study.

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

(XLSX)

Acknowledgments

We thank Dr. Mikko Kiljunen for help with SIA at the University of Jyväskylä, and Dr. David Hochman, a professional English editor, for English editing. We also thank the two anonymous reviewers and the editor for helpful comments and suggestions.

References

  1. 1. Vander Zanden MJ, Casselman JM, Rasmussen JB. Stable isotope evidence for the food web consequences of species invasions in lakes. Nature. 1999; 401: 464–467.
  2. 2. Clarke LR, Vidergar DT, Bennett DH. Stable isotopes and gut content show diet overlap among native and introduced piscivores in a large oligotrophic lake. Ecol Freshw Fish. 2005; 14: 267–277.
  3. 3. Hayden B, Harrod C, Kahilainen KK. Lake morphometry and resource polymorphism determine niche segregation between cool- and cold-water-adapted fish. Ecology. 2014; 95: 538–552. pmid:24669746
  4. 4. McIntyre JK, Beauchamp DA, Mazur MM, Overman NC. Ontogenetic trophic interactions and benthopelagic coupling in Lake Washington: evidence from stable isotopes and diet analysis. Trans Am Fish. 2006; 135: 1312–1328.
  5. 5. Xu J, Zhang M. Primary consumers as bioindicator of nitrogen pollution in lake planktonic and benthic food webs. Ecol Indic. 2012; 14: 189–196.
  6. 6. Vašek M, Brabec M, Blabolil P, Čech M, Draštík V, Jůza T, et al. Fish scale stable isotopes as potential indicators of nutrient pollution: Exploring the response of roach (Rutilus rutilus) scale δ15N and δ13C to a gradient of land use disturbance. Sci Total Environ. 2023; 865: 161198.
  7. 7. Kidd KA, Bootsma HA, Hesslein RH, Muir DCG, Hecky RE. Biomagnification of DDT through the benthic and pelagic food webs of Lake Malawi, East Africa: importance of trophic level and carbon source. Environ Sci Technol. 2001; 35: 14–20. pmid:11351999
  8. 8. Semerád J, Horká P, Filipová A, Kukla J, Holubová K, Musilová Z, et al. The driving factors of per- and polyfluorinated alkyl substance (PFAS) accumulation in selected fish species: The influence of position in river continuum, fish feed composition, and pollutant properties. Sci Total Environ. 2022; 816: 151662. pmid:34780822
  9. 9. Vašek M, Souza AT, Říha M, Kubečka J, Znachor P, Hejzlar J. Stable isotope evidence from archived fish scales indicates carbon cycle changes over a four-decade history of the Římov Reservoir (Czechia). Sci Total Environ. 2021; 755: 142550.
  10. 10. Vejříková I, Vejřík L, Syväranta J, Kiljunen M, Čech M, Blabolil P, et. al. Distribution of herbivorous fish is frozen by low temperature. Sci Rep. 2016; 6: 39600. pmid:28004804
  11. 11. Vejříková I, Eloranta AP, Vejřík L, Šmejkal M, Čech M, Sajdlová Z, et al. Macrophytes shape trophic niche variation among generalist fishes. PlosONE. 2017; 12: e0177114. pmid:28486550
  12. 12. Vejřík L, Vejříková I, Blabolil P, Eloranta AP, Kočvara L, Peterka J, et al. European catfish (Silurus glanis) as a freshwater apex predator drives ecosystem via its diet adaptability. Sci Rep. 2017; 7: 15970. pmid:29162872
  13. 13. Vejřík L, Vejříková I, Blabolil P, Sajdlová Z, Kočvara L, Kolařík T, et al. Trophic position of the species and site trophic state affect diet niche and individual specialization: from apex predator to herbivore. Biology. 2023; 12: 1113. pmid:37626997
  14. 14. Boecklen WJ, Yarnes CT, Cook BA, James AC. On the use of stable isotopes in trophic ecology. Annu Rev Ecol Evol Syst. 2011; 42: 411–440.
  15. 15. Pinnegar JK, Polunin NVC. Differential fractionation of d13C and d15N among fish tissues: implications for the study of trophic interactions. Funct Ecol. 1999; 13: 225–231.
  16. 16. Tronquart NH, Mazeas L, Reuilly-Manenti L, Zahm A, Belliard J. Fish fins as non-lethal surrogates for muscle tissues in freshwater food web studies using stable isotopes. Rapid Commun Mass Spectrom. 2012; 26: 1603–1608. pmid:22693116
  17. 17. McIntosh LM, Reid MA. Fish fins as a non-lethal alternative to muscle tissue in stable isotope studies of food webs in an Australian river. Mar Freshw Res. 2021; 72: 838–847.
  18. 18. Trebilco R, Baum JK, Salomon AK, Dulvy NK. Ecosystem ecology: size-based constraints on the pyramids of life. Trends Ecol. 2013; 28: 423–431. pmid:23623003
  19. 19. Estes JA, Terborgh J, Brashares JS, Power ME, Berger J, Bond WJ, et al. Trophic downgrading of planet. Earth Sci. 2011; 333: 301–306.
  20. 20. Vejřík L, Vejříková I, Kočvara L, Blabolil P, Peterka J, Sajdlová Z, et al. The pros and cons of the invasive freshwater apex predator, European catfish Silurus glanis, and powerful angling technique for its population control. J Environ Manage. 2019; 241: 374–382.
  21. 21. Kelly MH, Hagar WG, Jardine TD, Cunjak RA. Nonlethal sampling of sunfish and slimy sculpin for stable isotope analysis: how scale and fin tissue compare with muscle tissue. N Am J Fish Manag. 2006; 26: 921–925.
  22. 22. Vašek M, Vejřík L, Vejříková I, Šmejkal M, Baran R, Muška M, et al. Development of non-lethal monitoring of stable isotopes in asp (Leuciscus aspius): a comparison of muscle, fin and scale tissues. Hydrobiologia. 2017; 785: 327–335.
  23. 23. Church M, Ebersole J, Rensmeyer K, Couture R, Barrows F, Noakes D. Mucus: a new tissue fraction for rapid determination of fish diet switching using stable isotope analysis. Can J Fish Aquat Sci. 2009; 66: 1–5.
  24. 24. Busst GMA, Britton JR. Tissue-specific turnover rates of the nitrogen stable isotope as functions of time and growth in a cyprinid fish. Hydrobiologia. 2017; 805: 49–60.
  25. 25. Nolan ET, Britton JR. Diet of invasive pikeperch Sander lucioperca: developing non-destructive tissue sampling for stable isotope analysis with comparisons to stomach contents analysis. Knowl Manag Aquat. 2018; 419: 49.
  26. 26. Wills AA, Kidd AR 3rd, Lepilina A, Poss KD. Fgfs control homeostatic regeneration in adult zebrafish fins. Development. 2008; 135: 3063–3070. pmid:18701543
  27. 27. German DP, Miles RD. Stable carbon and nitrogen incorporation in blood and fin tissue of the catfish Pterygoplichthys disjunctivus (Siluriformes, Loricariidae). Environ Biol Fishes. 2010; 89: 117–133.
  28. 28. Logan JM, Lutcavage ME. Stable isotope dynamics in elasmobranch fishes. Hydrobiologia. 2010; 644: 231–244.
  29. 29. Bugoni L, McGill RAR, Furness RW. Effects of preservation methods on stable isotope signatures in bird tissues. Rapid Commun Mass Spectrom. 2008; 22: 2457–2462. pmid:18642324
  30. 30. Vander Zanden MJ, Clayton MK, Moody EK, Solomon CT, Weidel BC. Stable isotope turnover and half-life in animal tissues: a literature synthesis. PLoS One. 2015; 10: e0116182. pmid:25635686
  31. 31. Henderson CJ, Stevens TF, Lee SY. Assessing the suitability of a non-lethal biopsy punch for sampling fish muscle tissue. Fish Physiol Biochem. 2016; 42: 1521–1526. pmid:27178186
  32. 32. George SE, Baker TR, Baker BB. Nonlethal detection of PFAS bioaccumulation and biomagnification within fishes in an urban- and wastewater-dominant Great Lakes watershed. Environ Pollut. 2023; 321: 121123. pmid:36681373
  33. 33. Hanisch JR, Tonn WM, Paszkowski CA, Scrimgeour GJ. δ13C and δ15N signatures in muscle and fin tissues: Nonlethal sampling methods for stable isotope analysis of salmonids. N Am J Fish Manag. 2010; 30: 1.
  34. 34. Pollard S, Anderson JC, Bah F, Mateus M, Sidhu M, Simmons DBD. Nonlethal blood sampling of fish in the lab and field with methods for dried blood plasma spot omic analyses. Front Genet. 2022; 13: 795348. pmid:35401689
  35. 35. Sanderson BL, Tran CD, Coe HJ, Pelekis V, Steel EA, Reichert WL. Nonlethal sampling of fish caudal fins yields valuable stable isotope data for threatened and endangered fishes. Trans Am Fish. 2009; 138: 1166–1177.
  36. 36. Fincel MJ, VanDeHey JA, Chipps SROV. Non-lethal sampling of walleye for stable isotope analysis: a comparison of three tissues. Fish Manag Ecol. 2012; 19: 283–292.
  37. 37. Robbins CT, Felicetti LA, Sponheimer M. The effect of dietary protein quality on nitrogen isotope discrimination in mammals and birds. Oecologia. 2005; 144: 534–540. pmid:15800751
  38. 38. Thomas SM, Crowther TW. Predicting rates of isotopic turnover across the animal kingdom: a synthesis of existing data. J Anim Ecol. 2015; 84: 861–870. pmid:25482029
  39. 39. Graham CT, Harrison SSC, Harrod C. Development of non-lethal sampling of carbon and nitrogen stable isotope ratios in salmonids: effects of lipid and inorganic components of fins. Isotopes Environ Health Stud. 2013; 49: 555–566. pmid:23937861
  40. 40. Ankjærø T, Christensen JT, Grønkjær P. Tissue-specific turnover rates and trophic enrichment of stable N and C isotopes in juvenile Atlantic cod Gadus morhua fed three different diets. Mar Ecol Prog Ser. 2012; 461: 197–209.
  41. 41. Buchheister A, Latour RJ. Turnover and fractionation of carbon and nitrogen stable isotopes in tissues of a migratory coastal predator, summer flounder (Paralichthys dentatus). Can J Fish Aquat Sci. 2010; 67: 445–461.
  42. 42. Říha M, Gjelland KØ, Děd V, Eloranta AP, Rabaneda-Bueno R, Baktoft H, et al. Contrasting structural complexity differentiate hunting strategy in an ambush apex predator. Sci Rep. 2021; 11: 17472. pmid:34471177
  43. 43. Brooks JL, Dodson SI. Predation, body size, and composition of plankton. Science. 1965; 150: 28–35. pmid:17829740
  44. 44. Vejříková I, Vejřík L, Čech M, Říha M, Peterka J. Succession of submerged vegetation in a hydrologically reclaimed opencast mine during first ten years. Restor Ecol. 2022; 30: e13489.
  45. 45. Bartoň D, Brabec M, Sajdlová Z, Souza AT, Duras J, Kortan D, et al. Hydropeaking causes spatial shifts in a reproducing rheophilic fish. Sci Total Environ. 2022; 806: 150649. pmid:34597557
  46. 46. Vejřík L, Vejříková I, Peterka J, Čech M. Methods for capturing catfish and potential regulation of catfish population. In: Advances in Animal Science and Zoology. Ed.: Jenkins OP, Nova Science Pub., 2019. ISBN: 978-1-53616-048-2, 135–142.
  47. 47. Vejřík, Vejříková I, Blabolil P, Bartoň D, Sajdlová Z, Kočvara L, et al. (submited) Long-lines for research monitoring or efficient population regulation of an invasive large-bodied apex predator, the European catfish (Silurus glanis). Heliyon.
  48. 48. Hoffman JC, Sierszen ME, Cotter AM. Fish tissue lipid-C:N relationships for correcting d13C values and estimating lipid content in aquatic food-web studies. Rapid Commun Mass Spectrom. 2015; 29: 2069–2077.
  49. 49. Stat-Soft, Inc., Tulsa, OK, USA, 2009.
  50. 50. Bates D, Maechler M, Bolker B, Walker S. Fitting Linear Mixed-Effects Models Using lme4. J Stat Soft. 2015; 67: 1–48.
  51. 51. Pinheiro J, Bates D, DebRoy S, Sarkar D. R Core Team_nlme: Linear and Nonlinear Mixed Effects Models_. R package version 3.1–153, 2021. <URL: https://CRAN.R-project.org/package=nlme>.
  52. 52. R Software. version 4.1.2 The R Foundation for Statistical Computing. 2021.
  53. 53. Cano-Rocabayera O, Maceda-Veiga A, de Sostoa A. Fish fins and scales as non-lethally sampled tissues for stable isotope analysis in five fish species of north–eastern Spain. Environ Biol Fish. 2015; 98: 925–932.
  54. 54. Lawrence MJ, Raby GD, Teffer AK, Jeffries KM, Danylchuk AJ, Eliason EJ, et al. Best practices for non-lethal blood sampling of fish via the caudal vasculature. J Fish Biol. 2020; 97: 4–15. pmid:32243570
  55. 55. Guy CS, Blankenship HL, Nielsen LA. Tagging and marking. In Murphy BR, Willis DW (eds), Fisheries Techniques, 2nd ed. American Fisheries Society, Bethesda 1996; 353–383.
  56. 56. Jardine TD, Hunt RJ, Pusey BJ, Bunn SE. A non-lethal sampling method for stable carbon and nitrogen isotope studies of tropical fishes. Mar Freshw Res. 2011; 62: 83–90.
  57. 57. Roberts KN, Lund T, Hayden B, Poesch M. Season and species influence stable isotope ratios between lethally and non-lethally sampled tissues in freshwater fish. J Fish Biol. 2022; 100: 229–241. pmid:34739138
  58. 58. Winter ER, Nyqvist M, Britton JR. Non-lethal sampling for stable isotope analysis of pike Esox lucius: how mucus, scale and fin tissue compare to muscle. J Fish Biol. 2019; 95: 956–958.
  59. 59. Syväranta J, Cucherousset J, Kopp D, Crivelli A, Céréghino R, Santoul F. Dietary breadth and trophic position of introduced European catfish Silurus glanis in the River Tarn (Garonne River basin), southwest France. Aquat Biol. 2010; 8: 137–144.
  60. 60. Fredrickson A, Rieucau G, Fontenot QC, Lackmann A, Solomon RD. Non-lethal fin clip model validation for stable isotope analysis of spotted and alligator gar. Trans Am Fish. 2021; 151: 72–80.
  61. 61. Willis TJ, Sweeting CJ, Bury SJ, Handley SJ, Brown JCS, Freeman DJ, et al. Matching and mismatching stable isotope (d13C and d15N) ratios in fin and muscle tissue among fish species: a critical review. Mar Biol. 2013; 160: 1633–1644.
  62. 62. Martins EG, Minke-Martin V, Patterson DA, Power M. The effects of inter- and intra-population variability on the use of adipose fin-muscle tissue comparisons for stable isotope studies of sockeye salmon. Fish Res. 2019; 211: 148–157.
  63. 63. Říha M, Rabaneda Bueno R, Vejřík L,Jarić I, Prchalová M,Draštík V, et al. (in prep.) Hungry catfish–how the availability of prey influences the activity pattern of an apex predator.
  64. 64. Lv W, Ju T, Dong B, Yu B, Yin J. Effect of dietary stable isotopic ratios of carbon and nitrogen on the extent of their incorporation into tissues of rats. J Animal Sci Biotechnol. 2012, 3: 14. pmid:22958501