Determination of the Concentration of Heavy Metals in Artisanal Cheeses Produced in the Mexican States of Tabasco and Chiapas

Cheese consumption provides humans with minerals, proteins, carbohydrates, and vitamins. In Mexico, several cheese varieties are produced, each with its texture, scent, and flavor. The artisanal cheeses made in the states of Tabasco and Chiapas—including, among others, the varieties named crema (cream), doble crema (double cream), oaxaca, panela, fresco, bola, poro, cotija, and asadero—have a high demand in the domestic and foreign markets. The intensification of anthropic activity in these states causes an increased emission to the environment of contaminants like heavy metals, which could reach human foodstuffs through the food chains. In particular, heavy metal contents in cheeses consumed daily by these states’ local populations might represent a public health risk. Because of that, our objectives in this work were to determine the concentrations of lead, cadmium, nickel, copper, zinc, and iron in artisanal cheeses produced in the states of Tabasco and Chiapas and to determine the values of the hazard quotient (HQ), total hazard quotient (THQ), and cancer risk total (CRT) for adult and young men and women. The results of our analyses of cheese samples from the states of Tabasco and Chiapas showed that the average concentrations (mg kg−1) of cadmium (0.0023 ± 0.002, 0.0023 ± 0.002 mg kg−1, respectively, for each state), lead (0.0047 ± 0.00, 0.0051 ± 0.002), nickel (0.0039 ± 0.0046, 0.0031 ± 0.0039), copper (0.0199 ± 0.021, 0.0202 ± 0.022), zinc (0.1611 ± 0.18, 0.194 ± 0.21), and iron (61.84 ± 4.23, 65.76 ± 6.61 mg kg−1), the first three values lower than the limits established by the FAO/WHO and Codex Alimentarius. The value of THQ that we obtained was less than one, and that of CRT was within the limits established by the US-EPA, which means that the consumption of artisanal cheeses from Tabasco and Chiapas by humans does not imply a risk of disease or cancer.


Introduction
The large variety of cheeses produced in different geographic regions represents healthy foodstuffs for human consumption because they contain high concentrations of proteins, fats, vitamins, and minerals-including macro (Na, K, Ca) and micro (Zn, Cu, Fe) elements, each having their own organoleptic and physicochemical properties [1].In each geographic region, the natural grasslands, fodders, and grains fed to cattle might represent sources of contaminants of cheeses.In particular, heavy metals bioaccumulate in the muscles, organs, and bones of animals, thus contaminating dairy products like milk, whey, cream, butter, and cheeses whose consumption might be a risk to human health.The intake of foodstuffs containing heavy metals has become a public health issue because of their carcinogenic, mutagenic, and cytotoxic effects [2,3], because of which the Food and Agriculture Organization (FAO)/World Health Organization (WHO) [4], Codex Stan

Statistical Analysis
In this study, ANOVA with Tukey and paired Student's t-test were performed using the InfoStat/L version 2020 software to know if the differences between the moisture content and the concentrations of Pb, Cd, Ni, Cu, Zn, and Fe observed in the 88 cheese samples were significant (Tables 3-5).

Lead
The International Programme of Chemical Safety (IPCS) [11] and the International Agency for Research on Cancer (IARC) [12] classify lead as a probable human carcinogen (2A).The high toxicity of lead might cause diseases and, at high concentrations, several types of cancer [13,14].Our results showed that the concentration of lead in the 88 artisanal cheese samples we analyzed was below the 0.02 mg kg −1 and 2 mg kg −1 maximum limits for dairy products for human consumption established by the European Commission EC. (2001) [15] and FAO/WHO [4], respectively.The 44 analyzed samples from Tabasco had an average lead concentration of 0.0041, with a minimum of 0.0002 and a maximum of 0.0125 mg kg −1 .By artisanal cheese variety, the average, minimum, and maximum values of lead concentration in mg kg −1 were as follows: 0.0056, 0.0004, and 0.0125 for cream or soup and double crema; 0.0045, 0,0028, and 0.0084 for oaxaca; 0.0037, 0.0005, and 0.0068 for fresco; and 0.0041, 0.0031, and 0.0047 for panela.For the 44 samples from Chiapas, we found average, minimum, and maximum values of lead concentration in mg kg −1 of 0.0051, 0.0003, and 0.0089.For the latter state's artisanal cheese varieties, the average, minimum, and maximum values of lead concentration in mg kg −1 were as follows: 0.0062, 0.0029, and 0.0089 for crema; 0.0054, 0.0039, and 0.0088 for cotija; 0.0044, 0.0003, and 0.0083 for oaxaca; and 0.0028, 0.0005, and 0.0078 for panela.
Figure 1 shows the distribution of the lead concentration values that we observed in the samples of artisanal cheese varieties crema, oaxaca, panela, cotija, and fresco we collected in the states of Tabasco and Chiapas.The maximum lead concentration of 0.0125 mg kg −1 that we observed corresponded to a sample of crema artisanal cheese from the municipality of Huimanguillo, Tabasco could have been due to pollution during the manufacturing process [16].In general, the varieties cotija and crema had higher lead concentrations than the varieties oaxaca, fresh, and panela, which could be explained by the former varieties having a less content of water by weight (38.699% and 45.264%, respectively) than the latter (59.546, 60.799, and 61.325%) see Table 3.In addition, the Tukey analysis showed that there are significant differences in moisture content between cotija and cream cheeses with respect to panela, oaxaca, and fresh cheeses (Table 4).The results of the statistical analysis of the lead concentrations in the artisanal cheeses from the state of Tabasco and Chiapas gave a value of TC = −0.661for the Student's t-test for paired data, which is within the range of critical value (CV between −2.672 and +2.672), indicating that there are no significant differences in the lead content in the cheese samples, which confirms the results from the Tukey test (Table 5a), which could mean that in both processes the contamination sources were the same.and maximum values of lead concentration in mg kg −1 were as follows: 0.0062, 0.0029, and 0.0089 for crema; 0.0054, 0.0039, and 0.0088 for cotija; 0.0044, 0.0003, and 0.0083 for oaxaca; and 0.0028, 0.0005, and 0.0078 for panela.
Figure 1 shows the distribution of the lead concentration values that we observed in the samples of artisanal cheese varieties crema, oaxaca, panela, cotija, and fresco we collected in the states of Tabasco and Chiapas.The maximum lead concentration of 0.0125 mg kg −1 that we observed corresponded to a sample of crema artisanal cheese from the municipality of Huimanguillo, Tabasco could have been due to pollution during the manufacturing process [16].In general, the varieties cotija and crema had higher lead concentrations than the varieties oaxaca, fresh, and panela, which could be explained by the former varieties having a less content of water by weight (38.699% and 45.264%, respectively) than the la er (59.546, 60.799, and 61.325%) see Table 3.In addition, the Tukey analysis showed that there are significant differences in moisture content between cotija and cream cheeses with respect to panela, oaxaca, and fresh cheeses (Table 4).The results of the statistical analysis of the lead concentrations in the artisanal cheeses from the state of Tabasco and Chiapas gave a value of TC = −0.661for the Student's t-test for paired data, which is within the range of critical value (CV between −2.672 and +2.672), indicating that there are no significant differences in the lead content in the cheese samples, which confirms the results from the Tukey test (Table 5a), which could mean that in both processes the contamination sources were the same.In a study of fresco cheese from the Mexican state of Puebla, by Benítez found an average lead concentration of 2.96 mg kg −1 [16].Castro analyzed the lead concentration values in artisanal cheeses made in Santa Ana Xalmimilulco, Puebla, reporting averages of 0.11 ± 0.04 mg kg −1 for ranchero and 0.05 ± 0.03 mg kg −1 for Oaxaca cheese varieties, the authors attributed the presence of lead in the cheeses to the alfalfa fed to cattle being irrigation with polluted water from the Atoyac River [17].In studies made in Turkey, Yüzbaşi found lead concentrations between 0.0364 and 0.251 mg kg −1 in the artisanal Kasar cheese samples from Ankara [18], and Bakircioglu reported lead concentrations ranging from 0.60 ± 0.17 and 0.48 ± 17 mg kg −1 for cream and white cheese made in Edirne [19].
Moreno analyzed 57 varieties of cheese from different regions of Spain and found that their lead contents were between 5.0 ± 0.01 and 110.0 ± 8.2 µg kg −1 [1]; while in Italy, Lante reported a lead concentration of 0.06 mg kg −1 in crescenza and squacquerone cheeses [20].

Cadmium
The IARC and IPCS [11,12] determined that the ingestion of even a small amount of cadmium is highly toxic for humans; therefore, determining its presence in food is essential.The results of our determination of the content of Cd in cheese samples from Tabasco show average concentrations between 0.0023 ± 0.0025 and 0.0056 mg kg −1 .By cheese variety, the average, minimum, and maximum values of cadmium concentration in mg kg −1 were as follows: 0.0045, 0.0019, and 0.0156 for crema; 0.0009, 0.00, and 0.0031 for oaxaca; 0.0011, 0.0002, and 0.0037 for fresco; and 0.0019, 0.0001, and 0.0075 for panela.For the analyzed cheese samples from Chiapas, we found an average Cd concentration of 0.0023 ± 0.0025 with a minimum of 0.00 and a maximum of 0.0055 mg kg −1 .By cheese variety, the average, minimum, and maximum values of cadmium concentration in mg kg −1 were as follows: 0.0043, 0.0029, and 0.0077 for crema; 0.0042, 0.0032, and 0.0062 for cotija; 0.0019, 0.00, and 0.0046 for oaxaca; and 0.0014, 0.0002, and 0.0034 for panela.
The distribution of the cadmium concentrations we found in the crema, oaxaca, panela, cotija, and fresco cheese samples from Tabasco and Chiapas shows that the highest Cd concentrations correspond to crema and cotija cheeses, and the lowest Cd concentrations, to the panela, oaxaca, and fresco samples (Figure 2).The (Table 3) latter result may have been due to the lower water content of the crema (<45.264% by weight) and cotija samples (<38.699%)relative to that in the panela, oaxaca, and fresco (>59%).The results obtained from the statistical analysis of cadmium concentrations in Tabasco and Chiapas cheeses showed that there are no significant differences between cheese samples since the value of the paired t-Student test (TC = −1.318) is within the range of critical value (VC between −2.672 and +2.672).Furthermore, the results of the Tukey's test confirmed this result (Table 5b), which suggests that the possible sources of cadmium pollution were similar in both states.We found that the cadmium concentrations in all the samples we analyzed were under the maximum values established by the Codex Alimentarius 0.05 mg kg −1 [5] and the 1.0 mg kg −1 FAO/WHO [4].In a study made in the Mexican state of Puebla, Castro found the presence o cadmium in urine, blood, and milk of Holstein ca le in the Santa Ana Xalmimilulco region in the municipality of Huejo ingo and determined its concentration in milk samples to be 0.54 mg kg −1 , which the authors a ributed to the ca le's ingestion of alfalfa containing heavy metals [21].In their study of the heavy metal content in samples of fresco cheeses sold in markets in the city of Puebla, Benítez found heavy metals, including an average Cd concentration of 0.13 mg kg −1 [16].In Egypt, Meshref reported an average cadmium concentration of 0.09 mg kg −1 in Kareish in the Beni-Suef region [22].Elbarbary and Hamouda found 0.24 mg kg −1 of cadmium in feta cheese samples [23].In Turkey Çetinkaya reported an average cadmium concentration of 0.028 mg kg −1 in cami bogaz cheese samples from the Trabzon region [24].The concentrations of cadmium in cheeses reported in the literature that we found were higher than those we observed in our study

Nickel
The United States Department of Health and Human Services (DHHS) determined that metallic nickel and some nickel compounds are carcinogens, and the ingestion of high In a study made in the Mexican state of Puebla, Castro found the presence of cadmium in urine, blood, and milk of Holstein cattle in the Santa Ana Xalmimilulco region in the municipality of Huejotzingo and determined its concentration in milk samples to be 0.54 mg kg −1 , which the authors attributed to the cattle's ingestion of alfalfa containing heavy metals [21].In their study of the heavy metal content in samples of fresco cheeses sold in markets in the city of Puebla, Benítez found heavy metals, including an average Cd concentration of 0.13 mg kg −1 [16].In Egypt, Meshref reported an average cadmium concentration of 0.09 mg kg −1 in Kareish in the Beni-Suef region [22].Elbarbary and Hamouda found 0.24 mg kg −1 of cadmium in feta cheese samples [23].In Turkey, Çetinkaya reported an average cadmium concentration of 0.028 mg kg −1 in cami bogazi cheese samples from the Trabzon region [24].The concentrations of cadmium in cheeses reported in the literature that we found were higher than those we observed in our study.

Nickel
The United States Department of Health and Human Services (DHHS) determined that metallic nickel and some nickel compounds are carcinogens, and the ingestion of high Ni concentrations can be hazardous to human health [25,26].The FAO and European Commission (EC) have not established the permitted maximum intake limit of nickel for humans [4,15]; however, the Integrated Risk Information System of the United States Environmental Protection Agency (US-EPA-IRIS) established a 300 µg/kg limit [14,27], and the WHO established a tolerable daily intake of 11 µg kg −1 of body weight for children and 12 µg kg −1 of body weight for adults [28].
The results of our analysis of artisanal cheese samples from Tabasco showed an average nickel concentration of 0.0039 ± 0.0046 mg kg −1 , a maximum value of 0.0132 mg kg −1 , and a minimum of 0.0002 mg kg −1 .By cheese variety, the average, maximum and minimum values in mg kg −1 were as follows: 0.0059, 0.0002, and 0.0132 for crema; 0.0029, 0.0005, and 0.0049 for oaxaca; 0.0031, 0.0021, and 0.0042 for fresco; and 0.0029, 0.0004, and 0.0056 for panela.For the cheese samples from Chiapas, we found an average Ni concentration of 0.0031 ± 0.0039 mg kg −1 with a minimum of 0.0003 mg kg −1 and a maximum of 0.0091 mg kg −1 .For the same state by variety, we found average, minimum, and maximum Ni concentrations in mg kg −1 of 0.0050, 0.0015, and 0.0091 for crema; 0.0052, 0.0027, and 0.0079 for cotija; 0.0031, 0.0016, and 0.0052 for oaxaca; and 0.0017, 0.0006, and 0.0031 for panela.The distribution of nickel concentrations in cheese samples we analyzed from Tabasco and Chiapas (Figure 3) shows that matured cheeses with lower humidity content like crema and cotija have higher concentrations of nickel than the more humidity varieties like oaxaca, panela, and fresco (Tables 3 and 4).The statistical analysis of the nickel concentrations in the samples yielded a value of TC = +0.698,and Tukey's test indicated no significant differences (Table 5c), again showing similar contamination sources in both geographic areas.
Molecules 2023, 28, x FOR PEER REVIEW varieties like oaxaca, panela, and fresco (Tables 3 and 4).The statistical analys nickel concentrations in the samples yielded a value of TC = +0.698,and Tuk indicated no significant differences (Table 5c), again showing similar contam sources in both geographic areas.
Our results compare well with those of the analyses of Ni concentrations m Castro in samples of oaxaca (0.003 mg kg −1 ) and ranchero (0.01 mg kg −1 ) artisanal from the Mexican state of Puebla [17].However, the values of nickel concentr artisanal cheeses reported by Nöel for French samples (0.409 mg kg −1 ) [29], and by for milk (7.70 mg kg −1 ) and cheese (8.33 mg kg −1 ) samples from the Nigerian O Oyo states are much higher than those reported for Mexico [30]; these authors a the contamination with nickel of milk and cheese to oil industry activities and th agrochemical inputs.

Sources of Lead, Cadmium, and Nickel in Artisanal Cheeses from Tabasco and Chiap
The presence of lead, cadmium, and nickel in the artisanal cheeses made in and Chiapas might come from the natural and anthropic sources described below a.Oil industry.Large oil production facilities in Tabasco include Samaria, Golpe, Delta del Grijalva, Bellota, Ogarrio, Cinco Presidentes, and Cu Macuspana, and in Chiapas, Cactus.Southeastern México has 99 oil fields w active wells and 2360 km of oil ducts connecting wells and oil fields with p Our results compare well with those of the analyses of Ni concentrations made by Castro in samples of oaxaca (0.003 mg kg −1 ) and ranchero (0.01 mg kg −1 ) artisanal cheeses from the Mexican state of Puebla [17].However, the values of nickel concentration in artisanal cheeses reported by Nöel for French samples (0.409 mg kg −1 ) [29], and by Olujimi for milk (7.70 mg kg −1 ) and cheese (8.33 mg kg −1 ) samples from the Nigerian Ogun and Oyo states are much higher than those reported for Mexico [30]; these authors attributed the contamination with nickel of milk and cheese to oil industry activities and the use of agrochemical inputs.

Sources of Lead, Cadmium, and Nickel in Artisanal Cheeses from Tabasco and Chiapas
The presence of lead, cadmium, and nickel in the artisanal cheeses made in Tabasco and Chiapas might come from the natural and anthropic sources described below: a.
Oil industry.Large oil production facilities in Tabasco include Samaria, Jujo, El Golpe, Delta del Grijalva, Bellota, Ogarrio, Cinco Presidentes, and Cuenca de Macuspana, and in Chiapas, Cactus.Southeastern México has 99 oil fields with 935 active wells and 2360 km of oil ducts connecting wells and oil fields with pumping stations and oil processing centers [31].Manríquez determined heavy metals were present in the Mexican crude oils Maya, Istmo, and Olmeca, reporting concentrations between 8 and 277 mg L −1 of vanadium and from 2.5 to 52.0 mg L −1 of lead [32].Siebe reported the presence of lead (0.01-0.1 mg L −1 ), chromium (2.0 mg L −1 ), zinc (0.1-2.5 to 0.2 mg L −1 ), cadmium (0.02-0.04 mg L −1 ), and copper (0.2-0.4 mg L −1 ) [33].In the Cinco Presidentes oil field well in Tabasco, Fridler analyzed heavy metal concentrations average in groundwater reported for Nickel 9 ug L −1 , zinc 38 ug L −1 , copper 16 ug L −1 , cadmium 0.5 ug L −1 , and lead 8 ug L −1 ; in sediments, nickel (49-92 mg kg −1 ), chromium (73-138 mg kg −1 ), and lead (8-95 mg kg −1 ) [31].Villanueva and Botello determined the content of heavy metals in the sediments of the Laguna el Yucateco in the state of Tabasco, and in the muscles of the inhabiting aquatic organisms [34], reporting between 14.27 and 329.24 µg g −1 of lead, 0.76-5.52µg g −1 of cadmium, and 44.61 ± 5.38 µg g −1 of nickel in sediments, and in the muscle of several edible species, 0.33-4.30mg kg −1 of cadmium, 0.19-15.68mg kg −1 of lead, and 0.1-8.75mg kg −1 , of nickel-the latter concentrations being above the maximum limits established in the FAO/WHO and Codex Alimentarius [4,5].Therefore, the presence of heavy metals is an ecological risk factor because they can move from soils to groundwater, and from there, to aquatic ecosystems where edible species bioaccumulate them; b.
Urban areas and roads.The large towns and terrestrial communication networks in the regions of Tabasco and Chiapas where artisanal cheese samples were collected are possible sources of heavy metal pollution.Internal combustion engines used in vehicles and the oil industry might emit gasses containing high amounts of heavy metals [46].The solid microparticles generated during gasoline and diesel combustion travel long distances dispersing heavy metals and are deposited by gravity on croplands and grazelands [47]; besides, Akpoveta and Osakwe reported contents of 0.24 ppm of lead, 1.68 ppm of cadmium, 1.74 ppm of copper, and 1.43 ppm of zinc in gasolines, and of 1.01 ppm of lead, 1.50 ppm of cadmium, 1.77 ppm of copper, and 2.87 ppm of nickel in diesel [48]; d.
Processing tools and inputs.The use of metallic utensils and the addition of salt during the manufacturing of artisanal cheeses are possible sources of heavy metal pollution [49]; e.
Volcanic emissions.The eruption of the Chichonal volcano in 1982 emitted large quantities of ashes over Chiapas, Tabasco, Campeche, and southern Veracruz [50].
In this regard, Rincón analyzed the sediments from the Chichonal volcano finding 3.26-7.06mg kg −1 of cadmium, 3.2-4.5 mg kg −1 of lead, 3.2-4.3mg kg −1 of nickel, 0.83-2.76mg kg −1 of copper, 3.5-17 mg kg −1 of zinc, and 52-126 mg kg −1 of iron, thus showing that volcanic activity contributed to the distribution of heavy metals on the regions in Tabasco and Chiapas where we collected the samples of artisanal cheeses [51].

Copper
Copper is an essential microelement active in vital functions of the human body; however, the ingestion of high concentrations of copper might cause health issues [49].The results of our analyses of cheese samples made in Tabasco and Chiapas showed copper concentrations below the maximum limits established by the FAO/WHO [4].The average, minimum, and maximum copper concentrations in mg kg −1 that we observed in cheese samples from Tabasco were 0.0199 ± 0.021, 0.0059, and 0.0437.By cheese variety, the average, minimum, and maximum copper concentrations in mg kg −1 were 0.0285, 0.0129, and 0.0402 for crema, 0.0233, 0.0059, and 0.0884 for oaxaca, 0.0093, 0.0078, and 0.0108 for fresco, and 0.0148, 0.0059, and 0.0346 for panela.In the cheese samples from Chiapas, the average, minimum, and maximum copper concentrations in mg kg −1 that we recorded were 0.0202 ± 0.022, 0.0007, and 0.0354.By cheese variety, the average, minimum, and maximum copper concentrations in mg kg −1 were as follows: 0.0271, 0.0167, and 0.0354 for crema; 0.0141, 0.0067, and 0.0223 for oaxaca; 0.0267, 0.0113, and 0.0337 for cotija; and 0.0165, 0.0059, and 0.0285 for panela.No significant differences in copper concentrations were found between the artisanal cheeses from Tabasco and Chiapas.The paired data test was found to be TC = +0.3578, a value within the CV range (−2.672 to +2.672) that was confirmed by the Tukey's test (Table 5d).
The distribution of copper concentrations we observed in the cheese samples from Tabasco and Chiapas shown in Figure 4 shows a high variation in such concentrations.By cheese variety, our results determined the samples of crema contained copper at concentrations ranging between 0.0402 and 0.0129 mg kg −1 , of oaxaca between 0.0437 and 0.0007 mg kg −1 , of cotija between 0.0337 and 0.0113 mg kg −1 , of frescos between 0.0197 and 0.0078 mg kg −1 , and of panela between 0.0346 and 0.0039 mg kg −1 .We can attribute this variability and the low concentration of copper in the cheese samples from Tabasco and Chiapas to the states' soil type variation.In both states the soil types include Gleysol, Fluvisol, Histosol, Leptosol-Vertisol, Acrisol, Vertisol, Luvisol, and Ferralsol, whose differences in the values of pH, electric conductivity, cation exchange capacity, and organic matter content are factors determining the amount of extractable concentration of copper in each soil type, in studies carried out by Palma and the National Institute of Statistic and Geography [52,53].Studies carried out by De la Cruz report extractable copper concentrations of 0.99-2.45mg kg −1 in Vertisols and 0.98-2.45mg kg −1 in Fluvisoles, and Salgado reported extractable copper concentrations in subunits of Acrisols of 0.60 ± 0.30 mg kg −1 (subunit ACdyhfr) 0.80 ± 0.50 mg kg −1 (ACfrpl), 0.40 ± 0.30 mg kg −1 (ACfrum), 1.10 ± 0.20 mg kg −1 (ACglpl), 0.80 ± 0.50 mg kg −1 (AChupl), 0.60 ± 0.50 mg kg −1 (AChuum), 0.80 ± 0.20 mg kg −1 (ACumgl) and 0.70 ± 0.50 mg kg −1 (ACumpl) [35,36].The value of copper concentration that we found in the analyzed cheese samples was below the 0.40 mg L −1 maximum permitted limit in milk and dairy products by the European Commission (EC) [54].The copper contained in the artisanal cheese samples from Tabasco and Chiapas may contribute to the recommended daily ingestion of 3 mg for adults FAO/WHO [4].
copper concentration that we found in the analyzed cheese samples was below mg L −1 maximum permi ed limit in milk and dairy products by the E Commission (EC) [54].The copper contained in the artisanal cheese samples from and Chiapas may contribute to the recommended daily ingestion of 3 mg fo FAO/WHO [4].Oaxaca and ranchero cheeses made in Santa Ana Xalmimilulco, Puebla reported a 0.02 mg kg −1 copper content, and for cheese samples from markets in of Puebla [55].In Egypt, Meshref reported copper concentrations of between 0 0.53 mg kg −1 in samples of Kareish cheese made in the Beni-Suef region [22].In Elbarbary and Hamouda reported copper concentrations of 3.25 ± 1.06 mg kg − cheese [56], and Reinholds found an average copper concentration of 0.29 mg cheeses made in the Kvemo Kartli region of Georgia [14].Al Sidawi determ average copper concentration of 1.261 ± 0.739 mg kg −1 in imeruli cheese and 2.46 Oaxaca and ranchero cheeses made in Santa Ana Xalmimilulco, Puebla, Castro reported a 0.02 mg kg −1 copper content, and for cheese samples from markets in the city of Puebla [55].In Egypt, Meshref reported copper concentrations of between 0.002 and 0.53 mg kg −1 in samples of Kareish cheese made in the Beni-Suef region [22].In Europe, Elbarbary and Hamouda reported copper concentrations of 3.25 ± 1.06 mg kg −1 for feta cheese [56], and Reinholds found an average copper concentration of 0.29 mg kg −1 in cheeses made in the Kvemo Kartli region of Georgia [14].Al Sidawi determined an average copper concentration of 1.261 ± 0.739 mg kg −1 in imeruli cheese and 2.463 ± 2.314 mg kg −1 in sulguni cheese [49].Previous reports of copper content in cheese samples are similar to or higher than those we determined in our analysis of cheese varieties made in the states of Tabasco and Chiapas.

Zinc
Zinc is another essential element that when ingested at high concentrations might lead to neurological, hematological, immunological, renal, hepatic, cardiovascular, and genotoxic conditions [23].The average, minimum, and maximum concentrations of zinc in mg kg −1 that we determined in the cheese samples from the states of Tabasco were 0.161 ± 0.18, 0.0204, and 0.301 mg kg −1 .By cheese variety, the average, minimum, and maximum zinc concentrations in mg kg −1 we observed were as follows: 0.0874, 0.0204, and 0.2418 for crema; 0.1579, 0.0317, and 0.2645 for oaxaca; 0.2584, 0.2435, and 0.2785 for fresco; and 0.2437, 0.1436, and 0.2927 for panela.For the samples from Chiapas, the average, minimum, and maximum zinc concentrations in mg kg −1 that we observed were 0.194 ± 0.21, 0.0421, and 0.437 mg kg −1 .By cheese variety, the average, minimum, and maximum zinc concentrations in mg kg −1 we recorded were as follows: 0.1008, 0.0421, and 0.3699 for crema; 0.2437, 0.0677, and 0.4369 for oaxaca; 0.2216, 0.1664, and 0.3708 for cotija; and 0.2181, 0.1325, and 0.3308 for panela.It was found that there are no significant differences in the Zn concentrations in the Tabasco and Chiapas cheese samples.In Table 5e, it is observed that Tukey's test presents only one grouping (A) and the Student's t-test with TC (−1.376) is between the critical values VC (−2.672; +2.672).
As seen in Figure 5, we found an ample variation in the zinc concentration of the cheese samples that we analyzed, for example, these values (expressed in mg kg −1 ) were 0.0204-0.369for crema, 0.032-0.437for oaxaca, 0.188-0.278for cotija, 0.167-0.374for fresco, and 0.132-0.334for panela.Such variability in zinc content might have been due to the diversity of soil types in Tabasco and Chiapas, and to the physicochemical variables of these types [52], which agrees both with the similarities in soil types between both states and with the recognition that the content of zinc in the cheese samples was due to the soil types present in the geographic regions where they were collected, as previously reported in other studies.For example, De la Cruz determined Zn concentrations of 0.68-0.85mg kg −1 in Fluvisols, and of 0.66-0.87mg kg −1 in Vertisols, and in different Acrisol units [35], and Salgado reported Zn concentrations in mg kg −1 of 0.30 ± 0.20 (ACdyhfr), 0.50 ± 0.30 (ACfrpl), 0.40 ± 0.10 (ACfrum), 0.50 ± 0.20 (ACglpl), 0.60 ± 0.50 (AChupl), 0.30 ± 0.20 (AChuum), 0.30 ± 0.20 (ACumgl), and 0.51 ± 0.30 (ACumpl) [36].
As seen in Figure 5, we found an ample variation in the zinc conce cheese samples that we analyzed, for example, these values (expressed in 0.0204-0.369for crema, 0.032-0.437for oaxaca, 0.188-0.278for cotija, 0 fresco, and 0.132-0.334for panela.Such variability in zinc content might ha the diversity of soil types in Tabasco and Chiapas, and to the physicochem of these types [52], which agrees both with the similarities in soil types betw and with the recognition that the content of zinc in the cheese samples was types present in the geographic regions where they were collected, as previ in other studies.For example, De la Cruz determined Zn concentrations o kg −1 in Fluvisols, and of 0.66-0.87mg kg −1 in Vertisols, and in different Acr and Salgado reported Zn concentrations in mg kg −1 of 0.30 ± 0.20 (ACdyh (ACfrpl), 0.40 ± 0.10 (ACfrum), 0.50 ± 0.20 (ACglpl), 0.60 ± 0.50 (AChup (AChuum), 0.30 ± 0.20 (ACumgl), and 0.51 ± 0.30 (ACumpl) [36].The recommended daily intake of zinc for adult women and men is 8 mg and 11 mg kg −1 , respectively.The range of zinc concentration we found in cheese samples from Tabasco and Chiapas suggests that their consumption might contribute to satisfy the minimum requirement of zinc for humans.Our results agree with those of Castro for cheeses in the Mexican state of Puebla, who reported zinc concentrations of 0.18 ± 0.09 mg kg −1 for oaxaca and 0.74 ± 0.1 mg kg −1 for ranchero cheeses [17].In work carried out by Ghafari and Sobhanardakani reported 0.198 mg kg -1 of zinc in the cheeses produced from Hamedan, Iran [57].Al Sidawi reported zinc contents of 75.86 ± 52.528 mg kg −1 in imeruli and 124.8 ± 97.775 mg kg −1 in sulguni cheeses manufactured in the province Kartli, Georgia [49].For blue cheese sold in the European Union, Reinholds reported zinc concentrations between 10.4 and 39.5 mg kg −1 [14].For different cheeses made in Turkey, Mendil and Orak reported 12.0 mg kg −1 and 15.57mg kg −1 of zinc for white cheese [58,59]; Çetinkaya and collaborators found an average zinc concentration of 27.52 ± 1.85 mg kg −1 in cheese samples in the Cami Bogazi region and for white cheeses [24].The authors attributed the presence of zinc in cheeses to the use of contaminated machinery and containers during cheese making, and to the transport of zinc through the food web due to the environmental pollution of soil, water, and fodder.

Iron
Iron is essential for oxygen transportation and storage in the human body, but ingestion of large quantities of iron may lead to blood, heart, kidney, and endocrine system conditions, cellular damage, and mutations [25,60].Figure 6 shows the distribution of Fe in the cheese samples from Tabasco and Chiapas.
samples from Chiapas, the average, minimum, and maximum Fe concen kg −1 were 65.76 ± 6.61, 55.82, and 97.41, and by cheeses variety were as follow and 72.33 for crema; 68.09, 57.96, and 97.41 for oaxaca; 65.91, 62.19, and 7 and 64.72, 59.86, and 72.94 for panela.The results of the Fe concentrations o Chiapas cheeses showed significant differences since the paired Student' value of TC = −3.696,which is outside the range of CV (−2.672 to +2.672).T showed two groupings (A) and (B), thus confirming the previous results (T  The average, minimum, and maximum Fe concentrations in mg kg −1 that we found in cheese samples from Tabasco were 61.84 ± 4.23, 55.97, and 72.76, whereas for the samples from Chiapas were 62.58, 55.97, and 72.08 for crema, 62.15, 56.92, and 72.76 for oaxaca, 61.61, 58.38, and 68.29 for fresco, and 60.25, 56.07, and 67.77 for panela.For cheese samples from Chiapas, the average, minimum, and maximum Fe concentrations in mg kg −1 were 65.76 ± 6.61, 55.82, and 97.41, and by cheeses variety were as follows: 64.37, 55.82, and 72.33 for crema; 68.09, 57.96, and 97.41 for oaxaca; 65.91, 62.19, and 72.09 for cotija; and 64.72, 59.86, and 72.94 for panela.The results of the Fe concentrations of Tabasco and Chiapas cheeses showed significant differences since the paired Student's t-test gave a value of TC = −3.696,which is outside the range of CV (−2.672 to +2.672).The Tukey's test showed two groupings (A) and (B), thus confirming the previous results (Table 5f).
We attributed the presence of Fe in the analyzed cheese samples from Tabasco and Chiapas to soil types in the region, where Ferralsols rich in iron sesquioxides are dominant.For example, the Acrisols-soils which are easily identified by their yellowish to reddish dark color, strong acidity, and a B horizon accumulating alluvial clay-are rich in iron and aluminum sesquioxides.Salgado and Palma analyzed Fe concentrations in subunits of Acrisols finding 67.00 ± 52.00 (ACdyhfr), 76.00 ± 21.00 (Acfrpl), 51.00 ± 16.00 (Acfrum), 115.00 ± 58.00 (Acglpl), 85.00 ± 29.00 (Achupl), 94.00 ± 29.00 (Achuum), 72.00 ± 15.00 (Acumgl), and 64.00 ± 49.00 (Acumpl) mg kg −1 [36,52].Additionally, De la Cruz found Fe concentrations of 67.6 in Vertisols and 112.2 mg kg −1 in Fluvisols [35].The abovementioned soil types are rich in organic matter with high concentrations of humic and fulvic acids, which give soil acidic pH values between 4.6 and 5.7, chemical conditions that facilitate the biosorption of iron by grass rhizomes and its accumulation in stalks.
According to our results, the concentrations of iron that we observed in cheese samples from Tabasco and Chiapas are within the parameters for human consumption and may cover a large part of the recommended daily intake of 8 mg for adult males, 18 mg for adult females, 11 mg for young males, 15 mg for young females, and 27 mg for pregnant women.Our results were similar to those of Al Sidawi who analyzed cheese samples from the Kvemo Kartli (Georgia), finding Fe concentrations of 69.09 ± 64.918 mg kg −1 in imeruli and 101.1 ± 91.166 mg kg −1 in sulguni cheeses [49].Jalili also found Fe concentrations between 67.7 ± 4.5 and 71.3 ± 4.9 mg kg −1 in iron-fortified feta cheeses [61].In contrast, the values for Fe concentration we found exceed the values reported in previous studies made in Turkey by Kirdar in akcaka tik cheeses between 7.49 and 29.05 mg kg −1 [62].While Centinkaya in cami bogazi cheeses from Trabzon (Turkey) reported 0.371 ± 0.177 mg kg −1 Fe.Meshref found Fe concentrations between 1.763 and 17.739 ppm in kareish cheeses from the Beni-Suef region [22], and in Europe, Reinholds reported Fe concentrations between 1.57 and 12.4 mg kg −1 in blue cheeses [14].

Risk Analysis
The values of HQ we estimated for the concentrations of the heavy metals lead, cadmium, nickel, copper, and zinc that we determined in cheese samples from Tabasco are shown in Figure 7a,b for the samples from Chiapas.Zinc, copper, and cadmium were the metals that contributed the most to the HQ values of samples from Tabasco in the order young girls > adult woman > young boys > adult man (Figure 7a), while for cheese samples from Chiapas, zinc and copper were the metals having the highest contribution to the HQ values in the same order as for Tabasco.However, all the THQ values that we estimated were lower than one (Figure 7c,d), meaning that there is no risk of developing medical conditions from intake of heavy metals contained in the cheese varieties we analyzed.In previous studies, Castro reported values of THQ below one for arsenic in milk samples from the upper Balsas River region in the state of Puebla [17,21]; Reinholds found THQ between 0.05 and 0.14 for youngsters, and between 0.03 and 0.09 for adults for the intake of blue cheese consumed by the populations from Denmark, France, Italy, Spain, and the UK [14], and Zafarzadeh reported THQ values of 1.11, 1.33, and 5.42 for adults In previous studies, Castro reported values of THQ below one for arsenic in milk samples from the upper Balsas River region in the state of Puebla [17,21]; Reinholds found THQ between 0.05 and 0.14 for youngsters, and between 0.03 and 0.09 for adults for the intake of blue cheese consumed by the populations from Denmark, France, Italy, Spain, and the UK [14], and Zafarzadeh reported THQ values of 1.11, 1.33, and 5.42 for adults and youngsters due to the intake of cadmium-polluted butter from the Gorgan region in Iran [63], which poses a high health risk for the inhabitants of that geographic area.
Table 6 shows the results of the Tukey tests for CR and CRT considering the variables heavy metals (Pb, Cd, and Ni) and age group (youngsters and adults).The results revealed significant differences in the effect on CR of each metal in the order Cd > Pb > Ni, which means that the ingestion by youngsters and adults of even low concentrations of Cd could represent a high CR relative to the ingestion of Pb or Ni; however, the CR will depend on the consumed cheese variety.The significant differences in CR between age and weight groups (Table 6a,b) could have been due to adults having slower metabolic rates than youngsters [21].Table 6c shows no significant differences in CRT between adults and youngsters, possibly due to the low contents of Pb, Cd, and Ni of the artisanal cheeses made in Tabasco and Chiapas.The CR and CRT values we estimated for adults and youngsters are shown in Figure 8a,b.In our study, lead and cadmium were the heavy metals that contributed the most to cancer risk in both states, with a higher value for women than for men.The values of CRT were in the decreasing order of young girls > adult women > young boys > adult men.In general, as seen in Figure 8c,d, all the CRT values we estimated were within the range of 10 −4 -10 −6 established by the US-EPA [8].
Our CRT estimates were lower than the values of 0.0018-0.014reported by Castro for children and youngsters in the state of Puebla from the intake of milk containing heavy metals due to contamination of fodder irrigated with the polluted water from the Atoyac, Xochiac, and Xopanac rivers [21].
The CR and CRT values we estimated for adults and youngsters are shown in Figure 8a,b.In our study, lead and cadmium were the heavy metals that contributed the most to cancer risk in both states, with a higher value for women than for men.The values of CRT were in the decreasing order of young girls > adult women > young boys > adult men.In general, as seen in Figure 8c,d, all the CRT values we estimated were within the range of 10 −4 -10 −6 established by the US-EPA [8].

Artisanal Cheese Samples
For the analysis of heavy metals, we acquired 44 samples of artisanal cheeses from seven varieties (doble crema, crema, oaxaca, fresco, panela, mozzarella, and provolone) in four municipalities (Huimanguillo, Balancán, Tenosique, and Centro) of the state of Tabasco, and 44 samples from five varieties (doble crema, oaxaca, panela, cotija, and fresco) in four municipalities (Solosuchiapa, Juárez, Catazajá, and Reforma) of the state of Chiapas.The average weight of samples was 1 kg (Figures 9 and 10).

Sample Preparation
All the glassware used in the experimental procedures was washed in a 5% hydrochloric acid solution for 24 h and afterward rinsed with deionized water to remove possible contaminants that could interfere with the results of the analyses.Cheese samples were digested following the method of [19].Samples were first homogenized by being finely chopped (approximately 2 mm).A 1 g aliquot of each sample homogenate was weighted in an Ohaus analytical balance and placed in a porcelain crucible and calcined in a Thermo Scientific-BF51794C-1 furnace at 450-500 • C for 16 h and, after, slowly cooled to room temperature.The resulting ashes were treated with 1 mL of concentrated nitric acid and calcined once more at 450-500 • C for 6 h to destroy the organic matrix.After the second calcination, samples were dissolved in 2 mL of concentrated nitric acid and 2 mL of hydrogen peroxide, after which the resulting solutions were filtered with no.41 Whatman paper and the filtrates were diluted to 100 mL in a 100 mL volumetric flask.Each final solution was placed in a Teflon jar and stored at 10 • C until they were analyzed in a Thermo Scientific ICE 3000 Series atomic absorption spectrometer as described below.

Artisanal Cheese Samples
For the analysis of heavy metals, we acquired 44 samples of artisanal cheeses fr seven varieties (doble crema, crema, oaxaca, fresco, panela, mozzarella, and provolone four municipalities (Huimanguillo, Balancán, Tenosique, and Centro) of the state Tabasco, and 44 samples from five varieties (doble crema, oaxaca, panela, cotija, a fresco) in four municipalities (Solosuchiapa, Juárez, Catazajá, and Reforma) of the state Chiapas.The average weight of samples was 1 kg (Figures 9 and 10).

Sample Preparation
All the glassware used in the experimental procedures was washed in hydrochloric acid solution for 24 h and afterward rinsed with deionized water to r possible contaminants that could interfere with the results of the analyses.Cheese sa were digested following the method of [19].Samples were first homogenized by finely chopped (approximately 2 mm).A 1 g aliquot of each sample homogena weighted in an Ohaus analytical balance and placed in a porcelain crucible and ca in a Thermo Scientific-BF51794C-1 furnace at 450-500 °C for 16 h and, after, slowly to room temperature.The resulting ashes were treated with 1 mL of concentrated acid and calcined once more at 450-500 °C for 6 h to destroy the organic matrix.Af second calcination, samples were dissolved in 2 mL of concentrated nitric acid and of hydrogen peroxide, after which the resulting solutions were filtered with Whatman paper and the filtrates were diluted to 100 mL in a 100 mL volumetric Each final solution was placed in a Teflon jar and stored at 10 °C until they were an in a Thermo Scientific ICE 3000 Series atomic absorption spectrometer as described

Analytical Parameters for Atomic Absorption Spectrometry
To determine the concentration in samples of heavy metals (measured in m calibration curves were run at different wavelengths for each heavy metal analy follows: 217.0 nm for Pb, 228.8 nm for Cd, 232.0 nm for Ni, 324.8 nm for Cu, 213.9 Zn, and 248.3 nm for Fe.A mixture of air-acetylene combustion gas was used a wavelength was adjusted for each heavy metal.

Daily Intake of Metals
The health risks of human intake of artisanal cheeses made in Tabasco and C were estimated based on the results of heavy metal determinations in the an samples by assuming a per capita yearly intake of 6 kg-taking into account the intake value of between 4 and 8 kg established for developing countries by the OEC

Analytical Parameters for Atomic Absorption Spectrometry
To determine the concentration in samples of heavy metals (measured in mg L −1 ), calibration curves were run at different wavelengths for each heavy metal analyzed as follows: 217.0 nm for Pb, 228.8 nm for Cd, 232.0 nm for Ni, 324.8 nm for Cu, 213.9 nm for Zn, and 248.3 nm for Fe.A mixture of air-acetylene combustion gas was used and the wavelength was adjusted for each heavy metal.

Daily Intake of Metals
The health risks of human intake of artisanal cheeses made in Tabasco and Chiapas were estimated based on the results of heavy metal determinations in the analyzed samples by assuming a per capita yearly intake of 6 kg-taking into account the yearly intake value of between 4 and 8 kg established for developing countries by the OECD/FAO [12]-using Equation (1) as suggested by Castro [21].
Chronic Daily Intake CDI = (C metal)(D intake) B average weight where C metal is the metal concentration determined in cheese samples in mg kg −1 , D intake is the per capita consumption of cheese in kg yr −1 ), and B average weight is the body weight of individuals in kg.

Hazard Quotient
The hazard quotient was estimated by Equation (2).
where HQ is the hazard quotient, CDI is chronic daily intake, RfD is oral reference dose value of exposure to the chronic daily intake in mg kg −1 .Table 7 contains the values used for the calculation of HQ [21].

Total Hazard Quotient
The total hazard quotient is calculated by adding the HQ values for each metal (Equation ( 3)).THQ = ∑ (HQ Cd + HQ Pb + HQ Ni + HQ Cu + HQ Zn + HQ Fe ) (3 where THQ is the total hazard quotient, HQ Cd is the hazard quotient for cadmium, HQ Pb is the hazard quotient for lead, HQ Ni is the hazard quotient for nickel, HQ Cu is the hazard quotient for copper, HQ Zn is the hazard quotient for zinc, and HQ Fe is the hazard quotient for iron.According to the recommendations of the US-EPA [8], values of THQ smaller than one imply no risk, while values of THQ greater than one suggest a high risk for human health [21].

Cancer Risk and Total Cancer Risk
The total cancer risk was calculated by adding the cancer risk values for the carcinogenic metals cadmium, nickel, and lead (Equations ( 4) and ( 5)) recommended by Reinholds [14].
CRT = ∑ CR Pb+Cd+Ni (5) where CR is the cancer risk, CRT is the total cancer risk, CR Cd+Ni+Pb is the cancer risk for cadmium, nickel, lead, CDI is chronic daily intake, and slope factors (Sf ).The average body weights of the Mexican population were established by the OECD/FAOto be 74.8 kg for adult males, 68.7 kg for adult females, 62.9 kg for young females, and 70.4 kg for young males [12].

Conclusions
The contents of Pb, Cd, Ni, Cu, and Zn showed no significant differences between cheese samples from the states of Tabasco and Chiapas.The concentrations of Pb, Cd, and Ni that we determined in all cheese samples were below the values established by the Codex Alimentarius and FAO/WHO.The heavy metals Zn, Cu, and Ni contributed the most to the HQ values.All THQ values were smaller than one, which implies no hazard from the intake of the cheese varieties we analyzed and were in the descending order young females > adult females > young males > adult males.The heavy metals Pb and Cd contributed the most to the CR and CRT values-which followed the descending order young girls > adult females > young boys > adult males-that were within the values established by the US-EPA, therefore indicating no cancer risk due to the intake of the analyzed cheeses.Considering the population growth and increased anthropic activities in Tabasco and Chiapas, it is essential to determine the presence in locally made cheeses of other carcinogenic heavy metals like arsenic, mercury, chromium, and vanadium for their incorporation into the estimations of the THQ and CRT values.We attributed the concentrations of Cu, Zn, and Fe in the analyzed cheese samples to the soil types present in Tabasco and Chiapas.

Figure 1 .
Figure 1.Distribution of lead in the different varieties of artisanal from Tabasco and Chiapas.

Figure 2 .
Figure 2. Distribution of cadmium in the different varieties of artisanal cheeses from Tabasco and Chiapas.

Figure 2 .
Figure 2. Distribution of cadmium in the different varieties of artisanal cheeses from Tabasco and Chiapas.

Figure 3 .
Figure 3. Distribution of nickel in the different varieties of artisanal cheeses from Tab Chiapas.

Figure 3 .
Figure 3. Distribution of nickel in the different varieties of artisanal cheeses from Tabasco and Chiapas.

Figure 4 .
Figure 4. Distribution of copper in the different varieties of artisanal cheeses from Tab Chiapas.

Figure 4 .
Figure 4. Distribution of copper in the different varieties of artisanal cheeses from Tabasco and Chiapas.

Figure 5 .
Figure 5. Distribution of zinc in the different varieties of artisanal cheeses from Tabas

Figure 5 .
Figure 5. Distribution of zinc in the different varieties of artisanal cheeses from Tabasco and Chiapas.

Figure 6 .
Figure 6.Distribution of iron in the different varieties of artisanal cheeses from Tabas

Figure 6 .
Figure 6.Distribution of iron in the different varieties of artisanal cheeses from Tabasco and Chiapas.

23 Figure 7 .
Figure 7. Values of the risk index and contribution of each metal for adults and youngsters (a,b); values of the total risk index for adults and youngsters (c,d).

Figure 7 .
Figure 7. Values of the risk index and contribution of each metal for adults and youngsters (a,b); values of the total risk index for adults and youngsters (c,d).

Figure 8 .
Figure 8. Values of the cancer risk index and contribution of each metal, for adults and youngsters (a,b); values of the CRT index for adults and youngsters (c,d).

Figure 8 .
Figure 8. Values of the cancer risk index and contribution of each metal, for adults and youngsters (a,b); values of the CRT index for adults and youngsters (c,d).

Figure 9 .
Figure 9. Municipalities of Tabasco and Chiapas, where samples of artisanal cheese were collect

Figure 9 .
Figure 9. Municipalities of Tabasco and Chiapas, where samples of artisanal cheese were collected.

Figure 9 .
Figure 9. Municipalities of Tabasco and Chiapas, where samples of artisanal cheese were co

Table 1 .
Metal concentrations in mg kg −1 in samples of different varieties of artisan cheeses made in the state of Tabasco, Mexico.

Table 2 .
Metal concentrations in mg kg −1 in samples of different varieties of artisan cheeses made in the state of Chiapas, Mexico.

Table 3 .
Moisture content in artisanal cheeses from the states of Tabasco and Chiapas.

Table 4 .
Tukey test to determine significant differences between moisture content in artisanal cheeses from the states of Tabasco and Chiapas.Means with a common letter are not significantly different (p > 0.05) α = 0.05.

Table 5 .
Tukey test to determine significant differences between heavy metal content in artisanal cheeses from the states of Tabasco and Chiapas.Means with a common letter are not significantly different (p > 0.05) α = 0.05.

Table 6 .
Tukey's tests for significant differences in CR and CRT by heavy metal and age group.A, B and C letters indicates the significant different between CR and CRT.

Table 7 .
Oral reference dose (RfD) values and slope factors (Sf ) of metal considered to be carcinogenic.