Next Article in Journal
A Review of Environment Effects on Nitrate Accumulation in Leafy Vegetables Grown in Controlled Environments
Previous Article in Journal
Nutrient Composition of Demersal, Pelagic, and Mesopelagic Fish Species Sampled Off the Coast of Bangladesh and Their Potential Contribution to Food and Nutrition Security—The EAF-Nansen Programme
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

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

by
Vesna Vukašinović-Pešić
1,*,
Bogumiła Pilarczyk
2,
Tymoteusz Miller
3,
Monika Rajkowska-Myśliwiec
4,
Joanna Podlasińska
5,
Agnieszka Tomza-Marciniak
2,
Nada Blagojević
1,
Nevzeta Trubljanin
6,
Andrzej Zawal
3 and
Vladimir Pešić
6
1
Faculty of Metallurgy and Technology, University of Montenegro, 81000 Podgorica, Montenegro
2
Department of Animal Reproduction Biotechnology and Environmental Hygiene, Faculty of Biotechnology and Animal Husbundry, West Pomeranian University of Technology, 71-415 Szczecin, Poland
3
Institute of Marine and Environmental Sciences, University of Szczecin, 71-415 Szczecin, Poland
4
Department of Toxicology, Dairy Technology and Food Storage, Faculty of Food Science and Fisheries, West Pomeranian University of Technology, 71-415 Szczecin, Poland
5
Department of Environmental Management, Faculty of Environmental Management and Agriculture, West Pomeranian University of Technology, 71-415 Szczecin, Poland
6
Department of Biology, University of Montenegro, 81000 Podgorica, Montenegro
*
Author to whom correspondence should be addressed.
Foods 2020, 9(6), 731; https://doi.org/10.3390/foods9060731
Submission received: 22 April 2020 / Revised: 26 May 2020 / Accepted: 29 May 2020 / Published: 3 June 2020
(This article belongs to the Section Food Analytical Methods)

Abstract

:
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.

1. 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).

2. Materials and Methods

2.1. 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) Nikšić, Dragovoljići village (42°46′19.48″ N, 19°2′30.65″ E), (2) Malesija, Skorać village (42°22′2.83″ N, 19°23′20.30″ E) and (3) Biogradska Gora (42°53′39.54″ N, 19°36′17.20″ E). 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).

2.2. 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 HNO3 (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% HNO3 (Pb, Cd) or 15% HNO3 (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 r2 ≥ 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]. The calculated LODs were as follows (μg/L): Al 50.0, Fe 4.0, Mn 1.4, Zn 1.8, Cu 5.4, Cr 8.92, Ni 5.76, Li 2.0, Pb 1.04 and Cd 0.098.
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 HNO3 (Chempur, Poland) at 230 °C for 180 min and in HClO4 (Chempur, Poland) at 310 °C for 20 min. Then, selenate (Se6+) was reduced to selenite (Se4+) 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.

2.3. 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 ƛ-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.

3. Results and Discussion

3.1. 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ć.

3.2. 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 the studied metals, with the exception of Ni, Al, Li and Pb, were found in the shells. The concentration of metals in the shells varied: for Zn, from 3.7 ± 0.1 (Malesija and Nikšić) to 4.2 ± 0.3 (Biogradska Gora); Fe, from 2.8 ± 0.8 (Biogradska Gora) to 26.8 ± 4.7 (Nikšić); Mn, from 2.3 ± 1.0 (Malesija) to 6.9 ± 1.3a (Biogradska Gora); Cr, from 0.26 ± 0.15 (Biogradska Gora) to 0.50 ± 0.17 (Malesija); and for Cu, from 7.8 ± 1.8 (Biogradska Gora) to 12.9 ± 9.4 mg kg−1 d.w. (Nikšić).
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.

3.3. 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].

4. 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.

Author Contributions

Sample preparation, V.V.-P., N.B., N.T., V.P. and A.Z.; Chemical analysis, B.P., M.R.-M., J.P. and A.T.-M.; Statistical analysis, T.M.; Writing—original draft preparation, V.P.; Writing—review and editing, V.V.P. and V.P. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Acknowledgments

This paper was produced as part of the “Eco-DNA” project financed by the Ministry of Science of Montenegro.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Mooney, B.D.; Nichols, P.I.; Elliot, N.G. Seafood the Good Food II: Oil Profiles for Further Australian Seafoods and Influencing Factors; CSIRO Division of Marine Research and Deakin, Hobart A.C.T., Fisheries Research & Development Corporation (FRDC): Canberra, Australia, 2002; Available online: http://hdl.handle.net/102.100.100/199022?index=1 (accessed on 24 February 2020).
  2. Drozd, Ł.; Ziomek, M.; Szkucik, K.; Paszkiewicz, W.; Maćkowiak-Dryka, M.; Bełkot, Z.; Gondek, M. Selenium, copper, and zinc concentrations in the raw and processed meat of edible land snails harvested in Poland. J. Vet. Res. 2017, 61, 293–298. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  3. Ćirić, J.; Cerić, O.; Marković, R.; Janjić, J.; Spirić, D.; Popović, M.; Pećanac, B.; Baltić, B.; Baltić, M.Ž. Seasonal distributions of heavy metal concentrations in different snail (Helix pomatia) tissues from an urban environment in Serbia. Environ. Sci. Pollut. Res. Int. 2018, 25, 33415–33422. [Google Scholar] [CrossRef] [PubMed]
  4. Gomot, A.; Pihan, F. Growing snails used as sentinels to evaluate terrestrial environment contamination by trace elements. Chemosphere 2000, 40, 275–284. [Google Scholar] [CrossRef]
  5. Regoli, F.; Gorbi, S.; Fattorini, D.; Tedesco, S.; Notti, A.; Machella, N.; Bocchetti, R.; Benedetti, M.; Piva, F. Use of the land snail Helix aspersa as sentinel organism for monitoring ecotoxicologic effects of urban pollution: An integrated approach. Environ. Health Perspect 2006, 114, 63–69. [Google Scholar] [CrossRef] [Green Version]
  6. Laskowski, R.; Hopkin, S.P. Effect of Zn, Cu, Pb, and Cd on fitness in snails (Helix aspersa). Ecotoxicol. Environ. Saf. 1996, 34, 59–69. [Google Scholar] [CrossRef] [Green Version]
  7. Toader-Williams, A.; Golubkina, N. Investigation upon the edible snail’s potential as source of selenium for human health and nutrition observing its food chemical contaminant risk factor with heavy metals. Bull. Uasvm Agric. 2009, 66, 495–499. [Google Scholar]
  8. Gheoca, V. Can heliciculture act as a tool for edible land snails’s ACT AS A TOOL FOR natural populations’ management in Romania? Manag. Sustain. Dev. Sibiuromania 2013, 5, 21–25. [Google Scholar]
  9. Dallinger, R.; Wieser, W. Patterns of acumulation, distribution and liberation of Zn, Cu, Cd, and Pb in different organs of the land snail Helix pomatia L. Comp. Biochem. Phys. C 1984, 79, 117–124. [Google Scholar] [CrossRef]
  10. Gomot, A.; Pihan, F. Comparison of the bioaccumulation capacities of copper and zinc in two snail subspecies (Helix). Ecotoxicol. Environ. Saf. 1997, 38, 85–94. [Google Scholar] [CrossRef] [PubMed]
  11. Menta, C.; Parisi, V. Metal concentrations in Helix pomatia, Helix aspersa, and Arion rufus: A comparative study. Environ. Pollut. 2001, 115, 205–208. [Google Scholar] [CrossRef]
  12. Nowakowska, A.; Laciak, T.; Caputa, M. Heavy metals accumulation and antioxidant defence system in Helix pomatia (Pulmonata: Helicidae). Molluscan Res. 2012, 32, 16–20. [Google Scholar]
  13. Wenzl, T.; Haedrich, J.; Schaechtele, A.; Robouch, P.; Stroka, J. Guidance Document on the Estimation of LOD and LOQ for Measurements in the Field of Contaminants in Feed and Food; EUR 28099; Publications Office of the European Union: Luxembourg, 2016; ISBN 9789279617683. [Google Scholar] [CrossRef]
  14. Neubert, E. Revision of Helix Linnaeus, 1758 in its Eastern Mediterranean Distribution Area, and Reassignment of Helix Godetiana Kobelt, 1878 to Maltzanella Hesse, 1917 (Gastropoda, Pulmonata, Helicidae); Contributions to Natural History; Naturhistorisches Museum: Bern, Switzerland, 2014; Volume 26, pp. 1–200. [Google Scholar]
  15. Radwan, M.A.; El-Gendy, K.S.; Gad, A.F. Biomarkers of oxidative stress in the land snail, Theba pisana for assessing ecotoxicological effects of urban metal pollution. Chemosphere 2010, 79, 40–46. [Google Scholar] [CrossRef] [PubMed]
  16. Mikac, N.; Kwokal, Z.; Martincic, D.; Branica, M. Uptake of mercury species by transplanted mussels Mytilus galloprovincialis under estuarine conditions (Krka river estuary). Sci. Total Environ. 1996, 184, 173–182. [Google Scholar] [CrossRef]
  17. Truzzi, C.; Annibaldi, A.; Girolametti, F.; Giovannini, L.; Riolo, P.; Ruschioni, S.; Olivotto, I.; Illuminati, S. A Chemically Safe Way to Produce Insect Biomass for Possible Application in Feed and Food Production. Int. J. Environ. Res. Public Health. 2020, 17, 2121. [Google Scholar] [CrossRef] [Green Version]
  18. Pilarczyk, B.; Tomza-Marciniak, A.; Mituniewicz-Małek, A.; Wieczorek-Dąbrowska, M.; Pilarczyk, R.; Wójcik, J.; Balicka-Ramisz, A.; Bąkowska, M.; Dmytrów, I. Selenium content in selected products of animal origin and estimation of the degree of cover daily Se requirement in Poland. Int. J. Food Sci. Technol. 2010, 45, 186–191. [Google Scholar] [CrossRef]
  19. Pilarczyk, B.; Tomza-Marciniak, A.; Pilarczyk, R.; Hendzel, D.; Błaszczyk, B.; Bąkowska, M. Tissue distribution of selenium and effect of season and age on selenium content in roe deer from northwestern Poland. Biol. Trace Elem. Res. 2011, 140, 299–307. [Google Scholar] [CrossRef]
  20. Pilarczyk, B.; Hendzel, D.; Pilarczyk, R.; Tomza-Marciniak, A.; Błaszczyk, B.; Wieczorek-Dąbrowska, M.; Bąkowska, M.; Adamowicz, E.; Bujak, T. Liver and kidney concentrations of selenium in wild boars (Sus scrofa) from northwestern Poland. Eur. J. Wildl. Res. 2010, 56, 797–802. [Google Scholar] [CrossRef]
  21. Kaneko, J.J.; Ralston, N.V.C. Selenium and Mercury in Pelagic Fish in the Central North Pacific near Hawaii. Biol. Trace Elem. Res. 2007, 119, 242–254. [Google Scholar] [CrossRef]
  22. Olmedo, P.; Hernández, A.F.; Pla, A.; Femia, P.; Navas-Acien, A.; Gil, F. Determination of Essential Elements (Copper, Manganese, Selenium and Zinc) in Fish and Shellfish Samples. Risk and Nutritional Assessment and Mercury-Selenium Balance. Food Chem. Toxicol. 2013, 62, 299–307. [Google Scholar] [CrossRef]
  23. Cabañero, A.I.; Carvalho, C.; Madrid, Y.; Batoréu, C.; Cámara, C. Quantification and Speciation of Mercury and Selenium in Fish Samples of High Consumption in Spain and Portugal. Biol. Trace Elem. Res. 2005, 103, 17–36. [Google Scholar] [CrossRef]
  24. Burger, J.; Jeitner, C.; Gochfeld, M. Locational Differences in Mercury and Selenium Levels in 19 Species of Saltwater Fish from New Jersey. J. Toxicol. Environ. Health Part A 2011, 74, 863–874. [Google Scholar] [CrossRef] [PubMed]
  25. Annibaldi, A.; Truzzi, C.; Carnevali, O.; Pignalosa, P.; Api, M.; Scarponi, G.; Illuminati, S. Determination of Hg in farmed and wild Atlantic Bluefin Tuna (Thunnus thynnus L.) muscle. Molecules 2019, 24, 1273. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  26. Nica, D.V.; Bura, M.; Gergen, I.; Harmanescu, M.; Bordean, D.M. Bioaccumulative and conchological assessment of heavy metal transfer in a soil-plant-snail food chain. Chem. Cent. J. 2012, 6, 55. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  27. Vukašinović-Pešić, V.; Blagojević, N.; Vukanović, S.; Savić, A.; Pešić, V. Heavy metal concentrations in different tissues of the snail Viviparus mamillatus (Küster, 1852) from lacustrine and riverine environments in Montenegro. Turk. J. Fish. Aquat. Sci. 2017, 17, 557–563. [Google Scholar] [CrossRef]
  28. Blagojević, N.; Vukašinović-Pešić, V.; Grudić, V.; Pešić, V. The endemic freshwater snails as an environmental indicator of metal pollution of the Zeta river, Montenegro. J. Environ. Prot. Ecol. 2014, 15, 210–216. [Google Scholar]
  29. Vukašinović-Pešić, V.; Blagojević, N. Metal Pollution: Evidences from Plants, Aquatic Invertebrates and Fish from Lake Skadar. In Skadar/Shkodra Lake Environment; The Handbook of Environmental, Chemistry; Pešić, V., Karaman, G., Kostianoy, A., Eds.; Springer: Cham, Switzerland, 2018; Volume 80, pp. 141–152. [Google Scholar]
  30. Coughtrey, P.J.; Martin, M.H. The distribution of Pb, Zn, Cd, and Cu within the pulmonate mollusc Helix aspersa Müller. Oecologia 1976, 23, 315–322. [Google Scholar] [CrossRef]
Figure 1. Canonical discriminant score plot of (a) examined sampling localities, and (b) examined snail tissues.
Figure 1. Canonical discriminant score plot of (a) examined sampling localities, and (b) examined snail tissues.
Foods 09 00731 g001
Table 1. Mean values (mean ± standard error) of water contents and other allometric parameters of the studied snail population.
Table 1. Mean values (mean ± standard error) of water contents and other allometric parameters of the studied snail population.
SiteShell Height (cm)Shell Width (cm)Wet Weight (g)Dry Weight (g)Water Content (%)
Nikšić4.4 ± 0.43.9 ± 0.322.0 ± 4.09.2 ± 2.158.3 ± 8.4
Malesija4.7 ± 0.24.0 ± 0.129.8 ± 3.810.6 ± 1.664.2 ± 5.6
Biogradska Gora5.3 ± 0.24.6 ± 0.538.3 ± 5.117.4 ± 4.555.1 ± 6.4
Table 2. Operating parameters for the microwave oven (CEM, MDS 2000).
Table 2. Operating parameters for the microwave oven (CEM, MDS 2000).
ParameterStage
12345
power (%)100100100100100
pressure (PSI)204085135175
time (min)51010105
Table 3. Recovery (%) rate (mean ± SD) from reference materials MODAS-2 and MODAS-3.
Table 3. Recovery (%) rate (mean ± SD) from reference materials MODAS-2 and MODAS-3.
Reference MaterialZnNiFeMnCrCuAlLiPbCd
Modas 3 (M-3 Her Tis)99.6 ± 3.193.1 ± 3.695.2 ± 6.1101.0 ± 6.9106.2 ± 7.387.5 ± 5.4-108.6 ± 3.695.9 ± 6.293.8 ± 5.9
Modas 2 (M2 BotSed)113.0 ± 3.387.1 ± 2.898.2 ± 5.5107.0 ± 5,1108.2 ± 3.599.4 ± 4.487.0 ± 6.5102.6 ± 3.183.5 ± 4.297.2 ± 2.9
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.
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.
SiteSample Zn
mg kg−1
Ni
mg kg−1
Fe
mg kg−1
Mn
mg kg−1
Cr
mg kg−1
Cu
mg kg−1
Al
mg kg−1
Li
mg kg−1
Pb
mg kg−1
Cd
mg kg−1
Se
mg kg−1
Hg
µg kg−1
Nikšić (Helix secernenda)GrMean ± SD248.5aA * ± 22.745.41bbA * ± 1.7522532.9abA * ± 627.91074.9aA * ± 29.685.54aA * ± 2.7942.7 ± 5.132212.8A * ± 1483.635.08A * ± 1.0930.92A * ± 0.361.24a * ± 0.040.06abA * ± 0.01115.7ab * ± 6.0
HpMean ± SD387.0aA * ± 16.90.82b * ± 0.12205.7ab * ± 8.7405.1aA * ± 11.03.80a * ± 0.1945.7 ± 5.9214.9 ± 3.50.18 ± 0.041.96 ± 0.391.43a * ± 0.030.19ab * ± 0.0120.3a * ± 0.8
BCF1.56 ± 0.740.02 ± 0.070.01 ± 0.010.38 ± 0.370.04 ± 0.071.07 ± 1.160.00 ± 0.000.00 ± 0.040.06 ± 1.081.15 ± 0.753.35 ± 0.640.17 ± 0.13
RsMean ± SD41.4a ± 2.50.71b ± 0.3369.4 ± 12.614.2 ± 1.20.86 ± 0.1298.4A ± 4.362.0 ± 1.90.04 ± 0.000.90 ± 0.121.10 ± 0.150.18 ± 0.013.7 ± 0.3
BCF0.17 ± 0.110.02 ± 0.190.00 ± 0.020.01 ± 0.040.01 ± 0.042.30 ± 0.840.00 ± 0.000.00 ± 0.000.03 ± 0.330.89 ± 3.753.09 ± 1.360.03 ± 0.05
AgMean ± SD33.2aA * ± 3.51.21a * ± 0.0326.2a * ± 5.017.7a * ± 1.50.36a * ± 0.0570.4 ± 11.58.6 ± 1.20.12a * ± 0.000.25a * ± 0.060.57 ± 0.080.11 ± 0.004.0 ± 0.00
BCF0.13 ± 0.150.03 ± 0.020.00 ± 0.010.02 ± 0.050.00 ± 0.021.64 ± 2.250.00 ± 0.000.00 ± 0.000.01 ± 0.170.46 ± 21.88 ± 0.450.03 ± 0.00
FMean ± SD39.5 ± 2.40.66a * ± 0.1641.7a * ± 4.67.6a * ± 0.21.29 ± 0.1053.6a * ± 2.5138.4 ± 29.10.09a * ± 0.050.21 ± 0.130.48 ± 0.140.11 ± 0.005.7 ± 0.1
BCF0.16 ± 0.100.01 ± 0.090.00 ± 0.010.01 ± 0.010.01 ± 0.031.25 ± 0.490.00 ± 0.020.00 ± 0.040.01 ± 0.360.38 ± 3.51.95 ± 0.270.05 ± 0.02
DsMean ± SD63.6 ± 5.00.64 ± 0.17181.3 ± 13.214.7 ± 0.52.71 ± 0.4268.3 ± 16.481.2 ± 11.20.31 ± 0.270.94a ** ± 0.044.81a ** ± 1.030.19 ± 0.019.2a ** ± 2.0
BCF0.25 ± 0.220.01 ± 0.100.01 ± 0.020.01 ± 0.020.03 ± 0.151.6 ± 3.210.00 ± 0.010.01 ± 0.250.03 ± 0.113.88 ± 25.753.26 ± 1.270.08 ± 0.33
MMean ± SD78.2 ± 2.80.13 ± 0.0391.8 ± 25.949.6 ± 5.60.82 ± 0.2651.4 ± 0.2116.7 ± 0.40.31 ± 0.060.01 ± 0.002.72 ± 0.090.10 ± 0.005.4 ± 0. 3
BCF0.31 ± 0.120.00 ± 0.020.00 ± 0.040.05 ± 0.190.01 ± 0.091.20 ± 0.040.00 ± 0.000.01 ± 0.050.00 ± 0.002.19 ± 2.251.79 ± 0.540.05 ± 0.05
SMean ± SD3.7aA * ± 0.10.96a * ± 0.1726.8a * ± 4.75.7 ± 0.70.50a * ± 0.1712.9A * ± 9.4502.9 ± 29.60.41 ± 0.010.12 ± 0.110.08aA * ± 0.020.01aA * ± 0.000.8a * ± 0.00
BCF0.01 ± 0.000.02 ± 0.100.00 ± 0.000.00 ± 0.020.01 ± 0.060.30 ± 1.840.01 ± 0.020.01 ± 0.010.00 ± 0.300.06 ± 0.50.24 ± 0.180.00 ± 0.00
Malesija (Helix secernenda)GrMean ± SD108.8aA * ± 7.959.04A ± 5.4127068.7aA * ± 1977.4758.4A * ± 18.0110.56A * ± 11.4227.3a * ± 0.950288.1aA * ± 3197.351.47aA * ± 4.6028.32aA * ± 1.290.65 ± 0.020.06aA * ± 0.01121 ± 11
HpMean ± SD253.4A * ± 18.52.12 ± 0.36463.9 ± 23.832.4 ± 0.33.41 ± 0.5840.4 ± 0.5484.0 ± 4.50.72 ± 0.290.99 ± 0.2916.38A * ± 1.460.27A * ± 0.0224 ± 0.4
BCF2.33 ± 2.340.03 ± 0.070.02 ± 0.010.04 ± 0.020.03 ± 0.051.48 ± 0.550.01 ± 0.000.01 ± 0.060.03 ± 0.2225.2 ± 734.44 ± 2.220.19 ± 0.04
RsMean ± SD61.3a * ± 1.10.60b * ± 0.3289.5ab * ± 6.010.4ab * ± 0.90.92ab * ± 0.0167.1a * ± 4.296.0ab * ± 8.10.04 ± 0.000.07a * ± 0.042.99aA * ± 0.240.24abA * ± 0.025 ± 0.3
BCF0.56 ± 0.140.01 ± 0.060.00 ± 0.000.01 ± 0.050.01 ± 0.002.45 ± 4.670.00 ± 0.000.00 ± 0.000.00 ± 0.034.6 ± 123.98 ± 2.670.04 ± 0.03
AgMean ± SD45.7a * ± 1.60.24a * ± 0.3299.5a * ± 10.510.9a * ± 0.70.63a * ± 0.2141.1 ± 2.742.0 ± 11.80.08 ± 0.060.08a * ± 0.021.69A * ± 0.080.24abA * ± 0.013.5 ± 0.0
BCF0.42 ± 0.200.00 ± 0.060.00 ± 0.000.01 ± 0.040.00 ± 0.021.50 ± 30.00 ± 0.000.00 ± 0.010.00 ± 0.012.6 ± 43.93 ± 1.550.03 ± 0.00
FMean ± SD56.9 ± 0.40.28a * ± 0.06134.3a * ± 15.48.9a * ± 0.11.18a * ± 0.4475.8 ± 15.2130.9 ± 4.20.12 ± 0.000.18a * ± 0.030.91 ± 0.090.18 ± 0.016.7 ± 0.2
BCF0.52 ± 0.050.00 ± 0.010.00 ± 0.010.01 ± 0.000.01 ± 0.042.77 ± 16.90.00 ± 0.000.00 ± 0.000.01 ± 0.021.4 ± 4.52.96 ± 1.220.05 ± 0.01
DsMean ± SD41.1a * ± 2.51.44a * ± 0.13805.0a * ± 67.626.6a * ± 2.14.84 ± 0.4336.9 ± 0.71411.1a * ± 172.01.26a * ± 0.070.42 ± 0.065.24 ± 0.510.59A ± 0.0113.3 ± 0. 4
BCF0.38 ± 0.320.02 ± 0.020.03 ± 0.030.03 ± 0.120.04 ± 0.041.35 ± 0.770.03 ± 0.050.02 ± 0.010.01 ± 0.058.06 ± 25.59.62 ± 1.670.11 ± 0.04
MMean ± SD62.4 ± 12.70.42a * ± 0.25224.8 ± 14.529.0 ± 5.41.24a * ± 0.1064.3 ± 3.9432.2 ± 28.90.28a * ± 0.000.21 ± 0.042.16 ± 0.750.17 ± 0.026.6 ± 0.6
BCF0.57 ± 1.600.01 ± 0.050.01 ± 0.010.04 ± 0.30.01 ± 0.012.35 ± 4.330.01 ± 0.010.00 ± 0.000.01 ± 0.033.32 ± 37.52.83 ± 20.05 ± 0.05
SMean ± SD3.7aA ** ± 0.20.95a ** ± 0.4425.1a ** ± 10.42.3a ** ± 1.00.45a ** ± 0.099.0aA ** ± 4.6512.4a ** ± 46.80.43a ** ± 0.050.16a ** ± 0.090.10aA ** ± 0.020.01a ** ± 0.001.9A ** ± 0.1
BCF0.03 ± 0.020.02 ± 0.080.00 ± 0.000.00 ± 0.050.00 ± 0.010.33 ± 5.110.01 ± 0.010.01 ± 0.010.00 ± 0.070.15 ± 10.19 ± 0.440.01 ± 0.01
Biogradska Gora (Helix vladika)GrMean ± SD153.7A * ± 15.441.57A * ± 6.6829053.7A ** ± 1499.12895.7A ± 172.187.63A * ± 1.5936.7 ± 1.937653.1A * ± 3416.496.38A * ± 16.4949.72A * ± 5.530.99 ± 0.040.05 ± 0.01260.0 ± 1.1
HpMean ± SD151.6A * ± 6.90.95a * ± 0.07149.3 ± 16.4350.6A * ± 1.90.95 ± 0.1633.2 ± 11.5100.1 ± 2.60.10 ± 0.005.46 ± 0.1811.28 ± 0.170.11 ± 0.0253.9A * ± 2.1
BCF0.99 ± 0.440.02 ± 0.010.00 ± 0.010.12 ± 0.010.01 ± 0.100.90 ± 6.050.00 ± 0.000.00 ± 0.000.11 ± 0.0311.39 ± 4.252.35 ± 2.290.21 ± 1.91
RsMean ± SD61.0 ± 4.41.28 ± 0.3587.7 ± 3.213.6 ± 0.51.20 ± 0.2853.9 ± 6.568.4 ± 10.50.08 ± 0.000.12 ± 0.002.00 ± 0.070.13 ± 0.016.5A * ± 0.5
BCF0.39 ± 0.280.03 ± 0.050.00 ± 0.000.00 ± 0.000.01 ± 0.181.47 ± 3.420.00 ± 0.000.00 ± 0.000.00 ± 0.002.02 ± 1.752.73 ± 1.710.02 ± 0.45
AgMean ± SD51.7 ± 1.40.91 ± 0.2864.0 ± 0.323.8 ± 0.10.51 ± 0.3353.5 ± 1.616.6 ± 0.60.08 ± 0.002.78 ± 1.131.66 ± 0.260.10 ± 0.005.3 ± 0.3
BCF0.33 ± 0.090.02 ± 0.040.00 ± 0.000.01 ± 0.000.01 ± 0.211.46 ± 0.840.00 ± 0.000.00 ± 0.000.05 ± 0.201.68 ± 6.52.1 ± 0.710.02 ± 0.27
FMean ± SD52.9a * ± 1.80.52b * ± 0.1187.1 ± 7.716.0 ± 1.11.32ab * ± 0.2649.0 ± 1.3106.5 ± 12.00.03ab * ± 0.051.76 ± 0.360.65 ± 0.040.11 ± 0.015.4 ± 0.2
BCF0.34 ± 0.110.01 ± 0.010.00 ± 0.000.00 ± 0.010.01 ± 0.161.33 ± 0.680.00 ± 0.000.00 ± 0.000.03 ± 0.060.65 ± 12.2 ± 1.280.02 ± 0.18
DsMean ± SD63.7a * ± 2.70.20a * ± 0.10111.9a * ± 11.148.7a * ± 5.91.15a * ± 0.5842.4 ± 3.573.6 ± 8.00.03a * ± 0.053.09a * ± 0.036.10 ± 0.860.11 ± 0.0242.1 ± 33.0
BCF0.41 ± 0.170.00 ± 0.010.00 ± 0.010.02 ± 0.030.01 ± 0.361.15 ± 1.840.00 ± 0.000.00 ± 0.000.06 ± 0.006.16 ± 21.52.16 ± 2.850.16 ± 30
MMean ± SD44.5 ± 2.00.22 ± 0.1186.8 ± 3.114.9 ± 2.41.00 ± 0.1568.3 ± 3.6207.6 ± 21.60.11 ± 0.120.06 ± 0.002.61 ± 0.110.08 ± 0.016.8 ± 0.1
BCF0.29 ± 0.130.00 ± 0.020.00 ± 0.000.00 ± 0.010.01 ± 0.101.86 ± 1.890.00 ± 0.010.00 ± 0.010.00 ± 0.002.63 ± 2.751.55 ± 2.140.03 ± 0.09
SMean ± SD4.2A ** ± 0.31.15a ** ± 0.242.8a ** ± 0.86.9a ** ± 1.30.26 ± 0.157.8A ** ± 1.8445.6 ± 17.20.44 ± 0.120.16 ± 0.040.03aA ** ± 0.010.01A ** ± 0.001.5 ± 0.0
BCF0.03 ± 0.020.03 ± 0.030.00 ± 0.000.00 ± 0.010.00 ± 0.100.21 ± 1.950.01 ± 0.000.00 ± 0.000.00 ± 0.000.03 ± 0.250.14 ± 0.280.01 ± 0.00
Table 5. Pearson correlations between the metal concentrations and snail tissue and sampling localities, respectively. * p ≤ 0.05.
Table 5. Pearson correlations between the metal concentrations and snail tissue and sampling localities, respectively. * p ≤ 0.05.
VariableZnNiFeMnCrCuAlLiPbCdSeHg
Tissue−0.662 *−0.586 *−0.590 *−0.541 *−0.485 *−0.222−0.572 *−0.536 *−0.598 *−0.244−0.344 *−0.591 *
Locality−0.167−0.0270.0160.120−0.084−0.1810.0010.1070.0920.177−0.114−0.031
Table 6. Pearson correlation between the selected organs and sampling sites. Abbreviations: Gr—ground, Hp—hepatopancreas, Rs—reproductive system, Ag—albumen gland, F—foot, Ds—digestive system, M—mantle, S—shell. * p ≤ 0.05.
Table 6. Pearson correlation between the selected organs and sampling sites. Abbreviations: Gr—ground, Hp—hepatopancreas, Rs—reproductive system, Ag—albumen gland, F—foot, Ds—digestive system, M—mantle, S—shell. * p ≤ 0.05.
SiteHpRsAgFDsMS
Nikšić0.2330.3650.0010.798 *0.865 *0.865 *0.757 *
Malesija0.885 *0.739 *0.5550.859 *0.999 *0.891 *0.926 *
Biogradska Gora0.3650.821*0.3950.950*0.736 *0.902 *0.733 *
Table 7. Mean concentrations of metals in the studied snails from the three sampling sites of Montenegro and the legal limits (guidelines EU No 1881/2006), referred to as wet weight.
Table 7. Mean concentrations of metals in the studied snails from the three sampling sites of Montenegro and the legal limits (guidelines EU No 1881/2006), referred to as wet weight.
MetalsInvestigated SitesLegal Limit
(1881/2006/EU)
NikšićMalesijaBiogradska Gora
Pb (mg kg−1)0.26 ± 0.000.11 ± 0.010.86 ± 0.121.5 (molluscs)
Cd (mg kg−1)0.67 ± 0.051.49 ± 0.141.56 ± 0.091.0 (molluscs)
Hg (µg kg−1)2.94 ± 0.003.6 ± 0.09.00 ± 0.00500 (fish)

Share and Cite

MDPI and ACS Style

Vukašinović-Pešić, V.; Pilarczyk, B.; Miller, T.; Rajkowska-Myśliwiec, M.; Podlasińska, J.; Tomza-Marciniak, A.; Blagojević, N.; Trubljanin, N.; Zawal, A.; Pešić, V. Toxic Elements and Mineral Content of Different Tissues of Endemic Edible Snails (Helix vladika and H. secernenda) of Montenegro. Foods 2020, 9, 731. https://doi.org/10.3390/foods9060731

AMA Style

Vukašinović-Pešić V, Pilarczyk B, Miller T, Rajkowska-Myśliwiec M, Podlasińska J, Tomza-Marciniak A, Blagojević N, Trubljanin N, Zawal A, Pešić V. Toxic Elements and Mineral Content of Different Tissues of Endemic Edible Snails (Helix vladika and H. secernenda) of Montenegro. Foods. 2020; 9(6):731. https://doi.org/10.3390/foods9060731

Chicago/Turabian Style

Vukašinović-Pešić, Vesna, Bogumiła Pilarczyk, Tymoteusz Miller, Monika Rajkowska-Myśliwiec, Joanna Podlasińska, Agnieszka Tomza-Marciniak, Nada Blagojević, Nevzeta Trubljanin, Andrzej Zawal, and Vladimir Pešić. 2020. "Toxic Elements and Mineral Content of Different Tissues of Endemic Edible Snails (Helix vladika and H. secernenda) of Montenegro" Foods 9, no. 6: 731. https://doi.org/10.3390/foods9060731

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop