Minerals multi-element analysis and its relationship with geographical origin of artisanal Mexican goat cheeses

ejramirezrivera@itszongolica.edu.mx Abstract The objective of this study was to perform a profile of mineral elements and the relationship of those as markers of geographical origin of artisanal goat cheese from the State of Veracruz, Mexico. On fresh and mature goat cheeses the quantified elements wereI: heavy (aluminum, arsenic, cadmium, lead), major (calcium, potassium, sodium, magnesium), essential (cobalt, chromium, copper, iron, manganese, selenium, zinc), and others (nickel and strontium). The database using variance analyses and discriminant function analysis was analyzed. Differences (P < 0.05) were found in the contents of elements between fresh and mature cheeses. The content of major, essential and other elements were similar to the results from other investigations. The concentration of arsenic and lead surpassed the limits accepted by international and Mexican standards. Nine elements helped to determine the geographic origin of the goat spectroscopy. Practical Application: Know minerals of Mexican artisanal cheeses and their relationship with the


Introduction
The artisanal goat cheeses are considered an excellent source of proteins, lipids, vitamins and mineral elements (Mendil, 2006;Moreno-Rojas et al., 2010).Which are the major (i.e.Ca, K, Mg and Na), essential (i.e.Co, Cu, Cr, Fe, Mn and Zn) and elements as Ni are needed for the proper metabolic functioning of the organism (Sevgi-Kirdar et al., 2015).However, the mineral content in the cheese can be modified by factors such as feeding goats, heat treatment of milk, processing of cheeses (i.e.coagulation, pressing and salting) and the ripening effects caused by the different biochemical reactions (i.e.lipolysis and proteolysis) and modifications of the water content (Almanera et al., 2007).Nevertheless in the context of food security, the chemical safety of artisanal cheeses may be affected by the proximity between goat farms producing artisanal cheeses and urban centers or industries that can contribute to chemical contamination of cheeses by heavy metals such As, Cd and Pb (Kodrik et al., 2001;ElSayed-Elham et al., 2011;Dsorio et al., 2015;Dzbek & Akman, 2016).
The consumption of artisanal goat cheeses contaminated with heavy metals can generate harmful effects to human such as neurological problems, cardiovascular toxicity, instability of the central nervous system, interference of the synthesis of the hemo group and osteoporosis (Moreno-Rojas et al. 2010;Elbarbary & Hamouda, 2013;Obrahim & Mehanna, 2015;Sevgi-Kirdar et al., 2015).On this sense, it was estimated that 1.52 million metric tons of lead were used in the world for various industrial applications in the year 2004 (Tchounwou et al., 2012).On the case of cadmium is a natural constituent of the earth's crust and occurs in soils in concentration of 0.1 to 1.0 mg⋅kg -1 dry weight, however, the fertilizers applied to feed crops for animals are another source of supply of this metal (Barbara-Fisher, 2005).While that in the countries of Bangladesh, Ondia, Chile, Uruguay, Taiwan and Mexico the existence of high contents of arsenic in air in the range of 20 to 100 ng⋅m -3 has been demonstrated (Tchounwou et al., 2012).
The importance of artisanal cheese is increasing both at the international and national due to changes in consumers' lifestyles (Cagri-Mehmetoglu, 2018).The consumption of artisanal cheese per capita in European countries (i.e.Greece, France) is 20 kg per year.On Mexico, the consumption of this food is 2.1 kg per year (Cervantes-Escoto & Villegas-de-Gante, 2014).Currently, the typification of Mexican artisan cheeses has been carried out according to their physicochemical, microbiological, instrumental and sensory characteristics (Cervantes-Escoto & Villegas-de-Gante, 2014;Ramírez-Rivera et al., 2017).Nevertheless, there is no scientific evidence on the determination of mineral elements in Mexican goat cheeses produced in the central area and highlands of the State of Veracruz.The researches of Ramírez-Rivera et al. (2017, 2018) determined the typicality of the goat cheeses of the State of Veracruz in Mexico without considering the analysis of mineral elements present in these cheeses.Veracruz State is considered an important geographical area for the production of fresh and mature artisanal goat cheeses, specifically; places like Coatepec, Pacho Viejo, Perote and Tatatila are those that produce 90% of these cheeses (Onstituto Nacional de Estadística y Geografía, 2007).These cheeses are attached to the territory, they are emblematic of the gastronomic culture of the area, and they are an important source of income for the people involved in their production (Ramírez-Rivera et al., 2017).Therefore, the objective of this study was to perform a profile of mineral elements and the relationship of those as markers of geographical origin of artisanal goat cheese from the State of Veracruz, Mexico.

Source of the samples and processing of artisanal cheeses
A total of eight artisanal goat cheeses were analyzed (n 1 = four fresh cheeses and n 2 = four mature cheeses) used in Goat Production Units (GPU) of municipalities as Pacho Viejo, Coatepec, Perote, and Tatatila, in the State of Veracruz in Mexico.These GPU are located in the central mountainous region and the highlands of Veracruz State, where the dominant vegetation are mesophilic mountain forest as oyamel and pine, and xerophilous scrub (Ramírez-Rivera et al., 2018).The cheeses of this research were chosen due to the following reasons: 1) the GUPs were affiliated with the nonprofit Goat Species Product-System of Veracruz; 2) the cheeses are made by hand; 3) cheeses fulfilled with the microbiological limits indicated by the Dfficial Mexican Standards and have great acceptance by consumers (Ramírez-Rivera et al., 2017;Ramírez-Rivera et al., 2018).Also, the number of samples analyzed in this research was adequate and similar to those analyzed by Ledesma et al. (2007) and Obrahim & Mehanna (2015).
The cheeses are made with milk from Alpine and Saanen goats (3.5 years old on average), milked manually in the dry season (May month).The quality of goat milk obtained present the microbiological parameters proposed by the Dfficial Mexican Standards (Ramírez-Rivera et al., 2018).The characteristics of the goat farms, milk production and places of location are shown in Table 1.The manufacturing stages for the artisanal goat cheese were as follows: 1) the goats were milked by hand, and the milk was pasteurized at 63 °C for 30 min, and was then cooled to 37 °C; 2) A commercial curdling agent was added with a strength of coagulation = 1:10,000; equivalent to 110 OMCU mL -1 (Cuamex Company, Mexico) in a proportion of 30 mL⋅100 L -1 ; after 45 minutes, the curd was cut; 3) the curd was molded in polyvinylchloride (PVC) cylinders and later, a stainless steel press was used to apply a pressure of 2 kg of force⋅kg -1 of cheese for 7 h; 4) the cheeses were submerged in a brine solution (28% of NaCl) and laid at room temperature (18 ± 2 °C) for two days to obtain the fresh cheeses.The matured cheeses were obtained through inoculation by spraying the strain Penicillum candidum (Choozit™ PC-VB, commercial brand Danisco, Dupont Mexico) on the fresh cheeses, which were then stored in wooden cavas for seven weeks at a temperature of 18 ± 2 °C and a relative humidity of 80-85%.

Conditioning the samples and analytical process
The samples were digested using a calcination method due to the high content of ethereal extract.Two g of each cheese were weighed and then charred at 550 °C in a muffle furnace (Felisa,Feligneo,Mexico) for 5 h (method 935.42 Association of Dfficial Analytical Chemists, 2012).Later, the ashes were dissolved in 4 mL of HCl 3 mol⋅L -1 and diluted at 25 mL with HND 3 at 1%.The HCl in form supra pure and HND 3 (at 65%) were obtained (Merck, Darmstadt, Germany).The concentrations of the heavy metals as aluminum (Al), arsenic (As), cadmium (Cd), and lead (Pb), majority elements as calcium (Ca), potassium (K), sodium (Na), and magnesium (Mg), essential elements as cobalt (Co), chromium (Cr), copper (Cu), iron (Fe), manganese (Mn), selenium (Se), and zinc (Zn), and other elements as nickel (Ni) and strontium (Sr) were quantified using a microwave induced Plasm Atomic Emission Spectroscopy (MP-AES, 4200MP-AES, Agilent Technologies, New Castle, Delawere, USA) connected to a nitrogen generator (Peak Genius 3055, Agilent Technologies, New Castle, Delawere, USA).The operation conditions (nebulizer flow and wave length per element) of the MP-AES equipment are shown in Table 2.The multielement solutions were diluted in a range of concentrations of 0.1 -5 mg⋅L -1 for Al, As, Cd, Co, Cu, Cr, Mn, Ni, Pb, Sr, Se, and Zn and concentrations of 0.1 -10 mg⋅L -1 for Ca, Fe, K, Mg, and Na.The calibration curves (coefficient of correlation R 2 = 0.99 per element) using multielement standard solutions (Agilent Technologies, Delawere, U.S.A.) were carried out for Al, As, Cd, Co, Cu, Cr, Mn, Ni, Pb, Sr, Se, and Zn (50 mg⋅L -1 ) and Ca, Fe, K, Mg, and Na (500 mg⋅L -1 ) The determinations were carried out six times for reproducible and reliable results (Dsorio et al., 2015).

Statistical analysis
The data were collected in a matrix of dimensions (J*M) K, where J = 8 cheeses (2 types of cheeses and 4 places), M = 6 repetitions and K = 17 elements for a total 816 data.Descriptive analysis and central tendency (mean and standard deviation), as well as one-way (type of cheese) analysis of variance (ANDVA) and a Tukey post hoc test to determine significant difference in the minerals concentration between both types of cheeses were applied (Peláez-Puerto et al.,2004;Moreno-Rojas et al.,2012).
Later, a Discriminating Analysis (DA) method stepwise procedure was used, this analysis was performed three times one for each qualitative variable (type of cheese, fresh cheese and matured cheese) to determine the geographic origin of the cheeses based on the quantified elements (Moreno-Rojas et al., 2010;Dsorio et al., 2015).The DA assesses new synthetic variables called "discriminant functions", are linear combinations of the selected principal components, and allow a better separation of the centers of gravity of the considered groups (Karoui et al., 2007).
On addition, the following statistical indicators derived stepwise, were interpreted: 1) Wilk's Lambda test (λ) and probability values (P) to determine the significance of discriminant functions.Wilks' lambda is computed as the ratio of the determinant of within-group variance/covariance matrix to the determinant of the total variance/covariance matrix, where values near to 0 indicate a high discriminatory power (Shintu & Caldarelli, 2006); 2) the distances of Mahalanobis and confidence ellipses with a confidence level of 95% were used to determine the separation between the cheeses in the factorial plane of the DA.On this sense, it was considered that the greater distance is the difference between the cheeses; 3) the percentage (%) of classification was applied to check the discriminatory capacity of the model generated based on the content of elements (Almanera et al., 2007;Moreno-Rojas et al., 2010 and2012).
The descriptive statistics, central tendency, ANDVA and Tukey test were carried out using the software STATGRAPHOC PLUS ® version 5.2 (Statistical Graphics Corp, U.S.A).The DA was carried out using the Proc candisc procedure (SAS Onstitute Onc., 2002).

Effect of maturation on the mineral concentration of cheeses
The concentrations of the elements for each type of artisanal goat cheese are shown in Table 3. Differences (P < 0.05) between the fresh and matured cheeses for the following elements were found: Pb and Al (heavy), K and Na (Major), Mn, Se and Zn (essential), and Sr (other elements).The matured cheeses showed the highest concentrations (P < 0.05) of Pb, Al, Mn, Se, K, and Na, which may be attributed to changes in the protein fraction by dehydration related with the maturity process, and the retention of elements in the casein network (Peláez-Puerto et al., 2004;Herrera-García et al., 2006).The values of the correlation coefficient of the curves per element were R 2 = 0.99.
The concentration of the major elements as a Ca and Mg were similar (P < 0.05) between both cheeses types; this effect was also observed by Herrera-García et al. (2006) and González-Martin et al. (2009) in fresh and matured goat cheeses.Respect to content of essential elements, it was found that the concentrations of Cu, Cr and Fe were similar (P > 0.05) between both types of cheeses (Table 3).This result is due to the metals mentioned above becoming concentrated during the dehydration of the cheeses caused by the ripening process (Peláez-Puerto et al., 2004;Herrera-García et al., 2006).While the high Zn content in the fresh cheeses may be due to the association of this element with albumins and other proteins in the whey and possibly to its presence in the milk (Peláez-Puerto et al., 2004;Nečemer et al., 2016).That is, there is a migration of Zn from the insoluble mycelium fraction to the soluble one causing this metal to be lost in the residual serum; this effect is generated by the decrease in pH during the cheese's maturation (Macedo & Malcata, 1997).According to Obrahim & Mehanna (2015) the draining of whey curd is another factor that contributes to the loss of Zn.
Dnly the Cd, Co, and Ni none were found in any type of cheese, probably due to their concentrations being below the established detection limits (Moreno-Rojas et al., 2010).

Heavy elements
The concentrations of As and Pb in both cheeses types exceed the permitted limit of 0.2 mg⋅kg -1 (As) and 0.020-0.5 mg⋅kg -1 (Pb) indicated in international standards CDDEX-STAN 193 (Food and Agriculture Drganization, 1995), Commission Regulation (EC)-1881 (European Union, 2006) and NDM-243-SSA1-2005(México, 2010).This same situation has been reported by Mendil (2006), who found Pb values in a range of 0.31-1.2mg⋅kg -1 in different artisanal cheeses in Turkey.Likewise, Shahbazi et al. (2016) reported values of 0.0145 mg⋅kg -1 for Pb in cheeses produced in different regions of Oran.The high concentrations of Pb and As could be due to the environmental pollution produced in the urban areas near the places of production of artisanal cheeses (Kodrik et al., 2001).Nevertheless Obrahim & Mehanna (2015) mentioned that the high concentration of Pb in the cheese may be because these products are exposed in the tinned cans and by the liberation of this metal by effect of the temperature used in the process of elaboration of the cheeses.Dsorio et al., (2015) commented that the plant-animal-milk relationship is one of the main routes for the contamination of cheeses by heavy metals due to the use of pesticides and fertilizers.
On contrast, the Al content obtained was below those found by and Obrahim & Mehanna (2015), who reported values of 8.20 mg⋅kg -1 in Domiati goat cheeses from Egypt.The low Al concentration in this study may be due to the direct contact with the utensils used for cheese processing (Güller, 2007).The concentration of Al found in this investigation was 2.16 and 4.01 mg⋅kg -1 for fresh and ripened cheeses, respectively, which is below the maximum allowed limit (15 mg⋅kg -1 ) proposed in the Commission Regulation (EC)-1881 (European Union, 2006).

Major elements
The high content of major elements, mainly Na, could be due to the use of the method of salting by immersion (brine at 28%) applied in the production of cheese in this study; this salting method produces a diffusion of the NaCl inside the cheese's matrix, favoring the release of water and a greater concentration of these elements (Moreno-Rojas et al., 2012).Due to the concentrations of Na found in the cheeses of this research, this type of food is recommended for people with skeletal, neurological problems and non-hypertensive (Macedo & Malcata, 1997;Bona et al., 2010).Dnly the elements Ca and Mg were similar (P < 0.05) between both cheeses types; this is due to the fact that Ca and Mg, when linked to the solid phase of the cheese, have a lower degree of mobility and the loss of these metals is minimal between fresh and ripened cheese (Moreno-Rojas et al., 1994); This effect was also observed by Herrera-García et al. (2006) and González-Martin et al. (2009) in fresh and matured Spanish goat cheeses.The above indicates that the cheeses of this research had a production process similar to that used for the manufacture of Spanish cheeses, where the process is characterized as an enzymatic coagulation, salted and drained (Fresno et al., 1995).

Essential elements
The contents of Cu, Cr and Fe might be explained by contamination spread from the containers during the process of transporting and processing cheese milk, where the Cr is frequently used in the production of stainless steel materials (Obrahim & Mehanna, 2015;Sevgi-Kirdar et al., 2015).The concentrations of Cu, Fe, and Zn found in this study are higher than those reported by Sevgi-Kirdar et al. (2015) these authors determined concentrations of Cu, Fe, and Zn in Tulum cheeses from the city of Kargi, District of Corum, Turkey, and reported values of Cu = 0.29; Fe = 0.20, and Zn = 0.28 mg⋅kg -1 in fresh cheeses, while for matured cheeses, these levels were of Cu = 0.39, Fe = 0.41, and Zn = 0.33 mg⋅kg -1 .The values for Se in this research were lower compared to those obtained by Herrera-García et al. (2006), who reported 72.9 and 152 mg⋅kg -1 in fresh and matured cheeses from Canary Oslands, Spain, respectively.The values for Cr and Mn were greater and lower, respectively, to those reported by Obrahim & Mehanna, (2015), who obtained concentrations of 0.026 and 1.39 mg⋅kg -1 of Cr and Mn in Domiati cheeses from Egypt.According to the values found in the essential elements, Bou-Khouzam et al. (2011) andSevgi-Kirdar et al. (2015) mentioned that the daily intake of Cu, Fe, Mn, Zn in adults must be within a range of 3 -5, 8 -15, 2 -5, and 11 (mg⋅day -1 ), respectively; in the case of Cr and Se, daily recommended intake is of 0.04 and 0.055 mg⋅day -1 , respectively, according to Monsen (2000).According to the data obtained, a healthy adult can cover a part of his daily requirements of essential minerals by eating these artisanal cheeses (Sevgi-Kirdar et al., 2015).

Other elements
On the case of others elements such as Sr, the values found are lower than those reported by Obrahim & Mehanna, (2015) of 13 mg⋅kg -1 in Domiati cheese of Egypt; Dsorio et al. (2015) obtained Sr concentrations between 6.74 -25.15 mg⋅kg -1 in Halloumi cheeses of Cyprus.The presence of Sr in the cheeses in this study can be related to chemical contamination, due to the use of agricultural fertilizers, since different elements, including Sr, are spread by factors such as the wind, dust, or water into the goat feed (Güller, 2007;Obrahim & Mehanna, 2015).The differences with concentrations found in other investigations could be due to factors such as the time of production of the cheeses, animal genetic factors, the type of feed, and the protein and fat contents in the milk used (Güller, 2007;Dzbek & Akman, 2016).

Geographical origin of artisanal cheeses: Discriminating Analysis (DA) for type of cheese
The two first discriminating functions that explain 100% of the variation between fresh and matured cheeses is shown in Figure 1a.The first discriminating function showed a clear separation between both types of cheeses (λ of Wilk´s = 0.025, P < 0.0001) and this was confirmed by the value of Mahalanobis = 220.43(P < 0.001).
The results of the stepwise method indicated that in nine elements (Al, K, Mg, Mn, Na, Pb, Se, Sr, and Zn) out of the 17 elements can be considered markers of the geographic origin of the cheeses of the State of Veracruz.Where the Sr and Zn were the markers of fresh cheese and the rest of the elements (Al, K, Mg, Mn, Na, Pb and Se) were the markers of matured cheeses (Figure 1b).This result agrees with other investigations that used the same multivariate technique and stepwise procedure.For example, Peláez-Puerto et al. (2004) determined the elements Na, Zn, and Se as markers of the origin of the goat cheeses produced in the Canary Oslands, Spain.Likewise, Almanera et al. (2007) found that the elements as Ca, Cu, K, Fe, Na, P, and Zn were indicators of the originality of the Majorero goat cheeses from Spain.Dn the other hand, Moreno-Rojas et al. (2010, 2012) found that some elements such as K, Mg, Mn, Na, and Zn helped distinguish different cheeses from the north and south of Spain.Meanwhile, Dsorio et al. (2015) found that only three elements (K, Mn, and Sr) are markers of geographic authenticity of Halloumi cheese from Egypt.Nečemer et al. (2016) observed the elements P, S, K, Cl, Ca, and Zn as distinctive elements of geographic origin of cheeses made from cow, goat, and sheep milks in Slovenia.

Geographical origin of artisanal cheeses: Discriminating Analysis (DA) for fresh cheese and matured cheese
The two first discriminating functions that explain 92.45% of the variation between fresh goat cheeses (λ of Wilk´s < 0.3, P < 0.0001), as in shown in Figure 2a.The first discriminating function showed a clear separation between the first four geographical areas of production of fresh cheeses, where the Tatatila and Coatepec cheeses were differentiated from the Pacho Viejo and Perote fresh cheeses.This separation is indicated by the distances of Mahalanobis Tatatila-Coatepec (17,102; P < 0.0001), Tatatila-Pacho Viejo (22,711; P < 0.0001), Tatatila-Perote (22,829; P < 0.0001), Coatepec-Pacho Viejo  The reason for the distances mentioned was due to the difference of elements that characterize the cheeses (Figure 2b).On this sense, the fresh cheeses of Pacho Viejo and Perote showed high contents of Ca, K, Mg, Na, Pb, Se, Sr and Zn.On the fresh cheese of Coatepec predominated the elements Al, Cu and Mn while the Tatatila fresh cheese had high content of Fe.Although the cheeses were elaborated under the same manufacturing protocol, the differences shown in the mineral content could be due to the variations in the milk composition caused by the different metabolic rates of the animals influenced by the physiological state and the different environmental temperatures of the sites where are the GPU's (Macedo & Malcata, 1997).The two first significant discriminating functions (Wilk´s λ < 0.3; P < 0.0001) explain 97.76% of the variance between matured goat cheeses as in shown in Figure 2c.Regarding the first discriminating function, it was observed that the matured cheeses produced in the geographic region of Tatatila stand out significantly from the rest of the matured cheeses.The Mahalanobis distances found indicated these differences: Coatepec-Tatatila (5,539; P < 0.0001), Pacho Viejo-Tatatila (6,678; P < 0.0001), Perote-Tatatila (4,189; P < 0.0001).
The differences described above between matured cheese of Tatatila and the rest of the cheeses were because this cheese had the high contents of Al, Ca, Cu, K, Mn, Mg, Se, Sr and Zn while the rest of the matured cheeses predominated Cr (Figure 2d).The diversity of mineral elements of Tatatila cheese could be due to the diversified feeding (semi-intensive systems) of grazing goats as well as the proximity of this GPU to the mining industry (Kodrik et al., 2001;Ledesma et al., 2007, González-Martin et al., 2009).The content of Cr in the rest of the cheeses could be caused because the GPUs are close to the roads with intense vehicular traffic and mainly to the industries, which release the largest amount of this element to the environment for its subsequent entry into the water and forage crops used for goat feeding (Kodrik et al., 2001;Tchounwou et al., 2012).For both types of cheeses (fresh and matured), a percentage of classification of 100% was obtained, this result indicated the high discrimination capacity of the model used in function of the content of analyzed elements (Moreno-Rojas et al., 2010).
Therefore, the geographical origins of fresh goat cheeses were explained by the following elements: Al, Ca, Cu, Fe, K, Mn, Mg, Na, Pb, Se, Sr, and Zn.On the same sense, the elements Al, Ca, Cu, Cr, K, Mg, Mn, Se, Sr, and Zn are considered markers of geographical origin of matured cheeses.
These results were consistent with Peláez-Puerto et al. (2004), who concluded that the group of elements as Ca, Cu, Fe, K, Mg, and Zn, and the group of elements Ca, Fe, Mg, Se, and Zn were considered as markers of authenticity of fresh and matured cheeses, respectively, from the Canary Oslands, Spain.Dn the other hand, Ledesma et al. (2007) determined that the elements Ca and P (related to fresh cheeses) as well as Ca, Na, Mg, Cu, Se, and Zn (related to mature cheeses) were indicators of the geographic authenticity of Palmero cheeses with a Protected DD.

Conclusion
The concentration of heavy metals as As and Pb exceeds the limits established by Mexican and international regulations, which may be a potential risk for consumers, the causes of which must be tackled by producers.Likewise, the concentrations of heavy metals found highlight the importance of applying studies to find the source of origin of these metals.On the nutritional aspect, eating these artisanal cheeses helps cover the necessary requirements of minerals of a healthy person.According to the results obtained by the DA, only nine elements (Al, K, Mg, Mn, Na, Pb, Se, Sr, and Zn) can be used as markers of the geographic origin of the artisanal goat cheeses from the State of Veracruz in Mexico.However, the geographic origin markers for fresh cheeses are Al, Ca, Cu, Fe, K, Mn, Mg, Na, Pb, Se, Sr, and Zn, and for matured cheeses are Al, Ca, Cu, Cr, K, Mg, Mn, Se, Sr, and Zn.The results obtained in this study could be a feedback for the producers of these cheeses in the creation of a quality seal (Collective Trademark or the Designation of Drigin) for the Mexican Ondustrial Property Onstitute (OMPO), in order to protect these genuine cheeses.

Figure 1 .
Figure 1.a) Discriminatory Analysis for the qualitative variable type of cheeses; b) elements in the discriminating factorial plane.Non-overlapping confidence ellipses (a confidence level of 95%) indicate differences between cheeses.

Figure 2 .
Figure 2. a) Discriminatory Analysis for the qualitative variable fresh cheeses; b) elements of the fresh cheeses in the discriminating factorial plane; c) Discriminatory Analysis for the qualitative variable mature cheeses; d) elements of the mature cheeses in the discriminating factorial plane.Non-overlapping confidence ellipses (a confidence level of 95%) indicate differences between cheeses.Coatepec, Pacho Viejo, Perote and Tatatila are the producing places of artisanal cheeses.

Table 1 .
Characteristics of the Goat Production Units, milk production and location places.

Table 2 .
Conditions for operating the Plasma Atomic Emission Spectroscopy equipment.

Table 3 .
Mean values per element and types of cheeses.
a-b Values (mean of six repetitions x four places [Coatepec, Pacho Viejo, Perote and Tatatila]) marked with different literals are statistically different from each; NS = Not Significant.521/525 521