# Geographical Origin Identification of Tequila Based on Multielement and Stable Isotopes

Authors: Jian Zhang; Yiqian Ma; Feilong Shao; Yuanyu Lu; Libo Zhu; Xiansong Xu
Journal: Journal of Food Quality
 (2021)

---

## Abstract

To explore the regional characteristics of multielement and stable isotopes in tequila and the viability of tracing the geographical origin, both inductively coupled plasma mass spectrometry (ICP-MS) and stable isotope ratio mass spectrometry (IRMS) were applied in this study to determine the multielement content and stable isotope ratios of tequila. Then, principal component analysis (PCA) and partial least squares discriminant analysis (PLS-DA) as chemometrics were conducted in combination to establish a model for determining the geographical origin of Mexican tequila. According to the analytical results, the content of 22 elements varied, as did the ranges of stable isotopeδ13C and δ18O values for the tequila obtained from the 4 producing areas. Besides, there were regional characteristics shown to some extent. PCA and PLS-DA methods can be adopted to identify the geographical origin of tequila in Mexico. According to the discriminant verification of 16 samples, the accuracy of prediction reached 93.75%. It has been demonstrated that it is possible to distinguish tequila of different geographical origins in Mexico by combining the determination of the multielement content and stable isotope ratios with chemometric analysis, which lays a foundation for tracing the geographical origin of tequila.


---

## Body

## 1. Introduction

Tequila refers to a distilled alcoholic product made from blue agave grown across the 4 legal states of Mexico including Jalisco, Guanajuato, Michoacan, and Tamaulipas. In 2014, the Chinese market was opened up to Mexican tequila, which made the country of origin (geographical indication) trademark of Mexican tequila recognized. As the imports and consumption of tequila by Chinese consumers are increasing on a continued basis in China, there has been a phenomenon arising where illegal merchants bring fake or shoddy products on the market for sale. Consequently, the quality of tequila in the domestic market is made uneven, with a disconnect between quality and pricing emerging. Therefore, it is imperative to develop new technologies that are effective, fast, and accurate in identifying the geographical origin and quality of tequila.At present, such solutions as mineral element fingerprint method [1–4], aroma fingerprint method [5], and stable isotope mass spectrometry (IRMS) technology [6–8] have been commonly applied for food production. Among them, the stable isotope ratio method performs well in indicating regional and climatic characteristics. It is widely used to trace food geographical origin, for example, wine [9], beef [10], and tea [11, 12]. Besides, it can be applied for authentication, such as the brewing of soy sauce [13] and honey [14–16]. The content of mineral elements in tequila is largely determined by the environmental conditions of its growth, such as soil, water, and climate. Since it is difficult to metabolize or change them during the brewing and aging process, they can be taken as an important marker used to identify the geographical origin of tequila. It is possible to use the carbon 13 isotopic ratio (δ13CVPDB), oxygen 18 isotopic ratio (δ18OVSMOW), and the element content as complementary parameters of authenticity to identify tequila and tequila 100% agave (silver class) [17, 18]. Using the 49 mineral elements combined with the stable isotope ratios of C, N, O, S, and Sr, Bontempo et al. [19] identified the geographical origin of 10 cheese samples obtained from the Alps; the correct classification rate of the samples reached up to 94%. Shao Shengzhi et al. [20] combined principal component analysis (PCA) and linear discriminant analysis (LDA) as chemometrics. They used the contents of 23 mineral elements and the stable isotope ratios of C, N, O, and H to conduct discriminant verification on 23 rice samples obtained from Heilongjiang, Jiangsu, and Liaoning provinces. Ultimately, the accuracy reached as high as 91%. Combined with the stoichiometric PCA-LDA, Yuan Yuwei et al. [11] used the content of 27 mineral elements and the stable isotope ratios of C, N, O, and H to identify the tea obtained from Fujian, Shandong, and Zhejiang, with the accuracy rate reaching 99%. It can be judged from above that the use of stable isotopes and mineral elements in combination with stoichiometry can produce a better outcome in tracing the geographical origin of tea.Being one of the multivariate analysis techniques, PCA (principal component analysis) is frequently conducted as a dimension-reduction tool to reduce a large set of variables to a small set that remains capable of containing most of the information in the large set. Considering this, PCA was employed for the transformation of largely correlated variables of this work into smaller uncorrected variables [21]. PLS-DA (partial least squares discriminant analysis) is a linear classification method that combines the properties shown by partial least square regression with the discrimination power of a classification technique [22].In this study, mineral elements and stable isotope ratios were combined with principal component analysis (PCA) and partial least squares discriminant analysis (PLS-DA) as chemometrics. On this basis, a solution was developed to trace and identify the geographical origin of tequila obtained from the 4 legal regions of Mexico, which provides technical support for ensuring the authenticity of the Mexican tequila imported to China.


## 2. Materials and Methods

### 2.1. Instruments and Reagents

DRC-e Inductively Coupled Plasma Mass Spectrometer (ICP-MS) is from PerkinElmer, USA. DELTA V Plus Isotope Mass Spectrometer (IRMS) is from Thermo Fisher Company, USA, equipped with Trace GC gas chromatograph. The instrument control software used in this study is ISODATE 3.0. GasBench II preprocessing devices are GC PAL autosampler, PoraPlotQ column (25 m, 0.32 mm), screw cap vial, injection needle, and blowing needle. Gas chromatography column is RTX-WAX (30 m × 0.25 mm × 0.25 μm). Carbon isotope standard material in carbon black is δ13CVPDB = (−36.91 ± 0.10)%. Hydrogen and oxygen isotope water standard material GBW04458 is δ18OVSMOW= (−0.15 ± 0.07)‰ and GBW04459 is δ18OVSMOW = (−8.61 ± 0.08)%, with high purity helium (99.999%); standard reference gas CO2 was purchased from China Shengying Company. Multielement standard solution is from PerkinElmer, USA. Milli-Q Gradient pure water treatment system is from Millipore, France.


### 2.2. Materials

The importers were entrusted with the purchase of totally 60 original bottles of tequila as obtained from 4 legal production areas in Mexico, with 13 bottles from Jalisco, 15 bottles from Guanajuato, 15 bottles from Michoacan, and 17 bottles from Tamaulipas. Among them, 44 tequila wines were used for constructing the geographical origin traceability model, and the remaining 16 were used to carry out model verification.


### 2.3. Sample Preparation

Firstly, 25 mL of tequila sample was placed into a water bath at 100°C for the complete removal of alcohol through evaporation. Secondly, it was being diluted to 25 mL again by 1% nitric acid solution for later multielement analysis [23]. Thirdly, 20 μL of tequila sample was placed into the gas chromatography bottle and then was added with 1 mL of chromatographic pure acetone to analyze the stable carbon isotope ratio of ethanol. Besides, 0.5 mL of tequila sample was placed into a screw cap sample bottle, which was filled with a mixture of 0.3% CO2 and He for 300 s to remove the air present in the bottle. Then, it was equilibrated for 18–24 hours. Finally, the gaseous CO2 and oxygen present in the water were exchanged for an isotope, to facilitate the analysis of the stable oxygen isotope in the water.


### 2.4. Instrument Conditions

#### 2.4.1. ICP-MS Measurement Conditions [23]

RF power is 1100 W. Nebulizer flow rate is 0.88 L/min. Cooling air flow rate is 16 L/min. Auxiliary gas flow rate is 1.2 L/min. Ion lens voltage is 6 V. As for the determination of arsenic element, a dynamic reaction cell mode was adopted. Parameters are cell gas A: 0.3; RPq: 0.5.


#### 2.4.2. IRMS Measurement Conditions

(1) Stable Carbon Isotope Determination Conditions. Column is RTX-WAX (30 m × 0.25  mm, 0.25 μm). Injection volume is 1 μL. Carrier gas is high purity. Flow rate is 1.2 mL/min. Split ratio is 30 : 1. Inlet temperature is 200°C. The heating program is detailed as follows. The initial temperature is 40°C, which is kept for 5 min. Then, it raised to 80°C at a rate of 10°C/min, which is kept for 1 min. Finally, it raised to 200°C at a rate of 25°C/min, which was kept for 15 min. Combustion tube temperature is 1000°C. Ion source voltage is 2.97 KV. Vacuum degree is 1.6 × 10−6 mBar. Bombardment voltage is 123.8 eV.(2) Stable Oxygen Isotope Measurement Conditions. Carrier gas is helium. Carrier gas pressure is 0.9 bar. Reference gas is CO2. Reference gas pressure is 1.3 bar. Column temperature is 65°C. There is a Polar Plot Q column. Ion source voltage is 2.97 KV. Vacuum degree is 1.6 × 10−6 mBar. Bombardment voltage is 123.8 eV.


### 2.5. Data Processing

Thermo Electron Isodat version 3.0 software was applied to perform isotope ratio calculation and data acquisition. Principal component analysis (PCA) and partial least squares discriminant analysis (PLS-DA) were conducted usingR language. Origin 9.0 (Origin lab, USA) was adopted to carry out data mapping.


## 2.1. Instruments and Reagents

DRC-e Inductively Coupled Plasma Mass Spectrometer (ICP-MS) is from PerkinElmer, USA. DELTA V Plus Isotope Mass Spectrometer (IRMS) is from Thermo Fisher Company, USA, equipped with Trace GC gas chromatograph. The instrument control software used in this study is ISODATE 3.0. GasBench II preprocessing devices are GC PAL autosampler, PoraPlotQ column (25 m, 0.32 mm), screw cap vial, injection needle, and blowing needle. Gas chromatography column is RTX-WAX (30 m × 0.25 mm × 0.25 μm). Carbon isotope standard material in carbon black is δ13CVPDB = (−36.91 ± 0.10)%. Hydrogen and oxygen isotope water standard material GBW04458 is δ18OVSMOW= (−0.15 ± 0.07)‰ and GBW04459 is δ18OVSMOW = (−8.61 ± 0.08)%, with high purity helium (99.999%); standard reference gas CO2 was purchased from China Shengying Company. Multielement standard solution is from PerkinElmer, USA. Milli-Q Gradient pure water treatment system is from Millipore, France.


## 2.2. Materials

The importers were entrusted with the purchase of totally 60 original bottles of tequila as obtained from 4 legal production areas in Mexico, with 13 bottles from Jalisco, 15 bottles from Guanajuato, 15 bottles from Michoacan, and 17 bottles from Tamaulipas. Among them, 44 tequila wines were used for constructing the geographical origin traceability model, and the remaining 16 were used to carry out model verification.


## 2.3. Sample Preparation

Firstly, 25 mL of tequila sample was placed into a water bath at 100°C for the complete removal of alcohol through evaporation. Secondly, it was being diluted to 25 mL again by 1% nitric acid solution for later multielement analysis [23]. Thirdly, 20 μL of tequila sample was placed into the gas chromatography bottle and then was added with 1 mL of chromatographic pure acetone to analyze the stable carbon isotope ratio of ethanol. Besides, 0.5 mL of tequila sample was placed into a screw cap sample bottle, which was filled with a mixture of 0.3% CO2 and He for 300 s to remove the air present in the bottle. Then, it was equilibrated for 18–24 hours. Finally, the gaseous CO2 and oxygen present in the water were exchanged for an isotope, to facilitate the analysis of the stable oxygen isotope in the water.


## 2.4. Instrument Conditions

### 2.4.1. ICP-MS Measurement Conditions [23]

RF power is 1100 W. Nebulizer flow rate is 0.88 L/min. Cooling air flow rate is 16 L/min. Auxiliary gas flow rate is 1.2 L/min. Ion lens voltage is 6 V. As for the determination of arsenic element, a dynamic reaction cell mode was adopted. Parameters are cell gas A: 0.3; RPq: 0.5.


### 2.4.2. IRMS Measurement Conditions

(1) Stable Carbon Isotope Determination Conditions. Column is RTX-WAX (30 m × 0.25  mm, 0.25 μm). Injection volume is 1 μL. Carrier gas is high purity. Flow rate is 1.2 mL/min. Split ratio is 30 : 1. Inlet temperature is 200°C. The heating program is detailed as follows. The initial temperature is 40°C, which is kept for 5 min. Then, it raised to 80°C at a rate of 10°C/min, which is kept for 1 min. Finally, it raised to 200°C at a rate of 25°C/min, which was kept for 15 min. Combustion tube temperature is 1000°C. Ion source voltage is 2.97 KV. Vacuum degree is 1.6 × 10−6 mBar. Bombardment voltage is 123.8 eV.(2) Stable Oxygen Isotope Measurement Conditions. Carrier gas is helium. Carrier gas pressure is 0.9 bar. Reference gas is CO2. Reference gas pressure is 1.3 bar. Column temperature is 65°C. There is a Polar Plot Q column. Ion source voltage is 2.97 KV. Vacuum degree is 1.6 × 10−6 mBar. Bombardment voltage is 123.8 eV.


## 2.4.1. ICP-MS Measurement Conditions [23]

RF power is 1100 W. Nebulizer flow rate is 0.88 L/min. Cooling air flow rate is 16 L/min. Auxiliary gas flow rate is 1.2 L/min. Ion lens voltage is 6 V. As for the determination of arsenic element, a dynamic reaction cell mode was adopted. Parameters are cell gas A: 0.3; RPq: 0.5.


## 2.4.2. IRMS Measurement Conditions

(1) Stable Carbon Isotope Determination Conditions. Column is RTX-WAX (30 m × 0.25  mm, 0.25 μm). Injection volume is 1 μL. Carrier gas is high purity. Flow rate is 1.2 mL/min. Split ratio is 30 : 1. Inlet temperature is 200°C. The heating program is detailed as follows. The initial temperature is 40°C, which is kept for 5 min. Then, it raised to 80°C at a rate of 10°C/min, which is kept for 1 min. Finally, it raised to 200°C at a rate of 25°C/min, which was kept for 15 min. Combustion tube temperature is 1000°C. Ion source voltage is 2.97 KV. Vacuum degree is 1.6 × 10−6 mBar. Bombardment voltage is 123.8 eV.(2) Stable Oxygen Isotope Measurement Conditions. Carrier gas is helium. Carrier gas pressure is 0.9 bar. Reference gas is CO2. Reference gas pressure is 1.3 bar. Column temperature is 65°C. There is a Polar Plot Q column. Ion source voltage is 2.97 KV. Vacuum degree is 1.6 × 10−6 mBar. Bombardment voltage is 123.8 eV.


## 2.5. Data Processing

Thermo Electron Isodat version 3.0 software was applied to perform isotope ratio calculation and data acquisition. Principal component analysis (PCA) and partial least squares discriminant analysis (PLS-DA) were conducted usingR language. Origin 9.0 (Origin lab, USA) was adopted to carry out data mapping.


## 3. Results and Discussion

### 3.1. Precision and Accuracy

In our previous study [23], an investigation was conducted into the precision and accuracy of ICP-MS. On this basis, a total of 22 elements in tequila were determined through standard recovery experiments. For each element to be tested, the standard addition recovery rate falls within the range of 90∼110%, which is consistent with the results of the 8 elements found in the previous study of food carried out by this group [24]. It was finally demonstrated that ICP-MS can meet the detection requirements of various elements.According to the quality control requirements specified for international testing methods [25], the measurement deviation for the δ13C value of quality control standards and samples is required not to exceed 0.5% and 0.3%, respectively. The carbon isotope and oxygen isotope calibration of the reference gas CO2 was conducted using the carbon black standard material as well as the hydrogen and oxygen isotope water standard material. As shown in Table 1, the standard deviation shown by the δ13C value of the reference gas CO2 is 0.016%, and the standard deviation of the oxygen isotope ratio δ18O is less than 0.06%, which meets the testing requirements of the experimental requirements.Table 1
Precision ofδ13C values and δ18O values of the standard.
NumberCarbon black standardδ13C (%)GBW04458δ18O (%)GBW04459δ18O (%)1−21.73−30.20−30.342−21.71−30.09−30.383−21.74−30.15−30.334−21.70−30.17−30.305−21.73−30.12−30.356−21.75−30.13−30.377−21.73−30.23−30.378−21.74−30.10−30.41Means−21.73−30.15−30.36SD0.0160.050.03∗The data in the table has been automatically calibrated by the instrument to the value of the international standard V-PDB and V-SMOW.The same sample was used to repeatedly examine the accuracy of the experimental method. According to the abovementioned experimental conditions and procedures, the carbon isotope ratio of ethanol in the same sample was measured 8 times by GC-IRMS, while the oxygen isotope ratio of water in the same sample was measured 8 times by GasBench II-IRMS. Then, the standard deviation was calculated, and the results were obtained, as listed in Table2. The standard deviations shown by the δ13C value of ethanol and the δ18O value of water in tequila are 0.12% and 0.03%, respectively, which meets the accuracy requirements of the experiment.Table 2
Accuracy ofδ13C values and δ18O values of tequila by isotope ratio mass spectrometry.
Numberδ13C (%)δ18O (%)1−12.42−13.212−12.47−13.263−12.31−13.234−12.54−13.295−12.57−13.246−12.36−13.247−12.65−13.278−12.59−13.22Mean−12.49−13.24SD0.120.03∗The data shown in the table has been automatically calibrated by the instrument to the value of the international standard V-PDB and V-SMOW.


### 3.2. Geographical Origin Traceability Analysis

In order to assess the multielement content in tequila, carbon and oxygen stable isotope ratio analysis was conducted of their effectiveness in tracing the geographical origin of tequila obtained from 4 legal production areas in Mexico, including Jalisco, Guanajuato, Michoacan, and Tamaulipas. This method was used to determine the multielement content in tequila, the carbon isotope ratio of ethanol, and the oxygen isotope ratio of water. As indicated by the analytical results, there are variations in the multielement content of tequila obtained from different geographical origins as shown in Table3 [23]. Besides the carbon and oxygen isotope ratio shown Table 4, the δ13C experimental data are in a range of −14.29% and −12.60%. Similar results have been reported for a range of −12.1% to −13.2% [26] and a range of −11.0% to −14.0% [17], which suggests that the multielement content as well as carbon and oxygen isotope ratios of tequila varies between different geographical origins. Therefore, these can be used to assist geographical origin traceability analysis.Table 3
Element contents of tequila by ICP-MS.
ElementSample1234567891011Li0.7021.7661.9541.4493.1642.7663.0740.2710.4210.5230.572Na701.4512473.3422572.2312888.8094844.8754723.9824928.5223314.1011451.0762117.1612172.331Mg39.605126.94528.22996.34858.16825.74967.13737.358142.561108.0149.467Al15.0062.3130.8929.4227.78313.32616.81834.997.61740.57217.984K125.784221.6341244.744583.98363.914368.183322.903672.762144.387155.049267.881Ca106.771276.097177.873216.988195.95474.38193.409190.63475.065398.741242.484V0.1220.1610.080.1160.0560.0620.0920.530.1950.9750.214Cr0.6021.5790.8871.11811.73317.916.5791.0024.72416.33812.874Mn1.6714.520.7231.8021.9250.8422.683.6663.3187.2672.715Fe0.6815.7311.379.7913.9084.4235.64410.7086.67812.12913.901Co0.0690.0860.0480.0970.0840.0810.1260.1430.0950.5190.464Ni0.3440.2530.1490.5771.7310.2560.4660.4460.4321.6310.936Cu20.022420.7614.858978.1824.7284.5574.191328.7930.7591.4381.113Zn1.7957.4240.20226.8462.220.9372.0072.2192.1393.2363.143As1.2460.8530.4870.4830.3920.4130.2910.2190.810.9730.948Se0.2070.1450.1450.1040.2170.1630.09800.5820.6840.559Ag0.070.0330.0140.0110.0110.0160.0140.0130.00400Cd0.0160.02500.03400.01200.0070.00700.009Ba00.84901.6082.22600.5262.1092.93900Pb0.0310.2020.2770.8750.2270.1080.2880.020.2864.2424.431Bi0.4020.1050.0220.0080000000Sr0.3892.230.2611.2281.0760.3751.551.5186.4867.183.293Table 4
δ13C values and δ18O values of stable isotope ratio in tequila.
Geographical originEthanolδ13C valueWaterδ18O valueGuanajuato−13.67−13.24Guanajuato−13.56−13.11Guanajuato−13.62−13.29Guanajuato−13.5−13.38Guanajuato−13.71−13.44Michoacan−14.04−9.76Michoacan−14.13−9.62Michoacan−14.18−9.8Michoacan−14.3−9.72Michoacan−12.74−12.11Michoacan−12.6−11.97Jalisco−14.12−9.32Jalisco−14.29−9.41Jalisco−14.18−9.36Jalisco−14.09−9.33Jalisco−13.86−9.12Tamaulipas−14.09−9.26Tamaulipas−14.16−9.48Tamaulipas−13.75−13.14Tamaulipas−13.52−13.02Tamaulipas−13.83−13.13Tamaulipas−13.68−13.22


### 3.3. Model Building

There are significant differences exhibited by the results obtained from the single-factor analysis of variance for the stable carbon isotopeδ13C and oxygen isotope δ18O in tequila obtained from the four tequila-producing areas, the results of principal component analysis for the content of 22 metal elements, and the results of six principal components’ analysis for a metallic element. Table 5 indicates the cumulative contribution rate of the principal component analysis of 22 variances, respectively. Normally, parameters need to exceed 75%, and the characteristic root value is greater than or equal to 1 of the principal component, which will be deemed sufficient [20]. The biggest contributor is made up of 6 elements (Sr, Pb, Se, V, Ca, and Fe) with 91.59% contributions in this PCA model. Therefore, the geographical origin traceability model can be constructed to identify the exact geographical origin. Partial least squares discriminant analysis (PLS-DA) was performed on 44 samples obtained from 4 production areas, including Guanajuato, Jalisco, Michoacan, and Tamaulipas. The results are shown in Figures 1 and 2, respectively.Table 5
Component score coefficient matrix in tequila.
ElementsComponent123456Li−0.072−0.137−0.023−0.0820.3140.051Na−0.042−0.2090.0280.0300.1860.211Mg0.0650.0810.166-0.2130.1440.036Al0.080−0.038−0.0450.221−0.1720.383K−0.060−0.0800.0510.109−0.350−0.286Ca0.0960.0280.116−0.228−0.063−0.054V0.101−0.005−0.0020.131−0.1910.257Cr0.058−0.141−0.1190.0130.290−0.012Mn0.1050.0250.0370.002−0.0370.278Fe0.093−0.0430.1060.207−0.008−0.066Co0.111−0.030−0.0650.1580.041−0.128Ni0.074−0.083−0.0220.0130.2320.189Cu−0.0240.0710.2420.2230.0900.003Zn−0.0090.0680.2240.2040.224−0.118As0.0560.185−0.125−0.0420.109−0.150Se0.1080.021−0.063−0.1330.072−0.228Ag−0.0500.167−0.1400.0540.0620.256Cd−0.0270.1700.0920.1440.2370.029Ba0.000−0.0120.212−0.214−0.0750.249Pb0.103−0.018−0.0600.1690.091−0.271Bi−0.0330.200−0.1410.0180.0300.169Sr0.1150.0190.045−0.158−0.031−0.013∗Extraction method: principal component analysis.Figure 1
The two-dimensional figure of the PLS-DA model of identification of Mexican tequila.Figure 2
The three-dimensional figure of the PLS-DA model of identification of Mexican tequila.As shown in Figures1 and 2, the tequila in the four regions can be effectively distinguished by stable carbon, oxygen isotope, and multielement analysis. Tamaulipas is concentrated in the area formed by the negative semiaxis of PC1 and PC2. In comparison the remaining three regions (Guanajuato, Jalisco, and Michoacan) are concentrated in the area of the positive semiaxis of PC1. These can be differentiated more effectively.


### 3.4. Model Validation

According to Table6, one sample obtained from Guanajuato is misjudged as Tamaulipas, with the prediction accuracy reaching 75%. The other three production areas are all distinguished correctly, with the prediction accuracy reaching 100%. Therefore, the model is verified as reliable, with an overall prediction accuracy of 93.75%. Despite this, an effort will be continued in the future to increase the sample size and expand the existing database. Meanwhile, an attempt will be made to achieve technical cooperation with Mexican tequila related laboratories for acquiring more technical information about local tequila and tequila, such as local climatic data and local oxygen isotope data in surface water. In this way, a data base can be constructed that is more reliable, informative, and closer to the real-world local conditions.Table 6
PLS-DA model validation of tequila from four regions of Mexico.
Geographical originGuanajuatoJaliscoMichoacanTamaulipasGuanajuato4000Jalisco0300Michoacan1030Tamaulipas0005Prediction accuracy rate (%)75100100100


## 3.1. Precision and Accuracy

In our previous study [23], an investigation was conducted into the precision and accuracy of ICP-MS. On this basis, a total of 22 elements in tequila were determined through standard recovery experiments. For each element to be tested, the standard addition recovery rate falls within the range of 90∼110%, which is consistent with the results of the 8 elements found in the previous study of food carried out by this group [24]. It was finally demonstrated that ICP-MS can meet the detection requirements of various elements.According to the quality control requirements specified for international testing methods [25], the measurement deviation for the δ13C value of quality control standards and samples is required not to exceed 0.5% and 0.3%, respectively. The carbon isotope and oxygen isotope calibration of the reference gas CO2 was conducted using the carbon black standard material as well as the hydrogen and oxygen isotope water standard material. As shown in Table 1, the standard deviation shown by the δ13C value of the reference gas CO2 is 0.016%, and the standard deviation of the oxygen isotope ratio δ18O is less than 0.06%, which meets the testing requirements of the experimental requirements.Table 1
Precision ofδ13C values and δ18O values of the standard.
NumberCarbon black standardδ13C (%)GBW04458δ18O (%)GBW04459δ18O (%)1−21.73−30.20−30.342−21.71−30.09−30.383−21.74−30.15−30.334−21.70−30.17−30.305−21.73−30.12−30.356−21.75−30.13−30.377−21.73−30.23−30.378−21.74−30.10−30.41Means−21.73−30.15−30.36SD0.0160.050.03∗The data in the table has been automatically calibrated by the instrument to the value of the international standard V-PDB and V-SMOW.The same sample was used to repeatedly examine the accuracy of the experimental method. According to the abovementioned experimental conditions and procedures, the carbon isotope ratio of ethanol in the same sample was measured 8 times by GC-IRMS, while the oxygen isotope ratio of water in the same sample was measured 8 times by GasBench II-IRMS. Then, the standard deviation was calculated, and the results were obtained, as listed in Table2. The standard deviations shown by the δ13C value of ethanol and the δ18O value of water in tequila are 0.12% and 0.03%, respectively, which meets the accuracy requirements of the experiment.Table 2
Accuracy ofδ13C values and δ18O values of tequila by isotope ratio mass spectrometry.
Numberδ13C (%)δ18O (%)1−12.42−13.212−12.47−13.263−12.31−13.234−12.54−13.295−12.57−13.246−12.36−13.247−12.65−13.278−12.59−13.22Mean−12.49−13.24SD0.120.03∗The data shown in the table has been automatically calibrated by the instrument to the value of the international standard V-PDB and V-SMOW.


## 3.2. Geographical Origin Traceability Analysis

In order to assess the multielement content in tequila, carbon and oxygen stable isotope ratio analysis was conducted of their effectiveness in tracing the geographical origin of tequila obtained from 4 legal production areas in Mexico, including Jalisco, Guanajuato, Michoacan, and Tamaulipas. This method was used to determine the multielement content in tequila, the carbon isotope ratio of ethanol, and the oxygen isotope ratio of water. As indicated by the analytical results, there are variations in the multielement content of tequila obtained from different geographical origins as shown in Table3 [23]. Besides the carbon and oxygen isotope ratio shown Table 4, the δ13C experimental data are in a range of −14.29% and −12.60%. Similar results have been reported for a range of −12.1% to −13.2% [26] and a range of −11.0% to −14.0% [17], which suggests that the multielement content as well as carbon and oxygen isotope ratios of tequila varies between different geographical origins. Therefore, these can be used to assist geographical origin traceability analysis.Table 3
Element contents of tequila by ICP-MS.
ElementSample1234567891011Li0.7021.7661.9541.4493.1642.7663.0740.2710.4210.5230.572Na701.4512473.3422572.2312888.8094844.8754723.9824928.5223314.1011451.0762117.1612172.331Mg39.605126.94528.22996.34858.16825.74967.13737.358142.561108.0149.467Al15.0062.3130.8929.4227.78313.32616.81834.997.61740.57217.984K125.784221.6341244.744583.98363.914368.183322.903672.762144.387155.049267.881Ca106.771276.097177.873216.988195.95474.38193.409190.63475.065398.741242.484V0.1220.1610.080.1160.0560.0620.0920.530.1950.9750.214Cr0.6021.5790.8871.11811.73317.916.5791.0024.72416.33812.874Mn1.6714.520.7231.8021.9250.8422.683.6663.3187.2672.715Fe0.6815.7311.379.7913.9084.4235.64410.7086.67812.12913.901Co0.0690.0860.0480.0970.0840.0810.1260.1430.0950.5190.464Ni0.3440.2530.1490.5771.7310.2560.4660.4460.4321.6310.936Cu20.022420.7614.858978.1824.7284.5574.191328.7930.7591.4381.113Zn1.7957.4240.20226.8462.220.9372.0072.2192.1393.2363.143As1.2460.8530.4870.4830.3920.4130.2910.2190.810.9730.948Se0.2070.1450.1450.1040.2170.1630.09800.5820.6840.559Ag0.070.0330.0140.0110.0110.0160.0140.0130.00400Cd0.0160.02500.03400.01200.0070.00700.009Ba00.84901.6082.22600.5262.1092.93900Pb0.0310.2020.2770.8750.2270.1080.2880.020.2864.2424.431Bi0.4020.1050.0220.0080000000Sr0.3892.230.2611.2281.0760.3751.551.5186.4867.183.293Table 4
δ13C values and δ18O values of stable isotope ratio in tequila.
Geographical originEthanolδ13C valueWaterδ18O valueGuanajuato−13.67−13.24Guanajuato−13.56−13.11Guanajuato−13.62−13.29Guanajuato−13.5−13.38Guanajuato−13.71−13.44Michoacan−14.04−9.76Michoacan−14.13−9.62Michoacan−14.18−9.8Michoacan−14.3−9.72Michoacan−12.74−12.11Michoacan−12.6−11.97Jalisco−14.12−9.32Jalisco−14.29−9.41Jalisco−14.18−9.36Jalisco−14.09−9.33Jalisco−13.86−9.12Tamaulipas−14.09−9.26Tamaulipas−14.16−9.48Tamaulipas−13.75−13.14Tamaulipas−13.52−13.02Tamaulipas−13.83−13.13Tamaulipas−13.68−13.22


## 3.3. Model Building

There are significant differences exhibited by the results obtained from the single-factor analysis of variance for the stable carbon isotopeδ13C and oxygen isotope δ18O in tequila obtained from the four tequila-producing areas, the results of principal component analysis for the content of 22 metal elements, and the results of six principal components’ analysis for a metallic element. Table 5 indicates the cumulative contribution rate of the principal component analysis of 22 variances, respectively. Normally, parameters need to exceed 75%, and the characteristic root value is greater than or equal to 1 of the principal component, which will be deemed sufficient [20]. The biggest contributor is made up of 6 elements (Sr, Pb, Se, V, Ca, and Fe) with 91.59% contributions in this PCA model. Therefore, the geographical origin traceability model can be constructed to identify the exact geographical origin. Partial least squares discriminant analysis (PLS-DA) was performed on 44 samples obtained from 4 production areas, including Guanajuato, Jalisco, Michoacan, and Tamaulipas. The results are shown in Figures 1 and 2, respectively.Table 5
Component score coefficient matrix in tequila.
ElementsComponent123456Li−0.072−0.137−0.023−0.0820.3140.051Na−0.042−0.2090.0280.0300.1860.211Mg0.0650.0810.166-0.2130.1440.036Al0.080−0.038−0.0450.221−0.1720.383K−0.060−0.0800.0510.109−0.350−0.286Ca0.0960.0280.116−0.228−0.063−0.054V0.101−0.005−0.0020.131−0.1910.257Cr0.058−0.141−0.1190.0130.290−0.012Mn0.1050.0250.0370.002−0.0370.278Fe0.093−0.0430.1060.207−0.008−0.066Co0.111−0.030−0.0650.1580.041−0.128Ni0.074−0.083−0.0220.0130.2320.189Cu−0.0240.0710.2420.2230.0900.003Zn−0.0090.0680.2240.2040.224−0.118As0.0560.185−0.125−0.0420.109−0.150Se0.1080.021−0.063−0.1330.072−0.228Ag−0.0500.167−0.1400.0540.0620.256Cd−0.0270.1700.0920.1440.2370.029Ba0.000−0.0120.212−0.214−0.0750.249Pb0.103−0.018−0.0600.1690.091−0.271Bi−0.0330.200−0.1410.0180.0300.169Sr0.1150.0190.045−0.158−0.031−0.013∗Extraction method: principal component analysis.Figure 1
The two-dimensional figure of the PLS-DA model of identification of Mexican tequila.Figure 2
The three-dimensional figure of the PLS-DA model of identification of Mexican tequila.As shown in Figures1 and 2, the tequila in the four regions can be effectively distinguished by stable carbon, oxygen isotope, and multielement analysis. Tamaulipas is concentrated in the area formed by the negative semiaxis of PC1 and PC2. In comparison the remaining three regions (Guanajuato, Jalisco, and Michoacan) are concentrated in the area of the positive semiaxis of PC1. These can be differentiated more effectively.


## 3.4. Model Validation

According to Table6, one sample obtained from Guanajuato is misjudged as Tamaulipas, with the prediction accuracy reaching 75%. The other three production areas are all distinguished correctly, with the prediction accuracy reaching 100%. Therefore, the model is verified as reliable, with an overall prediction accuracy of 93.75%. Despite this, an effort will be continued in the future to increase the sample size and expand the existing database. Meanwhile, an attempt will be made to achieve technical cooperation with Mexican tequila related laboratories for acquiring more technical information about local tequila and tequila, such as local climatic data and local oxygen isotope data in surface water. In this way, a data base can be constructed that is more reliable, informative, and closer to the real-world local conditions.Table 6
PLS-DA model validation of tequila from four regions of Mexico.
Geographical originGuanajuatoJaliscoMichoacanTamaulipasGuanajuato4000Jalisco0300Michoacan1030Tamaulipas0005Prediction accuracy rate (%)75100100100


## 4. Conclusion

After the data was measured by the stable carbon isotopeδ13C of ethanol and the stable oxygen isotope δ18O of water in tequila, ICP-MS was applied to analyze the measurement results of 22 metal elements in tequila. On this base, the geographical origin model was constructed, the overall prediction accuracy of which reached 93.75%. Therefore, the model is deemed reliable. In the future, it is still necessary to further increase the sample size and improve the accuracy of the follow-up model in geographical origin prediction.



---

*Source: 6615264-2021-06-23.xml*