Metals Bioaccumulation in 15 Commonly Consumed Fishes from the Lower Meghna River and Adjacent Areas of Bangladesh and Associated Human Health Hazards

The lower Meghna River, the easternmost part of the Ganges Delta, faces severe anthropogenic perturbations as it receives a huge discharge and industrial effluents. To measure the metal concentrations and human health hazards, edible tissues of 15 commercially important fish species were collected from the local fish markets and the lower Meghna River, Bangladesh. Trace and heavy metals such as Pb, Cr, Cu, Zn, Mn, Fe, Hg, Ni, Ca, Co, Se, Rb, Sr, and As were detected using the Energy Dispersive X-ray Fluorescence (EDXRF) method. The hierarchy of mean metal concentrations obtained was: Fe (162.198 mg/kg) > Zn (113.326 mg/kg) > Ca (87.828 mg/kg) > Sr (75.139 mg/kg) > Cu (36.438 mg/kg) > Se (9.087 mg/kg) > Cr (7.336 mg/kg) > Mn (6.637 mg/kg) > Co (3.474 mg/kg) > Rb (1.912 mg/kg) > Hg (1.657 mg/kg) > Ni (1.467 mg/kg) > Pb (0.521 mg/kg) > As (BDL). Based on the metal concentration obtained, the carnivorous species contained more metals than omnivores and herbivores. Similarly, the euryhaline and benthic feeder fishes had more metals than the stenohalines and demersal fishes. The metal pollution index (MPI) suggested that the highly consumed fish species Tilapia (Oreochromis mossambicus) and Rui (Labeo rohita) accumulated higher metals than other fishes. Both the Targeted Hazard Quotient (THQ) and Hazard Index (HI) values for adult and child consumers were <1, indicating that consumers would not experience the non-carcinogenic health effects. Although children were more susceptible than adults, carcinogenic risk (CR) exposure of Cr for all the consumers was found in the acceptable range (10−6 to 10−4), but the CR exposure of Pb was negligible for all the consumers. The correlation, principal component analysis (PCA), and cluster analysis were conducted to identify the sources of metals identified from the fish tissue. The results indicated that the probable sources of the pollutants were anthropogenic, arising from agricultural activities, electroplating materials, and lubricants used near the study area. However, the present study showed a different metal concentration in the samples at different levels but within the threshold levels non-carcinogenic and carcinogenic health risks; hence, the fishes of the area, in general, are safe for human consumption.


Introduction
In recent decades, extensive population growth and subsequent economic development have aroused a global consciousness about heavy metals contamination owing to their metals in the most consumed fishes from Noakhali Fish Market and their potential health risks through their consumption. In addition, in the coastal rural area of Bangladesh, the general people and authorities are unaware of the health implication of heavy metals accumulation in fish. Therefore, this study aims to determine the heavy metal concentration in fish from the Noakhali fish market and evaluate the risk of heavy metals to human health, which will provide insights to local people and management authorities to take necessary steps regarding this issue.

Sample Collection, Preparation and Analysis
A total of 30 fish samples of 15 species (two samples per species) from different feeding habitats (Table 1) were collected from different fish markets of Noakhali and the lower Meghna River to evaluate heavy metal levels ( Figure 1). These species were selected because these are the mostly consumed throughout the country, available at almost all the fish markets and commercially important. After collection, all the fish samples were kept frozen (-20 • C) by storing them in an icebox and transported to the laboratory as soon as possible. After measuring the physical parameters (weight and length), each fish sample was cleaned and washed with deionized water. Then, it was chopped with a stainless-steel knife cleaned with acetone and hot distilled water before use. Fourteen trace and heavy metals (Pb, Cr, Cu, Zn, Mn, Fe, Hg, Ni, Ca, Co, Se, Rb, Sr, and As) were analyzed from the dissected fish muscle tissue. For analyzing the metals, fish flesh was taken into a beaker and kept in a muffle furnace at 300 • C for 3 h to make ash. The ash samples were ground for making powder using carbide mortar and pestle. For EDXRF analysis, each powdered piece was pressed into a pellet of 2.5 cm diameter with a hydraulic press pellet maker (Specac) using 7 (seven) tons pressure. The irradiation of all actual samples was performed by assigning a time-based program controlled by a software package provided with the EDXRF system. The standard materials were also irradiated under similar experimental conditions to construct the calibration curves for quantitative elemental determination in the respective samples. The generated X-ray spectra of the materials were stored in the computer.

Analytical Quality Control
Before using, all the glassware and plastics were washed in nitric acid solution for 15 min and then rinsed with deionized water. Analytical grade reagents (Merck; Germany) and de-ionized water were used for the analysis throughout the study. For the validation and accuracy, the analytical method certified reference material (CRM 320, Merck KGaA, Darmstadt, Germany) was used. The recovery rate of the of the selected metals in the applied analytical method ranged from 72% to 105%. Besides, during the analytical process, the contamination's influence was absent and the relative standard deviation (RSD) was ≤10% for all tests. To assess the metal pollution, the metal pollution index (MPI) was adopted as follows [33,34]: where CM 1 is the concentration of the first concerning metal, CM 2 is the concentration of the second concerning metal, CM 3 is the concentration of the third concerning metal, CM n is the concentration of the nth metal (mg/kg dry wt) in the tissue sample of a particular species.

Estimated Daily Intake (EDI)
Estimated daily intake (EDI) was calculated by the following equation [35,36]: where Cn is the concentration of metal in the selected fish muscles tissue (mg/kg dry wt); IGr is the acceptable ingestion rate, which is 55.5 g/day for adults and 52.5 g/day for children [37,38]; Bwt is the bodyweight: 70 kg for adults and 15 kg for children [37].

Target Hazard Quotient (THQ) for Non-Carcinogenic Risk Assessment
THQ was estimated by the ratio of EDI and oral reference dose (RfD). The ratio value <1 implies non-significant risk effects [39]. The THQ formula is expressed as follows [40,41].
where Ed is exposure duration (65 years) [37]; Ep is exposure frequency (365 days/year) [31]; At is the average time for the non-carcinogenic element (Ed × Ep).

Hazard Index (HI)
Hazard index (HI) was calculated for the multiple elements (Hg, As, Mn, and Cr) found in the fish samples, and the equation is as follows [8].
where THQ are the estimated risk value for individual metal [8]. When the HI value is higher than 10, the non-carcinogenic risk effect is considered high for exposed consumers [42,43].

Carcinogenic Risk (CR)
To assess the probability of developing cancer over a lifetime, the carcinogenic risk is evaluated for the consequence of exposure to the substantial carcinogens [44,45]. The acceptable range of the risk limit is 10 −6 to 10 −4 [46,47]. CRs higher than 10 −4 are likely to increase the probability of carcinogenic risk effect [48,49]. The established equation to assess the CR is as follows [38,41].
where CSF is the oral slope factor of a particular carcinogen (mg/kg-day) [47].

Statistical Analysis
The correlation matrix (CM), principal component analysis (PCA), and Hierarchical cluster analysis were performed through PAST (version 3). CM and PCA help to determine the correlation between heavy metals in fish tissue [50]. Hierarchical cluster analysis is one of the most widely used hierarchical algorithms, which results in clusters in which variables or individuals are added in sequence considering the hierarchy of the cluster [35]. Clustering of metals concentrations in fish muscles based on Bray-Curtis similarity was performed to plot elements in a separate cluster, thus differentiating the samples' contamination status [51].

The Concentration of Heavy Metals in Fish Muscle Tissue
Heavy metal contamination in fish is one of the severe threats to humans and aquatic animals. Determination of heavy metal concentration is the first step to evaluating the extent of pollution in fish. The concentration (wet weight) of selected metals was in the following descending order: Fe (162.198) > Zn (113.326) > Cu (36.438) > Cr (7.336) > Mn (6.637) > Hg (1.657) > Ni (1.467) > Pb (0.521) ( Table 2). Among the examined fishes, the average metal concentration was maximum in Tilapia and maintained the following descending trend: Tilapia > Rui > Gulia > Poa > Loitta > Chiring > Ricksha > Bata > Catla > Grass Carp > Kalibaus > Koral > Pabda > Koi > Big head carp. However, considering the feeding guild of the sampled fishes, carnivorous species had the highest metal concentrations followed by herbivores and omnivores. Besides, the euryhaline fishes possessed higher metal concentrations then the stenohaline fishes. Moreover, the average concentration of metals in demersal fishes was lower than the benthic fishes. However, the concentrations of metals in sampled fishes differed largely, which might be a result of different ecological needs, metabolism, and feeding patterns of the examined fishes [52,53]. Many studies reported the metal concentrations in fishes depend mostly on their habitat type [52,54]. It is commonly observed that the sediment is the major uptake pathway for metal contamination and plays a critical role in the heavy metal uptake for fish [55]. Fish living near the sediments of the waterbody and feeding on humic substances and benthic invertebrates accumulate and transfer heavy metals from sediments to fishes [55]. Therefore, benthic and benthopelagic fishes generally exhibit higher concentration of metals than demersal fishes [36,54]. Besides, a previous study reported that the metal concentrations in piscivorous species in the higher trophic level tend to accumulate more metals than omnivorous and herbivorous species [55], which supports our findings. However, this finding suggests that the metal concentrations in fishes are not only influenced by the habitat but also bio-accumulation through the food chain [14,[55][56][57][58]. In the present study, the mean concentration of Cu in the fish was 36.44 ± 9.19 (mg/kg wet weight). The concentration of Cu was found to vary from 30.29 ± 3.934 to 48.59 ± 9.935 (mg/kg) among all the fishes. The highest concentration of Cu was found in Grass carp (48.59 ± 9.935 mg/kg), whereas the lowest concentration was found in Poa (30.29 ± 3.934 mg/kg). However, the Cu concentration was higher compared to the national and international guideline values and the previous studies of the same geographic region (Table 3). Previously, the maximum concentration of Cu in Bangladesh was recovered from the Bangshi river [32]. Even most of the international guidelines restricted the Cu concentrations within 30 mg/kg. The concentration of Cu we obtained surpassed all the previous findings in Bangladesh and other international guidelines as well (Table 3). Besides, the range of Cu concentration found in the fishes from Asafo market, Ghana, ranged between 0.02-0.156 [59], and Pearl river, China, was within 1.17-6.72 [60]. All these studies reported lower concentration of Cu than our findings. Though the optimum concentration of Cu is important for the body as it produces hemoglobin and some other vital enzymes, the excess amount may lead to malfunction of liver and kidney [38]. Notably, for the trace element Cu, the maximum recommended level for 1-3 years old children is 1.0 mg/day, and for 19-70 years old males/females, it is 10 mg/day. Therefore, an excess amount of Cu over the recommended levels may lead to organ damage (kidney, liver) [61]. In terms of Pb, the concentrations in fish muscles ranged from 0.202 to 0.68 (mg/kg wet weight). The highest Pb concentration was found in Koi (0.68 mg/kg), whereas the lowest Pb concentration was found in Loitta (0.202 mg/kg). Pb concentrations found in the examined fishes maintained the following decreasing order: Koi > Chiring > Bata > Bighead carp > Rui > Catla > Gulia > Kalibaus > Grass carp > Pabda > Ricksha > Koral > Tilapia > Poa > Loitta. Based on the FAO [62] and WHO [77], maximum permissible concentrations for Pb are 2 and 0.5 mg/kg, respectively. Based on the concentrations obtained, Pb concentrations in the muscle of all fishes were below the threshold limit from the WHO [77]. However, the mean Pb concentrations in the present study were lower than the data reported earlier for the Koral and Poa in the same geographic region [75]. Besides, a more or less similar Pb concentration were obtained in fish tissue by Staniskiene et al. [78] and Copat et al. [79].
Optimum Cr concentration in the diet has an important role in lipid and glucose metabolism [31,80]. However, the excess Cr consumption may lead to acute pulmonary disorders and organ damage like lungs, kidney, and liver [18,81]. The recommended maximum permissible concentration for Cr is 50 mg/kg from the WHO [77]. In our study, the mean Cr concentration in the muscle of fish ranged from BDL to 9.685 mg/kg and did not exceed the proposed limit from the WHO [77]. Cr concentrations among the fish species maintained the following descending order: Poa > Koral > Koi > Grass carp > Chiring > Ricksha > Bata > Bighead carp. However, the mean Cr concentration in the muscle tissues of Poa, Koral, Koi, Grass carp, Chiring, Ricksha, Bata, and Bighead carp in the present study was found to be higher than the data reported for eight species from the Meghna river estuary [75] and was found to be lower than the concentration found in C. carpio and S. lucioperca in the Beysehir Lake [82].
The Fe concentration was the maximum obtained compared to all other elements analyzed in the different species of fishes. The maximum permissible concentration for Fe is 100 mg/kg [83]. In the present study, the mean concentration of Fe was 162.198 mg/kg, and ranged largely among the species. The highest concentration of Fe recorded in Bighead carp (208.51 mg/kg) and the lowest value was in Koral (135.58 mg/kg). Fe concentrations in the muscles of fifteen fish species were in the following decreasing sequence: Bighead carp > Grass carp > Ricksha > Poa > Catla > Bata > Chiring > Gulia > Tilapia > Kalibaus > Loitta > Pabda > Rui > Koi > Koral. Based on the values recovered, Fe concentrations found in the muscles of all fishes exceeded the permissible limit by the WHO [83]. Besides, the concentration of Fe was higher than the earlier study of Bhuyan et al. [28], where the range of Fe reported was 7.85 to 147.77 mg/kg. However, mean concentrations of Fe in all species were lower than fishes from Gorgan Bay [76]. Fe is an essential micronutrient for the fishes, as a vital component regarding cellular respiration and oxygen transfer [84]. However, acute Fe overdose is potentially life threatening and also slowly developing damages to organs like heart and liver [18]. Besides, the excess amount of Fe acts as a catalyst in Fenton reaction, responsible for generating free radicals which is toxic [85].
In the present study, the mean Zn concentrations in the muscle of fish species were 113.326 mg/kg. Zn concentrations were found in fish in the following sequence: Bata > Rui > Gulia > Tilapia > Grass carp > Chiring > Kalibaus > Poa > Koi > Loitta > Catla > Koral > Bighead carp > Pabda > Ricksha. According to the FAO/WHO [86], the maximum permissible amount of Zn for human consumption is 30 mg/kg. Zn concentrations found in the muscles of all fishes exceeded the guideline value [86]. The concentration level of Zn in the fishes is almost alike to the reported value of the Bangshi river [32]. However, the mean concentrations of Zn in all species were higher than the other international reports [60,79,87,88]. Zn has a tendency to be accumulated in the fatty tissues of fishes and other aquatic organisms, and likely to affect the reproductive physiology in fishes [89]. Besides, the chronic exposure to Cu and Zn is reported to be associated with Parkinson's disease [90].
According to the FAO/WHO [68], the maximum permissible concentration for Hg is 0.5 mg/kg for human. There is no known physiological requirement for Hg in animal metabolism, and high Hg exposures can result in severe toxicity [91]. In the present study, the mean Hg concentration in the muscle of fish species was 1.657 mg/kg. Hg concentrations decreased in the following order: Loitta > Catla > Koral > Poa > Ricksha > Bata > Rui > Gulia > Grass carp > Bighead carp > Tilapia > Chiring > Koi > Pabda > Kalibaus. From the hierarchy, the highest value of Hg was 2.899 mg/kg in Loitta, and the lowest value was 0.72 in Kalibaus. Hg concentrations found in muscles of all fishes were above the proposed limit by the FAO/WHO [68], which may pose a threatening consequence. However, the concentration of Hg was in line with the results of Ullah et al. [92], where the concentration range was from 0.021 to 0.121 mg/kg in the highly consumed cultured fish in Bangladesh. However, the mean concentration of Hg in all species were lower than the fish of the Pearl river and marine fish in Malaysia [60,93].

Metal Pollution Index (MPI)
The MPI was considered using heavy metal concentrations in the fish species and used to compare the total metal contents of the muscle of the examined fishes. The MPI is generally used to define the polluted degree of heavy metals in tissues of fish. It is considered that the higher value of estimated MPI describes the higher degree of contamination in fish [94]. The highest MPI value was obtained for Tilapia and the lowest for the Pabda. The distribution pattern of total concentrations of heavy metal accumulations in the studied fish species follow the order: Tilapia > Rui > Poa > Loitta > Ricksha > Gulia > Bighead carp > Catla > Koral > Grass carp > Koi > Kalibaus > Chiring > Bata > Pabda (Table 4). In recent years, the Tilapia has been the most consumed and cultured fish species in Bangladesh. Therefore, the high MPI value of Tilapia is a matter of metal contaminationrelated health hazards to local people.

Human Health Risk Assessment
Fish constitute a significant part of the diet of Bangladeshi people. Herein, we anticipated that the local population consumes fish and, since muscle is the most edible part of fish for humans, its intake risks must be taken into account. However, the risk assessment results are summarized in Tables 4-6 for EDI (estimated daily intake), THQ (target hazard quotient), and CR (carcinogenic risk), respectively.

Estimated Daily Intake (EDI)
Heavy metals tend to accumulate in various organs of aquatic organisms, especially in fish, which may enter into the human metabolism through consumption, causing severe health hazards [95]. Thus, the daily intake of some selected trace metals was estimated and compared with the recommended values to assess whether the metal levels found in fish samples from the Noakhali fish market were safe for human consumption (adults & children) (Table 5). This study considered only the fish muscle, as humans mostly consume this portion. The highest recorded EDI values were 0.1285 and 0.5677 (mg/day/person) found in Fe for adults and children, whereas the lowest recorded EDI values were found in Pb (0.0004 and 0.0018 mg/day/person for adults and children). Children's EDI values were higher than the EDI values of adults for all the metals. Ingestion of the metals through the intake of aquatic foods was the primary exposure path instead of possible risk effect from inhalation and direct dermal contact [96]. For ingestion, the results in the study area for adults and children were below the recommended daily allowance (RDA), presented in the following descending order: Fe > Zn > Cu > Cr > Mn > Hg > Ni > Pb. Therefore, EDIs, lower than RDA, indicated a possible lower health effect for the targeted groups of people (adults and children). However, it was not a permanent measurement process to conclude 'acceptable limit' and 'unacceptable limit', based on doses lower than RDA/Rfd [38,40].

THQ and HI
THQ and HI proposed by USEPA [100] are the parameters for risk assessment that compare the ingestion amount of a pollutant with a standard reference dose and have been widely used in the risk assessment of metals in contaminated foods [76]. In addition, the THQ value has been recognized as one of the reasonable parameters for the risk assessment of metals associated with consuming contaminated fish [101]. The threshold limit for THQ is 1 suggested by USEPA [49]. The result described that the mean THQ of all the species was below 1 for both adults and children ( Table 6). None of the metals in all the species exceeded the threshold limit, which indicates that the intakes of metals by consuming these species do not result in an appreciable hazard on the human body. The highest THQ were 8.17 × 10 −3 and 3.61 × 10 −2 found in Cr, and the lowest THQ were 1.03 × 10 −4 and 4.06 × 10 −4 found in Cu for both adults and children, respectively. For all the cases, the THQ was higher in children than adults. The findings enhanced the necessity of evaluating hazard index (HI), where surpassed HI unit expositions determined the alarming concern of health risk for the local consumers [50,102]. The investigated HI did not surpass the suggested limit. Our findings revealed that the HI of metals for species maintained the descending order: Tilapia > Poa > Koral > Koi > Grass carp > Chiring > Ricksha > Bighead carp > Loitta > Bata > Gulia > Catla > Rui > Kalibaus > Pabda. HI exceeding 1 indicates that the metals are toxic and hazardous to human health [101]. In the present study, the average HI values for all fish species were below the threshold value, which indicates that the intakes of metals by consuming those fishes do not result in an appreciable hazard risk for the human body. However, due to the absence of a definite dose relationship, THQ and HI are not considered as a direct measurement of risk concern [103].

Carcinogenic Risk (CR) Assessment
Due to the unavailability of the carcinogenic slope factor for maximum metals, carcinogenic risk (CR) was calculated only for Pb and Cr ( Table 7). The range of the CR found in the selective organisms for Pb and Cr was 1.362 × 10 −9 to 4.583 × 10 −9 and 2.341 × 10 −6 to 8.629 × 10 −5 in adults, respectively, while 6.009 × 10 −9 − 2.023 × 10 −8 and 1.033 × 10 −5 − 3.809 × 10 −4 in children. Generally, a CR value above 10 −4 is unacceptable, whereas CR ranging from 10 −4 to 10 −6 is regarded as an acceptable carcinogenic risk, and below 10 −6 is negligible [40]. In our study, CR value of Pb was negligible for both adults and children, and the CR exposures of Cr was in the acceptable range for both adults and children. The results also specified that children were more susceptible to CR exposures than adults.

Source Identification
The strong and moderate correlation between elements indicates their sources are similar, especially from the point and non-point sources [28]. If no correlation exists among the elements, even a single factor does not control the metals [104]. In the correlation matrix, there was a strong positive correlation between Ca vs. Sr (0.80323) with a 99% confidence level (p < 0.01 significance), a moderate negative correlation between Hg vs. Pb (−0.6004), and a moderate positive correlation between Fe vs. Rb (0.5551) with 95% significance level (p < 0.05 significance) ( Table 8). Such correlation indicates that their origins are probably similar, and they might have a common anthropogenic source like industrial effluents, municipal wastes, and agricultural inputs. The study region is strongly polluted with Ca and Sr and moderately contaminated with Hg, Pb, Fe, and Rb. In PCA, the components were taken into account whose eigenvalues were greater than 0.5 ( Figure 2 and Table 9). PCA explained 99.99% of the data variation and a total of 8 significant PCs were extracted with an eigenvalue > 1. PC1 explained 30% of the total variances and exhibited an eigenvalue of 2.41. PC1 was dominated by Ca and Sr with the loadings of 0.54 and 0.52, respectively. The employed PCA revealed that the source of origin of the metal was anthropogenic. Ca, Sr, Pb, and Fe were the dominant compounds in PCA analysis due to their high loading scores in respective components. PC2 explained around 25% of the total variance and where Pb contained the highest loading scores (0.62). Besides, PC3 explained 17% of the total variance and was dominated by Fe (0.55), whereas PC4 explained 10% of the total variance with a maximum loadings of Zn (0.56) and Hg (0.57). The loadings of Zn and Hg are very close in PC4, which represented a similar source of these metals. However, PC5 and PC6 explained around 7.82% and 6% of total variance with maximum loading of Zn (0.55) and Fe (0.67), respectively. Besides, PC7 explained 3% of the total variance with a maximum loading of Hg (0.51) and Pb (0.68), and PC8 explained 0.37% of the total variance with moderately favorable loading of Sr (0.66). From the component seven, the importance of Hg and Pb are very close, which reflects the precise origin of the metals. Although Pb occurs naturally in the environment, anthropogenic activities such as fossil fuel burning, mining, and manufacturing around the area contribute to the release of high concentrations [105]. Hg is utilized in the electrical industry (switches, thermostats, batteries) and other numerous industrial processes, including the production of caustic soda, in nuclear reactors, as antifungal agents for wood processing, as a solvent for reactive and precious metal, and as a preservative of pharmaceutical products [106].  Cluster analysis classifies variables into homogenous clusters in the form of dendrogram with variables that show similarities in the same group and dissimilarities between different groups [107]. Hierarchical cluster analysis (HCA) was used to determine the relationship between the metal concentration and their possible source. The HCA was established at (Dlink/Dmax) × 100 < 0.2, the Euclidean distance of similarities in variables, which represented two distinct groups of clusters ( Figure 3). Cluster 1 consisted of Hg, Rb, Pb, and Cu that could have been arising from agricultural activities, electroplating materials and lubricants used near the study area. Besides, Fe, Ca, Sr, and Zn confined in cluster 2 that could have been attributed to natural (rock and soil weathering, etc.) or human activities like chemical and pharmaceutical industries, tanneries, industrial effluents, and others. Furthermore, correlations among the metals identified in the multivariate analyses also indicated the resemble accumulative characteristics for those metals presented in the aquatic organisms [108].

Conclusions
This study provides information on the levels of fourteen trace and heavy metals (Pb, Cr, Cu, Zn, Mn, Fe, Hg, Ni, Ca, Co, Se, Rb, Sr, and As) in the fifteen commercial fish species collected from the local fish market of Noakhali and the lower Meghna River. Iron (Fe) showed the highest accumulation level, whereas the Pb, Cr, Ni, and As levels in the muscle tissue of the studied fish species were lower than the permitted limits suggested by the WHO and FAO. However, the concentrations of Cu, Zn, Fe, Hg, and Ni in the fish muscle tissue exceeded the permissible limits suggested by the WHO and FAO. Among the organisms Tilapia was the most susceptible to metal accumulation and poses the maximum risks. Besides, the carnivorous, benthic, and the euryhaline species were the highest accumulator of metals on that area, which revealed the influence of habitat preferences and bio-magnification of metals through food cycle. EDI, THQ, and HI values for both adults and children were within the threshold limit and depicted that none experienced non-carcinogenic health risks. On the other hand, the carcinogenic health risks for Pb and Cr in all fish species were in the safe range (10 −6 to 10 −4 ) for adults and children. However, the correlation matrix, PCA, and hierarchical cluster dendrogram demonstrated that most of the elements in fishes arose from the anthropogenic sources. Hence, to protect the consumers from the derogative health risk effect, the release of toxic chemicals should be checked in a proper monitoring process.