Bioaccumulation and potential human health risks of metals in commercially important fishes and shellfishes from Hangzhou Bay, China

Hangzhou Bay is facing severe anthropogenic perturbation because of its geographic position. We studied species-specific bioaccumulation of metals in commercially important fishes and shellfishes, and calculated the potential human health hazards through their consumption, which has not been reported earlier from this area. The hierarchy of metal concentration in organisms was in the decreasing order of Zn (10.32 ± 7.13) > Cu (2.40 ± 2.66) > As (0.42 ± 0.26) > Cr (0.11 ± 0.08) > Cd (0.07 ± 0.07) > Pb (0.05 ± 0.02) > Hg (0.012 ± 0.009). Except for Cd and As concentrations in fishes, metal concentrations have not exceeded the national and international guideline values. P. laevis and P. trituberculatus were the most bioaccumulative of the species studied. According to the non-carcinogenic risk assessment, children were more susceptible to metal contamination than adults. The carcinogenic risk (CR) values indicated that children were likely to experience carcinogenic threats for taking cancer-causing agents As and Cd through fish consumption. In terms of organisms, intake of two crab species, P. trituberculatus and E. sinensis, as well as the oyster species P. laevis, could be detrimental to consumers.

. Sampling locations and the study area of the Hangzhou Bay. The development zones (DZ), bridges and ports were graphically presented following previous literature 39 .

Materials and methods
Study area. The Hangzhou Bay, a typical funnel-shaped macrotidal estuary 47,49 , is located in the northern Zhejiang province, covering about 8500 km 240 . This area is surrounded by six megacities and nine industrial parks (Fig. 1). Generally, industrials parks are the major contributors to the country's economy 38 ; hence hundreds of millions of m 3 wastewater are discharged into this area from those industries 38,50 . Moreover, except for the megacity Shanghai, there are 22 economic-technological development zones (ETDZ) and high-tech industry development zones (HIDZ) in this area ( Fig. 1) 39 . Among these development zones (DZ), different areas are famous for distinguished industrial and economic activities. For example, Ningbo is a port city that handles millions of tons of cargo and ranks the world's first 39 . The logistics industry is famous in the City of Zhoushan, and the textile industry in the City of Shaoxing supports a considerable contribution to the regional economy. Besides, Shanghai city, located at the north bank of HB (Fig. 1), is one of the largest cities in the world, covering an area of 6430 km 2 and the most prominent economic hub in China 51 . Because of this geographical setting, the HB faces severe anthropogenic disturbance and acts as the ultimate sink of pollutants from these DZ.

Sampling parameters, methods and analysis.
This report is the sequel of our recently published paper from this area 48 . Therefore, for the detailed sampling methods and analysis, the paper mentioned above is referred to. In brief, samples were collected from 33 sites during four seasons (Spring-May 2018, Summer-July 2018, Autumn-October 2018, Winter-January 2019) to analyze Cd, Cu, Pb, Zn, Hg, As and Cr in the organism's muscles. The bottom trawling collected four types of organisms, including six fish species, two crabs, two prawns, and an oyster species (Table 1). The sampled organisms were identified following existing literature [52][53][54] based on their key characteristics. The edible muscles from these organisms we cut and stored in -20 °C. All procedures strictly followed the specification of marine organism's analysis in China 55 . In brief, the muscle tissue samples were freeze-dried in the laboratory. An amount of 0.5 g samples were then placed in nitric acid (10 ml) and perchloric acid (1 ml). The solution was then put into a microwave digester for 2 h. After that for each samples the total concentrations of metals were analyzed using the Flameless Graphite Furnace Atomic Absorption Spectrometry, Flame Atomic Absorption Spectrometry and Atomic Fluorescence Spectrometry (Supplementary  Table 1). Moreover, the Chinese National certified reference material (BW-HZ001 quality control sample of heavy metals in organisms) was used for the validation and accuracy. The detection level, recovery range of all Bioaccumulation factor. Bioaccumulation factor (BAF) is generally calculated as the ratio between metal concentrations in muscles and water 56,57 and considered as the degree of metal concentrations in the organism 57,58 .
Here C m is the metal concentration in the organism's muscle, C w is the metal concentration in water. Further, BAF was classified as less bioaccumulative (BAF < 1000), bioaccumulative (1000 < BAF < 5000) and highly bioaccumulative (BAF > 5000) 59 .
Human health risk assessment. Estimated daily intake (EDI). EDI was assessed based on the metal concentrations in food and their daily consumption amount and calculated by the following equation:  60,61 .
Targeted hazard quotient (THQ). Targeted Hazard Quotient (THQ) is the way of determining non-carcinogenic risk to the local inhabitants provided by USEPA 62,63 . THQ was calculated as the ratio of EDI and oral reference dose (R f d) 64,65 following the equation below: The R f Dvalues for all the metals are given in Supplementary Table 2. The ratio < 1 reveals no non-carcinogenic risk effects, and the ratio > 1 implies that the community is likely to have a non-carcinogenic risk.
Hazard index (HI). HI is calculated as the sum of individual non-carcinogenic risk (THQ) for all the metals (Cd, Cu, Cr, Pb. As, Hg, Zn) following the equation: If the HI value exceeds the threshold value (10), the exposed consumer will face significant non-carcinogenic health risk 23,66,67 .
Carcinogenic risk (CR). The carcinogenic risk (CR) was assessed to evaluate the possibility of cancer in an individual over the lifetime for the exposure of cancer-causing agents 28,57,68 . The acceptable range of the carcinogenic risk is 10 -4 to 10 -6 , and the CR values higher than 10 -4 will probably build the likelihood of cancer-causing hazard impact [69][70][71] . CR was calculated as follows; Here, CSF is the oral slope factor of cancer causing agents (mg/kg/day) 69 , which is only available for Pb (0.0085) , Cd (6.3), As (1.5) and Cr (0.5) 69 .
Statistical analysis. The mean metal concentrations and their standard deviation (SD) in the organism's muscles were calculated using Microsoft excel. The Shapiro-Wilk tests were conducted using PAST (version 3.0) normality test, and dataset were square root transformed before further multivariate analysis. Analysis of Variance (ANOVA) and Kruskal-Wallis tests (when ANOVA did not appear) were applied with the Origin pro for the non-parametric test. Levene's test of homogeneity in terms of ANOVA was adopted to calculate the homogeneity of variance. Pearson correlation through PAST was analysed to identify the relationship between heavy metals in the organism's muscle. Principle component analysis (PCA) was implemented (by Origin pro) to show the association of heavy metals in the organism's muscles. The hierarchical cluster analysis based on Euclidean distance and the Ward-Linkage method was investigated to determine the connection between metal concentration and potential sources.
Ethics approval and consent to participate. This study involves fish, crab, shrimp, and oyster animal testing approved by the Laboratory Animal Ethics Committee of the Second Institute of Oceanography of the Chinese Ministry of Natural Resources, and all methods were carried out in accordance with relevant guidelines and regulations.

Results and discussion
Metal concentrations in fishes and shellfishes. The concentrations of selected seven metals from 11 species are presented in Table 1. The mean metal concentrations (wet weight) of Cu, Pb, Zn, Cd, Hg, As and Cr in selected fishes varied as 2.40 ± 2.66, 0.05 ± 0.02, 10.32 ± 7.13, 0.07 ± 0.07, 0.012 ± 0.009, 0.42 ± 0.26 and 0.11 ± 0.08 mg/kg, respectively. The hierarchy of mean metal concentrations was Zn > Cu > As > Cr > Cd > Pb > H g. Among different species groups, all the metal concentrations were higher in crab, prawn and oyster than the fish species but not exceeded the national and international guideline values. Some species of fish exceeded the guideline values of Cd and As concentrations (Table 1).
Overall the concentrations of Zn were highest compared to all metals in all species. The maximum concentrations of Zn were recovered from P laevis. They maintained the following decreasing trend: P. laevis > E. sinensis > P. trituberculatus > E. annandalei > E. carinicauda > C. nasus > C. joyneri > L. ocellicauda > C. lucidus > H. nehereus > M. miiuy (Table 1). In our study, the concentrations of Zn ranged from 2.66 to 41.5 mg/kg (mean 7.15 mg/kg), and the average concentration of Zn was 4.3, 187, 133, 860, 24 and 93 folds higher than the average concentrations of Cu, Pb, Cd, Hg, As and Cr respectively. The mean concentrations of Zn in fish, crab and Scientific Reports | (2022) 12:4634 | https://doi.org/10.1038/s41598-022-08471-y www.nature.com/scientificreports/ prawn were within the national and FAO/WHO guideline values. But the Zn concentrations in oyster exceeded the national first limit standard. The concentrations of Zn in oysters varied from 41.5 to 16.52 mg/kg (average 28.56 ± 11.85 mg/kg). Besides, the Zn concentrations in all fish species were lower than the related studies at the Pearl River 72 , the Taihu Lake 73 , the Daya Bay 74 and the Palk Bay 75 ; but higher than the concentration recovered from the Santa Maria Bay, Mexico 76 . Similarly, in crabs, prawns and oysters, the Zn concentrations were lower than the previous study at the Xiangshan Bay 11 and the Taihu Lake 77 (Table 1). Invertebrates like oysters and barnacles accumulate specific metals (Cu and Zn) to phenomenally high concentrations 78 . Along with the metal characteristics and organisms bioaccumulation capability, the various environmental processes might be essential to regulate the bioavailability 79 . Such as, significant correlations were observed between the sediment and oysters heavy metal concentrations 78 . In the previous report, it was discerned that most of the heavy metal concentrations in the sediment of this area were within the national and international standard 48 . Besides, high metal concentrations are sometimes found in the sediment but not in the oyster's tissue 79 . The high particulate matter in the sediment and body size of the oysters may cause a tremendous difference between the metal concentrations 78,79 . Moreover, a high amount of metal can be found in the gills, mantle and viscera 78 , which was not focused in this study. According to the national guideline value of China, the permissible limit for Cu concentrations is 20, 100 and 10 mg/kg, respectively, in fish, crustaceans and shellfish. FAO assigned that 30 mg/kg is the permissible limit of Cu concentrations for fish and shellfish 82 . The concentrations of Cu in our study varied from 0.07 to 14 mg/kg (mean 2.40 mg/kg) among all species. The maximum concentrations of Cu recovered from the collected oyster species P. laevis. Among the fish species, the maximum concentration of Cu found in L. ocellicauda and minimum in H. nehereus (Table 1). Comparing with the fish species, crab, prawn and oyster comprised considerably higher Cu concentrations. Cu is a vital element for the body as it forms hemoglobin and other essential enzymes, but excess Cu consumption is responsible for the malfunction of the liver and kidney 23,84,85 . However, in our study, none of the species contained a higher concentration of Cu than the national and international guideline values. The mean Cu concentrations in fishes (0.45 mg/kg) were higher than the concentrations recovered from the Meiliang bay 83 and Santa Maria Bay 76 but lower than most of the other related studies 66,[72][73][74][75] . However, the concentrations of Cu in crustaceans was lower than the reported value in the Xiangshan Bay 11 , whereas, in oysters it was higher than the previous report 77, 83 (Table 1).
In comparison with Zn and Cu, the concentrations of Pb, Cd, Hg and Cr concentrations were very lower in all the species. Among them, the concentrations of Pb ranged from 0.006 to 0.3 mg/kg (mean 0.059 mg/kg); and the maximum concentrations of Pb found in P. laevis (mean 0.08 ± 0.01 mg/kg) and minimum in E. carinicauda (0.01 ± 0.001 mg/kg) ( Table 1). The national permissible limit for Pb concentrations is 2, 2 and 0.1, respectively, for fish, crustaceans and shellfish. FAO set up the guideline values as 1 and 0.5 mg/kg for fish and shellfish, respectively. Therefore, the Pb concentrations in our study were not only within the guideline values, but also lower than the national and international reports 6,11,66,72,[74][75][76][77]85 . The concentrations of Hg did not vary largely among the tested species, though the maximum concentrations recovered from P. laevis. Similar to other metals, the concentration of Cr was maximum in P. laevis followed by P. trituberculatus and E. sinensis. Both of the elements' concentrations were within the national and international standard values set up by FAO 81,82 (Table 1).
Cd concentration was maximum in P. laevis and minimum in C. joyneri. According to the Chinese national standard, the guideline value for Cu concentration in fish is 0.06 mg/kg, 2 mg/kg for both crustacean and shellfish. According to FAO, the permissible limit for Cd concentrations is 0.2 and 2 mg/kg for fish and shellfish. In our study, H. nehereus, L. ocellicauda and M. miiuy crossed the FAO permissible limit, whereas other species were in line with the standard criteria. The mean Cd concentration in fish (0.023 mg/kg) was almost similar to the values reported from the Taihu lake (0.03 mg/kg) 73 , higher than the Daya bay (0.011); but much lower than Meiliang Bay (0.173 mg/kg) 83 , Pearl River (4.61 g/kg) 72 , and Karnafuli estuary (0.39 mg/kg) 66 . Besides, the Cd concentrations in shellfish exceeded the value reported earlier from the Meiliang bay (0.19 mg/kg) 83 but were lower than the Taihu lake (1.7 mg/kg) 77 . In terms of crustacean, both crab species, especially P. trituberculatus had higher concentrations of Cd than the similar study in the Xiangshan bay 11 ; and both of the prawn species had higher Cd concentrations of that report (Table 1).
In terms of As, the concentrations ranged from 0.04 to 1.48 mg/kg (mean 0.42 ± 0.26 mg/kg), and the mean concentration was higher in crab and prawn than in other groups. P. trituberculatus comprised the maximum As concentrations, and the hierarchy of spices in terms of As concentration was as follows: P. trituberculatus > E. sinensis > P. laevis > E. carinicauda > E. annandalei > C. joyneri > C. nasus > C. lucidus > L. ocellicauda > M. miiyu > H. nehereus. FAO set up the maximum permissible limit of As is 0.26 mg/kg in fish and 1 mg/kg in shellfish 81,82 . In our study, C. joyneri exceeded the permissible limit of As concentration, but other species was within the permissible limit. The inorganic As is more deadly than the organic form and can act as a cancer-causing agent for human beings if they consume the inorganic arsenic for a longer period 84 . But, it is a matter of satisfaction that most of the As concentrations in seafood are in organic form and can be directly excreted through the urine 66,84,86 . Our finding is comparable with the similar study in the bay and coastal region in China and worldwide ( Table 1). The mean As concentration in Daya Bay was reported as 0.10 mg/kg in fish 74 , which is lower than our study. In Karnafuli estuary Bangladesh, the reported As concentration was 4.89 mg/kg in fish 66 , and in our study, the concentration was very lower than their findings. Besides, the concentration of As in shellfish, crab and prawn in Taihu lake 77 and Xiangshan Bay 11 were higher than the concentrations recorded in our study (Table 1).
Three species (C. nasus, C. lucidus, P. laevis) were sampled in four seasons for comparing the seasonal variation of heavy metal concentration. On their muscles, the concentrations of Zn were higher than other metals throughout the year. In C. nasus, the concentrations of Zn were maximum in summer (11.84 mg/kg) followed by spring (8.66 mg/kg), autumn (4.8 mg/kg) and winter (3.50 mg/kg). Cu concentrations were maximum in spring (1.09 mg/kg) followed by winter (0.47 mg/kg), summer (0.32 mg/kg) and autumn (0.21 m/kg) (Supplementary Fig. 1 Fig. 1). But, the concentrations of Cu were maximum in autumn (14.03 mg/kg), followed by spring, summer and winter. The As concentrations were considerably higher in autumn (1.48 mg/kg) than in the other three seasons. Other metal concentrations not varied largely on P. laevis. As a result, the average metal concentrations in P. laevis maintained the seasonality as follows: spring > summer > autumn > winter. The findings of our study are comparable with similar regional studies on seasonal pollution in the Meiliang Bay of Taihu Lake 6 . They reported that higher concentrations of metals were found in summer and winter, but both the fish and oysters showed great seasonality in terms of metal types 6 . In fact, various species take up various types of metals from both water and foodstuffs, and their physical properties are additionally influenced by ingestion and filtration rates, food quality, physiological states and ecological variables 87,88 . Therefore, it is not astounding that the seasonality of metals in different species was unique, regardless of having been sampled in the same period 87 . In our study, the demersal species comprised the highest metal concentrations, followed by benthopelagic and pelagic species. Demersal species lives close to the bottom 89 , and many of them possess special features (such as modified fins for fish) to crawl over the bottom (crab, oyster). Those species are bottom feeders, mostly carnivorous and primarily feed on fish, benthic macroinvertebrates and zooplankton, in contrast with the pelagic fishes. Source identification of metals. Two-way ANOVA revealed that most metal concentrations were significantly different among the sites and seasons at a 95% confidence level (P < 0.05). In contrast to these findings, Cr concentrations in fish muscles were not significantly different among the sites and seasons. Similarly, the concentrations of Pb and As were not significant among the sites, and Cd was not significantly different among the seasons. One-way ANOVA tested the seasonality of metal concentrations, and all metals in C. nasus and C. lucidus showed significant variation in four different seasons ( Supplementary Fig. 1). The Kruskal-Walis test was implemented for P. laevis as the sample size was smaller. It revealed a significant difference in metal concentrations among the seasons (P < 0.01). Levene's test of ANOVA revealed that the metals were not homogenously distributed among the sites and seasons. As Hangzhou Bay is a complex nature of ecosystems, we have applied both univariate and multivarite approaches i.e., pearson correlation, principal component analysis (PCA) and cluster analysis (CA) to identify the sources of metals. Pearson correlation analysis of metal concentrations in organism's muscles is shown graphically in Supplementary Fig. 2. According to the correlation analysis, strong positive relationships were found between the metals in the organism's muscle. Such correlations between particular elements in nature may depict a similar degree of contamination and discharge from a similar source of contamination 23,91 . Among them, the strongest relationships were found between Cu-Zn (r = 0.73), Cu-As (r = 0.72), Cu-Cd (r = 0.87) at 95% confidence level. Zn, As, Cr and Cd showed significant positive relationships with all other metals, but Pb    Fig. 2). Therefore, this correlation reveals that most of the metals which are significantly correlated with each other may originate from the same source, either anthropogenic or natural, whereas Pb might be released from different source. Among the multivariate analyses, PCA is based on eigenvalue analysis of the correlation matrix which estimates the correlation structure of the variables in a multidimensional data set 92 . Each variable has a loading which shows how much each variable contributes to the meaningful variation in the data and to interpret variables relationship. Practically, higher number PC components that explain only a small proportion of variance is ignored 93 . In this study, the association of metals in organisms' muscles were interpreted through PCA bi-plot (Fig. 2). The first two components of the PCA have explained 77.54% of total variation in metal concentrations and the rest components showed small amount of variation. PC1 was dominated by Cu, Zn, Cd and As with the loadings of 0.41, 0.45, 0.41 and 0.41, respectively. On the other hand, PC2 was highly associated with Pb (loading 0.75) and moderately with Cr (loading 0.42). In addition, partial representation was found for Hg in PC1 (loading 0. 37) and PC3 (loading 0.71). Therefore, analysis of the PCA results interprets that Cu, Zn, Cd and As were originated from similar source i.e. anthropogenic origin as their levels exceed (in some cases) or close to the legal limits, and on the other hand, Pb being the far below of the permissible limit was originated from natural source. Again partial association of Hg in two PCs suggests mixed sources of origin.
Besides, associations for interdependence among the different metals in the study area were also established by means of the hierarchical cluster analysis based on Euclidean distance. The results demonstrated that among the seven metals, six metals (Cu, As, and Zn in cluster 1; Cd, Cr and Hg in cluster 2) were grouped in two major clusters whereas, Pb showed a clear distinction forming an individual group. This findings also indicated that most of the metals were originated from the same sources excepting Pb. It is reported earlier that the heavy metals in this area originates from similar sources 40,94 and are related to anthropogenic activities 40 . Besides, the transport of terrigenous clastic particles was an unequivocal definitive of the formation and dissemination. At the same time, the transition mechanism of trace elements of various substrates and contamination of marine aquaculture were likewise significant contributors to trace metals' sources in this region 94 (Fig. 3).
Bioaccumulation factor (BAF) of aquatic species. The BAFs of the heavy metals in selected organism's muscles are graphically represented in Fig. 4. BAFs associate the number recovered from the portion of collected concentration in a specific organ of a species and the habitat quickly and effectively 23,85,95 . In our study, the highest mean BAF was encountered for Cu, lowest for Cr; and followed the decreasing trend of Cu (1640) > Zn (736) > Hg (387) > Cd (366) > As (288) > Pb (173) > Cr (82) in the organism's muscle. Actually, Cu is www.nature.com/scientificreports/ effectively persistent in an organism's muscles due to being an essential component of living tissue 66,96,97 . In terms of organisms, all species were bioaccumulative, but P. trituberculatus and P. laevis were extremely bioaccumulative. Among the selected species, the oyster, crab and prawns were more bioaccumulative than the fishes. Among all the species, the BAF value was maximum in P. laevis (13998) and minimum in C. joyneri (1098) (Fig. 4). In fact, the bioaccumulation of an aquatic species relies on its characteristics, metabolism of the inspected tissue, invasion pathways, and habitat condition 66,98 . Therefore, each species has a specific bioaccumulative receptivity regarding various metals 23,69,99,100 .

Potential human health risk assessment. Non-carcinogenic risk (EDI, THQ and HI).
Heavy metals poisonous quality could incline people to health hazards through consumption of defiled aquatic food; therefore, it is crucial for assessment 87 . Based on the oral reference dose R f D, the EDI was calculated to measure both the non-carcinogenic and carcinogenic risk of metal consumption through seafood 101 . Besides, the EDI value denotes the exposure of heavy metals, which is performed to evade any detrimental impact on human health 102 . The hazard analysis was conducted for two age groups, children and adults. The calculated EDI of the examined metals associated with fish, crab, prawn and oyster are presented in Supplementary Table 2. For both groups, the EDI of metals through the organism muscles were organised as follows Zn > Cu > As > Cr > Cd > Pb > Hg (Supplementary Table 2). The EDI was higher for the children than the adult, especially for children the EDI value of As, Cd and Zn; and for adults only As had higher EDI values than the Recommended Daily Allowance (RDA) provided by WHO 103 (Supplementary Table 2).Therefore, the higher EDI than the RDA guidelines, revealed that here is a possibility of higher health impact associated with As, Cd and Zn to the consumers, mostly children would be more susceptible. The estimated THQ values (Fig. 5) were higher in the children group than the adult. Particularly the THQ of As was significantly higher in both groups than other metals. For the adult people, the THQ of As for the consumption of E. sinensis, P. trituberculatus, and P. laevis exceeded the threshold value (> 1) suggested by USEPA 104 . Whereas, for the children, the THQ of As was higher (> 1) for all the species but maximum in three species similar to the adult. For both groups, except for As, all the other metals in all the species possessed lower THQ values, which were far lower than the threshold value (< 1). The accumulated THQ of all metals in the concern of HI, also much higher for the children. Though the HI was within the threshold limit (both for adults and children) for all the organisms, it was much higher for oysters, prawns and crabs (Fig. 5). Therefore, the local consumer can suffer from As contamination for the long term seafood consumption. Moreover, an organism's age and size could be important factors for metal accumulation and health risks 28,105 . Though no definite trend of metal accumulation was observed in Capoeta umbla fish species, higher concentrations were found in medium-aged fish compared to the small-sized or older fish 106 . However, metal concentration is determined by the organism's feeding rate with the developmental stage rather than age 105 . Besides, the larger size organisms tend to accumulate more metals 105 ; consequently, they may pose more threat to human health. A previous study reported that the larger-sized Labeo rohita accumulated more metals and had higher hazard potential than smaller fishes like Glossogobius giuris, Puntius sophore and Puntius chola. Our study did not focus on size-based metal accumulation; hence all the fishes sampled were within 200-230 gm and posed little threat to humans (low THQ). However, our results emphasize that the organism's trophic guild and habitat preferences play significant role in the metal accumulation and consequent health hazard. Many earlier studies reported that the contamination level differ largely on the food habitat of the organisms [107][108][109] . Therefore, though the two crab species P. trituberculatus and E. sinensis weighed lower than the fish (200 ± 40 g and 158 ± 42 g, respectively), they accumulated higher metals and hazard potential. Similarly, the oyster species P. laevis was weighed 178 ± 70 g, also showed high bioaccumulation capacity and THQ value for both adults and children, whereas fishes showed low THQ values although weighed more. But, suppose the organisms sampled from the same habitat type and feeding guild, then the organism's size is an important factor regarding metal accumulation and the hazard potential to their www.nature.com/scientificreports/ consumers. However, the THQ and HI have no particular dose relationship that's why they are not considered as an immediate estimation of hazard concern 110 . Carcinogenic risk (CR). The estimated CR of several elements (Pb, Cd, As, Cr) for adults and children are addressed in Table 2. CR values lower than 10 -6 denote the metals' negligible exposure, whereas 10 -6 to 10 -4 means the acceptable range, and higher than 10 -4 indicates the terrible exposure 84,111 . For the adults, all the metals possessed the CR values within the limit of 10 -4 -10 -6 , though the As showed high carcinogenic risk than other elements. In line with the previous studies, children were more susceptible for CR exposure 23,101,112 .
For the children, both the As and Cd showed high CR values than the permissible limit for several species. As showed terrible CR exposure for the children through the consumption C. joyneri, E. sinensis, P. trituberculatus, E. carinicauda and P. laevis. Whereas, besides those species, E. annandalei showed terrible CR exposure for children in terms of Cd (Table 2). For both groups, CR values of Pb were negligible, and the metal exposure for the seafood consumption was organised as As > Cd > Cr > Pb. Therefore, the local consumers are at carcinogenic risk associate with the toxic As, and for children, Cd is an additional matter of carcinogenic concern. It was reported in the Persian Gulf that carcinogenic risk in terms of As could be detrimental to local consumers, though it did not show any non-carcinogenic risk. Actually, 90% of the carcinogenic risk has been found in the As polluted seafood 66 . However, the inorganic form of As is more deadly than the organic form 57,113 , and around 90% of the aggregated As can be evaluated as organic form 84 . Chronic exposure of inorganic As to human lead to various malfunctions 114 including the organ failure such as respiratory tract, circulatory system, digestive system, nervous system and liver 115 . Besides, Cd is also responsible for endocrine malfunctioning, which can cause the failure of the essential organ such as the kidney and brain 116 . Moreover, long term contamination of Cd may cause the dysfunction of the blood circulatory system 117 , bone softening and prostate cancer 118 . Therefore, the USEPA categorized Cd as a priority contaminant and considered as "carcinogenic to humans" (Carcinogenic classification-B) 23 . However, the tissue samples in this study were digested by nitric acid and perchloric acid system, and total heavy metal concentrations were used in the health risk assessment. The use of total heavy metal in health risk assessment is seen to overestimate the risks calculated and cause uncertainties 119 . Besides, the total metal consumption depends on the gross quantity of dietary intake besides the aquatic foods and their metal concentration. Therefore, the health risk of heavy metal should be considered the bioavailability of metal concentration and the total dietary intake along with seafood in future.

Conclusions
The concentrations of seven heavy metals in the 11 aquatic organism's muscles were analysed from the HB area.
Overall, the hierarchy of mean metal concentrations was Zn > Cu > As > Cr > Cd > Pb > Hg. Among different species groups, metal concentrations were higher in shellfishes (crab, prawn, and oyster) than finfish. Notably, the demersal species P. laevis, P. trituberculatus and E. sinensis possessed the maximum metal concentrations, while the pelagic species had several degrees of lower metal concentrations. Though most of the species contained the less metal concentrations than national and international guideline values, As concentrations in C. joyneri and Cd concentrations in H. nehereus, L. ocellicauda and M. miiuy exceeded the FAO permissible limit. Seasonally, the metal concentrations varied with both metal and species types. Most of the metals (except Pb) showed a positive relationship, and Cu-Zn, Cu-As and Cu-Cd showed the strongest correlation. All the species possessed bioaccumulative capability, but P. laevis and P. trituberculatus were highly bioaccumulative. In terms of health risks, both adults and children may experience a high risk for As contamination through seafood consumption. All the fish species asserts non-carcinogenic risks (EDI, THQ and HI) to children, whereas for the adult people E. sinensis, P. laevis and P. trituberculatus consumption would be harmful. However, the CR index values depicted that children are susceptible to carcinogenic risk of AS and Cd contamination, hence unsafe for consumption. Mostly the demersal species of crab, prawn and oyster (especially P. laevis, P. trituberculatus and E. sinensis) may pose a high risk to the consumer for both carcinogenic and non-carcinogenic risk. Therefore, it is highly recommended to thoroughly examine all the aquatic organisms consumed by human including the metal concentration in various organs. The management authorities and policymakers should take in consideration the current metal and metalloids status in the HB area to provide a healthful environment.

Data availability
The datasets generated during and/or analysed during the current study are available from the corresponding author on reasonable request.