# Multielement Analysis of Pakchoi (Brassica rapa L. ssp. chinensis) by ICP-MS and Their Classification according to Different Small Geographical Origins

Authors: Trung Nguyen-Quang; Giang Do-Hoang; Minh Truong-Ngoc
Journal: Journal of Analytical Methods in Chemistry
 (2021)

---

## Abstract

Statistical interpretation of the concentrations of 42 elements, determined using solution-based inductively coupled plasma mass spectrometry (ICP-MS) analysis and multivariate statistical methods, such as principal components analysis (PCA), was used to establish the provenance of pakchoi (Brassica rapa L. ssp. chinensis) from 6 areas in Ha Noi, Vietnam. Although pakchoi is widely cultivated and manufactured, no universal method is used to discriminate the origin of this vegetable. Our study introduced for the first time a method to classify pakchoi in small geographical areas. 42 metallic elements of pakchoi were detected by ICP-MS, which were further analyzed using multivariate statistical analysis to perform clusters based on the geographical locations. Eleven elements, i.e., 28Si; 56Fe; 59Co; 63Cu; 69Ga; 75As; 85Rb; 93Nb; 107Ag; 118Sn, and 137Ba, were identified as discriminators to distinguish pakchoi from those areas. Results from this study presented a new method to discriminant the geographical origins of pakchoi, which could apply to other types of vegetables on the food market.


---

## Body

## 1. Introduction

Discriminating the geographical origins of various food and drink products is growing considerably in recent years. In general, inductively coupled plasma mass spectrometry (ICP-MS) is one of the most useful approaches for multielement analysis of foods and drinks. It has been shown that ICP-MS data could be used successfully for discriminating wines from different geographical areas in Slovenia, Romania, and Croatia [1–3]. Chuzinska and Baralkiewicz determined 13 mineral elements for the classification of 55 honey samples in Poland [4]. Whereas, Chung et al. distinguished verified cultivation regions of potato by the stable isotope composition analysis of bioelements [5]. The geographical origin of vegetables, such as cabbages, was also discriminated by 21 inorganic elements through ICP-MS analysis [6]. The data of elements compositions are often combined with multivariate statistical analysis, such as principal component analysis (PCA), to reduce the dimension of the data matrix, visualize the similarities and differences among samples, and identify the key discriminants distinguishing the samples [7].Pakchoi (Brassica rapa L.) is a major vegetable in Vietnamese diets. The consumption of pakchoi may provide many health benefits. Previous studies have shown that pakchoi potentially exerts inhibitory activity against cancer [8]. Phytochemical investigations have shown that flavonoids, hydroxycinnamic acids, and glucosinolate are the major bioactive ingredients of pakchoi, which possess high antioxidant activity [9, 10].The chemical compositions of pakchoi are varied with cultivating conditions, such as soils, weather, and water resources [11, 12]. The growth and quality of pakchoi are significantly affected by inorganic elements, for example, selenate at low concentrations (0.5–1.0 mg/kg) could promote the growth of pakchoi and reduce cadmium contents in this vegetable [13]. The use of organic fertilizers and soil improvers, such as leonardite, might enrich the contents of macronutrients (Mg, Ca, K, and S) and micronutrients (Fe, Cu, Mn, and Zn).The classification of different pakchoi cultivars has been conducted by different methods. Wiesner reported that pakchoi origin can be classified through their glucosinolate profile using HPLC-DA-MS analysis. In the glucosinolate profile, levels of 3-butenyl glucosinolate and 2-hydroxy-3-butenyl glucosinolate were considered important discriminators. By contrast, 5-methylsulfinylpentyl glucosinolate was found to be a characteristic ontogenetic variation of mature leaves of pakchoi by using PCA as the statistical model [14]. Another research also assessed the morphological and genetic diversity of pakchoi based on PCA and clustering analysis [15]. On the other hand, no previous approach categorized the geographical origins based on the characterization of elemental compositions of pakchoi. Thus, this study presents a promising method to discriminate against the geographical origins of pakchoi by using ICP-MS elemental analysis coupled with statistical principal component analysis. The proposed method may be utilized for discriminating against the origins of other vegetables.


## 2. Materials and Methods

### 2.1. Samples

Totally 60 pakchoi (Brassica rapa L. ssp. chinensis) samples were obtained from 6 different locations in Hanoi, Vietnam (Figure 1), which are Linh Nam (Hoang Mai district) 20°58′33.7″N 105°53′21.9″E, Tien Phong (Me Linh district) 21°09′27.2″N 105°44′19.4″E, Thanh Da (Phuc Tho district) 21°06′32.6″N 105°36′43.0″E, Thanh Xuan (Soc Son district) 21°14′04.6″N 105°45′59.3″E, Van Duc (Gia Lam district) 20°56′23.0″N 105°53′17.7″E, and Van Noi (Dong Anh district) 21°08′49.0″N 105°48′50.2″E. The samples were preliminarily sampled and prepared followed by Vietnam Standard TCVN 9016 : 2011 and TCVN 8551 : 2010 and then washed several times with distilled water.Figure 1
The geographical locations of the pakchoi cultivation regions.


### 2.2. Chemicals and Reagents

Nitric acid 65% (HNO3) and hydrogen peroxide 30% (H2O2) solutions were purchased from Merck, USA. Ultrapure deionized water with the resistivity of 18.2 MΩ cm was obtained from a Milli-Q Plus water puriﬁcation system (Millipore, Bedford, MA, USA). Twenty-six multielement standard solutions including 11B, 23Na, 24Mg, 27Al, 28Si, 39K, 43Ca, 51V, 52Cr, 55Mn, 56Fe, 59Co, 60Ni, 63Cu, 66Zn, 69Ga, 75As, 85Rb, 88Sr, 107Ag, 111Cd, 121Sb, 133Cs, 137Ba, 202Hg, and 208Pb (TraceCERT, periodic table mix 1 for ICP, product no. 92091, lot: BCBW5563) and nine rare earth elements (45Sc, 89Y, 139La, 140Ce, 141Pr, 146Nd, 147Sm, 153Eu, and 157Gd) 10 mg/L each element were provided by Sigma-Aldrich Company. A standard solution containing 50 µg/L of 47Ti, 90Zr, 93Nb, 95Mo, 105Pd, 118Sn, and 182W in 6% ethanol/0.14 M HNO3 was used to determine the sensitivity factors for all elements across the entire mass range for the measurement of diluted samples made in the semiquantitative mode. In case of measuring digested samples, ethanol was omitted in the calibration solution. High-purity ethanol was used for preparing matrix-matched standards. Internal standards for quantitative analysis were prepared in 6% ethanol/0.14 M HNO3 for diluted samples and in 0.14 M HNO3 for digested samples. Blanks for the measurement of diluted and digested samples were a 6% ethanol/0.14 M HNO3 solution containing 50 µg/L of the internal standard and a 0.14 M HNO3 procedure blank submitted to the microwave treatment and including the internal standard, respectively. Standards and the internal standard were prepared by appropriate dilution from 1000 mg/L standard stock solutions.


### 2.3. Sample Preparation and ICP-MS Measurements

All pakchoi samples were lyophilized at −45°C for 3 days and then ground with a pulverizer to obtain find powder (<400μm particle size) before storing in clean plastic bags. Samples (0.25 g) were weighed into a polytetrafluoroethylene (PTFE) vessel, and then, 4 mL H2O, 2 mL HNO3, and 2 mL H2O2 were added. After 30 min, the sample was digested in a Mars X-press plus microwave digestion system (CEM, NC, USA). The digestion program was as follows: the sample was heated to 12°C for 15 min, held for 10 min, before heating to 160°C for 10 min, held for 10 min, and finally heated to 180°C within 10 min and held for 30 min. After cooling, the solution was diluted to 25 mL in a volumetric flask with ultrapure water. The contents were then diluted with deionized water and analyzed [16]. An Agilent 7900 ICP-MS instrument (Agilent Technologies, Tokyo, Japan) was utilized for the measurement of 42 elements in the pakchoi samples, which were 11B, 23Na, 24Mg, 27Al, 28Si, 39K, 43Ca, 45Sc, 47Ti, 51V, 52Cr, 55Mn, 56Fe, 59Co, 60Ni, 63Cu, 66Zn, 69Ga, 75As, 85Rb, 88Sr, 89Y, 90Zr, 93Nb, 95Mo, 105Pd, 107Ag, 111Cd, 118Sn, 121Sb, 133Cs, 137Ba, 139La, 140Ce, 141Pr, 146Nd, 147Sm, 153Eu, 157Gd, 182W, 202Hg, and 208Pb. The analytical parameters of the ICP-MS were RF power at 1550 W, RF matching at 2.0 V, cell entrance at −40 V, cell exit of −60 V, cell energy discrimination at 5.0 V, spray chamber temperature at 2°C, argon was used as carrier gas at flowrate 1.09 L/min, and helium was used as auxiliary gas at 4.3 L/min. Data quantitation was achieved regarding matrix-matched multielement standards that had been prepared in 1% HNO3.In this study, instrument detection limits were calculated using the raw intensity data from the standard and the blank (using ultrapure 2% nitric acid matrix) as per the following equation: IDL = 3SDblankx Cx/(Sx–Sblank), where SDblank is the standard deviation of the intensities of the multiple blank measurements, Cx is the mean signal for the standard, Sx is the signal for Cx, and Sblank is the signal for blank. Method detection limits (MDLs) were calculated as follows: MDL = IDL x constant volume/sample weight.Continuing calibration verification standards were prepared from single element ICP standards (Merck) consisting of 20 mg/L Ca and Mg for the high standard series, and 250µg/L Al, B, Cu, Rb, Sr, and Zn and 25 µg/L Cd, Co, Cs, Ni, Tl, and V for the low standard series. The calibration verifications were measured after every 10 samples.Duplicates of two pakchoi samples were prepared. Possible matrix effects were checked by running an interference check sample consisting of Ca (50 mg/L), Na (100 mg/L), Mg (150 mg/L), Fe (200 mg/L), Cu (250 mg/L), and Zn (500mg/L). In addition, spike recovery tests and serial dilutions were performed on pakchoi samples. The spiked samples were prepared at concentration levels of 20 and 100µg/L for the elements Al, Cu, and Sr and of 100 and 500 µg/L for the elements B, Mn, Rb, and Zn. A serial dilution check (1 : 10 followed by 1 : 3, thus 1 : 30 final dilution) was performed on one pakchoi.


### 2.4. Statistical Analysis

The statistical analysis of the data was performed using STATISTICA 12 (Dell Software, United States). The principal component analysis (PCA) was applied to the data acquired from the ICP-MS analysis to evaluate the discriminants among 6 groups of pakchoi samples from 6 geographical locations. The outputs of multivariate statistical analysis included the scree plot to show contributions of principal components to the PCA model, the score scatterplot illustrates the separation of 6 different groups, the loading scatterplot explains the influences of the elements to the clustering, and the moving range charts illustrate the means and ranges of common distributions of variables on the cases.


## 2.1. Samples

Totally 60 pakchoi (Brassica rapa L. ssp. chinensis) samples were obtained from 6 different locations in Hanoi, Vietnam (Figure 1), which are Linh Nam (Hoang Mai district) 20°58′33.7″N 105°53′21.9″E, Tien Phong (Me Linh district) 21°09′27.2″N 105°44′19.4″E, Thanh Da (Phuc Tho district) 21°06′32.6″N 105°36′43.0″E, Thanh Xuan (Soc Son district) 21°14′04.6″N 105°45′59.3″E, Van Duc (Gia Lam district) 20°56′23.0″N 105°53′17.7″E, and Van Noi (Dong Anh district) 21°08′49.0″N 105°48′50.2″E. The samples were preliminarily sampled and prepared followed by Vietnam Standard TCVN 9016 : 2011 and TCVN 8551 : 2010 and then washed several times with distilled water.Figure 1
The geographical locations of the pakchoi cultivation regions.


## 2.2. Chemicals and Reagents

Nitric acid 65% (HNO3) and hydrogen peroxide 30% (H2O2) solutions were purchased from Merck, USA. Ultrapure deionized water with the resistivity of 18.2 MΩ cm was obtained from a Milli-Q Plus water puriﬁcation system (Millipore, Bedford, MA, USA). Twenty-six multielement standard solutions including 11B, 23Na, 24Mg, 27Al, 28Si, 39K, 43Ca, 51V, 52Cr, 55Mn, 56Fe, 59Co, 60Ni, 63Cu, 66Zn, 69Ga, 75As, 85Rb, 88Sr, 107Ag, 111Cd, 121Sb, 133Cs, 137Ba, 202Hg, and 208Pb (TraceCERT, periodic table mix 1 for ICP, product no. 92091, lot: BCBW5563) and nine rare earth elements (45Sc, 89Y, 139La, 140Ce, 141Pr, 146Nd, 147Sm, 153Eu, and 157Gd) 10 mg/L each element were provided by Sigma-Aldrich Company. A standard solution containing 50 µg/L of 47Ti, 90Zr, 93Nb, 95Mo, 105Pd, 118Sn, and 182W in 6% ethanol/0.14 M HNO3 was used to determine the sensitivity factors for all elements across the entire mass range for the measurement of diluted samples made in the semiquantitative mode. In case of measuring digested samples, ethanol was omitted in the calibration solution. High-purity ethanol was used for preparing matrix-matched standards. Internal standards for quantitative analysis were prepared in 6% ethanol/0.14 M HNO3 for diluted samples and in 0.14 M HNO3 for digested samples. Blanks for the measurement of diluted and digested samples were a 6% ethanol/0.14 M HNO3 solution containing 50 µg/L of the internal standard and a 0.14 M HNO3 procedure blank submitted to the microwave treatment and including the internal standard, respectively. Standards and the internal standard were prepared by appropriate dilution from 1000 mg/L standard stock solutions.


## 2.3. Sample Preparation and ICP-MS Measurements

All pakchoi samples were lyophilized at −45°C for 3 days and then ground with a pulverizer to obtain find powder (<400μm particle size) before storing in clean plastic bags. Samples (0.25 g) were weighed into a polytetrafluoroethylene (PTFE) vessel, and then, 4 mL H2O, 2 mL HNO3, and 2 mL H2O2 were added. After 30 min, the sample was digested in a Mars X-press plus microwave digestion system (CEM, NC, USA). The digestion program was as follows: the sample was heated to 12°C for 15 min, held for 10 min, before heating to 160°C for 10 min, held for 10 min, and finally heated to 180°C within 10 min and held for 30 min. After cooling, the solution was diluted to 25 mL in a volumetric flask with ultrapure water. The contents were then diluted with deionized water and analyzed [16]. An Agilent 7900 ICP-MS instrument (Agilent Technologies, Tokyo, Japan) was utilized for the measurement of 42 elements in the pakchoi samples, which were 11B, 23Na, 24Mg, 27Al, 28Si, 39K, 43Ca, 45Sc, 47Ti, 51V, 52Cr, 55Mn, 56Fe, 59Co, 60Ni, 63Cu, 66Zn, 69Ga, 75As, 85Rb, 88Sr, 89Y, 90Zr, 93Nb, 95Mo, 105Pd, 107Ag, 111Cd, 118Sn, 121Sb, 133Cs, 137Ba, 139La, 140Ce, 141Pr, 146Nd, 147Sm, 153Eu, 157Gd, 182W, 202Hg, and 208Pb. The analytical parameters of the ICP-MS were RF power at 1550 W, RF matching at 2.0 V, cell entrance at −40 V, cell exit of −60 V, cell energy discrimination at 5.0 V, spray chamber temperature at 2°C, argon was used as carrier gas at flowrate 1.09 L/min, and helium was used as auxiliary gas at 4.3 L/min. Data quantitation was achieved regarding matrix-matched multielement standards that had been prepared in 1% HNO3.In this study, instrument detection limits were calculated using the raw intensity data from the standard and the blank (using ultrapure 2% nitric acid matrix) as per the following equation: IDL = 3SDblankx Cx/(Sx–Sblank), where SDblank is the standard deviation of the intensities of the multiple blank measurements, Cx is the mean signal for the standard, Sx is the signal for Cx, and Sblank is the signal for blank. Method detection limits (MDLs) were calculated as follows: MDL = IDL x constant volume/sample weight.Continuing calibration verification standards were prepared from single element ICP standards (Merck) consisting of 20 mg/L Ca and Mg for the high standard series, and 250µg/L Al, B, Cu, Rb, Sr, and Zn and 25 µg/L Cd, Co, Cs, Ni, Tl, and V for the low standard series. The calibration verifications were measured after every 10 samples.Duplicates of two pakchoi samples were prepared. Possible matrix effects were checked by running an interference check sample consisting of Ca (50 mg/L), Na (100 mg/L), Mg (150 mg/L), Fe (200 mg/L), Cu (250 mg/L), and Zn (500mg/L). In addition, spike recovery tests and serial dilutions were performed on pakchoi samples. The spiked samples were prepared at concentration levels of 20 and 100µg/L for the elements Al, Cu, and Sr and of 100 and 500 µg/L for the elements B, Mn, Rb, and Zn. A serial dilution check (1 : 10 followed by 1 : 3, thus 1 : 30 final dilution) was performed on one pakchoi.


## 2.4. Statistical Analysis

The statistical analysis of the data was performed using STATISTICA 12 (Dell Software, United States). The principal component analysis (PCA) was applied to the data acquired from the ICP-MS analysis to evaluate the discriminants among 6 groups of pakchoi samples from 6 geographical locations. The outputs of multivariate statistical analysis included the scree plot to show contributions of principal components to the PCA model, the score scatterplot illustrates the separation of 6 different groups, the loading scatterplot explains the influences of the elements to the clustering, and the moving range charts illustrate the means and ranges of common distributions of variables on the cases.


## 3. Results and Discussion

### 3.1. Elemental Analysis of Pakchoi by ICP-MS

Table1 shows the contents of 42 elements in pakchoi obtained from 6 areas in Hanoi, including Tien Phong, Thanh Da, Linh Nam, Thanh Xuan, Van Duc, and Van Noi. K was the most abundant element in pakchoi samples which accounted for 52–62% of the total elements. Furthermore, Ca was the second major element present in pakchoi which accounted for 30–40% of the total elements. The highest total element contents (about 21 mg/g of fresh material) were found in the samples from Van Duc and Thanh Xuan, whereas the samples from Linh Nam and Van Noi showed the lowest contents of elements at about 16 mg/g of fresh material. The concentration of each element in the samples varied significantly due to the differences in geographical locations.Table 1
Contents (µg/g, fresh weight base) of the 42 elements in pakchoi from 6 locations in Hanoi.
ElementLODLOQCultivation regionTien PhongThanh DaLinh NamThanh XuanVan DucVan Noi(n = 10)(n = 10)(n = 10)(n = 10)(n = 10)(n = 10)11B0.031 ± 0.0020.103 ± 0.006<MQL2.225 ± 0.0012.296 ± 0.1491.679 ± 0.2541.103 ± 0.0281.367 ± 0.12923Na2.220 ± 0.0547.399 ± 0.179740.7 ± 3.2194.33 ± 3.43343 ± 271056 ± 7352.6 ± 3.6401 ± 824Mg9.090 ± 0.08130.29 ± 0.269528.8 ± 5.1543.4 ± 5.9536.5 ± 49.1543.1 ± 42629.6 ± 10.5518.6 ± 6.727Al1.100 ± 0.0073.666 ± 0.02358.85 ± 5.6242.1 ± 7.7130.5 ± 9.6102.8 ± 1.2305.8 ± 20.272.63 ± 4.1628Si28.07 ± 0.10393.55 ± 0.34245.77 ± 3.237.53 ± 1.83212.4 ± 10.8302.5 ± 33.1126.7 ± 1.8163.9 ± 9.939K12.01 ± 0.05140.02 ± 0.16911369 ± 4410464 ± 158683 ± 8511915 ± 2012113 ± 1339220 ± 3643Ca4.003 ± 0.01113.34 ± 0.0365422 ± 416425 ± 626626 ± 617517 ± 568268 ± 525983 ± 2445Sc0.667 ± 0.0042.223 ± 0.0130.135 ± 0.0110.272 ± 0.0260.215 ± 0.1210.213 ± 0.0760.241 ± 0.0120.207 ± 0.01147Ti0.429 ± 0.0031.429 ± 0.0091.383 ± 0.0122.726 ± 0.1211.246 ± 0.1222.062 ± 0.1691.608 ± 0.0731.255 ± 0.19151V0.056 ± 0.0060.186 ± 0.0190.028 ± 0.0020.065 ± 0.060.042 ± 0.0030.034 ± 0.0030.054 ± 0.0020.031 ± 0.00252Cr0.072 ± 0.0010.239 ± 0.0030.123 ± 0.0030.121 ± 0.0070.111 ± 0.0040.192 ± 0.0010.157 ± 0.0060.214 ± 0.00655Mn0.114 ± 0.0020.379 ± 0.00620.34 ± 1.079.976 ± 0.31110.98 ± 0.1122.9 ± 0.1320.44 ± 0.1814.11 ± 0.3756Fe2.030 ± 0.0076.765 ± 0.02320.39 ± 1.0550.89 ± 0.1831.59 ± 2.3526.17 ± 2.2121.12 ± 1.9426.97 ± 0.4359Co0.005 ± 0.0010.016 ± 0.0030.018 ± 0.0010.027 ± 0.0020.041 ± 0.0440.023 ± 0.0020.026 ± 0.0020.017 ± 0.00160Ni0.023 ± 0.0030.076 ± 0.0090.117 ± 0.0010.065 ± 0.0030.061 ± 0.0050.279 ± 0.0020.274 ± 0.0920.364 ± 0.01563Cu0.027 ± 0.0020.089 ± 0.0060.35 ± 0.020.219 ± 0.0110.441 ± 0.0030.979 ± 0.0190.702 ± 0.0740.786 ± 0.03666Zn0.703 ± 0.0262.343 ± 0.0867.148 ± 0.099.204 ± 0.8118.674 ± 0.57215.73 ± 1.2116.74 ± 0.30212.73 ± 0.15669Ga0.070 ± 0.0010.233 ± 0.0030.399 ± 0.0170.561 ± 0.0050.636 ± 0.0031.079 ± 0.0251.835 ± 0.0370.525 ± 0.03575As0.005 ± 0.0010.016 ± 0.003<MQL0.022 ± 0.0020.041 ± 0.0020.073 ± 0.190.044 ± 0.0040.032 ± 0.00385Rb0.004 ± 0.0010.013 ± 0.0033.111 ± 0.0442.474 ± 0.2184.783 ± 0.3218.548 ± 0.5816.638 ± 0.0584.587 ± 0.13988Sr0.091 ± 0.0160.303 ± 0.05317.77 ± 0.3215.57 ± 0.5713.6 ± 1.218.5 ± 0.817.18 ± 0.0115.81 ± 0.0289Y0.018 ± 0.0010.059 ± 0.0060.028 ± 0.0040.046 ± 0.0040.032 ± 0.0010.044 ± 0.0030.047 ± 0.0020.0302 ± 0.002290Zr0.0021 ± 0.00010.006 ± 0.00030.006 ± 0.0010.004 ± 0.0010.017 ± 0.0010.0048 ± 0.00030.021 ± 0.0010.0062 ± 0.000493Nb0.0014 ± 0.00010.004 ± 0.00030.004 ± 0.0010.004 ± 0.0020.003 ± 0.0030.0058 ± 0.0040.016 ± 0.0010.045 ± 0.00295Mo0.032 ± 0.0020.106 ± 0.0060.075 ± 0.0030.032 ± 0.0030.026 ± 0.0020.155 ± 0.0020.070 ± 0.0310.093 ± 0.003105Pd0.0091 ± 0.00020.030 ± 0.00060.034 ± 0.0020.056 ± 0.0020.038 ± 0.0030.039 ± 0.0070.045 ± 0.0030.035 ± 0.002107Ag0.0041 ± 0.00010.013 ± 0.00030.106 ± 0.0070.005 ± 0.0010.005 ± 0.0010.042 ± 0.0032.687 ± 0.4490.018 ± 0.001111Cd0.0032 ± 0.00010.010 ± 0.00030.011 ± 0.0010.023 ± 0.0020.036 ± 0.0030.027 ± 0.0020.040 ± 0.0030.033 ± 0.003118Sn0.061 ± 0.0020.203 ± 0.0060.093 ± 0.0040.015 ± 0.0010.255 ± 0.0130.121 ± 0.0010.155 ± 0.0090.123 ± 0.007121Sb0.0035 ± 0.00010.011 ± 0.00030.042 ± 0.0030.059 ± 0.0040.041 ± 0.0310.122 ± 0.0050.131 ± 0.0270.075 ± 0.005133Cs0.0027 ± 0.00020.008 ± 0.00060.012 ± 0.0010.039 ± 0.0030.041 ± 0.0020.035 ± 0.0030.044 ± 0.0010.020 ± 0.001137Ba0.523 ± 0.0481.743 ± 0.1591.318 ± 0.0941.31 ± 0.092.616 ± 0.1112.763 ± 0.5475.953 ± 0.4791.64 ± 0.11139La0.009 ± 0.0010.029 ± 0.0030.016 ± 0.0020.042 ± 0.0020.044 ± 0.0020.036 ± 0.0020.023 ± 0.0020.020 ± 0.002140Ce0.0052 ± 0.00010.017 ± 0.00030.018 ± 0.0010.055 ± 0.0020.031 ± 0.0020.045 ± 0.0020.057 ± 0.0030.028 ± 0.007141Pr0.0017 ± 0.00010.005 ± 0.00030.003 ± 0.0000.011 ± 0.0010.0043 ± 0.00030.006 ± 0.0010.0057 ± 0.0020.0028 ± 0.0001146Nd0.010 ± 0.0010.033 ± 0.0030.012 ± 0.0010.033 ± 0.0020.015 ± 0.0010.032 ± 0.0130.024 ± 0.0020.014 ± 0.001147Sm0.0025 ± 0.00010.008 ± 0.0003<MQL0.008 ± 0.0010.008 ± 0.0010.014 ± 0.0010.0032 ± 0.00020.0053 ± 0.0003153Eu0.0028 ± 0.00020.009 ± 0.00060.003 ± 0.0010.004 ± 0.0010.006 ± 0.0010.007 ± 0.4980.014 ± 0.0010.0038 ± 0.0001157Gd0.0014 ± 0.00020.004 ± 0.00060.003 ± 0.0010.011 ± 0.0010.006 ± 0.0010.013 ± 0.0010.015 ± 0.0010.0038 ± 0.0002182W0.0011 ± 0.00010.003 ± 0.0060.003 ± 0.0010.002 ± 0.0010.023 ± 0.0020.014 ± 0.0010.368 ± 0.0070.012 ± 0.001202Hg0.0005 ± 0.00000.001 ± 0.00010.003 ± 0.0010.0011 ± 0.00010.005 ± 0.0010.0058 ± 0.00030.0045 ± 0.0030.0072 ± 0.0001208Pb0.014 ± 0.0010.046 ± 0.0030.023 ± 0.0020.0010 ± 0.00010.301 ± 0.0020.489 ± 0.0090.354 ± 0.0020.512 ± 0.034MQL, method quantification limit.


### 3.2. Principal Components Analysis (PCA)

An unsupervised principal component analysis (PCA) was conducted to visualize the effects of geographical locations on elements and the discriminant of 60 pakchoi samples. Figure2 shows the scree plot which is used to illustrate the contributions of principal components (PCs) in a PCA model. As can be seen that the first three PCs accounted for more than 60% of the total variation of samples (Table 2), we can suppose that these PCs carried major information of variables. The quality of the PCA models was represented by R2 and Q2 values. The R2X(cum) value at 0.499 and the Q2(cum) value at 0.617 obtained from the PCA model were established by PC1 and PC2. The result indicated that 49.9% and 61.7% of the total variation could be explained and predicted, respectively, by the first two PCs.Figure 2
Eigenvalues scree plot of the principal components.Table 2
Principal component analysis summary.
PCR2X(Cum)Q2X(Cum)Eigenvalues10.28830.258212.108820.49930.61728.861530.61220.83714.741340.71370.94984.265650.77260.98752.472960.81610.99751.826070.84860.99961.364280.87540.99991.125390.89781.00000.9402100.91121.00000.5660In the two-dimensional score scatterplot based on PC1 and PC2 (Figure3), pakchoi samples from 6 areas were sharply separated, which highlighted the possibility to distinguish the origins of this vegetable merely by the metal elements distribution. In addition, the contribution of variables to classify was determined by the loading scatterplot (Figure 4). The contents of  69Ga, 157Gd, and 153Eu had the highest weight on the first PC, while 56Fe and 64Cu showed the largest contribution to the separation on PC2.Figure 3
PCA scores scatterplot of principal components 1 and 2.Figure 4
PCA loading scatterplots.Moving range charts illustrated the means and ranges of common distributions of variables on the cases. As can be seen, Thanh Da samples had significantly higher contents of56Fe than the samples from other areas that could be used as a marker for distinguishing (Figure 5(b)). Meanwhile, in the samples from Linh Nam, the elemental isotopes 59Co and 118Sn were dramatically increased as compared to samples from other sites that could be found as two important discriminating elements for this area (Figures 5(c) and 5(j)). Whereas, a unique combination of three elements (28Si, 63Cu, and 85Rb) was identified as the main discriminator of pakchoi from Thanh Xuan (Figure 5(a), 5(d), and 5(g)). Similarly, a group of three elements, namely, 69Ga, 107Ag, and 137Ba, contributed strongly to the discrimination of samples from Van Duc (Figure 5(e), 5(i), and 5(k)). The levels of 107Ag and 137Ba increased sharply in Van Duc samples, compared to significantly low contents in the samples of other areas. For the Van Noi samples, the content of 93Nb was three-fold to tenth-fold higher than the levels of this element in the other samples (Figure 5(h)). Thus, 93Nb was responsible for the unique classification of samples from Van Noi. Interestingly, no elemental isotope was identified that contained at a dramatically high level in Tien Phong samples; however, it was possible to discriminate this area’s samples based on the significantly low content of 69Ga, which was also identified as a discriminator for Van Duc (Figure 5(a)). This contrast was also seen on scores scatterplots where samples from Tien Phong and Van Duc distributed on two contrary ranges of PC1.Figure 5
Moving range chart of (a)28Si, (b) 56Fe, (c) 59Co, (d) 63Cu, (e) 69Ga, (f) 75As, (g) 85Rb, (h) 93Nb, (i) 107Ag, (j) 118Sn, and (k) 137Ba. The vertical axis illustrates the content (µg/g, fresh weight base) of each element in the sample. The horizontal axis shows names of the areas.
(a)(b)(c)(d)(e)(f)(g)(h)(i)(j)(k)To the best of our knowledge, this is the first time the mineral elements are used for discrimination of pakchoi geographical origins. The chemometric-based ICP-MS approach was applied to distinguish many types of foods and plants, such as cabbage [6], chili [17], and tea [18]. PCA was used widely in those reports that could find clearly and preciously the most important elements for the classification of the analytes.


## 3.1. Elemental Analysis of Pakchoi by ICP-MS

Table1 shows the contents of 42 elements in pakchoi obtained from 6 areas in Hanoi, including Tien Phong, Thanh Da, Linh Nam, Thanh Xuan, Van Duc, and Van Noi. K was the most abundant element in pakchoi samples which accounted for 52–62% of the total elements. Furthermore, Ca was the second major element present in pakchoi which accounted for 30–40% of the total elements. The highest total element contents (about 21 mg/g of fresh material) were found in the samples from Van Duc and Thanh Xuan, whereas the samples from Linh Nam and Van Noi showed the lowest contents of elements at about 16 mg/g of fresh material. The concentration of each element in the samples varied significantly due to the differences in geographical locations.Table 1
Contents (µg/g, fresh weight base) of the 42 elements in pakchoi from 6 locations in Hanoi.
ElementLODLOQCultivation regionTien PhongThanh DaLinh NamThanh XuanVan DucVan Noi(n = 10)(n = 10)(n = 10)(n = 10)(n = 10)(n = 10)11B0.031 ± 0.0020.103 ± 0.006<MQL2.225 ± 0.0012.296 ± 0.1491.679 ± 0.2541.103 ± 0.0281.367 ± 0.12923Na2.220 ± 0.0547.399 ± 0.179740.7 ± 3.2194.33 ± 3.43343 ± 271056 ± 7352.6 ± 3.6401 ± 824Mg9.090 ± 0.08130.29 ± 0.269528.8 ± 5.1543.4 ± 5.9536.5 ± 49.1543.1 ± 42629.6 ± 10.5518.6 ± 6.727Al1.100 ± 0.0073.666 ± 0.02358.85 ± 5.6242.1 ± 7.7130.5 ± 9.6102.8 ± 1.2305.8 ± 20.272.63 ± 4.1628Si28.07 ± 0.10393.55 ± 0.34245.77 ± 3.237.53 ± 1.83212.4 ± 10.8302.5 ± 33.1126.7 ± 1.8163.9 ± 9.939K12.01 ± 0.05140.02 ± 0.16911369 ± 4410464 ± 158683 ± 8511915 ± 2012113 ± 1339220 ± 3643Ca4.003 ± 0.01113.34 ± 0.0365422 ± 416425 ± 626626 ± 617517 ± 568268 ± 525983 ± 2445Sc0.667 ± 0.0042.223 ± 0.0130.135 ± 0.0110.272 ± 0.0260.215 ± 0.1210.213 ± 0.0760.241 ± 0.0120.207 ± 0.01147Ti0.429 ± 0.0031.429 ± 0.0091.383 ± 0.0122.726 ± 0.1211.246 ± 0.1222.062 ± 0.1691.608 ± 0.0731.255 ± 0.19151V0.056 ± 0.0060.186 ± 0.0190.028 ± 0.0020.065 ± 0.060.042 ± 0.0030.034 ± 0.0030.054 ± 0.0020.031 ± 0.00252Cr0.072 ± 0.0010.239 ± 0.0030.123 ± 0.0030.121 ± 0.0070.111 ± 0.0040.192 ± 0.0010.157 ± 0.0060.214 ± 0.00655Mn0.114 ± 0.0020.379 ± 0.00620.34 ± 1.079.976 ± 0.31110.98 ± 0.1122.9 ± 0.1320.44 ± 0.1814.11 ± 0.3756Fe2.030 ± 0.0076.765 ± 0.02320.39 ± 1.0550.89 ± 0.1831.59 ± 2.3526.17 ± 2.2121.12 ± 1.9426.97 ± 0.4359Co0.005 ± 0.0010.016 ± 0.0030.018 ± 0.0010.027 ± 0.0020.041 ± 0.0440.023 ± 0.0020.026 ± 0.0020.017 ± 0.00160Ni0.023 ± 0.0030.076 ± 0.0090.117 ± 0.0010.065 ± 0.0030.061 ± 0.0050.279 ± 0.0020.274 ± 0.0920.364 ± 0.01563Cu0.027 ± 0.0020.089 ± 0.0060.35 ± 0.020.219 ± 0.0110.441 ± 0.0030.979 ± 0.0190.702 ± 0.0740.786 ± 0.03666Zn0.703 ± 0.0262.343 ± 0.0867.148 ± 0.099.204 ± 0.8118.674 ± 0.57215.73 ± 1.2116.74 ± 0.30212.73 ± 0.15669Ga0.070 ± 0.0010.233 ± 0.0030.399 ± 0.0170.561 ± 0.0050.636 ± 0.0031.079 ± 0.0251.835 ± 0.0370.525 ± 0.03575As0.005 ± 0.0010.016 ± 0.003<MQL0.022 ± 0.0020.041 ± 0.0020.073 ± 0.190.044 ± 0.0040.032 ± 0.00385Rb0.004 ± 0.0010.013 ± 0.0033.111 ± 0.0442.474 ± 0.2184.783 ± 0.3218.548 ± 0.5816.638 ± 0.0584.587 ± 0.13988Sr0.091 ± 0.0160.303 ± 0.05317.77 ± 0.3215.57 ± 0.5713.6 ± 1.218.5 ± 0.817.18 ± 0.0115.81 ± 0.0289Y0.018 ± 0.0010.059 ± 0.0060.028 ± 0.0040.046 ± 0.0040.032 ± 0.0010.044 ± 0.0030.047 ± 0.0020.0302 ± 0.002290Zr0.0021 ± 0.00010.006 ± 0.00030.006 ± 0.0010.004 ± 0.0010.017 ± 0.0010.0048 ± 0.00030.021 ± 0.0010.0062 ± 0.000493Nb0.0014 ± 0.00010.004 ± 0.00030.004 ± 0.0010.004 ± 0.0020.003 ± 0.0030.0058 ± 0.0040.016 ± 0.0010.045 ± 0.00295Mo0.032 ± 0.0020.106 ± 0.0060.075 ± 0.0030.032 ± 0.0030.026 ± 0.0020.155 ± 0.0020.070 ± 0.0310.093 ± 0.003105Pd0.0091 ± 0.00020.030 ± 0.00060.034 ± 0.0020.056 ± 0.0020.038 ± 0.0030.039 ± 0.0070.045 ± 0.0030.035 ± 0.002107Ag0.0041 ± 0.00010.013 ± 0.00030.106 ± 0.0070.005 ± 0.0010.005 ± 0.0010.042 ± 0.0032.687 ± 0.4490.018 ± 0.001111Cd0.0032 ± 0.00010.010 ± 0.00030.011 ± 0.0010.023 ± 0.0020.036 ± 0.0030.027 ± 0.0020.040 ± 0.0030.033 ± 0.003118Sn0.061 ± 0.0020.203 ± 0.0060.093 ± 0.0040.015 ± 0.0010.255 ± 0.0130.121 ± 0.0010.155 ± 0.0090.123 ± 0.007121Sb0.0035 ± 0.00010.011 ± 0.00030.042 ± 0.0030.059 ± 0.0040.041 ± 0.0310.122 ± 0.0050.131 ± 0.0270.075 ± 0.005133Cs0.0027 ± 0.00020.008 ± 0.00060.012 ± 0.0010.039 ± 0.0030.041 ± 0.0020.035 ± 0.0030.044 ± 0.0010.020 ± 0.001137Ba0.523 ± 0.0481.743 ± 0.1591.318 ± 0.0941.31 ± 0.092.616 ± 0.1112.763 ± 0.5475.953 ± 0.4791.64 ± 0.11139La0.009 ± 0.0010.029 ± 0.0030.016 ± 0.0020.042 ± 0.0020.044 ± 0.0020.036 ± 0.0020.023 ± 0.0020.020 ± 0.002140Ce0.0052 ± 0.00010.017 ± 0.00030.018 ± 0.0010.055 ± 0.0020.031 ± 0.0020.045 ± 0.0020.057 ± 0.0030.028 ± 0.007141Pr0.0017 ± 0.00010.005 ± 0.00030.003 ± 0.0000.011 ± 0.0010.0043 ± 0.00030.006 ± 0.0010.0057 ± 0.0020.0028 ± 0.0001146Nd0.010 ± 0.0010.033 ± 0.0030.012 ± 0.0010.033 ± 0.0020.015 ± 0.0010.032 ± 0.0130.024 ± 0.0020.014 ± 0.001147Sm0.0025 ± 0.00010.008 ± 0.0003<MQL0.008 ± 0.0010.008 ± 0.0010.014 ± 0.0010.0032 ± 0.00020.0053 ± 0.0003153Eu0.0028 ± 0.00020.009 ± 0.00060.003 ± 0.0010.004 ± 0.0010.006 ± 0.0010.007 ± 0.4980.014 ± 0.0010.0038 ± 0.0001157Gd0.0014 ± 0.00020.004 ± 0.00060.003 ± 0.0010.011 ± 0.0010.006 ± 0.0010.013 ± 0.0010.015 ± 0.0010.0038 ± 0.0002182W0.0011 ± 0.00010.003 ± 0.0060.003 ± 0.0010.002 ± 0.0010.023 ± 0.0020.014 ± 0.0010.368 ± 0.0070.012 ± 0.001202Hg0.0005 ± 0.00000.001 ± 0.00010.003 ± 0.0010.0011 ± 0.00010.005 ± 0.0010.0058 ± 0.00030.0045 ± 0.0030.0072 ± 0.0001208Pb0.014 ± 0.0010.046 ± 0.0030.023 ± 0.0020.0010 ± 0.00010.301 ± 0.0020.489 ± 0.0090.354 ± 0.0020.512 ± 0.034MQL, method quantification limit.


## 3.2. Principal Components Analysis (PCA)

An unsupervised principal component analysis (PCA) was conducted to visualize the effects of geographical locations on elements and the discriminant of 60 pakchoi samples. Figure2 shows the scree plot which is used to illustrate the contributions of principal components (PCs) in a PCA model. As can be seen that the first three PCs accounted for more than 60% of the total variation of samples (Table 2), we can suppose that these PCs carried major information of variables. The quality of the PCA models was represented by R2 and Q2 values. The R2X(cum) value at 0.499 and the Q2(cum) value at 0.617 obtained from the PCA model were established by PC1 and PC2. The result indicated that 49.9% and 61.7% of the total variation could be explained and predicted, respectively, by the first two PCs.Figure 2
Eigenvalues scree plot of the principal components.Table 2
Principal component analysis summary.
PCR2X(Cum)Q2X(Cum)Eigenvalues10.28830.258212.108820.49930.61728.861530.61220.83714.741340.71370.94984.265650.77260.98752.472960.81610.99751.826070.84860.99961.364280.87540.99991.125390.89781.00000.9402100.91121.00000.5660In the two-dimensional score scatterplot based on PC1 and PC2 (Figure3), pakchoi samples from 6 areas were sharply separated, which highlighted the possibility to distinguish the origins of this vegetable merely by the metal elements distribution. In addition, the contribution of variables to classify was determined by the loading scatterplot (Figure 4). The contents of  69Ga, 157Gd, and 153Eu had the highest weight on the first PC, while 56Fe and 64Cu showed the largest contribution to the separation on PC2.Figure 3
PCA scores scatterplot of principal components 1 and 2.Figure 4
PCA loading scatterplots.Moving range charts illustrated the means and ranges of common distributions of variables on the cases. As can be seen, Thanh Da samples had significantly higher contents of56Fe than the samples from other areas that could be used as a marker for distinguishing (Figure 5(b)). Meanwhile, in the samples from Linh Nam, the elemental isotopes 59Co and 118Sn were dramatically increased as compared to samples from other sites that could be found as two important discriminating elements for this area (Figures 5(c) and 5(j)). Whereas, a unique combination of three elements (28Si, 63Cu, and 85Rb) was identified as the main discriminator of pakchoi from Thanh Xuan (Figure 5(a), 5(d), and 5(g)). Similarly, a group of three elements, namely, 69Ga, 107Ag, and 137Ba, contributed strongly to the discrimination of samples from Van Duc (Figure 5(e), 5(i), and 5(k)). The levels of 107Ag and 137Ba increased sharply in Van Duc samples, compared to significantly low contents in the samples of other areas. For the Van Noi samples, the content of 93Nb was three-fold to tenth-fold higher than the levels of this element in the other samples (Figure 5(h)). Thus, 93Nb was responsible for the unique classification of samples from Van Noi. Interestingly, no elemental isotope was identified that contained at a dramatically high level in Tien Phong samples; however, it was possible to discriminate this area’s samples based on the significantly low content of 69Ga, which was also identified as a discriminator for Van Duc (Figure 5(a)). This contrast was also seen on scores scatterplots where samples from Tien Phong and Van Duc distributed on two contrary ranges of PC1.Figure 5
Moving range chart of (a)28Si, (b) 56Fe, (c) 59Co, (d) 63Cu, (e) 69Ga, (f) 75As, (g) 85Rb, (h) 93Nb, (i) 107Ag, (j) 118Sn, and (k) 137Ba. The vertical axis illustrates the content (µg/g, fresh weight base) of each element in the sample. The horizontal axis shows names of the areas.
(a)(b)(c)(d)(e)(f)(g)(h)(i)(j)(k)To the best of our knowledge, this is the first time the mineral elements are used for discrimination of pakchoi geographical origins. The chemometric-based ICP-MS approach was applied to distinguish many types of foods and plants, such as cabbage [6], chili [17], and tea [18]. PCA was used widely in those reports that could find clearly and preciously the most important elements for the classification of the analytes.


## 4. Conclusions

This study is the first report of the method to determine the geographical origin of various pakchoi samples based on ICP-MS analysis and multivariate statistical analysis (PCA). The results illustrated that eleven elemental discriminators for vegetables were identified, namely,28Si, 56Fe, 59Co, 63Cu, 69Ga, 75As, 85Rb, 93Nb, 107Ag, 118Sn, and 137Ba. By using the PCA model, 60 samples of 6 areas in Hanoi were classified clearly with one-half of the variation that could be explained. The findings of this study imply that the discrimination of pakchoi based on their geographical locations could be well achieved by the combination of elemental profiling and multivariate statistical analysis. The reported method is convenient, fast, and environmentally friendly, with potential applications in distinguishing origins of a larger amount, ranges, and types of vegetables.



---

*Source: 8860852-2021-02-09.xml*