Toxic Elements and Mineral Content of Different Tissues of Endemic Edible Snails (Helix vladika and H. secernenda) of Montenegro

The objective of the present study is to determine the differences between the mineral content of various organs of Helix vladika and H. secernenda, the two most abundant edible snail species in Montenegro. The bioaccumulation of 12 examined elements (zinc, manganese, copper, aluminum, cadmium, lead, nickel, iron, chromium, lithium, selenium and mercury) was determined in the hepatopancreas, albumen gland, digestive tract, reproductive system, mantle, foot and shell from three sampling sites (Biogradska Gora, Nikšić and Malesija). The examined populations of H. vladika and H. secernenda showed a significant difference in their lithium and selenium contents. The levels of the most examined metals (Zn, Mn, Cu, Al, Cd, Pb, Se and Hg) varied significantly among organs. The digestive tract and hepatopancreas tend to bioaccumulate selenium and cadmium. The general mean concentration of cadmium in the examined snail tissues exceeded the maximum allowable level at the Biogradska Gora and Malesija sites. Therefore, the use of the Montenegrin edible snails collected from the wilderness for human consumption seems to be limited by their higher bioaccumulation capacity for toxic elements such as cadmium.


Introduction
Snail meat has long been recognized as a high-quality food. Snail meat contains high amounts of protein, a relatively low lipid content [1] and a balanced intake of omega-3 and omega-6 fatty acids that have led to its recognition, due to its nutritional and therapeutic benefits [2,3]. On the other hand, land snails are considered as appropriate sentinel species [4,5] because trace metals tend to accumulate in their tissues, particularly in the digestive glands [6]. The edible snails of the genus Helix could be a powerful tool for environmental pollution assessment [3]. Depending on the soil mineral content and level of contamination, the snails can accumulate large quantities of heavy metals such as lead and cadmium, therefore making the snail meat a very high-risk food product [7].
In Montenegro, the harvesting of edible snails has boomed in the first decade of the 20th century. Based on the data from the International Heliciculture Institute from Cherasco (Italy), Montenegro has 22 snail farms with 7 hectares [8] established before 2010. All of these farms use the extensive (the Italian) method with outdoor breeding, but in practice, they are the most commonly worked with populations of snails collected in the wild that are bred for a while and then exported. At that time, many farmers were attracted by the idea of getting almost immediate benefits. However, after 2010, inadequate legislation as well as poor management in snail farming led to almost all farmers having given up on snail farming, what led to a quick failure of this snail farming practice and the disappearance of almost all farms in Montenegro.
Several species of the genus Helix are present in Montenegro, among them the most widespread and abundant are Helix vladika Kobelt, 1898, and H. secernenda Rossmässler, 1847. The latter species is with wide distribution from the Croatian coastland to Northern Greece, while H. vladika is endemic to Montenegro, but also found in the neighboring parts of Serbia and Albania. Due to its larger size and morphological resemblance to H. pomatia, which is the most widely used edible snail in Europe, H. vladika is more intensely exploited in Montenegro.
It is known that edible snails can accumulate higher concentrations of some essential elements such as selenium (Se), copper (Cu) and zinc (Zn), and therefore can be beneficial for human nutrition. Many studies have mentioned snail meat as a source of selenium for human health and nutrition [2,7]. Most of these studies (e.g., [2,3,7,[9][10][11][12]) have been done on the Roman snail (Helix pomatia), which is the most widely used edible snail species in Europe, while other species are very poorly researched and little is known (or unknown as is the case with our species) about the mineral content of their meat.
This study aims to determine the contents of zinc, manganese (Mn), copper, aluminum (Al), cadmium (Cd), lead (Pb), nickel (Ni), iron (Fe), chromium (Cr), lithium (Li), selenium and mercury (Hg) in different snail tissue of the two most abundant edible snails of the genus Helix in Montenegro. Statistical analyses were conducted to establish possible correlations in the element content distribution with different snail tissues between species and sampling localities. The obtained data are allowed to identify different snail tissues as a potential source of some essential elements for human nutrition. On the other hand, this study was aimed also to evaluate the risk assessment of the accumulation of toxic elements such as Pb, Cd and Hg which are considered dangerous and priority pollutants by the EU legislation (European Parliament and Council of European Union, Directive 2000/60/EC).

Material
The study was performed on two species of edible snails of the genus Helix, i.e., H. vladika and H. secernenda, that inhabit the territory of Montenegro. The snails were collected from three sampling sites: (1)  The Biogradska Gora site where H. vladika was collected is located in a protected area, the Biogradska Gora National Park. Two other sites, Malesija (Skorać) and Nikšić (Dragovoljići), are inhabited by H. secernenda and are located outside the motorized traffic, and were not affected by any source of pollution. Samples were taken in May (at Biogradska Gora and Nikšić) and in September 2018 (Malesija). From each site, 30 snail specimens with almost a similar size were randomly taken. In addition, three subsamples of soils were randomly taken (at a depth of 10 cm) from each site. Snails were put in plastic bags, transferred in the laboratory and thawed at room temperature on a clean tissue paper. The soft tissues were then separated from the shell, and dissected and pooled into six different parts, namely the mantle, hepatopancreas, digestive tract, albumen glad, reproductive system and foot.
The ethical guidelines have been strictly followed in dealing with experiments on animals using the minimum number to obtain a significant result of the experimental work. For the allometric studies, five specimens from each site were selected and measured for the following parameters: shell height, shell weight, tissue wet weight, tissue dry weight and water content (Table 1).

Chemical Analysis
After dissection, the tissue samples were dried at 60 • C for 72 h in sterile Petri dishes to constant dry weights. The soil sediment samples were dried at 60 • C for at least 72 h to a constant dry weight.
The concentrations of Pb and Cd were analyzed by means of graphite furnace atomic absorption spectrometry with Zeeman background correction (GF-AAS, 4110 ZL Perkin Elmer) and other elements (Al, Fe, Mn, Zn, Cu, Cr, Ni, Li) using inductively coupled plasma atomic emission spectrometry equipped with the Meinhard TR 50-C1 nebulizer (ICP-AES, Yobin YvonJY-24).
The samples were weighed in an amount of 0.2 g ± 0,001 g into the Teflon ® vessels of a microwave oven (CEM MDS 2000). As the oxidizing agent, 5 mL of 65%, analytical grade HNO 3 (Merck KGaA, Darmstadt, Germany) was added per sample. The mineralization process proceeded in 5 steps under the conditions shown in Table 2. The digested samples were filtered into volumetric flasks and topped up to 25 mL using bidistilled water. The calibration curve technique was used for the quantification of toxic elements. Single-element calibration standards (Merck, Germany) were used as the stock standards for preparing the working standards. Working standards were prepared by serial volume/volume dilution in glass volumetric flasks with 0.2% HNO 3 (Pb, Cd) or 15% HNO 3 (others elements) in Milli-Q water (Merck, Darmstadt, Germany). The linearity was considered satisfactory and acceptable if the correlation coefficient, r, exceeded 0.999, and in the case of Pb and Cd, we found a correlation coefficient of r 2 ≥ 0.9997. For quality control, blanks, calibration standards of each analyte and standard reference materials were tested during each batch.
The limits of detection (LOD) were determined and calculated based on the data from the calibration curve [13]. In order to evaluate the accuracy of the methods, the following certified materials were used: MODAS-2, Bottom Sediment (M-2 BotSed) or MODAS-3, Herring Tissue (M-3 HerTis) and 2 blank samples. The results obtained by means of the GF-AAS and ICP-AES methods were compared with the certified values, and recoveries were calculated (Table 3). Total Hg (THg) concentrations were determined using atomic absorption spectroscopy (AAS). The assays were run in an AMA 254 mercury analyzer (Altach Ltd., Czech Republic). The AMA 254 analyzer allowed the determination of Hg in samples without having to perform prior mineralization in wet conditions. The detection limit for this device is 0.01 ng THg. The limit of detection was calculated based on three times the standard deviation of the blank. A blank (i.e., an empty sample nickel boat) was analyzed periodically to verify that mercury was not being carried over between samples.
The analytical procedure was checked by determining THg concentrations in samples of one reference material CRM015-050 (metals on sediment analyte concentrations; n = 3). The recovery rate was 99.4%.
The Se concentrations in samples were determined using spectrofluorometric methods (Shimadzu RF-5001 PC, Shimadzu Corporation, Tokyo, Japan). The samples were digested in HNO 3 (Chempur, Poland) at 230 • C for 180 min and in HClO 4 (Chempur, Poland) at 310 • C for 20 min. Then, selenate (Se 6+ ) was reduced to selenite (Se 4+ ) using 3 mL of 9 % HCl (Chempur, Poland). EDTA and hydroxylamine hydrochlorine were used as masking agents. Then, selenium was derivatized with 2,3-diaminonaphthalene (Sigma Aldrich, USA), under the conditions of controlled pH (pH 1-2) with the formation of a selenodiazole complex. Fluorescence of the organic (cyclohexane, Chempur, Poland) layer was measured using an emission wavelength of 518 nm and an excitation wavelength of 378 nm. Blank samples were determined concurrently with the samples proper. The LOD was 0.003 µg/g. The accuracy of the analyses was verified using certified reference material NIST SRM 1946 (Lake Superior Fish Tissue). Recovery ranged from 93% to 98% of the reference value. The precision (RSD %) of the analysis was 2.9%. All samples were analyzed in triplicate and their average value was assessed.

Statistical Analysis
All statistical analyses were performed using the Statistica 13.0 PL software. Data were expressed as mean ± standard deviation. The bioaccumulation factor (BAF) was calculated on a dry weight basis, as BAF = concentration of the target element in the snail tissue/concentration of the same element in the ground.
The normal distribution in each group was tested with the Kolmogorov-Smirnov test, applicable for small samples. Analyses of variance (ANOVA) were computed to test if the observation depended on one or several factors acting simultaneously.
When the assumption of homogenous variances was fulfilled, the post-hoc Tukey's t-test was applied to compare the groups for statistically significant differences. In the case when variances were heterogenous, a Kruskal-Wallis test was used. Differences were considered as significant at the level of p < 0.05. The Mann-Whitney test was used to investigate whether the mineral contents varied significantly between the studied snail species and between sampling seasons.
Pearson correlation analysis was applied to examine the relationship between the concentrations of the studied metals in the selected organs and ground.
In calculating the canonical discriminant functions, the forward method was performed in order to examine the discrimination power properties of variables in discrimination of the investigated tissue. Discriminant analysis (DA) was run on the so-called raw data, introducing appropriate organ grouping. When assessing the received orthogonal canonical functions, the comparison of our own values and the 1 ƛ -Wilks test were used. On the basis of the first two discriminatory functions, the coordinate system of the canonical values was plotted to illustrate the grouping of these values.

Distribution of Element Concentrations in Studied Edible Snail Species of Montenegro
The investigated species are the most widespread edible snails in Montenegro [14], and therefore important for human consumption. Regarding the growing interest in edible snails as a potential group of alternative foods, the determination of the factors that affect the content of metals is gaining in importance. The allometric characteristics of the investigated populations of the two analyzed edible snail species, Helix vladika and H. secernenda, are given in Table 1.
The mean concentrations of the analyzed minerals and toxic elements in the tissues of the studied snail populations from three sampling sites in Montenegro, together with bioaccumulation factor (BAF), are reported in Table 4. The Mann-Whitney test revealed a significant difference between Helix vladika (Biogradska Gora) and H. secernenda (Nikšić and Malesija) only in terms of the content of Li (p = 0.006) and Se (p = 0.002).
In general, the concentration of the most tested metals was lower at the Biogradska Gora site inhabited by Helix vladika and highest at the Malesia site inhabited by H. secernenda (see Table 4). However, for the general mean metal concentration in all snails, no significant difference between the studied localities were found (p = 0.7563). The reason for this may be that the examined localities in our study were not affected by any kind of pollution. In the study ofĆirić et al. [3], the metal content in the examined tissues of H. pomatia from four different environment locations in Pančevo city in Serbia were significantly higher at the polluted sites.
It is well known that metal content variability is not only affected by the possible level of pollution, but also by the geochemical characteristics of the habitat that snails inhabit as well as the season influence [2,3]. In our study, the canonical analysis performed on the metals variability at the investigated sampling sites separated the Malesia locality (see Figure 1a). On the other hand, the resulting scatter plot of the first and second canonical scores shows a part overlap of the Biogradska Gora and Nikšić localities, suggesting that these two sites share similarities in metal content variability. The first discriminant function explains 40% of the total variance and it is dominated by the Fe content (R = 0.99). The second discriminant function explains 25% of the total variance and it is dominated by the Zn content (R = 0.98). The variable with the greatest ability for discrimination along the first discriminant was Fe, while Zn plays an important role in the discrimination along the second discriminant function.
Samples from Biogradska Gora (H. vladika) and Nikšić (H. secernenda) were taken in May while snails from Malesija (H. secernenda) were collected in September, which may have an impact on the metal content in the snail issue. A number of studies have shown that the maximum metal content in snail tissue is reached in autumn and winter and the lowest levels in summer [3,15,16]. In the study ofĆirić et al. [3], the highest concentrations of Cu, Cd, Zn, Fe and Mn in snail tissues were found in October and the lowest in June. In our study, however, the Mann-Whitney test did not reveal a significant difference in seasonal sampling between the studied populations, neither between H. secernenda populations from Malesija and Nikšić. Table 4. Minerals and toxic elements content of the different snail tissues from the three sampling sites. The concentrations of elements in the ground and snail tissues were expressed as mean ± standard deviation calculated on a dry weight (d.w.) basis. a,b-the mean values marked with small letters differ statistically significantly at p ≤ 0.05 horizontally; A-the mean values marked with capital letters differ statistically significantly at p ≤ 0.05 vertically. * Tukey's t-test statistics, ** Kruskal-Wallis statistics. Abbreviations: Gr-ground, Hp-hepatopancreas, Rs-reproductive system, Ag-albumen gland, F-foot, Ds-digestive system, M-mantle, S-shell.   1.5 ± 0.0 BCF 0.03 ± 0.02 0.03 ± 0.03 0.00 ± 0.00 0.00 ± 0.01 0.00 ± 0.10 0.21 ± 1.95 0.01 ± 0.00 0.00 ± 0.00 0.00 ± 0.00 0.03 ± 0.25 0.14 ± 0.28 0.01 ± 0.00 Foods 2020, 9, x FOR PEER REVIEW 9 of 15

Distribution of Element Concentrations in Different Snail Tissues
The metal accumulation in the different tissues of the snails of the genus Helix depends on the investigated snail tissue [3,9]. In our study, we found significant differences in the concentrations of the studied elements between the ground and examined tissues (Table 4). With the exception of Ni, Fe, Cr and Li, all other studied elements showed significant differences in the metal content between

Distribution of Element Concentrations in Different Snail Tissues
The metal accumulation in the different tissues of the snails of the genus Helix depends on the investigated snail tissue [3,9]. In our study, we found significant differences in the concentrations of the studied elements between the ground and examined tissues (Table 4). With the exception of Ni, Fe, Cr and Li, all other studied elements showed significant differences in the metal content between the analyzed snail tissues. (Table 4). The comparison of each metal concentration in tissue at subsequent sampling sites revealed significant differences in the Zn (p = 0.032), Cu (p = 0.017) and Cd (p = 0.039) contents. The BAF values for Zn from 1 to 2.3, indicate a bioaccumulation of this metal in the hepatopancreas.
The highest concentrations of Zn, Ni, Mn, Pb, Cd and Hg were found in the hepatopancreas. The concentration of Zn in the hepatopancreas was highest in snails from Nikšić (387.0 ± 16.9) and lowest in snails from Biogradska Gora (151.6 ± 6.9 mg kg −1 d.w.). The highest concentration of Zn in the hepatopancreas is in agreement with similar studies (e.g., [2,11]) on the metal content in snail tissues of Helix pomatia.
The concentration of Ni in the hepatopancreas ranged from 0.82 ± 0.12 in snails from Nikšić to 2.12 ± 0.36 mg kg −1 d.w. in snails from Malesija. On the other hand, the concentration of Mn in the latter tissue was highest in snails from Nikšić (405.12 ± 11.0) and lowest in snails from Malesija (32.4 ± 0.3 mg kg −1 d.w.). The BAF values for Ni and Mn below 1 ( Table 4) indicate that they do not bioaccumulate in snail tissues.
The Pb concentration in the hepatopancreas ranged between 0.99 ± 0.29 (Malesija) to 5.46 ± 0.18 mg kg −1 d.w. (Biogradska Gora). Further, the Cd concentration in the latter tissue was highest in snails from Malesija (16.38 ± 1.46 mg kg −1 d.w.) and lowest in snails from Nikšić (1.43 ± 0.03 mg kg −1 d.w.). The BAF values for Pb were below 1, indicating that it did not bioaccumulate in the snail tissues. On the other hand, the high BAF values for Cd, ranging from 1.2 (Nikšić) up to 25.2 (Malesija), indicate the bioaccumulation of this element in the hepatopancreas.
The Hg concentration varied from 0.02 ± 0.0008 in snails from Nikšić to 0.054 ± 0.0021 mg kg −1 d.w. in snails from Biogradska Gora. The Fe, Cr, Al, Li and Se contents were highest in the digestive tissue. The concentration of Fe in the latter tissue was highest in snails from Malesija (805.0 ± 67.6) and lowest in snails from Biogradska Gora (111.9 ± 11.1 mg kg −1 d.w.). The Cr concentration was lowest in snails from Biogradska Gora (1.15 ± 0.58) and highest in the digestive tissue of snails from Malesija (4.84 ± 0.43 mg kg −1 d.w). Similarly, concentrations of Al and Li were highest in Malesija (1411.1 ± 172.0 and 1.26 ± 0.07 for Al and Li, respectively) and lowest in Biogradska Gora (73.6 ± 8.0, 0.03 ± 0.05 mg kg −1 d.w. for Al and Li, respectively). As can be seen from Table 4, the BAF values for Hg, Fe, Cr, Al and Li below 1 indicate that they did not bioaccumulate in the snail tissues.
The selenium concentration in the digestive tissues varies from 0.106 ± 0.02 in snails from Biogradska Gora to 0.587 ± 0.015 mg kg −1 d.w. in snails from Malesija. The high BAF values for Se from 2.2 (Biogradska Gora) to 9.6 (Malesija) tend for a bioaccumulation of this element in the digestive tract. The general mean value of selenium in the examined snail tissues of H. vladika (Biogradska Gora) was 0.04 mg kg −1 wet weight (w.w.) and from 0.05 to 0.08 mg kg −1 w.w. for H. secernenda (Nikšić and Malesija). These values were lower in comparison with those reported for the raw meat of Roman snails from Poland (0.09 mg kg −1 ) [2] as well as from those reported for the meat of H. pomatia snails from Moldavia, Ukraine and Russia, which varied in the range between 0.130 and 0.423 mg kg −1 [7]. Nevertheless, in our study, the Se concentration varies among snail tissues as well as between the investigated sites, reaching the highest value of 0.587 ± 0.015 mg kg −1 d.w. in the digestive tissue of snails from Malesia. Based on the available literature, the highest Se content in snail meat was found in H. pomatia from Moldavia and was 0.423 mg kg −1 [7].
The general mean values of selenium in our study were comparable with the content of the latter element in the mealworm (Tenebrio molitor) larvae, which ranged from 0.057 to 0.085 mg kg −1 w.w. [17]. The values of the Se content in the snails of our study (0.04-0.08 mg kg −1 w.w.) were similar with those reported for pork (0.078 mg kg −1 w.w.) and beef (0.064 mg kg −1 w.w.) [18], as well as for the liver (0.06 mg kg −1 w.w.), kidney (0.41 µg/g w.w.), heart and lungs (0.05 mg kg −1 w.w.) of roe deer [19].
On the other hand, our values for selenium were lower in comparison with those reported for wild boars and fishes. The concentration of selenium was 0.19 mg kg −1 w.w. in the liver and 1.20 mg kg −1 w.w in the kidney of wild boars from the northwest part of Poland [20]. In fish, the concentration of selenium varies normally between 0.2 and 0.9 mg kg −1 [21][22][23], reaching the highest values in species that are at the top of the food chain such as swordfish and tuna [24,25]. The level of selenium in the bluefin tuna (Thunnus thynnus) caught in the Mediterranean area varies from 0.6 ± 0.3 mg kg −1 (wild population) to 1.1 ± 0.9 mg kg −1 w.w. in farmed specimens [25]. Therefore, the studied populations of edible snails in Montenegro can be considered as a poor source of selenium for human consumption.
The most Cu was accumulated in the reproductive organ tissue. The highest Cu amount (98.4 ± 4.3) was found in snails from Nikšić and the lowest in snails from Biogradska Gora (53.9 ± 6.5 mg kg −1 d.w.). The BAF values, from 1.5 to 2.5, indicate a bioaccumulation of this metal. Dallinger and Wiesser [9] showed that the distribution of Cu in a snail tissue of H. pomatia is relatively even, without a clear preference for any organ.
The lowest concentration of Ni and Pb were found in the mantle tissue. The concentration of Ni and Pb were higher in the mantle tissue of snails from Malesija (0.42 ± 0.25 and 0.21 ± 0.04 for Ni and Pb, respectively) in comparison with snails from Nikšić (0.13 ± 0.03 and 0.01 ± 0.00 mg kg −1 d.w., for Ni and Pb, respectively). The lowest concentration of Al was found in the albumen gland tissue of snails from Nikšić (8.6 ± 1.2), and the lowest concentration of Li was found in the foot tissue of snails from Biogradska Gora (0.03 ± 0.05 mg kg −1 d.w.).
Our study confirmed previous studies that species of the genus Helix can be used as a bioindicator of a heavy metal accumulation, since they concentrate metals in their tissues [26]. The comparative analysis of the mean values shows that the accumulation of metals in the studied snail populations has the following order: Al > Fe > Mn > Zn > Cu > Cr > Li > Ni > Pb > Cd > Se > Hg. The accumulation of metals by the albumen gland followed the following order: Fe > Cu > Zn > Al > Mn > Cd > Pb > Ni > Cr > Se > Li > Hg, Al > Fe > Zn > Cu > Mn > Cd > Cr > Pb > Ni > Li > Se >Hg for the digestive organ tissue, Al > Fe > Cu > Zn > Mn > Cr > Pb > Cd > Ni > Se > Li >Hg for the foot tissue, Fe > Al > Zn > Mn > Cu > Cd > Pb > Cr > Ni > Li > Se > Hg for the hepatopancreas, Al > Fe > Zn > Cu > Mn > Cd > Cr > Ni > Li > Se > Pb > Hg for the mantle, Fe > Al > Cu > Zn > Mn > Cd > Cr > Ni > Pb > Se > Li > Hg for the reproductive organ tissue and Al > Fe > Cu > Mn > Zn > Ni > Li > Cr > Pb > Cd > Se for the shells. The conducted Pearson correlation analysis in our study showed statistically significant negative correlations between the Zn, Ni, Fe, Mn, Cr, Al, Li, Pb, Se and Hg concentrations and the examined tissue (Table 5). In the ground, the concentration of heavy metals follows a similar general pattern as in the snail organisms: Al>Fe>Mn>Zn>Cr> Li> Ni> Pb> Cu> Cd >Hg>Se. Positive correlations were found between the metal concentrations in the ground and the metal concentrations in the examined snail tissues, in particular (Table 6), (1) Nikšić-foot, digestive system, mantle and shell, (2) Malesija-hepatopancreas, reproductive system, foot, digestive system and mantle and (3) Biogradska Gora-reproductive system, foot, digestive system, mantle and shell. The canonical analysis performed on the metals variability in the examined tissue revealed three distinct clusters (Figure 1b). Most organs such as the reproductive system, albumen gland and foot were clustered together. On the other hand, the hepatopancreas and shell form separate clusters. The first discriminant function explains 40% of the total variance and it is dominated by the Fe content (R = 0.9998). The second discriminant function explains 25% of the total variance and it is dominated by the Zn content (R = 0.9878). The variables with the greatest ability for discrimination along the first discriminant function are Cd, Al and Fe, whereas Se and Ni are those that mostly discriminate snail tissues along the second discriminant function.
The results of our study showed that the shell accumulated the lowest metal amounts. This is in agreement with similar studies conducted on freshwater snails [27][28][29]. Some studies showed that shells play an important role as a sink for some heavy metals [12]. On the other hand, the hepatopancreas accumulates the highest metal content. This is consistent with other studies which revealed that the digestive gland is the tissue which accumulates the highest metal content [2,3,30]. The hepatopancreas and kidney have crucial roles in the detoxification process [12]. As can be seen from Table 4, the accumulation of Zn and Cd in the hepatopancreas are at least 3-fold higher than that in the snail foot. Nowakowaska et al. [12] showed that the hepatopancreatic accumulation of Zn, Cd and Mg was many times higher than that in the kidney and foot. Further, the concentration of Fe and Mn in the hepatopancreas are, respectively, higher that in the foot tissue. In the present study, the Fe concentration of the hepatopancreas was higher by 3.5-4.9-fold in H. secernenda and 1.7-fold in H. vladika than that in the foot tissue.

Toxic Elements Content in Edible Snails of Montenegro and Comparision with Legal Limit for Food
In our study, the values of toxic elements were compared with those established by the national and EU regulations (EU commission regulation No 1881/2006 of 19 December, 2006) ( Table 7). The general mean Pb concentration of the examined snail samples varied from 0.11 ± 0.01 (Malesija) to 0.86 ± 0.1 mg kg −1 w.w. (Biogradska Gora), which is below the allowable level for mollusc meat (1.5 mg kg −1 w.w.). The mean Hg concentration of the examined snail samples at each site was lower than the legal limit (0.5 mg kg −1 w.w.). The general mean Cd concentration of the examined snail tissues was 1.49 ± 0.14 mg kg −1 w.w. for the population of H. secernenda from Malesija and 1.565 ± 0.0 mg kg −1 w.w. for the population of H. vladika from Biogradska Gora. These values significantly exceeded the maximum allowable level established by EU regulations (1 mg kg −1 w.w.). The obtained values in our study were lower than the values for the Helix pomatia populations collected from the wilderness in Moldova, which for Pb ranged between 1.2018 and 48.0840, and between 2.0273 and 3.1281 mg kg −1 for Cd [7]. Hg (µg kg −1 ) 2.94 ± 0.00 3.6 ± 0.0 9.00 ± 0.00 500 (fish)

Conclusions
Our results showed that the concentration of the metals content varied among the snail tissues, which affected their nutritional benefit for human consumption. The hepatopancreas is often considered the tastiest part of the snail meat and, as shown in this study, tends to bioaccumulate selenium and to a lesser extent also zinc. The digestive tract tissue is a target organ for the bioaccumulation of selenium, whose level at the Malesija site reached a concentration of 0.587 mg k −1 d.w. Nevertheless, both the hepatopancreas and digestive tract tissues accumulate also toxic elements such as cadmium, whose concentrations in the present study exceeded the permissible limits established by the national and EU regulations. Therefore, the increased content of Cd in our study indicates a possible risk for human health and limits the use of the populations of Helix vladika and H. secernenda collected from the wilderness as a food source in human nutrition.