Health Risk Assessment of Heavy Metals and Lipid Quality Indexes in Freshwater Fish from Lakes of Warmia and Mazury Region, Poland

The objectives of study were to determine heavy metals content (Zn, Cu, Mn, Fe and Hg) and fatty acids in selected organs of roach, Rutilus rutilus (L.); bream, Abramis brama (L.); pike, Esox lucius (L.); Eurasian perch, Perca fluviatilis (L.) collected from reservoirs of Warmia and Mazury region (northeastern Poland). Heavy metals were determined with atomic absorption spectrometry AAS. The fatty acids were analyzed using gas chromatography. In a few cases, differences in the content of heavy metals and fatty acids were not significant between species. The muscles of fish characterized significantly higher values of mercury than other organs (p ≤ 0.05), except for bream. The reverse regularity was observed in the case of content of Cu, Zn, Mn and Fe. Fatty acids having a desirable dietary effect in humans (DFA-Hypocholesterolaemic fatty acids) (74.00–74.84) were more than OFA (hypercholesterolaemic fatty acids), i.e., those undesirable (24.03–24.79). The lipid quality indexes AI (index of atherogenicity) (0.40–0.44) and TI (index of thrombogenicity) (0.18–0.24) in muscles of fish were low, which means that the meat of the fish may be recommended for human health. THQ (target hazard quotient) and HI (hazard index) as individual foodstuff were below 1, whereas HI for a specific receptor/pathway combination exceeded 1. This may suggest that eating meat from a given species is safe from a health point of view.


Introduction
A properly balanced diet, including the amount of fish consumed, is a way to maintain the health and balance of the human body, as well as physiological functions [1]. Fish and fish products play a major role from a nutritional point of view because they are a rich source of nutrients, and provide a relatively low caloric content [2]. Moreover, fish flesh has a good digestibility [3]. According to Pal et al. [2], fish also provide a good balance of protein, and their availability and affordability is better than other animal protein sources [4]. It is also known that fish and fish oils (along with some vegetable oils) are a good source for essential polyunsaturated fatty acids in human nutrition [3,5]. Many clinical trials have shown the relationship between the intake of n-3 PUFA (polyunsaturated fatty acids) and beneficial effects in different diseases (optimal cardiovascular system, brain and vision functioning, cancer prevention, arthritis, hypertension and diabetes mellitus) [2,[5][6][7][8]. Since both linoleic acid (LA) and α-linolenic acid (ALA) cannot be synthesized in the body, they are therefore essential [9][10][11]. LA can be converted via γ-linolenic acid (18:3n-6) and dihomo-γ-linolenic acid (20:3n-6; DGLA) to arachidonic acid (20:4n-6; AA). ALA is the precursor for synthesis of eicosapentaenoic acid (20:5n-3; EPA) and further to docosapentaenoic acid (22:5n-3; DPAn-3) and docosahexaenoic acid (22:6n-3; DHA) [12]. However, among fatty acids, there are some saturated fatty acids (SFA) that have different

Mercury
Duplicate samples of organs were weighed into quartz boats (<270 mg ± 0.0001 g). The total mercury was processed with atomic absorption spectrometry thermal decomposition using Milestone DMA-80 with dual-cell and UV enhanced photodiodes (Italy). Conditions for the determination of mercury is described in an earlier publication [24].

Zinc, Copper, Iron and Manganese
For the determination of zinc, copper, iron and manganese in muscles of fish examined, samples were ashed at 450 °C using laboratory furnaces (Nabertherm, Germany). After obtaining the white ash, it was dissolved in 1 M HNO3 (Suprapur-Merck, Darmstadt, Germany) and transferred with deionized water (Merck-Millipore Elix Advantage 3, USA) into a volumetric flask of volume 25 mL. After weighing, samples of liver and gills were placed in a boro-silica glass tube and mineralized with a nitric and perchloric acid mixture (Merck, Darmstadt, Germany, 3:1, v/v) at 190 °C. For this purpose, a heating block with a programmable temperature and digestion time was used (DK 20, VELP Scientifica, Italy). The resulting solution was transferred to flasks with a volume of 25 mL, using deionized water.

Instrumental Analysis and Quality Control
Samples were prepared in two parallel wells. The total mercury was determined using Milestone DMA-80 with dual-cell and UV enhanced photodiodes (Italy). The detection limit (LOD) was 0.02 μg/kg. The quality control of methods was tested using the reference material: BCR CRM Ethical permit: fish were bought at the fish farm and were already dead. According to European and Polish Law, the research done on the tissue of commercially caught fish is free from obtaining permission from the Local Ethical Commission.

Mercury
Duplicate samples of organs were weighed into quartz boats (<270 mg ± 0.0001 g). The total mercury was processed with atomic absorption spectrometry thermal decomposition using Milestone DMA-80 with dual-cell and UV enhanced photodiodes (Italy). Conditions for the determination of mercury is described in an earlier publication [24].

Zinc, Copper, Iron and Manganese
For the determination of zinc, copper, iron and manganese in muscles of fish examined, samples were ashed at 450 • C using laboratory furnaces (Nabertherm, Germany). After obtaining the white ash, it was dissolved in 1 M HNO 3 (Suprapur-Merck, Darmstadt, Germany) and transferred with deionized water (Merck-Millipore Elix Advantage 3, USA) into a volumetric flask of volume 25 mL. After weighing, samples of liver and gills were placed in a boro-silica glass tube and mineralized with a nitric and perchloric acid mixture (Merck, Darmstadt, Germany, 3:1, v/v) at 190 • C. For this purpose, a heating block with a programmable temperature and digestion time was used (DK 20, VELP Scientifica, Italy). The resulting solution was transferred to flasks with a volume of 25 mL, using deionized water.

Instrumental Analysis and Quality Control
Samples were prepared in two parallel wells. The total mercury was determined using Milestone DMA-80 with dual-cell and UV enhanced photodiodes (Italy). The detection limit (LOD) was 0.02 µg/kg. The quality control of methods was tested using the reference material: BCR CRM 422 (muscles of cod, Gadus morhua (L.)) with a certified value of mercury. The recovery rate of Hg was 100.2% (n = 4).
Concentrations of zinc, copper, iron and manganese were measured using the flame atomic absorption spectrometry (iCE 3000 Series AAS, Thermo Scientific, Loughborough, England). Four blanks and four standards were analyzed with each batch of samples. The absorptions wavelength was as follows: 213.9 nm for zinc and 324.8 nm for copper. The calibration curves were prepared using four solution standards (1000 µg/L) with 0.1 M HNO 3 supplied by J.T.Baker ® (Deventer, Netherlands). The calibration curves were linear within the range of metal concentration (regression coefficients R 2 ≥ 0.999). The absorption wavelengths were as follows: 248.3 nm for iron, 213.9 nm for zinc, 324.8 nm for copper, 279.5 nm for manganese. The detection limits (LOD) were 0.5 mg/kg for Fe, 0.1 mg/kg for Zn, 0.05 mg/kg for Cu and 0.05 mg/kg for Mn. Sensitivity was as follows: 0.05, 0.05, 0.02 and 0.02 mg/L. The quality control of methods was tested using the elements in reference material: BCR CRM 422 (muscles of cod Gadus morhua (L.)) with a certified value of zinc, copper and mercury. The recovery rates were: 105.0% Zn, 103.0% Cu. 96% Fe and 103% Mn, respectively.

Fatty Acids Analysis
The lipids were extracted according to the Folch's procedure [25]. The fatty acid methyl esters were prepared from total lipids with the Peisker method with chloroform, methanol and sulphuric acid (100:100:1 v/v) [26].
The fatty acids of methyl esters of each sample were analyzed using 7890A Agilent Technologies chromatograph with a flame-ionization detector (FID) (Agilent Technologies, INC., Santa Clara, California, USA) under the following conditions: capillary column (dimension 30 m × 0.25 µm with a 0.32 mm internal diameter, liquid phase StabilwaxR), temperature: flame-ionization detector −250 • C, injector −230 • C, column −190 • C, carrier gas-helium with a flow rate 1.5 mL/min. Individual fatty acids were identified by comparing the relative retention time peaks to the known Supelco's standards.

The Lipid Quality Indexes (AI and TI)
They were calculated using the following pattern by Ulbricht and Southgate [27] and Garaffo et al. [28]: Flesh-lipid quality (FLQ) were calculated presented by Abrami et al. [29] and Senso et al. [30]: Hypercholesterolaemic fatty acids (OFA): Hypocholesterolaemic fatty acids (DFA): EDI is the estimated daily intake (µg/kg body weight/day), and is obtained as follows: where C is the average concentration of mercury in food stuffs (µg/g wet weight), IR is the daily ingestion rate (g/daily) and BW is the average body weight (60 kg) [31].

Target Hazard Quotient (THQ)
According to Ahmed et al. [32] and US EPA [33], THQ estimated the non-carcinogenic health risk of consumers due to intake of heavy metal contaminated fish use an oral reference dose of Hg, Zn, Cu, Fe and Mn (RfD = 3.00 × 10 −4 , 3.00 × 10 −1 , 7.00 × 10 −1 , 1.4 × 10 −1 ). When THQ < 1, there is health benefit from fish consumption and the consumers are safe, whereas THQ > 1 suggests a high probability of adverse risk of human health.
where Efr is the exposure frequency (365 days/year), ED is the exposure duration (70 years), FiR is the fish ingestion rate (g/person/day), C is the average concentration of mercury in food stuffs (µg/g wet weight), RfD is the oral reference dose (mg/kg/day) (US EPA 2017), BW is the average body weight of local residents (60 kg) [31], TA is the average exposure time (365 days/year × ED).

The Combined Risk of Many Heavy Metals
The TTHQ of heavy metals for individual foodstuff was treated as the mathematical sum of each individual metal THQ value [34]: HI for a specific receptor/pathway combination (e.g., diet) was calculated using the following pattern [34]: When the HI exceeds 1, there may be concern for potential health risks.

Statistical Analysis
Significant differences in the content of fatty acids and lipid quality indexes in the muscle tissue of fish studied were estimated using a one-way analysis of variance (ANOVA). Similar, differences in the concentration of heavy metals, between species and organs of the same species were calculated using STATISTICA 13.1 program (StatSoft, Kraków, Poland). Differences were found to be significant at p < 0.05. When the Bartlett's test showed that the variances were heterogeneous, values in particular groups were transformed (log x). The correlation coefficients between the content of heavy metals in organs of fish and their size (body weight and total length) were calculated using a STATISTICA 13.1 software. The significance level of p < 0.05 was used.

Differences in the Content of Heavy Metals
The muscles of fish examined had significantly higher values of mercury than other organs (p ≤ 0.05), except for bream, which had no statistically significant differences between the content of mercury in the muscles and liver of this species (p > 0.05) ( Table 1). Reverse regularity was observed for the content of Cu, Zn, Mn and Fe. The liver of the examined fish contained more Cu than other organs (p ≤ 0.05). The liver of roach and perch had a significantly higher content of Fe, whereas bream liver had more Zn (p ≤ 0.05). There were no significant differences between the content of Fe in the liver and gills of two species (bream and pike) and between the concentration of Zn in the liver and gills of perch. In other cases, higher contents of marked metals were observed in the gills of the examined fish (p ≤ 0.05). n-number of fish; SD-standard deviation; a, b, c -significant difference (p ≤ 0.05). The same letter indicates the absence of significant differences between organs of the same fish studied.
In most cases, it was not found that the feeding type of fish had a significant effect on the content of such elements as Cu, Zn, Fe and Mn ( Figure 2). Predatory fish (perch and pike) which occupy the last link in the trophic chain of the aquatic ecosystem contained more mercury in the muscles and liver than those of the studied fish belonging to the lower trophic levels (bream and roach) (p ≤ 0.05). This regularity was not found in gills of the examined fish (p > 0.05). Generally, the concentration of 18 elements, including Hg, Fe, Cu, Mn and Zn significantly differed between species such as sichel (Pelecus cultratus), ruffe (Gymnocephalus cernua) and European perch (Perca fluviatilis), caught at a polluted segment of the Danube River near Belgrade (p < 0.0001) [35]. Similarly, Sayegh Petkovšek et al. [36] found that fish (Abramis brama danubii, Alburnus alburnus alburnus, Barbus meridionalis petenyi, Carassius auratius gibelio, Cyprinus carpio, Lepomis gibossus, Leuciscius cephalus cephalus, Perca fluviatilis, Rutilus rutilus and Scardinus Generally, the concentration of 18 elements, including Hg, Fe, Cu, Mn and Zn significantly differed between species such as sichel (Pelecus cultratus), ruffe (Gymnocephalus cernua) and European perch (Perca fluviatilis), caught at a polluted segment of the Danube River near Belgrade (p < 0.0001) [35]. Similarly, Sayegh Petkovšek et al. [36] found that fish (Abramis brama danubii, Alburnus alburnus alburnus, Barbus meridionalis petenyi, Carassius auratius gibelio, Cyprinus carpio, Lepomis gibossus, Leuciscius cephalus cephalus, Perca fluviatilis, Rutilus rutilus and Scardinus erythrophtlalmus erythrophtlalmus) collected in the Šalek lakes (Slovenia) could be used to assess the bioavailability of metals in freshwater biota, although the concentration of metals (e.g., Hg and Zn) varied significantly between fish species. To compare the current results with the research of other authors, Table 2 presents the content of heavy metals of selected fish organs from various aquatic ecosystems located in Europe [35][36][37][38][39][40][41][42][43][44]. Based on these data, it was found that in most cases, these results did not coincide with the results covered by the study of this work. First of all, the mercury content does not accumulate along the trophic chain, except for Hg in the muscles of fish caught from the Kirchera River [37]. In addition, the fish muscles (Table 2) did not contain more mercury than the other examined organs presented in the table. The authors also did not find clear relationships between the content of Cu, Mn, Zn and Fe and the length of the trophic chain of a given water reservoir. Some relationships, where the amount of Zn, Cu and Mn was higher in fish organs inhabiting lower levels of the chain, were observed in fish from the West Morava River Basin [38].  References [35,38]-results expressed in mg/kg dry weight.  References [35], [38]-results expressed in mg/kg dry weight.     Significant positive correlations were found for Zn-Fe, Zn-Mn, while significant negative correlations were observed for Zn-Cu in the muscles of fish (Scardinus knezevici, Alburnus scoranza, Cyprinus carpio, Rutilus prespensis, Anguilla anguilla and Perca fluviatilis) collected from Skadar Lake (Montenegro) [39]. The results obtained by these authors are consistent with the results of the current study, but only for roach muscles. However, the current results are in contradiction with previous studies received by Łuczyńska and Tońska [45]. Figure 4 and 5 only shows the significant correlation between content of metals and size of fish. In other cases, no significant correlations were observed. The content of Hg in muscles of all fish species was positively correlated with body weight and total length (p < 0.018) (Figures 4 and  5). Only the concentration of Hg in gills and liver of perch increased along with increased body weight (p = 0.000 and p = 0.004) ( Figure 4) and length (p = 0.008 and p = 0.002), respectively ( Figure  5). Significant positive correlation were also found between the content of Cu in the muscle tissue of perch and body weight (p = 0.001) and length (p = 0.010), whereas the content of Cu in gills of the same species decreased as body weight (p = 0.05) and total length (p = 0.003) increased. Significant positive correlations between body weight (p = 0.014) ( Figure 4) and total length (p = 0.020) ( Figure  5) and the content of Zn were found in muscles of pike. There is a negative correlation between the content of Zn, Mn and Fe in muscles of roach and total length (p = 0.037, p = 0.031, p = 0.011, respectively) and between the level of Fe and their body weight (p = 0.043) (Figure 4). Negative correlation factors between Cu content in gills of roach and body weight (p = 0.048) and length (p = 0.019) were noted. In the case of bream, the content of zinc in gills decreased as body weight increased (p = 0.008) (Figure 4). Similarly, Zn in gills and Fe in liver of bream decreased as total length increased (p = 0.017, p = 0.027, respectively). The obtained results could have been influenced by the variety of conditions of the aquatic environment in which the fish lived, related to the water reaction and the availability of metals. Significant positive correlations were found for Zn-Fe, Zn-Mn, while significant negative correlations were observed for Zn-Cu in the muscles of fish (Scardinus knezevici, Alburnus scoranza, Cyprinus carpio, Rutilus prespensis, Anguilla anguilla and Perca fluviatilis) collected from Skadar Lake (Montenegro) [39]. The results obtained by these authors are consistent with the results of the current study, but only for roach muscles. However, the current results are in contradiction with previous studies received by Łuczyńska and Tońska [45].

Correlation between Metal Pairs and Size
Figures 4 and 5 only shows the significant correlation between content of metals and size of fish. In other cases, no significant correlations were observed. The content of Hg in muscles of all fish species was positively correlated with body weight and total length (p < 0.018) (Figures 4 and 5). Only the concentration of Hg in gills and liver of perch increased along with increased body weight (p = 0.000 and p = 0.004) ( Figure 4) and length (p = 0.008 and p = 0.002), respectively ( Figure 5). Significant positive correlation were also found between the content of Cu in the muscle tissue of perch and body weight (p = 0.001) and length (p = 0.010), whereas the content of Cu in gills of the same species decreased as body weight (p = 0.05) and total length (p = 0.003) increased. Significant positive correlations between body weight (p = 0.014) ( Figure 4) and total length (p = 0.020) ( Figure 5) and the content of Zn were found in muscles of pike. There is a negative correlation between the content of Zn, Mn and Fe in muscles of roach and total length (p = 0.037, p = 0.031, p = 0.011, respectively) and between the level of Fe and their body weight (p = 0.043) (Figure 4). Negative correlation factors between Cu content in gills of roach and body weight (p = 0.048) and length (p = 0.019) were noted. In the case of bream, the content of zinc in gills decreased as body weight increased (p = 0.008) (Figure 4). Similarly, Zn in gills and Fe in liver of bream decreased as total length increased (p = 0.017, p = 0.027, respectively). The obtained results could have been influenced by the variety of conditions of the aquatic environment in which the fish lived, related to the water reaction and the availability of metals.     There were no significant correlations between Fe and Hg in muscles of perch (Skandar Lake, Montenegro) and their age or size (length and weight) [39]. However, positive correlation between Hg in muscles of pike, perch, roach, bream and their size were found in previous studies [24,46,47]. Nozari et al. [48] found a positive correlation between these parameters in muscles of pike from Anzali International Wetland (Iran), while Miller et al. [49] reported a similar correlation in muscles of perch from Swedish and Finnish aquatic environments. This is accordance with the current results (Figures 4  and 5). According to Łuczyńska et al. [47] zinc in the gills of roach was negatively correlated with their weight (r = −0.693, p = 0.026) and length (r = −0.668, p = 0.035). During the study, the correlation was also negative, but not statistically significant (Figures 4 and 5). The same authors found a positive correlation between the Zn content in gills of perch and their weight (r = 0.634, p = 0.049) and no statistically significant correlation between the copper content in the roach and perch organs and their size [47]. This is in contradiction with the results being studied in the current study (Figures 4 and 5).

Fatty Acids and Lipid Quality Indexes
Among the species under study, pike had a significantly higher n-3/n-6 ratio (3.56) and Σ n-3 PUFA (42.61%) than the other fish species (p ≤ 0.05) ( Table 3). There were no significant differences between the values of Σ SFA (31.73-32.94%) and Σ n-6 PUFA (12.2-14.98%) in muscles of the examined fish. Σ MUFA (24.12%) in muscles of bream was significantly higher (p ≤ 0.05) than other fish species, with the exception of roach. Among the SFA and MUFA, the predominant fatty acids were palmitic acid C16:0 and oleic acid C18:1, whereas arachidonic C20:4 n-6, eicosapentaenoic EPA n-3 and docosahexaenoic DHA n-3 acids were the most abundant of PUFA.
Atherogenic index (AI), hypocholesterolemic (OFA) and hypercholesterolemic fatty acids (DFA) in muscles of the examined fish did not differ significantly (p > 0.05) ( Figure 6). The flesh-lipid quality index (FLQ) in the muscles of fish included in this study gave rise to the following sequence: pike > perch ≈ roach > bream (p ≤ 0.05). No difference was found in terms of thrombogenicity index (TI) (Figure 6) between the species within the following groups (bream, roach and perch) and (roach, pike and perch) (p > 0.05). The low values of AI and TI indicates that the tissue of the studied fish is beneficial from a health point of view. Fatty acids with undesirable dietary effect in humans (OFA) were less than DFA, i.e., those desirable. As with heavy metals, there was no clear relationship between the tested fatty acids or selected quality indicators and the fish species or their occurrence in organs and biometric data.  There were no significant correlations between Fe and Hg in muscles of perch (Skandar Lake, Montenegro) and their age or size (length and weight) [39]. However, positive correlation between Hg in muscles of pike, perch, roach, bream and their size were found in previous studies [24,46,47]. Nozari et al. [48] found a positive correlation between these parameters in muscles of pike from Anzali International Wetland (Iran), while Miller et al. [49] reported a similar correlation in muscles of perch from Swedish and Finnish aquatic environments. This is accordance with the current results (Figures 4 and 5). According to Łuczyńska et al. [47] zinc in the gills of roach was negatively correlated with their weight (r = −0.693, p = 0.026) and length (r = −0.668, p = 0.035). During the study, the correlation was also negative, but not statistically significant (Figures 4 and 5). The same authors found a positive correlation between the Zn content in gills of perch and their weight (r = 0.634, p = 0.049) and no statistically significant correlation between the copper content in the roach and perch organs and their size [47]. This is in contradiction with the results being studied in the current study (Figure 4 and 5).

Fatty Acids and Lipid Quality Indexes
Among the species under study, pike had a significantly higher n-3/n-6 ratio (3.56) and Σ n-3 PUFA (42.61%) than the other fish species (p ≤ 0.05) ( Table 3). There were no significant differences between the values of Σ SFA (31.73-32.94%) and Σ n-6 PUFA (12.2-14.98%) in muscles of the examined fish. Σ MUFA (24.12%) in muscles of bream was significantly higher (p ≤ 0.05) than other fish species, with the exception of roach. Among the SFA and MUFA, the predominant fatty acids were palmitic acid C16:0 and oleic acid C18:1, whereas arachidonic C20:4 n-6, eicosapentaenoic EPA n-3 and docosahexaenoic DHA n-3 acids were the most abundant of PUFA.
Atherogenic index (AI), hypocholesterolemic (OFA) and hypercholesterolemic fatty acids (DFA) in muscles of the examined fish did not differ significantly (p > 0.05) ( Figure 6). The fleshlipid quality index (FLQ) in the muscles of fish included in this study gave rise to the following sequence: pike > perch ≈ roach > bream (p ≤ 0.05). No difference was found in terms of thrombogenicity index (TI) (Figure 6) between the species within the following groups (bream, roach and perch) and (roach, pike and perch) (p > 0.05). The low values of AI and TI indicates that the tissue of the studied fish is beneficial from a health point of view. Fatty acids with undesirable dietary effect in humans (OFA) were less than DFA, i.e., those desirable. As with heavy metals, there was no clear relationship between the tested fatty acids or selected quality indicators and the fish species or their occurrence in organs and biometric data.   Many studies indicate that the fatty acids that dominated the muscles of the fish in the current study are the same groups of fatty acids found in the muscles of the species of fish in other studies [50][51][52][53][54][55]. Kainz et al. [56] found that n-3 and n-6 PUFA in fish of pre-alpine Lake Lunz (Austria) decreased with increasing trophic position, demonstrating that these essential fatty acids did not biomagnify with increasing trophic level. These observations are not consistent with the results of the current study (Table 3). Previous studies conducted by Łuczyńska et al. [51] also found that n-3 and n-6 PUFA did not differ between predatory and non-predatory fish. Based on the results, it can be stated that the values of these indices were not affected by the level of the trophic chain and, thus, the diet of the fish. Łuczyńska et al. [57], did not find any significant differences either between these indices for fish such as bream and perch, which occupy two different links in the aquatic ecosystem food chain. The results of the current research on AI and TI ( Figure 6) are consistent with studies conducted by Linhartová et al. [58] which found that the indicators in all analyzed fish, except for Nile tilapia, were below 0.5, which indicates significant benefit to human health if these fish are included in the human diet. AI and TI for seven freshwater fish species, including perch and bream, from the Czech Republic were close to the values given for the Eskimo diet, which indicates the high nutritional value of these fish and, therefore, not a threat [59]. According to Linhartová et al. [58] the higher values of these coefficients, the higher risk of developing cardiovascular diseases, because AI indicates the risk of diseases such as atherosclerosis (deposition of fat in the walls of the arteries) and TI determines the possibility of blood clots. Tilami et al. [59] showed that AI and TI for bream and perch were 0.30 and 0.38 or 0.25 and 0.22, respectively. Similar values were found in these species of fish subject to this study ( Figure 6).

Human Health Risk
Fish consumption in 2013 year (Poland) was 12.2 kg per capita/year (for an adult with a body weight of 60 kg) [60]. Therefore, the estimated daily intake (EDI) of Cu, Zn, Mn, Fe and Hg from the 33.42 g/person/day were: 0.076-0.117, 2.330-5.582, 0.027-0.068, 0.489-0.729 and 0.036-0.193, respectively ( Table 4). The THQ and HI as individual foodstuff (TTHQ) were below 1 (Table 5), whereas HI for a specific receptor/pathway combination (e.g., diet) exceeded 1. This may suggest that eating meat from a given species is a non-carcinogenic health risk for consumers.  Łuczyńska et al. [57], studied European perch and bream from the Polish market, and found that the THQ values were also below 1 and there was no carcinogenic health risk to the consumer by consuming fish. According to these authors the daily intake (EDI) of mercury was: 0.009 µg/body weight (bream) and 0.076 (perch). EDI of mercury from the portion of fish examined was higher (Table 4). THQ and HI of Cu, Zn and Hg in European perch and roach from Lake Pluszne were lower than 1, which is important for the local community and other people who use the lake for recreational purposes [42], although the EDI of zinc and copper from the 33.698 g portions of fish was higher and the EDI of mercury was lower than the examined fish (Table 4). According to Addo-Bediako et al. [61], when THQ < 1, these adverse health effects are unlikely, whereas THQ > 1 suggests a high probability of adverse health effects. A THQ value (for Mn and Zn) below 1 in the roach muscles from the Miankaleh international wetland was also obtained by the Alipour et al. [42]. There was also no threat from metals, including Cu, Fe, Mn and Zn, in the consumption of fish living in lakes, rivers, streams, and the sea in Sakarya (Turkey) [62]. This is accordance with the results obtained in the current study. On the basis of the 18 elements under study, the authors observed that both bream and zander are among the best fish species to monitor the quality of fish tissue, but fisheries monitoring programs should also include other economically important fish species (e.g., common carp, pike and barbell) and some exotic species, such as Prussian carp [63].

Conclusions
Based on the hazard factor for individual metals THQ < 1, it was found that the consumption of the studied fish does not constitute a carcinogenic health risk. Similarly, when the mathematical sum of each THQ value for a given species of fish is considered, their consumption does not raise concerns related to adverse health risks. However, when a specific receptor/pathway combination (e.g., diet) was calculated, which posed a threat from the consumption of all fish species simultaneously contaminated with test elements, there could be concerns about potential health risk because HI was below 1. Looking at the low daily intake of fish, it appears that these fish are safe from a nutritional and health point of view. It was also observed that these fish are a rich source of fatty acids and the proportions between fatty acids with desirable dietary effects (DFA) and those with undesirable effects (OFA) were in favor of the former. In addition, AI and TI lipid quality indices in fish muscles were low, which means that fish tissue is beneficial to consumers. There was no clear relationship between the tested heavy metals, fatty acids or selected quality indicators and the fish species or their occurrence in organs and biometric data. Such a relationship could be observed only in relation to mercury, which accumulated in fish occupying individual levels of the trophic chain, with humans at the end. However, its levels did not exceed acceptable standards.
Author Contributions: J.Ł. had the original idea for the study and wrote the manuscript and conducted the statistical analysis; B.P. and J.Ł. collected and analyzed the data; B.P. contributed in the writing of the section "Analytical methods". All authors read and approved the final manuscript.

Funding:
The study was financed by the statutory theme number 17.610.008-300 and the project was financially supported by the Minister of Science and Higher Education in the range of the program entitled "Regional Initiative of Excellence" for the years 2019-2022, project number 010/RID/2018/19, and the amount of funding was 12.000.000 PLN.

Conflicts of Interest:
The authors declare that there is no conflict of interests.