Geographical Origin Identification of Tequila Based on Multielement and Stable Isotopes

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 δ 13 C and δ 18 O 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 veriﬁcation of 16 samples, the accuracy of prediction reached 93.75%. It has been demonstrated that it is possible to distinguish tequila of diﬀerent 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.


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. erefore, 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][2][3][4], aroma fingerprint method [5], and stable isotope mass spectrometry (IRMS) technology [6][7][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][15][16]. e 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 (δ 13 C VPDB ), oxygen 18 isotopic ratio (δ 18 O VSMOW ), 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.
ey 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.

Materials.
e 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.

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]. irdly, 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% CO 2 and He for 300 s to remove the air present in the bottle. en, it was equilibrated for 18-24 hours. Finally, the gaseous CO 2 and oxygen present in the water were exchanged for an isotope, to facilitate the analysis of the stable oxygen isotope in the water. [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.

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. e heating program is detailed as follows. e initial temperature is 40°C, which is kept for 5 min. en, 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 CO 2 . Reference gas pressure is 1.3 bar. Column temperature is 65°C. ere 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.

Data Processing.
ermo 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 using R language. Origin 9.0 (Origin lab, USA) was adopted to carry out data mapping.

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 δ 13 C value of quality control standards and samples is required not to exceed 0.5% and 0.3%, respectively. e carbon isotope and oxygen isotope calibration of the reference gas CO 2 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 δ 13 C value of the reference gas CO 2 is 0.016%, and the standard deviation of the oxygen isotope ratio δ 18 O is less than 0.06%, which meets the testing requirements of the experimental requirements. e 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. en, the standard deviation was calculated, and the results were obtained, as listed in Table 2. e standard deviations shown by the δ 13 C value of ethanol and the δ 18 O value of water in tequila are 0.12% and 0.03%, respectively, which meets the accuracy requirements of the experiment.

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. is 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 Table 3 [23]. Besides the carbon and oxygen isotope ratio shown Table 4, the δ 13 C 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. erefore, these can be used to assist geographical origin traceability analysis.

Model Building.
ere are significant differences exhibited by the results obtained from the single-factor analysis of variance for the stable carbon isotope δ 13 C and oxygen isotope δ 18 O 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]. e biggest contributor is made up of 6 elements (Sr, Pb, Se, V, Ca, and Fe) with 91.59% contributions in this PCA model. erefore, 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. e results are shown in Figures 1 and 2, respectively.
As shown in Figures 1 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. ese can be differentiated more effectively. Table 6, one sample obtained from Guanajuato is misjudged as Tamaulipas, with the prediction accuracy reaching 75%. e other three production areas are all distinguished correctly, with the prediction accuracy reaching 100%. erefore, 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.

Conclusion
After the data was measured by the stable carbon isotope δ 13 C of ethanol and the stable oxygen isotope δ 18 O 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%. erefore, 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.
Data Availability e precision and accuracy of ICP-MS in determining 22 elements in tequila wine through standard recovery experiments are from ref. [18].

Conflicts of Interest
e authors declare that they have no conflicts of interest.