Next Article in Journal
Environmentally Friendly Approach to Pectin Extraction from Grapefruit Peel: Microwave-Assisted High-Pressure CO2/H2O
Next Article in Special Issue
Exploring the Influence of Soil Types on the Mineral Profile of Honey: Implications for Geographical Origin Prediction
Previous Article in Journal
Active Low-Density Polyethylene-Based Films by Incorporating α-Tocopherol in the Free State and Loaded in PLA Nanoparticles: A Comparative Study
Previous Article in Special Issue
Traceability Research on Dendrobium devonianum Based on SWATHtoMRM
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Effects of Geographical Origin and Tree Age on the Stable Isotopes and Multi-Elements of Pu-erh Tea

1
Institute of Urban Agriculture, Chinese Academy of Agricultural Sciences, Chengdu National Agricultural Science & Technology Center, Chengdu 610213, China
2
College of Food Science, Heilongjiang Bayi Agricultural University, Daqing 163319, China
3
China Food Flavor and Nutrition Health Innovation Center, Beijing Technology and Business University, Beijing 102488, China
*
Authors to whom correspondence should be addressed.
These authors contributed equally to this work.
Foods 2024, 13(3), 473; https://doi.org/10.3390/foods13030473
Submission received: 15 December 2023 / Revised: 24 January 2024 / Accepted: 26 January 2024 / Published: 2 February 2024

Abstract

:
Pu-erh tea is a famous tea worldwide, and identification of the geographical origin of Pu-erh tea can not only protect manufacture’s interests, but also boost consumers’ confidence. However, tree age may also influence the fingerprints of Pu-erh tea. In order to study the effects of the geographical origin and tree age on the interactions of stable isotopes and multi-elements of Pu-erh tea, 53 Pu-erh tea leaves with three different age stages from three different areas in Yunnan were collected in 2023. The δ13C, δ15N values and 25 elements were determined and analyzed. The results showed that δ13C, δ15N, Mg, Mn, Fe, Cu, Zn, Rb, Sr, Y, La, Pr, Nd, Sm, Eu, Gd, Tb, Dy, Ho, Er, Tm, Yb, and Lu had significant differences among different geographical origins (p < 0.05). Mn content was significantly influenced by region and tree age interaction. Based on multi-way analysis of variance, principal component analysis and step-wised discriminant analysis, 24 parameters were found to be closely related to the geographical origin rather than tree age, and the geographical origin of Pu-erh tea can be 100.0% discriminated in cross-validation with six parameters (δ13C, δ15N, Mn, Mg, La, and Tb). The study could provide references for the establishment of a database for the traceability of Pu-erh tea, and even the identification of tea sample regions with different tree ages.

1. Introduction

Pu-erh tea, made by large-leaf tea species (Camellia sinensis var.assamica), is one of the top ten famous teas in China and is a geographical landmark product of Yunnan [1]. Pu-erh tea is made from Yunnan large-leaf sun-blue maocha. Under the protection of geographical indications, the tea is made by special processing technology, including primary processing (picking, fixing, rolling, and sun-drying) and the blending and pressing of Pu-erh raw tea, the mentation and finishing of Pu-erh ripe tea, as well as the post processing [2]. Exceptionally, plenty of nutrients such as protein, amino acids, carbohydrates, tea polyphenols, and tea pigments have been reported in Pu-erh tea [3]. Pu-erh has multiple health benefits such as anti-cancer [4], antioxidant [5], anti-hypertensive [6], and hypolipidemic properties [7,8,9]. In recent years, the demand for Pu-erh tea has increased rapidly, and consumers have increasingly higher requirements for the quality of Pu-erh tea. The government and enterprises are committed to promoting the development of Pu-erh tea production in the direction of intensification, continuity, technology and digitization.
At present, Pu-erh tea has been listed as a geographical indication product. However, due to its high quality and price, especially raw Pu-erh tea, there are frequently phenomena of substandard quality and fake origin in the market. As a result, there is an urgent need to establish a stable and reliable traceability technology of raw Pu-erh tea’s origin, and to obtain a fresh and authentic identification model by sampling from the original place of origin. Currently, the common traceability methods applied to tea, herbs and spices include chromatographic, spectroscopic and electrochemical methods such as high-performance liquid chromatography [10], gas chromatography mass spectrometry [11], near-infrared spectroscopy [12,13], electronic nose and electronic tongue [14]. In addition, mineral elements [15,16,17] and stable isotopic ratios [18] were also effective fingerprints for geographical traceability. Liu et al. detected four stable isotopes (C, N, H, and O) and 20 elements in green tea samples from different provinces of China, and the correct discrimination rate of green tea samples from different origins was 92.30% [19]. Ni et al. examined the ability to discriminate the geographic origin of green tea by determining the multi-element content and stable isotope signature of flat green tea samples from different origins, combined with the decision tree (DT) method. Under the validation of cross-validation and the “blind” dataset, the prediction accuracy was more than 70.00%, and the discrimination accuracy of green tea from different origins was 90.00% [20].
Among them, the multi-element and stable isotopic ratios were closely related with local geological background (soil [21], water [22], etc.) and environment (temperature [23], precipitation [24], etc.), which proved to be effective tools for geographical traceability. As for the geographical traceability of Pu-erh tea, Zhang et al. analyzed 41 elements in 98 Pu-erh tea samples in Yunnan Province. The results showed that the average concentrations of Fe and Pb in tea samples from Pu-erh were significantly higher than those in other production areas, Mn and Cr were generally higher in Xishuangbanna production area, while Ba and rare earth elements showed higher concentrations in Lancang Pu-erh tea samples. As a result, the geographical origins of Pu-erh tea could be distinguished based on the concentrations of 12 elements combined with chemometric analysis [25]. Li et al. conducted an exploration of the content of eight microelements in raw Pu-erh tea to assess the safety risk related to the storage year [26]. The above literature indicates that researchers mostly focus on mineral element and isotope analysis of tea leaves from different regions, and it is feasible to identify the origin of Pu-erh tea by mass-spectrometric techniques. However, the mineral element and isotopic ratios of tea leaves are not only related to the origin (soil, water, climate) [27,28], but the above parameters may also be influenced by factors such as variety and age of the tree [29,30]. However, there are few reports on the influence of tree age on the fingerprint analysis of fresh Pu ‘er leaves from different geographical origin. Therefore, we will focus on the of different origin, tree age and other factors as a new research idea in this paper.
In this study, the stable isotopic ratios and mineral elemental contents in fresh leaves of Pu-erh tea from different regions and tree ages were collected. The C and N stable isotopes and the elemental contents in tea leaves were analyzed. The indicators mainly influenced by the region and less affected by tree age were screened to establish the robust discriminative model. Based on the above results, the study could lay the theoretical research foundation for the systematic research and database construction work on the origin traceability technology of Pu-erh tea.

2. Materials and Methods

2.1. Sample Cultivation and Collection

The main production areas (Jinggu County and Ning’er County in Yunnan Pu’er tea City, and Bangdong Township in Lincang City) were selected within the limited area of Geographical Indication Product Pu-erh tea [31]. Two tea gardens were selected in each county. About 200 g of young ‘bud’ leaves (one bud and two leaf) from large-leaf tea species was sampled from tea trees at each sampling site from 26 to 29 March in 2023. Furthermore, specific information including the geographical location information (longitude, latitude, and altitude information) and tree age of Pu-erh tea trees was also recorded. The detailed information is shown in Table 1.

2.2. Sample Pretreatment

The leaves were cleaned before drying so as to remove dust and dirt. All tea samples were put into a dryer at 40 °C to obtain a constant weight. Subsequently, the dried samples were finely ground into a uniform powder using a plant crusher. The resulting powder was then passed through a 100-mesh sieve to achieve a homogeneous particle size. All samples were uniformly stored at 4 °C for further analysis.

2.3. Multi-Element Analysis

The digestion process for each sample closely followed the methodology outlined in our prior research [1]. Approximately 0.25 g of the homogenized sample was subjected to a 2-h treatment with 6 mL of concentrated HNO3 in Teflon digestion vessels. Subsequently, 2 mL of BV-III grade H2O2 was added to each vessel and allowed to react for 30 min. After the release of nitrogen oxides, the digestion vessels were introduced into a microwave digestion instrument (CEM MARS Xpress, Charlotte, NC, USA) and heated gradually to 180 °C for 40 min.
The ICP-MS operational parameters were as follows: radio frequency power at 1280 W, atomizing chamber temperature at 2 °C. The cooling gas, carrier gas, and auxiliary gas were set at flow rates of 1.47 L min−1, 1 L min−1, and 1 L min−1, respectively. To ensure accuracy, the CRM of tea flour (GBW10016) underwent digestion and determination using the same procedure. All sample determinations were performed in triplicate, with a re-measurement undertaken if the relative standard deviation of internal standard concentration exceeded 5%. Element concentration data were corrected based on dry matter after being adjusted for water content measured before digestion. The quality control (LOD, LOQ, recovery, etc.) of the instrument for the mineral element determination is shown in Table S1.

2.4. Stable Carbon and Nitrogen Isotope Analysis

Dry tea samples (0.5–0.6 mg) were carefully weighed into 6.0 mm × 4.0 mm tin capsules and introduced into an elemental analyzer (vario PYRO cube, Elementar Company, Langenselbold, Germany) equipped with an autosampler. Carbon and nitrogen elements within the samples were combusted at 1020 °C, converting them into CO2 and NOx gases. Subsequently, the NOx was reduced to N2 through a copper wire at 600 °C before entering an isotope ratio mass spectrometer (IsoPrime100, Isoprime Company, Stockport, UK) via a Conflo III dilutor.
The final stable isotope ratios are expressed as δ notation relative to international standard (Vienna Pee Dee Belimnite (VPDB) for carbon, atmospheric nitrogen (AIR) for nitrogen), according to the following equation:
δ (‰) = (Rsample/Rstandard − 1) × 1000,
where δ (‰) represents the δ13C and δ15N values, and R is the ratio of 13C/12C or 15N/14N.
For scale normalization and quality assurance, tea samples were analyzed together with reference materials including USGS40 (L-glutamic acid; δ13CVPDB = −26.389‰, δ15Nair = −4.5‰) and urea (δ13CVPDB = −43.26‰, δ15Nair = −0.56‰) for δ13C and δ15N values. Each sample was analyzed three times. The instrumental precision for stable isotope ratio measurements based on the reference materials was ≤0.2‰ for δ13C values, ≤0.2‰ for δ15N values, ≤3‰ values, respectively.

2.5. Statistical Analysis

The statistical analyses of the data, including one-way analysis of variance (one-way ANOVA), multiway analysis of variance (multiway ANOVA), principal component analysis (PCA) and linear discriminant analysis (LDA), were carried out with SPSS for Windows version 22.0 (SPSS Inc., Chicago, IL, USA).
One-way ANOVA was applied to elements to test whether the differences in average elemental values are related to considered geographical origins. With post-hoc analysis conducted using either Dunnett’s or Tukey’s test for multiple comparisons according to the result of Bartlett’s test for equal variances. Multiway ANOVA was applied was to quantify the contributions of geographical origin, tree age and their interactions (three factors) to the total variance in element levels. A factor with a larger ratio of relative variance indicates the greater influence relative to the other factors.
Principal component analysis (PCA) is used to transform a set of correlated variables into a set of uncorrelated principal components (PCs) that explain the greatest possible amount of variation in the data, and to provide a comprehensive data visualization [32]. Upon applying PCA to the analytical data for three geographical origins, tea samples could be preliminarily clustered (the first four PCs). In addition, we used Fisher’s linear discrimination analysis (LDA) to assess the effectiveness of the elements for the identification of tea origin traceability. Linear discriminant analysis (LDA) is a supervised procedure that maximizes the variances between categories and minimizes the variances within categories by creating new variables (discriminant functions). The reliability of the discriminant model was also verified by the cross-validation method (leave-one-out method).

3. Results

3.1. Comparison of Isotopic Ratios and Mineral Contents from Different Regions

The mean values and standard deviations of mineral element contents in Pu-erh tea samples from different regions are shown in Table 2. The mineral elements (Mg, Mn, Fe, Cu, Zn, Rb, Sr, Y, La, Pr, Nd, Sm, Eu, Gd, Tb, Dy, Ho, Er, Tm, Yb and Lu) and isotopes (δ15N and δ13C) had significant differences among different geographical origins (p < 0.05). The mineral elements (Mg, Mn, Fe, Cu, Rb, Sr, La, Pr, Nd, Sm, Eu, Gd, Tb, Dy, Ho, Er, Tm, Yb and Lu) and isotopes (δ15N and δ13C) had significant differences between the two geographical origins of Jinggu and Bangdong (p < 0.01). The mineral elements (Mn, Y, Gd, Tb, Dy, Ho, Er, Tm, Yb and Lu) had significant differences between the two geographical origins of Jinggu and Ning’er (p < 0.01). The mineral elements (Mg, Mn, Rb, Sr, Y, La, Pr, Nd, Sm, Eu, Gd, Tb, Dy, Ho, Er, Tm, Yb and Lu) and isotopes (δ15N and δ13C) had significant differences between the two geographical origins of Bangdong and Ning’er (p < 0.01). Specifically, the highest elemental content of K was found in Ning’er tea samples, and the Fe, Sr, Tb contents and δ15N were significantly higher in Pu-erh tea samples from Jinggu than in other regions, the elemental contents of Mg, K and Ca were higher in Pu-erh tea samples from Bangdong, while the δ13C value and the elemental contents of Mg, Mn, Y, Sm, Eu, Gd, Dy, Ho, Er, Tm, Yb, and Lu were significantly higher in Pu-erh tea samples from Ning’er than in other regions. Box plots of stable isotope ratios and mineral contents in Pu-erh tea in different regions are shown in Figure 1.
The mean values and standard deviations of mineral element contents in deep soil (30–60 cm) from different regions are shown in Table 3. The mineral elements (Mg, K, Ca, Mn, Fe, Cu, Zn, Rb, Sr, Y, La, Pr, Nd, Sm, Eu, Gd, Tb, Dy, Ho, Er, Tm, Yb and Lu) had significant differences among different geographical origins (p < 0.05). The mineral elements (Mg, K, Mn, Zn, Rb, Y, Pr, Nd, Sm, Eu, Gd, Tb, Dy, Ho and Er) had significant differences between Jinggu and Bangdong (p < 0.01). The mineral elements (Mg, K, Mn, Fe, Cu, Zn, Rb, Sr, Y, Eu, Tb, Dy, Ho, Er, Tm, Yb and Lu) had significant differences between Jinggu and Ning’er (p < 0.01). The mineral elements (Mg, Fe, Cu, Zn, Rb and Sr) had significant differences between Bangdong and Ning’er (p < 0.01). In summary, the trends of the elements (Mn, Rb, Tb and Dy) in the soil in the three geographic regions were consistent with those in Pu-erh tea.
As can be seen from Table 4, the canonical correlation analysis (CCA) extracted a total of 9 groups of typical variables, of which 7 groups of typical variables had a correlation coefficient of 0.317, and 8 groups of typical variables had a correlation coefficient of 0.111. As can be seen in Figure 2, when CCA1 was taken as a benchmark, the contents of Nd, La, Pr, Sm, Gd, Mn, Yb, Dy, Fe, Tm, and Ce were positively correlated with soil elemental content, and Eu, Cu, Rb, K, Se, Ca, Tb, and Y were negatively correlated with soil elemental content, of which the content of Ho in tea was weakly correlated with soil elemental content; when CCA2 was taken as a benchmark, the content of Zn, and Sr in tea was positively correlated with soil elemental content, Eu, Cu, Rb, K, Se, Ca, Tb, and Y were positively correlated with soil elemental content, and Eu, Cu, Rb, K, Se, Ca, Tb, and Y were positively correlated with soil elemental content. The contents of Eu, Cu, Rb, K, Se, Ca, Tb, and Y were negatively correlated with the soil element contents, and the content of Mg in tea was weakly negatively correlated with the soil element contents. At the same time, the soil characteristics of the different regions (pH, EC, etc.) are shown in Table 5. Among them, the highest pH value was found in Bangdong, which was significantly different from the other two regions (p < 0.01). The lowest EC value was found in Ning’er, which was highly significantly different from the other two regions (p < 0.01).

3.2. Comparison of Isotopic Ratios and Mineral Contents from Different Tree Ages

The mean values and standard deviations of mineral element contents in Pu-erh tea samples of different ages are shown in Table 6. In addition to Zn, other mineral elements (Mg, K, Ca, Mn, Fe, Cu, Rb, Sr, Y, La, Pr, Nd, Sm, Eu, Gd, Tb, Dy, Ho, Er, Tm, Yb, and Lu) and isotopes (δ13C and δ15N) did not show any significant differences (p > 0.05) among tree ages. Among them, the elemental content of Zn was higher in Pu-erh tea samples aged 100~200 than in other ages. Specifically, the elemental differences among tree ages for each region are shown in Table 7. In the region of Bangdong, significant differences (p < 0.05) were found for δ13C and the contents of Ca, Zn, Ho, and Yb in tea leaves among different tree ages. However, the contents of Ca, Pr, Nd, Sm, Eu and Gd in tea leaves had significant differences in the region of Ning’er (p < 0.05).

3.3. Multi-Way Analysis of Variance for Stable Isotopic Ratios and Elements

A combined analysis of variance across three regions and three tree ages was performed using the general linear model (GLM) procedure of SPSS (Table 8). Regions and tree age were considered as fixed factors, and the effects were portioned into different sources, such as region (R), age (A), and region × age (R × A). In total, the stable isotope values (δ15N and δ13C) and the mineral contents (Mg, Mn, Rb, Sr, Y La, Pr, Nd, Sm, Eu, Gd, Tb, Dy, Ho, Er, Tm, Yb, and Lu) in Pu-erh tea were highly significantly influenced by region (p < 0.01), the contents of Fe and Zn in Pu-erh tea were significantly influenced by region (p < 0.05), whereas R × A had significant effects on Mn content (p < 0.01).

3.4. Principal Component Analysis of Isotope Ratios and Mineral Content of Pu-erh Tea from Different Regions

Through the above effects of region, tree age and their interaction on the isotope ratios and mineral contents of Pu-erh tea, 22 characteristic mineral elements related to the regions, including Mg, Mn, Fe, Zn, Rb, Sr, Y, La, Pr, Nd, Sm, Eu, Gd, Tb, Dy, Ho, Er, Tm, Yb, Lu, δ13C and δ15N were screened. The screened 22 characteristic mineral elements were subjected to principal component analysis of different regions, and the results are shown in Table S2. Four principal components were obtained with a cumulative contribution of 82.54%. Principal component 1 mainly contains δ13C, δ15N, Fe, Cu, Sr, Y, La, Pr, Nd, Sm, Eu, Gd, Tb, Ho, Er, Tm, Yb and Lu, and the contribution rate is 57.76%, principal component 2 mainly contains δ15N, Cu, Sr, Y, Pr, Nd, Sm and Dy element information, and the contribution rate is 13.43%, principal component 3 mainly contains the information of δ15N and Mg, Cu, Rb, Pr, Nd, Sm, Eu, Tb and Dy elements, and the contribution rate is 6.76%, principal component 4 mainly contains δ13C, Fe, Cu, Rb, Sr and La information with a contribution rate of 4.60%. The four most important variables were Nd, Sm, Eu, Gd, Tb, Ho, Er, Tm, Yb and Lu in PC1, and δ15N, Sr, Y and Dy in PC2, and Mg in PC3, respectively (Figure 3a). Meanwhile, three origin samples (Ning’er, Jinggu, Bangdong) were correctly distinguished. And they are distributed in different spatial regions. It shows that the selected mineral elements and stable isotope origin traceability fingerprint can be used to distinguish Pu-erh tea from different origins (Figure 3b). The Pu-erh tea samples from three geographical origins were distributed in different spatial distributions. The results indicated that these elements could be used for the identification of the geographical origin of Pu-erh tea.

3.5. Discriminant Analysis of Isotope Ratio and Mineral Element Content of Pu-erh Tea from Different Regions

Stepwise discriminant analysis was conducted on Mg, Mn, Fe, Zn, Rb, Sr, Y, La, Pr, Nd, Sm, Eu, Gd, Tb, Dy, Ho, Er, Tm, Yb, Lu, δ13C and δ15N, and seven indicators were selected to establish the discriminant model. The cross-validation model was also used to obtain the classification results of Pu-erh tea samples from three regions, as shown in Table S3. Based on the discriminant model using the screened elements, the classification results showed that both original and cross-validation of correct discriminant rates could reach 100%. The result indicated that a discriminative model consisting of these elements was fully capable of correctly discriminating samples.
The discrimination model formula is as follows:
YBangdong = −44.943 δ13C + 0.081 δ15N + 0.42 Mg + 2.863 × 10−5 Mn + 8.186 × 10−6 Rb + 0.175 La + 0.101 Tb − 675.771
YNing’er = −48.217 δ13C − 3.003 δ15N + 0.054 Mg + 4.555 × 10−5 Mn − 3.021 × 10−4 Rb + 0.129 La + 0.068 Tb − 739.961
YJinggu = −45.218 δ13C − 3.404 δ15N + 0.031 Mg + 4.504 × 10−5 Mn − 1.039 × 10−4 Rb + 0.220 La + 0.034 Tb − 621.148
As shown in Figure 4, the distribution of Pu-erh tea from different regions was obtained. As was seen from the figure, the Pu-erh tea samples from the three regions were completely distinguished and located in different spaces, and there was a certain spatial range between the regions, indicating that the selected indicators related to the region were accurate and effective.
As can be seen from Table 9, the model was validated by determining the indicators (δ13C, δ15N, Mg, Mn, Rb, La, and Tb) of nine Pu-erh tea samples from Bangdong. All the samples were correctly discriminated, and a correct discrimination rate of 100.00% was obtained.

4. Discussion

In our study, there were significant differences in stable isotopes among the three regions. The lowest δ13C value was found in Pu-erh tea leaves from the Bangdong region, while no significant difference was found between the other two regions. The δ13C value in tea leaves was shown as follows: Ning’er > Jinggu > Bangdong, while the tendency of latitude for the three regions was the same, which indicated that the δ13C values in tea leaves increased with the higher latitude, which is consistent with the previous results [33]. Generally, tea areas at low latitudes have high average annual temperatures and receive more light radiation on the surface of the ground. Pu-erh tea, as a C3 plant, utilizes the Calvin photosynthetic pathway for CO2 assimilation. The δ13C values in plants have been found to correlate with the ratio of intercellular CO2 and CO2 from the surrounding environment (Pi/Pa) [34]. This correlation reflects the variability in CO2 sources (both concentrations and δ13C values) across different regions. Specifically, higher latitude is often associated with higher atmospheric CO2 concentrations, which in turn results in variations in physiological processes and δ13C values within Pu-erh tea. In addition, the δ13C values of Pu-erh tea samples with different geographical origins ranged from −24.58‰ to −27.15‰, falling in the range of C3 plants (−34‰ to −24‰) [35], and the result was agreed with a previous report in which δ13C values for Chinese green tea varied between −28.5‰ and −24.5‰ [36]. The observed variation in δ15N value was found in Pu-erh tea from the Bangdong region, while there was no significant difference in values. Numerous investigations have showcased the correlation between fertilizer application in agricultural techniques and the δ15N levels in plants. Traditionally, synthetic nitrogen-based fertilizers exhibit nitrogen isotope values ranging from −4‰ to 4‰, while organic fertilizers generally manifest higher δ15N values and exhibit a considerably broader spectrum (2–30‰) than their synthetic counterparts [37]. It is worth mentioning that distinct synthetic nitrogen fertilizers may also display varying δ15N values.
In addition, the contents of mineral elements such as Rb, La, Pr, and Nd were lower in Pu-erh tea leaves from the Bangdong region, while no significant difference was found between the other two regions. There were significant differences in the contents of mineral elements (Fe) in Pu-erh tea leaves in the Jinggu region compared with other regions. This might be because the mineral elements in plants were closely related to the mineral contents in the regional soil [38]. The mineral contents in soils were mainly influenced by conditions such as soil-forming parent material, soil pH, climate, and precipitation, which formed the specific elemental fingerprints among different regions. The characteristics of the “soil-plant” system led to variations in the composition of trace elements in plant tissues. At the same time, compared with the previous mineral elements of Pu-erh tea, especially the Mg values of Pu-erh tea samples with different geographical origins ranged from 1797.00 to 2273.50 mg/kg, agreeing with a previous report in which Mn values for Pu-erh tea varied between 381.65 to 878.15 mg/kg [39]. For Pu-erh tea from Yunnan Province, the total amount of Eu (0.010 ± 0.002 μg/g) was close to a previous report in which the content of Eu was 6.49 ± 2.30 μg/kg, and the Yb content (0.019 ± 0.006 μg/g) was also similar with the previous report, with the Yb content was 12.25 ± 4.66 μg/kg [40]. Furthermore, the mineral element content (e.g., Zn) in Pu-erh tea is also influenced by the age of the tea trees. Differences in Zn content can arise as a result of variations in the growing environment of tea trees of different ages, given that Zn constitutes one of the vital trace elements needed for plant growth and development. However, additional comprehensive investigations are required to establish the pattern of change in the content of other mineral elements across different ages of Pu-erh tea trees.
Although previous studies have obtained good discrimination for the geographical origins of tea by using the stable isotope techniques or mineral element techniques [41,42], the combination of two techniques was used for the first time to effectively identify the origin of fresh Pu-erh tea, obtaining a high overall correct classification rate (100.0%) and cross-validation rate (100.0%).

5. Conclusions

In this study, δ13C and δ15N values and 24 elemental contents of Pu-erh tea samples from different regions and tree ages were comprehensively analyzed. Based on the multi-way analysis of variance, only those significantly influenced by region were screened to better distinguish between samples of different geographical origins. Using the screened elements, step-wised discriminant analysis with 100.0% correct discrimination and 100.0% cross-validation was obtained, and a discriminant model was established with only six parameters (δ13C, δ15N, Mn, Mg, La, and Tb), indicating that the elements screened in relation to origin were accurate and effective. Therefore, it is technically feasible to screen the effective stable isotope and mineral elements to discriminate the origin of Pu-erh tea. In addition, the study could also provide a reference for the establishment of a database for the traceability of Pu-erh tea’s geographical origin. Because the tree age of 100 years is a variable, the data should be available for identification in the second year. However, it requires resampling for validation. Two tea plantations were chosen from each region, which were representative to a certain extent but still this is not many. Therefore, further expansion of the production areas under study is needed.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/foods13030473/s1, Table S1: Stable isotopes and mineral contents in Pu-erh tea from different tree ages; Table S2: Principal component analysis table of characteristic mineral elements; Table S3: Classification results of Pu-erh tea from different regions; Table S4: Stable isotope and mineral element content in Pu-erh tea samples from different regions (Tree age, variety, etc.).

Author Contributions

H.-Y.L. and L.-L.Q. conceived the idea and designed the experiments. M.-M.C. and Q.-H.L. analyzed the data and wrote the manuscript draft. H.-D.Z., Y.-L.L. and Y.S. performed the experiments. Q.-H.L., Y.X., Y.L. and H.-Y.L. edited and revised the manuscript, H.-Y.L. and Z.-L.L. provided the funding. The final version of the manuscript was approved by all authors. All authors have read and agreed to the published version of the manuscript.

Funding

The National Natural Science Foundation of China (No. 32202169), the Local Financial Funds of National Agricultural Science and Technology Center, Chengdu (No. NASC2022KR06), the Agricultural Science and Technology Innovation Program (No. ASTIP-IUA-2023003), and the Open Project of China Food Flavor and Nutrition Health Innovation Center (No. CFC2023B-027).

Institutional Review Board Statement

Not Applicable.

Informed Consent Statement

Not Applicable.

Data Availability Statement

Data is contained within the article or supplementary material.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Liu, Z.; Yuan, Y.W.; Zhang, Y.Z.; Shi, Y.Z.; Hu, G.Z.; Zhu, J.H.; Rogers, K.M. Geographical traceability of Chinese green tea using stable isotope and multi-element chemometrics. Rapid Commun. Mass Spectrom. 2019, 33, 778–788. [Google Scholar] [CrossRef]
  2. Chen, W.; Keita, Y.; Osamu, S.K.; Akiko, M.; Alexis, G.; Naohiro, Y. Development of a methodology using gas chromatography-combustion-isotope ratio mass spectrometry for the determination of the carbon isotope ratio of caffeine extracted from tea leaves (Camellia sinensis). Rapid Commun. Mass Spectrom. 2012, 26, 978–982. [Google Scholar] [CrossRef]
  3. Lin, F.J.; Wei, X.L.; Liu, H.Y.; Li, H.; Xia, Y.; Wu, D.T.; Zhang, P.Z.; Gandhi, G.R.; Li, H.B.; Gan, R.Y. State-of-the-art review of dark tea: From chemistry to health benefits. Trends Food Sci. Technol. 2021, 109, 126–138. [Google Scholar] [CrossRef]
  4. Huyen, T.V.; Fu, V.S.; Kun, V.T.; Haibin, S.; Gregory, A.C. Systematic characterization of the structure and radical scavenging potency of Pu-erh tea polyphenol theaflavin. Org. Biomol. Chem. 2019, 17, 9942–9950. [Google Scholar] [CrossRef]
  5. Su, J.J.; Wang, X.Q.; Song, W.J.; Bai, X.L.; Li, C.W. Reducing oxidative stress and hepatoprotective effect of the water extracts from Pu-erh tea on rats fed with high-fat diet. Food Sci. Hum. Wellness 2016, 5, 199–206. [Google Scholar] [CrossRef]
  6. Lv, H.P.; Zhu, Y.; Tan, J.F.; Guo, L.; Dai, W.D.; Lin, Z. Bioactive compounds from Pu-erh tea with therapy for hyperlipidaemia. J. Funct. Foods 2016, 19, 194–203. [Google Scholar] [CrossRef]
  7. Huang, F.J.; Wang, S.L.; Zhao, A.H.; Zheng, X.J.; Zhang, Y.J.; Lei, S.; Ge, K.; Qu, C.; Zhao, Q.; Yan, C.; et al. Pu-erh tea regulates fatty acid metabolism in mice under high-fat diet. Front. Pharmacol. 2019, 10, 63. [Google Scholar] [CrossRef]
  8. Luo, D.; Chen, X.J.; Zhu, X.; Liu, S.; Li, J.; Xu, J.P.; Zhao, J.H.; Ji, X. Pu-erh tea relaxes the thoracic aorta of rats by reducing intracellular calcium. Front. Pharmacol. 2019, 10, 1430. [Google Scholar] [CrossRef]
  9. Gabriella, R.; Cristina, M.; Anita, G.; Claudia, P.; Giancarlo, A.; Marina, C.; Luca, R. Ripe and raw Pu-erh tea: LCMS profiling, antioxidant capacity and enzyme inhibition activities of aqueous and Hydro-alcoholic extracts. Molecules 2004, 24, 473–475. [Google Scholar] [CrossRef]
  10. Navratilova, K.; Hrbek, V.; Kratky, F. Green tea: Authentication of geographic origin based on UHPLC-HRMS fingerprints. J. Food Compos. Anal. 2019, 78, 121–128. [Google Scholar] [CrossRef]
  11. Lin, J.; Zhang, P.; Pan, Z.Q. Discrimination of oolong tea (Camellia sinensis) varieties based on feature extraction and selection from aromatic profiles analysed by HS-SPME/GC-MS. Food Chem. 2013, 141, 259–265. [Google Scholar] [CrossRef]
  12. Kumar, V.; Roy, B.K. Population authentication of the traditional medicinal plant Cassia tora L. based on ISSR markers and FTIR analysis. Sci. Rep. 2018, 8, 10714. [Google Scholar] [CrossRef] [PubMed]
  13. Kovcs, Z.; Dalmadi, I.; Lukcs, L. Geographical origin identification of pure Sri Lanka tea infusions with electronic nose, electronic tongue and sensory profile analysis. J. Chemom. 2009, 24, 121–130. [Google Scholar] [CrossRef]
  14. Zhang, Y.H.; Jiang, L.; Feng, N.C.; Gao, Y.Y.; Zhu, M.L. Content Comparison of Seven Harmful Elements in Lycium barbarum Berry from the Geographical Origin Fields and Supermarkets. Int. J. Agric. Biol. 2020, 6, 1565–1572. [Google Scholar] [CrossRef]
  15. Pilgrim, T.S.; Watling, R.J.; Grice, K. Application of trace element and stable isotope signatures to determine the provenance of tea (Camellia sinensis) samples. Food Chem. 2010, 118, 921–926. [Google Scholar] [CrossRef]
  16. Fern, N.C.; Mart, N.M.J.; Pablos, F. Differentiation of tea (Camellia sinensis) varieties and their geographical origin according to their metal content. J. Agric. Food Chem. 2001, 49, 4775–4779. [Google Scholar] [CrossRef]
  17. Budnov, G.; Vláčil, D.; Mestek, O.; Volka, K. Application of infrared spectroscopy to the assessment of authenticity of tea. Talanta 1998, 47, 255–260. [Google Scholar] [CrossRef] [PubMed]
  18. Zhang, M.L.; Huang, C.W.; Zhang, J.Y. Accurate discrimination of tea from multiple geographical regions by combining multi-elements with multivariate statistical analysis. J. Food Meas. Charact. 2020, 14, 3361–3370. [Google Scholar] [CrossRef]
  19. Liu, W.W.; Chen, Y.; Liao, R.X.; Zhao, J.; Yang, H.; Wang, F.H. Authentication of the geographical origin of Guizhou green tea using stable isotope and mineral element signatures combined with chemometric analysis. Food Control 2021, 125, 3–8. [Google Scholar] [CrossRef]
  20. Ni, K.; Wang, J.; Zhang, Q.F.; Yi, X.Y.; Ma, L.F.; Shi, Y.Z.; Ruan, J.Y. Multi-element composition and isotopic signatures for the geographical origin discrimination of green tea in China: A case study of Xihu Longjing. J. Food Compos. Anal. 2018, 67, 104–109. [Google Scholar] [CrossRef]
  21. Rashid, M.H.; Fardous, Z.; Chowdhury, M.A.Z.; Alam, M.K.; Bari, M.L.; Moniruzzaman, M.L.; Gan, S.H. Determination of heavy metals in the soils of tea plantations and in fresh and processed tea leaves: An evaluation of six digestion methods. Chem. Cent. J. 2016, 10, 7. [Google Scholar] [CrossRef]
  22. Han, W.; Shi, Y.; Ma, L.; Ruan, J.; Zhao, F. Effect of liming and seasonal variation on lead concentration of tea plant (Camellia sinensis (L.) O. Kuntze). Chemosphere 2007, 66, 84–90. [Google Scholar] [CrossRef]
  23. Du, Y.; Shin, S.; Wang, K.; Lu, J.; Liang, Y. Effect of temperature on the expression of genes related to the accumulation of chlorophylls and carotenoids in albino tea. J. Hortic. Sci. Biotechnol. 2015, 84, 365–369. [Google Scholar] [CrossRef]
  24. Jayasekera, S.; Kaur, L.; Molan, A.; Garg, M.L.; Moughan, P. Effects of season and plantation on phenolic content of unfermented and fermented Sri Lankan tea. Food Chem. 2014, 152, 546–551. [Google Scholar] [CrossRef] [PubMed]
  25. Zhang, J.Y.; Ma, G.C.; Chen, L.Y.; Liu, T.; Liu, X.; Lu, C.Y. Profiling elements in Pu-erh tea from Yunnan province, China. Food Addit. Contam. Part B 2017, 10, 155–164. [Google Scholar] [CrossRef]
  26. Li, F.; Lu, Q.H.; Li, M.; Yang, X.M.; Xiong, C.Y.; Yang, B. Comparison and risk assessment for trace heavy metals in raw Pu-erh tea with different storage years. Biol. Trace Elem. Res. 2020, 195, 696–706. [Google Scholar] [CrossRef]
  27. Zhao, H.Y.; Yang, Q.L. The suitability of rare earth elements for geographical traceability of tea leaves. J. Sci. Food Agric. 2019, 99, 6509–6514. [Google Scholar] [CrossRef]
  28. Zhao, H.Y.; Yu, C.D.; Li, M. Effects of geographical origin, variety, season and their interactions on minerals in tea for traceability. J. Food Compos. Anal. 2017, 63, 15–20. [Google Scholar] [CrossRef]
  29. Ding, B.; Zeng, G.F.; Wang, Z.Y.; Xie, J.J.; Wang, L.; Chen, W.R. Authenticity determination of tea drinks in the Chinese market by liquid chromatography coupled to isotope ratio mass spectrometry. Microchem. J. 2019, 144, 139–143. [Google Scholar] [CrossRef]
  30. Liu, H.Y.; Wei, Y.M.; Zhang, Y.Q.; Wei, S.; Zhang, S.S.; Guo, B.L. The effectiveness of multi-element fingerprints for identifying the geographical origin of wheat. Int. J. Food Sci. Technol. 2017, 52, 1018–1025. [Google Scholar] [CrossRef]
  31. GB/T 22111-2008; Geographical Indication Products Pu-erh Tea. Available online: https://www.chinesestandard.net/PDF/English.aspx/GBT22111-2008 (accessed on 8 August 2023).
  32. Wold, S.; Esbensen, K.; Geladi, P. Principal component analysis. Chemom. Intell. Lab. Syst. 1987, 2, 37–52. [Google Scholar] [CrossRef]
  33. Wang, J.; Li, X.; Wu, Y.; Qu, F.; Liu, L.; Wang, B.; Wang, P.; Zhang, X. HS−SPME/GC−MS Reveals the Season Effects on Volatile Compounds of Green Tea in High−Latitude Region. Foods 2022, 11, 3016. [Google Scholar] [CrossRef] [PubMed]
  34. Farquhar, G.D.; O’Leary, M.H.; Berry, J.A. On the Relationship Between Carbon Isotope Discrimination and the intercellular carbon dioxide concentration in leaves. Aust. J. Plant Physiol. 1982, 9, 121–137. [Google Scholar] [CrossRef]
  35. Vogel, J.C. Fractionation of the Carbon Isotopes during Photosynthesis; Springer: Berlin/Heidelberg, Germany, 1980; pp. 111–135. [Google Scholar] [CrossRef]
  36. Deng, X.F.; Liu, Z.; Zhan, Y.; Ni, K.; Zhang, Y.Z.; Ma, W.Z.; Shao, S.Z.; Lv, X.N.; Yuan, Y.W.; Rogers, K.M. Predictive geographical authentication of green tea with protected designation of origin using a random forest model. Food Control 2020, 107, 106807. [Google Scholar] [CrossRef]
  37. Kendall, C. Tracing nitrogen sources and cycling in catchments. In Isotope Tracers in Catchment Hydrology; Kendall, C., McDonnell, J.J., Eds.; Elsevier: Amsterdam, The Netherlands, 1998; pp. 519–576. [Google Scholar] [CrossRef]
  38. Zhang, J.Y.; Nie, J.Y.; Zhang, L.B.; Xu, G.F.; Zheng, H.D.; Shen, Y.M.; Kuang, L.X.; Gao, X.Q.; Zhang, H. Multielement authentication of apples from the cold highlands in southwest China. J. Sci. Food Agric. 2022, 102, 241–249. [Google Scholar] [CrossRef]
  39. Lv, H.P.; Lin, Z.; Zhang, Y.; Liang, Y.R. Study on the content of the major mineral elements in Pu-erh tea. J. Tea Sci. 2013, 33, 411–419. [Google Scholar] [CrossRef]
  40. Zhu, J.Y.; Chen, L.; Chen, Y.; Rong, Y.T.; Jiang, Y.W.; Liu, F.Q.; Zhou, Q.H.; Wei, X.H.; Yuan, H.B.; Zhang, J.J.; et al. Effect of geographical origins and pile-fermentation on the multi-element profiles of ripen Pu-erh tea revealed by comprehensive elemental fingerprinting. Food Control 2013, 105, 109978. [Google Scholar] [CrossRef]
  41. Lagad, R.A.; Alamelu, D.; Laskar, A.H.; Rai, V.K.; Singh, S.K.; Aggarwal, S.K. Isotope signature study of the tea samples produced at four different regions in India. Anal. Methods 2013, 5, 1604. [Google Scholar] [CrossRef]
  42. Lou, Y.X.; Fu, X.S.; Yu, X.P.; Ye, Z.H.; Cui, H.F.; Zhang, Y.F. Stable isotope ratio and elemental profile combined with support vector machine for provenance discrimination of Oolong tea (Wuyi-Rock Tea). J. Anal. Methods Chem. 2017, 2017, 5454231. [Google Scholar] [CrossRef]
Figure 1. Box plots of stable isotope ratio and mineral content in Pu-erh tea in different regions. a–c in the same row indicated that there are significant differences among regions at p < 0.05 level. Note: ◦ indicates outliers, and the number in the upper right corner represents the number of in-dividuals that are outliers; * denotes outliers, the number in the upper right corner represents the number of individuals that are outliers.
Figure 1. Box plots of stable isotope ratio and mineral content in Pu-erh tea in different regions. a–c in the same row indicated that there are significant differences among regions at p < 0.05 level. Note: ◦ indicates outliers, and the number in the upper right corner represents the number of in-dividuals that are outliers; * denotes outliers, the number in the upper right corner represents the number of individuals that are outliers.
Foods 13 00473 g001
Figure 2. Canonical correlation analysis plot.
Figure 2. Canonical correlation analysis plot.
Foods 13 00473 g002
Figure 3. Discriminant function scores of Pu-erh tea of different regions. (a) The Pu-erh tea samples from different origin sources were plotted using principal component scores; (b) The scores of Pu-erh tea from different origin sources were plotted using the scores of PC1, PC2, and PC3.
Figure 3. Discriminant function scores of Pu-erh tea of different regions. (a) The Pu-erh tea samples from different origin sources were plotted using principal component scores; (b) The scores of Pu-erh tea from different origin sources were plotted using the scores of PC1, PC2, and PC3.
Foods 13 00473 g003
Figure 4. Discriminant function score map of Pu-erh tea from different regions.
Figure 4. Discriminant function score map of Pu-erh tea from different regions.
Foods 13 00473 g004
Table 1. The sample numbers, location, and tree age of tea geographical origins.
Table 1. The sample numbers, location, and tree age of tea geographical origins.
RegionNumber of SamplesN Latitude (deg)E Longitude (deg)Altitude (m)Tree Age (Year)
Jinggu2423.7227100.68771842–190120~100 (12),
100~200 (6),
>200 (6)
Bangdong1423.9374–23.9397100.3532–100.35621633–173920~100 (8),
100~200 (3),
>200 (3)
Ning’er1523.2548101.0822161420~100 (5),
100~200 (5),
>200 (5)
Table 2. Stable isotope ratios and mineral contents in Pu-erh tea from different regions.
Table 2. Stable isotope ratios and mineral contents in Pu-erh tea from different regions.
ElementJingguBangdongNing’er
δ13C (‰) *−25.16 ± 0.83 a−26.77 ± 0.61 b−24.82 ± 1.16 a
δ15N (‰) **4.06 ± 3.35 a−0.11 ± 0.66 c2.39 ± 1.30 b
Mg (mg/kg) **1776.176 ± 184.688 b2077.253 ± 152.236 c1601.886 ± 201.680 a
K (mg/kg)19,617.81 ± 1591.58 a19,701.25 ± 1035.70 a20,368.84 ± 1833.42 a
Ca (mg/kg)4193.70 ± 820.56 a4184.00 ± 486.66 a3864.72 ± 469.04 a
Mn (mg/kg) **376 ± 139 c790 ± 219 b1159 ± 359 a
Fe (mg/kg) *155 ± 164 a70 ± 22 b84 ± 20 b
Cu (mg/kg) *12.7 ± 2.5 b15.2 ± 1.1 a14.3 ± 2.5 a
Zn (mg/kg) *34.7 ± 5.9 ab37.3 ± 4.8 a32.9 ± 4.2 b
Rb (mg/kg) *75 ± 21 a28 ± 6 b65 ± 17 a
Sr (mg/kg) **14 ± 5 a7 ± 2 c11 ± 4 b
Y (μg/kg) *876.96 ± 79.77 b843.16 ± 120.60 b1020.33 ± 125.51 a
La (μg/kg) *119.86 ± 38.80 a25.08 ± 4.42 b123.80 ± 68.58 a
Pr (μg/kg) *24 ± 8 a6 ± 2 b29 ± 14 a
Nd (μg/kg) *88.38 ± 30.35 a24.53 ± 6.16 b115.51 ± 54.43 a
Sm (μg/kg) **18 ± 6 b6 ± 2 c25 ± 12 a
Eu (μg/kg) **5.0 ± 1.2 b2.7 ± 0.8 c6.5 ± 2.3 a
Gd (μg/kg) **16 ± 5 b7 ± 3 c26 ± 11 a
Tb (μg/kg) **1122.8 ± 238.5 a630.6 ± 124.7 b219.6 ± 110.1 c
Dy (μg/kg) **13 ± 4 b6 ± 2 c21 ± 9 a
Ho (μg/kg) **2.7 ± 0.7 b1.2 ± 0.5 c4.3 ± 1.6 a
Er (μg/kg) **9 ± 3 b3 ± 1 c13 ± 5 a
Tm (μg/kg) **1.2 ± 0.5 b0.4 ± 0.2 c1.9 ± 0.8 a
Yb (μg/kg) **8 ± 3 b3 ± 1 c13 ± 5 a
Lu (μg/kg) **3.6 ± 1.8 c8.5 ± 2.9 b11.8 ± 3.9 a
Data are shown as the mean ± standard deviation. a–c in the same row indicated that there are significant differences among regions at p < 0.05 level. * means significant difference (p < 0.05), ** means highly significant difference (p < 0.01).
Table 3. Mineral contents in deep soil (30–60 cm) from different regions.
Table 3. Mineral contents in deep soil (30–60 cm) from different regions.
ElementJingguBangdongNing’er
Mg (μg/kg) *4973.044 ± 1006.566 a2381.777 ± 168.522 b1535.505 ± 76.886 b
K (mg/kg) *27.28 ± 2.41 a10.85 ± 0.08 11.55 ± 0.27 b
Ca (μg/kg) *66.26 ± 22.42 a51.38 ± 5.59 ab37.26 ± 1.24 b
Mn (μg/kg) **1264 ± 46 b1470 ± 40 a341 ± 26 c
Fe (mg/kg) *38 ± 5 a32 ± 1 a17 ± 1 b
Cu (μg/kg) *57.2 ± 7.8 a49.3 ± 1.6 a18.4 ± 0.3 b
Zn (μg/kg) **156.1 ± 14.7 a96.9 ± 3.1 b42.9 ± 2.1 c
Rb (μg/kg) **134 ± 13 a49 ± 4 c82.36 ± 2.04 b
Sr (μg/kg) **28 ± 4 b36 ± 0 a20 ± 1 c
Y (μg/kg) *18.18 ± 1.35 a11.66 ± 0.85 b9.68 ± 1.80 b
La (μg/kg) *40.03 ± 2.05 a27.05 ± 1.64 b31.82 ± 7.64 ab
Pr (μg/kg) *9 ± 0 a5 ± 0 b6 ± 2 b
Nd (μg/kg) *34.81 ± 1.35 a18.09 ± 1.27 b20.85 ± 7.20 b
Sm (μg/kg) *7 ± 0 a3 ± 0 b4 ± 2 b
Eu (μg/kg) **0.8 ± 0.1 a0.5 ± 0.0 b0.3 ± 0.0 c
Gd (μg/kg) *6 ± 0 a2 ± 0 b3 ± 1 b
Tb (μg/kg) *0.7 ± 0.1 a0.3 ± 0.0 b0.4 ± 0.1 b
Dy (μg/kg) *4 ± 0 a2 ± 0 b2 ± 1 b
Ho (μg/kg) *0.7 ± 0.1 a0.4 ± 0.0 b0.3 ± 0.1 b
Er (μg/kg) *2 ± 0 a1 ± 0 b1 ± 0 b
Tm (μg/kg) **0.2 ± 0.0 a0.2 ± 0.0 b0.1 ± 0.0 c
Yb (μg/kg) **2 ± 0 a1 ± 0 b1 ± 0 c
Lu (μg/kg) **0.2 ± 0.0 a0.1 ± 0.0 b0.1 ± 0.0 c
Data are shown as the mean ± standard deviation. a–c in the same row indicated that there are significant differences among regions at p < 0.05 level. * means significant difference (p < 0.05), ** means highly significant difference (p < 0.01).
Table 4. The canonical correlation coefficients.
Table 4. The canonical correlation coefficients.
Canonical VariableCorrelation Coefficient
1−0.077
2−0.150
3−0.069
4−0.255
5−0.246
6−0.215
70.317
80.111
90.371
Table 5. The soil characteristics of the different regions (pH, EC, etc.).
Table 5. The soil characteristics of the different regions (pH, EC, etc.).
RegionBangdongJingguNing’er
pH5.48 ± 0.15 a5.07 ± 0.17 b5.03 ± 0.06 b
EC (μs/cm)46.98 ± 0.93 a34.99 ± 5.97 b24.21 ± 5.64 c
a–c in the same row indicated that there are significant differences among regions at p < 0.05 level.
Table 6. Stable isotopes and mineral contents in Pu-erh tea from different tree ages.
Table 6. Stable isotopes and mineral contents in Pu-erh tea from different tree ages.
Element20~100100~200>200
δ13C (‰)−25.19 ± 1.27 a−25.56 ± 1.07 a−25.86 ± 1.37 a
δ15N (‰)1.82 ± 2.30 a2.43 ± 2.93 a2.36 ± 2.43 a
Mg (mg/kg)1758.798 ± 257.603 a1799.159 ± 335.310 a1808.026 ± 233.497 a
K (mg/kg)20,348.89 ± 1777.85 a19,778.20 ± 1233.41 a19,528.80 ± 1504.73 a
Ca (mg/kg)4109.37 ± 567.86 a3826.25 ± 757.32 a4137.37 ± 431.45 a
Mn (μg/kg)821 ± 481 a889 ± 363 a855 ± 400 a
Fe (mg/kg)1171 ± 124 a86 ± 50 a80 ± 20 a
Cu (mg/kg)13.9 ± 2.5 a15.2 ± 2.5 a13.4 ± 1.7 a
Zn (mg/kg) *34.4 ± 4.7 ab37.1 ± 5.9 a32.4 ± 4.0 b
Rb (mg/kg)62 ± 25 a51 ± 23 a55 ± 25 a
Sr (mg/kg)10 ± 4 a112 ± 6 a12 ± 5 a
Y (μg/kg)935.15 ± 137.55 a942.68 ± 126.28 a916.89 ± 159.38 a
La (μg/kg)109.72 ± 68.93 a68.80 ± 46.51 a94.24 ± 74.79 a
Pr (μg/kg)25 ± 16 a16 ± 9 a20 ± 14 a
Nd (μg/kg)95.70 ± 61.48 a61.80 ± 35.75 a79.98 ± 55.41 a
Sm (μg/kg)21 ± 13 a13 ± 7 a17 ± 11 a
Eu (μg/kg)5.4 ± 2.6 a4.3 ± 1.8 a5.1 ± 2.3 a
Gd (μg/kg)20 ± 13 a14 ± 7 a18 ± 11 a
Tb (μg/kg)638.3 ± 478.3 a520.7 ± 306.6 a560.2 ± 351.2 a
Dy (μg/kg)16 ± 10 a12 ± 6 a146 ± 9 a
Ho (μg/kg)3.3 ± 2.0 a2.5 ± 1.2 a3.0 ± 1.8 a
Er (μg/kg)10 ± 6 a8 ± 4 a9 ± 5 a
Tm (μg/kg)1.4 ± 1.0 a1.1 ± 0.7 a1.3 ± 0.9 a
Yb (μg/kg)9 ± 6 a7 ± 4 a9 ± 5 a
Lu (μg/kg)9.9 ± 5.4 a7.9 ± 3.3 b9.2 ± 4.0 c
Data are shown as the mean ± standard deviation. a–c in the same row indicated that there are significant differences among tree ages at p < 0.05 level. * means significant difference (p < 0.05).
Table 7. Stable isotopes and mineral content of Pu-erh tea of different tree ages in each origin.
Table 7. Stable isotopes and mineral content of Pu-erh tea of different tree ages in each origin.
ElementJingguBangdongNing’er
20~100100~200>20020~100100~200>20020~100100~200>200
δ13C (‰) *−24.99 ± 0.90 a−25.46 ± 0.33 a−25.33 ± 1.11 a−26.97 ± 0.34 b−26.19 ± 0.53 a−27.15 ± 0.50 b−24.58 ± 1.04 a−25.09 ± 1.43 a−25.04 ± 1.22 a
δ15N (‰)2.53 ± 3.20 a6.82 ± 3.29 a5.37 ± 1.67 a−0.53 ± 0.50 a0.21 ± 0.88 a−0.01 ± 0.38 a2.33 ± 1.33 a2.08 ± 0.88 a2.83 ± 1.65 a
Mg (mg/kg)1801.02 ± 210.65 a1731.73 ± 184.49 a1754.36 ± 161.95 a2022.54 ± 196.35 a2146.70 ± 150.99 a2062.52 ± 98.59 a1620.76 ± 220.96 a1543.25 ± 241.60 a1622.79 ± 125.77 a
K(mg/kg)19,745.13 ± 1778.10 a19,765.25 ± 1961.40 a19,130.87 ± 1048.56 a19,640.30 ± 794.35 a20,228.30 ± 1420.24 a19,235.15 ± 694.88 a21,046.64 ± 1905.34 a19,409.59 ± 661.17 a19,972.49 ± 2154.85 a
Ca (mg/kg) *4573.75 ± 567.83 a3250.35 ± 1008.22 b4123.57 ± 558.13 a3886.28 ± 424.46 b4532.52 ± 350.73 a4133.20 ± 509.79 ab3892.74 ± 449.67 ab3525.65 ± 424.20 b4147.75 ± 388.87 a
Mn (mg/kg)345 ± 147,577 a390 ± 162,990 a443 ± 106,601 a2 ± 196 a2 ± 151 a2 ± 99 a1207 ± 362 a1130 ± 299 a1093 ± 451 a
Fe (mg/kg)182 ± 210 a138 ± 103 a100 ± 32 a76 ± 39 a66 ± 7 a68 ± 4 a90 ± 25 a76 ± 10 a80 ± 15 a
Cu (mg/kg)12.1 ± 2.6 a14.5 ± 3.1 a12.5 ± 0.9 a15.6 ± 1.1 a15.6 ± 0.8 a14.4 ± 0.9 a14.5 ± 2.3 a15.1 ± 3.3 a13.1 ± 2.1 a
Zn (mg/kg) *35.5 ± 4.8 a34.9 ± 11.0 a32.4 ± 3.7 a37.1 ± 3.6 ab41.7 ± 2.2 a33.0 ± 4.1 b32.6 ± 4.7 a34.5 ± 2.5 a31.9 ± 4.6 a
Rb (mg/kg)77 ± 21 a62 ± 24 a82 ± 18 a26 ± 5 a27 ± 8 a32 ± 5 a67 ± 17 a65.477 ± 15 a60 ± 21 a
Sr (mg/kg)12 ± 3 a20 ± 3 a15 ± 2 a7 ± 2 a6 ± 2 a8 ± 2 a11 ± 4 a10 ± 3 a13 ± 6 a
Y (μg/kg)882.84 ± 89.57 a871.83 ± 76.97 a866.40 ± 83.77 a820.49 ± 126.44 a891.47 ± 98.12 a817.53 ± 144.60 a1017.80 ± 122.24 a1020.78 ± 135.79 a1024.94 ± 145.27 a
La (μg/kg)104.81 ± 28.51 a118.79 ± 36.43 a161.04 ± 46.52 a22.65 ± 3.31 a26.91 ± 6.35 a25.70 ± 2.23 a149.27 ± 68.38 a78.72 ± 41.09 a117.96 ± 75.00 a
Pr (μg/kg) *21 ± 7 a27 ± 7 a29 ± 8 a6 ± 2 a7 ± 1 a6 ± 1 a36 ± 14 a18 ± 7 b28 ± 14 ab
Nd (μg/kg) *77.34 ± 30.11 a99.98 ± 26.08 a106.19 ± 31.11 a20.61 ± 4.17 a26.65 ± 7.40 a26.33 ± 5.63 a139.23 ± 52.80 a71.99 ± 27.79 b111.59 ± 55.50 ab
Sm (μg/kg) *16 ± 6 a20 ± 6 a21 ± 4 a5 ± 1 a7 ± 2 a7 ± 2 a30 ± 11 a16 ± 6 b24 ± 12 ab
Eu (μg/kg) *4.5 ± 0.8 a5.7 ± 2.0 a5.6 ± 0.8 a2.1 ± 0.4 a2.9 ± 0.6 a2.9 ± 1.0 a7.4 ± 2.2 a4.6 ± 1.7 b6.6 ± 2.3 ab
Gd (μg/kg) *14 ± 4 a8 ± 6 a17 ± 6 a5 ± 1 a8 ± 3 a8 ± 4 a31 ± 10 a18 ± 7 b26 ± 12 ab
Tb (μg/kg) *1252.0 ± 207.2 a864.3 ± 64.2 b1036.8 ± 192.7 a654.1 ± 115.1 a652.9 ± 127.3 a584.9 ± 145.2 a222.6 ± 66.5 a182.2 ± 100.2 a249.6 ± 196.4 a
Dy (μg/kg)12 ± 3 a14 ± 5 a14 ± 3 a4 ± 1 a6 ± 1 a7 ± 3 a24 ± 9 a15 ± 5 a21 ± 10 a
Ho (μg/kg) *2.5 ± 0.6 a3.1 ± 0.7 a3.1 ± 0.8 a0.9 ± 0.2 b1.2 ± 0.2 ab1.5 ± 0.7 a4.8 ± 1.6 a3.2 ± 1.1 a4.3 ± 1.9 a
Er (μg/kg)8 ± 2 a1 ± 4 a9 ± 3 a3 ± 1 a3 ± 1 a4 ± 2 a14 ± 5 a10 ± 3 a13 ± 5 a
Tm (μg/kg)1.0 ± 0.4 a1.7 ± 0.7 a1.3 ± 0.5 a0.3 ± 0.1 a0.4 ± 0.1 a0.5 ± 0.3 a2.2 ± 0.9 a1.4 ± 0.5 a2.0 ± 0.8 a
Yb (μg/kg) *7 ± 3 a10 ± 4 a10 ± 3 a2 ± 1 b3 ± 1 ab4 ± 2 a14 ± 4 a9 ± 3 a12 ± 5 ab
Lu (μg/kg)3.7 ± 2.7 a3.5 ± 1.1 a3.6 ± 1.7 a7.6 ± 2.2 a8.2 ± 3.4 a9.6 ± 3.1 a12.9 ± 4.8 a9.8 ± 1.7 a11.8 ± 2.6 a
Data are shown as the mean ± standard deviation. a,b in the same row indicated that there are significant differences among regions at p < 0.05 level. * means significant difference (p < 0.05).
Table 8. Mean square of each stable isotope and element by analysis of variance.
Table 8. Mean square of each stable isotope and element by analysis of variance.
Source of VariationRegion (R)Age (A)R × AError
DF22444
cc (‰)16.162 **0.4530.7830.918
δ15N (‰)89.32 **4.983.313.24
Mg (mg/kg)940,372.990 **214.53216,805.58135,097.374
K (mg/kg)320,175.524,407,718.361,902,670.372,576,350.83
Ca (mg/kg)343,041.76441,090.52973,706.97223,680.31
Mn (mg/kg)212,6241 **45,78238,559,311 **88,978,814
Fe (mg/kg)28,113 *12,67971248161
Cu (mg/kg)12,380.39732.21244.75492.3
Zn (mg/kg)77,533.4 *48,583.222,130.521,778.1
Rb (mg/kg)7,570,777 **208,510230,305211,491
Sr (mg/kg)231,854 **32,51634,62413,801
Y (mg/kg)178.84 **4.562.2515.11
La (mg/kg)51.79 **2.902.402.32
Pr (μg/kg)2357 **8215295
Nd (mg/kg)35.45 **1.232.201.38
Sm (μg/kg)1466 **569561
Eu (μg/kg)60.9 **1.34.92.7
Gd (μg/kg)1439 **289757
Tb (mg/kg)2068.2 **0.00.00.0
Dy (μg/kg)954 **114839
Ho (μg/kg)39.5 **0.71.61.4
Er (μg/kg)363 **21813
Tm (μg/kg)9.5 **0.00.50.4
Yb (μg/kg)364 **52213
Lu (μg/kg)215.4 **3.85.511.3
* means significant effect (p < 0.05), ** means highly significant effect (p < 0.01).
Table 9. The external validation based on the indicator of new Pu-erh tea samples.
Table 9. The external validation based on the indicator of new Pu-erh tea samples.
Regionδ13Cδ15NMgMnRbLaTb
Bangdong−26.011.402,724,050.971,235,282.2590,404.89523.638.20
−25.982.742,772,521.471,066,798.5186,209.27289.256.09
−27.433.192,632,641.731,043,034.7896,223.16214.164.73
−27.912.752,686,877.79916,569.28102,914.99367.906.83
−26.162.573,067,714.29969,110.57133,737.01375.745.66
−25.362.502,658,854.62737,281.79108,884.34242.075.16
−25.682.312,482,755.051,006,464.11102,004.56391.337.49
−25.321.932,760,370.631,260,159.98120,953.10295.247.58
−26.302.572,543,127.79947,274.9995,336.29369.316.56
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Chen, M.-M.; Liao, Q.-H.; Qian, L.-L.; Zou, H.-D.; Li, Y.-L.; Song, Y.; Xia, Y.; Liu, Y.; Liu, H.-Y.; Liu, Z.-L. Effects of Geographical Origin and Tree Age on the Stable Isotopes and Multi-Elements of Pu-erh Tea. Foods 2024, 13, 473. https://doi.org/10.3390/foods13030473

AMA Style

Chen M-M, Liao Q-H, Qian L-L, Zou H-D, Li Y-L, Song Y, Xia Y, Liu Y, Liu H-Y, Liu Z-L. Effects of Geographical Origin and Tree Age on the Stable Isotopes and Multi-Elements of Pu-erh Tea. Foods. 2024; 13(3):473. https://doi.org/10.3390/foods13030473

Chicago/Turabian Style

Chen, Ming-Ming, Qiu-Hong Liao, Li-Li Qian, Hai-Dan Zou, Yan-Long Li, Yan Song, Yu Xia, Yi Liu, Hong-Yan Liu, and Ze-Long Liu. 2024. "Effects of Geographical Origin and Tree Age on the Stable Isotopes and Multi-Elements of Pu-erh Tea" Foods 13, no. 3: 473. https://doi.org/10.3390/foods13030473

APA Style

Chen, M. -M., Liao, Q. -H., Qian, L. -L., Zou, H. -D., Li, Y. -L., Song, Y., Xia, Y., Liu, Y., Liu, H. -Y., & Liu, Z. -L. (2024). Effects of Geographical Origin and Tree Age on the Stable Isotopes and Multi-Elements of Pu-erh Tea. Foods, 13(3), 473. https://doi.org/10.3390/foods13030473

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

Article Metrics

Back to TopTop