CN114414653A - Cherry producing area tracing method based on fruit stones - Google Patents
Cherry producing area tracing method based on fruit stones Download PDFInfo
- Publication number
- CN114414653A CN114414653A CN202210082187.9A CN202210082187A CN114414653A CN 114414653 A CN114414653 A CN 114414653A CN 202210082187 A CN202210082187 A CN 202210082187A CN 114414653 A CN114414653 A CN 114414653A
- Authority
- CN
- China
- Prior art keywords
- tracing
- cherry
- sample
- voltage
- cherries
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 241000167854 Bourreria succulenta Species 0.000 title claims abstract description 73
- 235000019693 cherries Nutrition 0.000 title claims abstract description 73
- 238000000034 method Methods 0.000 title claims abstract description 42
- 235000013399 edible fruits Nutrition 0.000 title claims abstract description 28
- 239000004575 stone Substances 0.000 title claims abstract description 19
- 238000001514 detection method Methods 0.000 claims abstract description 30
- 239000000843 powder Substances 0.000 claims abstract description 25
- 238000004519 manufacturing process Methods 0.000 claims abstract description 24
- 229910052788 barium Inorganic materials 0.000 claims abstract description 20
- 229910052748 manganese Inorganic materials 0.000 claims abstract description 20
- 229910052712 strontium Inorganic materials 0.000 claims abstract description 20
- 229910052701 rubidium Inorganic materials 0.000 claims abstract description 17
- 238000001035 drying Methods 0.000 claims abstract description 5
- 238000000227 grinding Methods 0.000 claims abstract description 5
- 239000007789 gas Substances 0.000 claims description 13
- XKRFYHLGVUSROY-UHFFFAOYSA-N Argon Chemical compound [Ar] XKRFYHLGVUSROY-UHFFFAOYSA-N 0.000 claims description 12
- 238000009616 inductively coupled plasma Methods 0.000 claims description 9
- 229910052786 argon Inorganic materials 0.000 claims description 6
- 238000000605 extraction Methods 0.000 claims description 6
- 238000005259 measurement Methods 0.000 claims description 6
- 239000002253 acid Substances 0.000 claims description 5
- 238000005103 elemental analyser isotope ratio mass spectroscopy Methods 0.000 claims description 5
- 238000000120 microwave digestion Methods 0.000 claims description 5
- 238000002156 mixing Methods 0.000 claims description 5
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 5
- 230000001133 acceleration Effects 0.000 claims description 4
- 239000012159 carrier gas Substances 0.000 claims description 4
- 238000000889 atomisation Methods 0.000 claims description 3
- 230000029087 digestion Effects 0.000 claims description 3
- 239000001307 helium Substances 0.000 claims description 3
- 229910052734 helium Inorganic materials 0.000 claims description 3
- SWQJXJOGLNCZEY-UHFFFAOYSA-N helium atom Chemical compound [He] SWQJXJOGLNCZEY-UHFFFAOYSA-N 0.000 claims description 3
- 230000002572 peristaltic effect Effects 0.000 claims description 3
- 238000005070 sampling Methods 0.000 claims description 3
- 230000000087 stabilizing effect Effects 0.000 claims description 3
- 239000000126 substance Substances 0.000 claims description 3
- 238000012360 testing method Methods 0.000 claims description 3
- 239000000463 material Substances 0.000 claims description 2
- 238000002203 pretreatment Methods 0.000 claims 1
- 239000002366 mineral element Substances 0.000 abstract description 16
- 235000015203 fruit juice Nutrition 0.000 abstract description 3
- 238000004458 analytical method Methods 0.000 description 18
- IJGRMHOSHXDMSA-UHFFFAOYSA-N Atomic nitrogen Chemical compound N#N IJGRMHOSHXDMSA-UHFFFAOYSA-N 0.000 description 8
- 238000001095 inductively coupled plasma mass spectrometry Methods 0.000 description 7
- 238000000513 principal component analysis Methods 0.000 description 7
- 150000002500 ions Chemical class 0.000 description 6
- 229910052791 calcium Inorganic materials 0.000 description 5
- 238000007621 cluster analysis Methods 0.000 description 5
- 229910052802 copper Inorganic materials 0.000 description 5
- 229910052759 nickel Inorganic materials 0.000 description 5
- 229910052700 potassium Inorganic materials 0.000 description 5
- OKTJSMMVPCPJKN-UHFFFAOYSA-N Carbon Chemical compound [C] OKTJSMMVPCPJKN-UHFFFAOYSA-N 0.000 description 4
- XEEYBQQBJWHFJM-UHFFFAOYSA-N Iron Chemical compound [Fe] XEEYBQQBJWHFJM-UHFFFAOYSA-N 0.000 description 4
- QVGXLLKOCUKJST-UHFFFAOYSA-N atomic oxygen Chemical compound [O] QVGXLLKOCUKJST-UHFFFAOYSA-N 0.000 description 4
- 229910052799 carbon Inorganic materials 0.000 description 4
- 230000000052 comparative effect Effects 0.000 description 4
- 235000013305 food Nutrition 0.000 description 4
- 229910052757 nitrogen Inorganic materials 0.000 description 4
- 239000001301 oxygen Substances 0.000 description 4
- 229910052760 oxygen Inorganic materials 0.000 description 4
- 229910052718 tin Inorganic materials 0.000 description 4
- UFHFLCQGNIYNRP-UHFFFAOYSA-N Hydrogen Chemical compound [H][H] UFHFLCQGNIYNRP-UHFFFAOYSA-N 0.000 description 3
- 239000001257 hydrogen Substances 0.000 description 3
- 229910052739 hydrogen Inorganic materials 0.000 description 3
- 229910052742 iron Inorganic materials 0.000 description 3
- ATJFFYVFTNAWJD-UHFFFAOYSA-N Tin Chemical compound [Sn] ATJFFYVFTNAWJD-UHFFFAOYSA-N 0.000 description 2
- 229910052782 aluminium Inorganic materials 0.000 description 2
- 229910052787 antimony Inorganic materials 0.000 description 2
- 229910052785 arsenic Inorganic materials 0.000 description 2
- 238000004364 calculation method Methods 0.000 description 2
- 235000015120 cherry juice Nutrition 0.000 description 2
- 229910052804 chromium Inorganic materials 0.000 description 2
- 235000011389 fruit/vegetable juice Nutrition 0.000 description 2
- 229910052745 lead Inorganic materials 0.000 description 2
- 229910052749 magnesium Inorganic materials 0.000 description 2
- 230000008569 process Effects 0.000 description 2
- 229910052708 sodium Inorganic materials 0.000 description 2
- 229910052719 titanium Inorganic materials 0.000 description 2
- 229910052720 vanadium Inorganic materials 0.000 description 2
- 229910052725 zinc Inorganic materials 0.000 description 2
- 229910052779 Neodymium Inorganic materials 0.000 description 1
- GRYLNZFGIOXLOG-UHFFFAOYSA-N Nitric acid Chemical compound O[N+]([O-])=O GRYLNZFGIOXLOG-UHFFFAOYSA-N 0.000 description 1
- 229910052777 Praseodymium Inorganic materials 0.000 description 1
- 229910052772 Samarium Inorganic materials 0.000 description 1
- BQCADISMDOOEFD-UHFFFAOYSA-N Silver Chemical compound [Ag] BQCADISMDOOEFD-UHFFFAOYSA-N 0.000 description 1
- 241001465385 Solandra grandiflora Species 0.000 description 1
- 239000000654 additive Substances 0.000 description 1
- 230000000996 additive effect Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000008901 benefit Effects 0.000 description 1
- 229910052797 bismuth Inorganic materials 0.000 description 1
- 229910052793 cadmium Inorganic materials 0.000 description 1
- 229910052792 caesium Inorganic materials 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 230000001186 cumulative effect Effects 0.000 description 1
- 238000007405 data analysis Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000004993 emission spectroscopy Methods 0.000 description 1
- 238000000556 factor analysis Methods 0.000 description 1
- 229910052733 gallium Inorganic materials 0.000 description 1
- 238000000265 homogenisation Methods 0.000 description 1
- 229910052738 indium Inorganic materials 0.000 description 1
- 229910052746 lanthanum Inorganic materials 0.000 description 1
- 229910052744 lithium Inorganic materials 0.000 description 1
- 238000004949 mass spectrometry Methods 0.000 description 1
- 229910052750 molybdenum Inorganic materials 0.000 description 1
- 229910017604 nitric acid Inorganic materials 0.000 description 1
- 238000010606 normalization Methods 0.000 description 1
- 238000001543 one-way ANOVA Methods 0.000 description 1
- 230000003647 oxidation Effects 0.000 description 1
- 238000007254 oxidation reaction Methods 0.000 description 1
- 238000002360 preparation method Methods 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 229910052702 rhenium Inorganic materials 0.000 description 1
- 229910052707 ruthenium Inorganic materials 0.000 description 1
- 229910052706 scandium Inorganic materials 0.000 description 1
- 229910052709 silver Inorganic materials 0.000 description 1
- 239000004332 silver Substances 0.000 description 1
- 238000007619 statistical method Methods 0.000 description 1
- 230000009897 systematic effect Effects 0.000 description 1
- 229910052716 thallium Inorganic materials 0.000 description 1
- 229940088594 vitamin Drugs 0.000 description 1
- 229930003231 vitamin Natural products 0.000 description 1
- 235000013343 vitamin Nutrition 0.000 description 1
- 239000011782 vitamin Substances 0.000 description 1
- 150000003722 vitamin derivatives Chemical class 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N27/00—Investigating or analysing materials by the use of electric, electrochemical, or magnetic means
- G01N27/62—Investigating or analysing materials by the use of electric, electrochemical, or magnetic means by investigating the ionisation of gases, e.g. aerosols; by investigating electric discharges, e.g. emission of cathode
- G01N27/626—Investigating or analysing materials by the use of electric, electrochemical, or magnetic means by investigating the ionisation of gases, e.g. aerosols; by investigating electric discharges, e.g. emission of cathode using heat to ionise a gas
Landscapes
- Chemical & Material Sciences (AREA)
- Chemical Kinetics & Catalysis (AREA)
- Electrochemistry (AREA)
- Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- General Physics & Mathematics (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Other Investigation Or Analysis Of Materials By Electrical Means (AREA)
Abstract
The invention discloses a cherry producing area tracing method based on fruit stones, which comprises the following steps: (1) taking kernels of a sample to be detected, drying and grinding to obtain a powder sample; (2) determination of Mn, Co, Rb, Sr, Ba and delta in powder samples2The content of H element; (3) and substituting the measured value into a discrimination function, calculating a Y value and judging the origin of the sample to be measured according to the size of the Y value. The cherry producing area tracing method takes the fruit stone as a detection object, and compared with fruit juice as the detection object, the fruit stone penetrates through the whole growth cycle of the fruit and is more representative of the producing area; meanwhile, the index for indicating the place of production in the present invention contains both mineral elements (Mn, Co, Rb, Sr, and Ba) and stable isotopes (δ)2H) The discrimination rate of the cherry producing areas after the two are combined is greatly improved; the origin tracing method of the invention has the whole discrimination rate of the cherry origin up to 92.1%.
Description
Technical Field
The invention belongs to the technical field of traceability identification, and particularly relates to a cherry producing area traceability method based on fruit stones.
Background
The cherry is popular with consumers because of its thick and fresh pulp, firmness and juiciness, and rich vitamin and iron. The cherry is a large cherry with thick peel produced in America, Chilean and other countries, and in recent years, China has introduced cherry trees in Shandong, Liaoning, Sichuan and other places to form a large cherry with taste and quality similar to those of China abroad. In the marketing link, however, prices of imported cherries and local cherries are far different, and the cherries with the same specification may have several times or even dozens of times of price difference due to different production places, so that the production place attribute of the cherries has very important significance for determining the price authenticity of the cherries.
In order to determine the production area of food, traceability systems are established successively in countries and regions such as European Union, America, Japan and the like, and food safety traceability systems are also established legally in the food safety law of China, so that the traceability of food is ensured.
Currently, researchers have conducted studies on the production area of cherries by analyzing the elements or stable isotopes contained in cherry juice. For example, chinese patent application publication No. CN109752441A discloses a cherry/cherry origin tracing method based on multiple mineral elements, which comprises: (1) collecting cherry/cherry samples from different producing areas; (2) performing denucleation and homogenization pretreatment on a sample; (3) performing super microwave pretreatment on the homogenized sample; (4) performing multi-element analysis of ICP-MS (inductively coupled plasma mass spectrometry) and ICP-OES (inductively coupled plasma emission spectroscopy) on the sample subjected to the super microwave treatment to obtain the multi-element content in the cherries/cherries sample; (5) and (3) analyzing the multi-element content data of the known cherry/cherry producing samples obtained in the step (4) by adopting a mode of combining a PCA method and a PLS-DA method, establishing a producing area tracing model, further analyzing mineral elements of the samples to be detected, and analyzing and judging the producing area of the samples to be detected according to the producing area tracing model.
The method for tracing the origin of the producing area has the following defects: (1) cherry juice is taken as a detection object, but the juice contains excessive sugar and water, which bring great influence on detection and needs more complicated pretreatment steps, and on the other hand, partial information contained in the pulp does not exist in the whole growth cycle of the fruit along with the change of growth and water content, so that the basis of tracing the production place by taking the juice as the source is deviated; (2) the origin tracing model needs to detect 27 main contribution rate elements, the detection cost is high, and the detection process is complicated.
Disclosure of Invention
The invention aims to provide a cherry producing area tracing method based on fruit stones, the method is simple and efficient, the fruit stones are used as producing area tracing bases, and the detection result is more accurate.
In order to achieve the purpose, the technical scheme of the invention is as follows:
a cherry producing area tracing method based on fruit stones comprises the following steps:
(1) taking kernels of a sample to be detected, drying and grinding to obtain a powder sample;
(2) determination of Mn, Co, Rb, Sr, Ba and delta in powder samples2The content of H element;
(3) substituting the measured values into the following discriminant functions, calculating a Y value and judging the origin of the sample to be measured according to the Y value:
Y(Australia)＝–185.044+4.177Mn–5.773Co+0.248Rb+5.761Sr+3.935Ba–4.586δ2H；
Y(USA)＝–469.467+6.053Mn–20.203Co–1.736Rb+12.548Sr–17.171Ba–7.743δ2H；
Y(New Zealand)＝–277.490+4.829Mn–7.884Co–1.164Rb+18.949Sr–17.518Ba–5.679δ2H；
Y(Chilean)＝–246.049+4.451Mn–13.999Co+1.653Rb+11.061Sr–21.685Ba–5.415δ2H；
Y(China)＝–328.687+7.972Mn–2.906Co–1.312Rb+9.055Sr–30.025Ba–6.153δ2H。
The cherry producing area tracing method takes the fruit stones as the detection objects, and compared with the traditional method that fruit juice is taken as the detection objects, the fruit stones penetrate through the whole growth cycle of fruits and are more representative of producing areas; meanwhile, the index for indicating the place of production in the present invention contains both mineral elements (Mn, Co, Rb, Sr, and Ba) and stable isotopes (δ)2H) The discrimination rate of the cherry producing areas after the two are combined is greatly improved; the origin tracing method of the invention has the whole discrimination rate of the cherry origin up to 92.1%.
Preferably, in the above method for tracing the production place of cherries based on kernels, in the step (1), the kernels do not contain nuts.
Preferably, in the method for tracing the production place of cherries based on the kernels, in the step (2), the powder sample is pretreated and then the contents of Mn, Co, Rb, Sr and Ba elements are measured by using an inductively coupled plasma mass spectrometer;
in addition, elemental analysis-isotope ratio mass spectrometry is used to measure delta in powder samples2The content of H was tested.
Before the inductively coupled plasma mass spectrometer is used for measuring the contents of Mn, Co, Rb, Sr and Ba elements, a powder sample needs to be pretreated so that each mineral element is dissolved in a solution for detection.
Preferably, the pretreatment comprises the following steps in sequence:
1) powder samples were mixed with acid according to 1 g: (30-35) uniformly mixing the materials in a mass-volume ratio of mL, and performing microwave digestion;
2) after complete digestion, acid is removed to 1-2 mL;
3) adding water to desired volume, and mixing.
Preferably, in the above method for tracing the production place of cherries based on the cherries, the detection parameters of the inductively coupled plasma mass spectrometer are as follows:
radio frequency power: 1600 w;
sampling depth: 10.0 mm;
argon carrier gas flow: 0.7L/min;
temperature of the atomization chamber: 2 ℃;
flow rate of argon diluted gas: 0.35L/min;
extraction of lens 1 voltage: 0V;
extraction of lens 2 voltage: -195V;
omega deflection voltage: -80V;
omega lens voltage: 8.6V;
collision cell entrance voltage: -40V;
collision cell exit voltage: -60V;
collision cell gas: helium gas;
collision cell gas flow rate: 5.0 mL/min;
eight-pole deflection voltage: -18.0V;
sample introduction time: 30 s;
peristaltic pump speed: 0.30 rps;
the stabilizing time is as follows: 35 s.
At detection of delta2And during the content of H, the pretreatment is not required, the powder sample and the standard substance are put into a chamber to be tested together for balancing for at least 48H, and then the elemental analyzer and the stable isotope mass spectrometer are sequentially adopted for testing.
Wherein, the detection parameters of the element analyzer are as follows:
H2pressure: 25 psi;
CO pressure: 13 psi;
he pressure: 2 psi.
The detection parameters of the stable isotope mass spectrometer are as follows:
H2acceleration voltage: 4184.67V;
CO acceleration voltage: 4492.39V;
an ionization mode: an EI ion source;
ion source voltage: 3.06 kV;
vacuum degree: 1.2X 10-6mBar；
Current: 1.5 mA;
δ2accuracy of continuous measurement of H<0.4 per mill, and the detection frequency n is 10.
Preferably, in the method for tracing the production place of cherries based on the kernels, if the number of the samples to be detected is greater than 5, the equipment is checked by using the standard sample after detecting at least 5 samples to be detected.
Compared with the prior art, the invention has the beneficial effects that:
the cherry producing area tracing method takes the fruit stones as the detection objects, and compared with the traditional method that fruit juice is taken as the detection objects, the fruit stones penetrate through the whole growth cycle of fruits and are more representative of producing areas; meanwhile, the index for indicating the place of production in the present invention contains both mineral elements (Mn, Co, Rb, Sr, and Ba) and stable isotopes (δ)2H) The discrimination rate of the cherry producing areas after the two are combined is greatly improved; the origin tracing method of the invention has the whole discrimination rate of the cherry origin up to 92.1%.
Drawings
FIG. 1 is a graph showing the results of principal component analysis in example 1;
wherein PC1 denotes the first principal component, PC2 denotes the second principal component, and PC3 denotes the third principal component;
FIG. 2 is a graph showing the results of cluster analysis in example 1;
wherein, A: australia, B: U.S., C: new zealand, D: chile, E: china.
Detailed Description
The technical solution of the present invention will be further described in detail with reference to the accompanying drawings and the detailed description.
The apparatus used in this example is as follows: model 7900 inductively coupled plasma mass spectrometer (Agilent); microwave digestion apparatus type CEM MARS 6 (CEM corporation); flash 2000 model elemental analyzer-Delta V advanced model stable isotope mass spectrometer (Thermo Fisher, inc., for carbon and nitrogen isotope determination); isoprime model 100 stable isotope mass spectrometer (Elementar for hydrogen and oxygen isotope determination), 2.5Triad model lyophilizer (LabConco);
CPA225D type analytical balance (accurate to 0.01mg, sydows scientific instruments (beijing) ltd); tin (9 mm. times.5 mm, Thermo Fisher Co., Ltd.), silver (8 mm. times.5 mm, Elementar Co., Ltd.).
Example 1
In order to trace the origin of the cherry production place, the fruit stone is taken as a detection object in the embodiment, and a proper tracing method is explored. The exploration process is as follows:
1. sample preparation
The powder samples were obtained by removing the fruit pulp and the fruit nut from 7 cherry samples from Australia, 5 cherry samples from the United states, 7 cherry samples from New Zealand, 10 cherry samples from Chile and 9 cherry samples from China (38 cherry samples in total), respectively, drying the fruit nut pieces, and grinding to powder.
2. Determination of mineral element content
(1) Pretreatment: putting 0.25g of powder sample into a microwave digestion instrument, adding 8mL of nitric acid, uniformly mixing, and performing microwave digestion; after digestion, acid is removed to about 2mL, then water is added to the solution to a constant volume of 25mL, and the solution is mixed uniformly for later use;
(2) inductively coupled plasma mass spectrometry (ICP-MS) measurements: detecting the mineral elements (Li, Na, Mg, Al, K, Ca, Sc, Ti, V, Cr, Mn, Fe, Co, Ni, Cu, Zn, Ga, As, Se, Rb, Sr, Mo, Ru, Cd, In, Sn, Sb, Cs, Ba, La, Pr, Nd, Sm, Re, Tl, Pb and Bi) In the sample 37;
the detection parameters of the inductively coupled plasma mass spectrometer are as follows:
radio frequency power: 1600 w;
sampling depth: 10.0 mm;
argon carrier gas flow: 0.7L/min;
temperature of the atomization chamber: 2 ℃;
flow rate of argon diluted gas: 0.35L/min;
extraction of lens 1 voltage: 0V;
extraction of lens 2 voltage: -195V;
omega deflection voltage: -80V;
omega lens voltage: 8.6V;
collision cell entrance voltage: -40V;
collision cell exit voltage: -60V;
collision cell gas: helium gas;
collision cell gas flow rate: 5.0 mL/min;
eight-pole deflection voltage: -18.0V;
sample introduction time: 30 s;
peristaltic pump speed: 0.30 rps;
the stabilizing time is as follows: 35 s.
3. Determination of the content of stable isotopes
(1) Determination of content of carbon and nitrogen stable isotopes
0.3mg (carbon stable isotope delta) was taken13C determination)/0.8 mg (nitrogen stable isotope determination) of the powder sample was placed in a tin cup, sealed, and transferred by an autosampler into a Flash 2000 type elemental analyzer and a Delta V additive type stable isotope mass spectrometer for determination of the contents of carbon and nitrogen stable isotopes, respectively.
The detection parameters of the Flash 2000 type element analyzer are as follows:
he carrier gas flow rate: 200 mL/min;
CO2reference gas flow rate: 90 mL/min;
oxygen flow rate: 180 mL/min;
the temperature of the oxidation furnace: 960 ℃;
the temperature of the reduction furnace is 640 ℃;
the trap temperature was 50 ℃.
The detection parameters of the Delta V Advantage type stable isotope mass spectrometer are as follows:
H2acceleration voltage: 4184.67V;
CO acceleration voltage: 4492.39V;
an ionization mode: an EI ion source;
ion source voltage: 3.06 kV;
vacuum degree: 1.2X 10-6mBar；
Current: 1.5 mA.
(2) Stable isotopes of hydrogen and oxygen (delta)2H、δ18O) content determination
And (3) putting 0.5mg of powder sample into a silver cup, putting the powder sample and the standard substance into an instrument to be tested together for balancing for at least 48h, and then respectively measuring the content of the hydrogen stable isotope and the content of the oxygen stable isotope by adopting an Isoprime 100 type stable isotope mass spectrometer.
Wherein, the detection parameters of the element analyzer are as follows:
H2pressure: 25 psi;
CO pressure: 13 psi;
he pressure: 2 psi.
H2Pressure: 25 psi;
CO pressure: 13 psi;
he pressure: 2 psi.
The detection parameters of the stable isotope mass spectrometer are as follows:
H2acceleration voltage: 4184.67V;
CO acceleration voltage: 4492.39V;
an ionization mode: an EI ion source;
ion source voltage: 3.06 kV;
vacuum degree: 1.2X 10-6mBar；
Current: 1.5 mA;
δ2accuracy of continuous measurement of H<0.4 per mill, and the detection frequency n is 10.
4. Data analysis
And (3) carrying out one-factor variance analysis, Duncan multiple comparison analysis, principal component analysis, clustering analysis and discriminant analysis on the data obtained by detection by adopting SPSS 20.0.
(1) One-way ANOVA and Duncan multiple comparison analysis
When the mineral element and stable isotope measurement is carried out, only 24 mineral elements (Na, Mg, Al, K, Ca, Ti, V, Cr, Mn, Fe, Co, Ni, Cu, Zn, As, Se, Rb, Sr, Mo, Cd, Sn, Sb, Ba, Pb) and 3 stable isotopes (delta, and delta) are found13C、δ2H、δ18O) is detected in the cherry pitTo; after the results of single-factor analysis of variance and Duncan multiple comparison analysis (results shown in Table 1) on the results of detection of the 24 mineral elements and the 3 stable isotopes, only 9 mineral elements (K, Ca, Mn, Co, Ni, Cu, Rb, Sr and Ba) and 2 stable isotopes (delta2H、δ18O) content there was a significant difference (p) between the cherry samples from different sources<0.05)。
TABLE 1
Note: the values in the table are mean ± sd, and different letters in the same column indicate significant differences (p < 0.05).
(2) Principal component analysis
For the above 11 parameters (K, Ca, Mn, Co, Ni, Cu, Rb, Sr, Ba, delta)2H and delta18O) main component analysis (R type: centralization and normalization of the data) the cumulative variance contribution of the 4 principal components was found to reach 78.22%. Wherein the 1 st principal component (PC1) mainly integrates 4 parameter information (K, Co, Rb and delta)2H) The 2 nd principal component (PC2) mainly integrates 3 parameter information (Ni, Cu and delta)18O), the 3 rd principal component (PC3) mainly integrates 2 parameter information (Mn and Ba), and the 4 th principal component (PC4) mainly integrates 2 parameter information (Ca and Sr).
The information of a plurality of parameters in the sample can be comprehensively embodied through principal component analysis, and the standardized scores of PC1, PC2 and PC3 are used for drawing to obtain the graph shown in FIG. 1. As can be seen from FIG. 1, a certain productive division can be made among Chinese cherries in Australia, the United states, New Zealand, Chile and China by the 1 st, 2 nd and 3 rd principal components (PC1, PC2 and PC 3).
(3) Cluster analysis
The normalized scores of the first 4 principal components were used for cluster analysis (systematic clustering, euclidean distance was used for clustering distance, Ward method), and the results are shown in fig. 2.
As can be seen from fig. 2, the tree was cut at a clustering distance of 4.50, and the 38 car centimeter samples were classified into 5 classes, with class 1 being primarily chinese car centimeter samples, class 2 being primarily american car centimeter samples, class 3 being primarily chile centimeter samples (also containing 2 australian car centimeter samples), class 4 being primarily australian car centimeter samples (also containing 2 new zealand car centimeter samples), class 5 being primarily new zealand car centimeter samples (also containing 1 chile car centimeter sample). Therefore, multi-parameter principal component analysis using multi-element combinations and stable isotope ratios can be used to classify cherries of different origins.
(4) Discriminant analysis
The results of the difference analysis, the principal component analysis and the cluster analysis of 9 characteristic mineral element indexes and 2 stable isotope indexes in the cherry samples in different producing areas show that the discrimination of the cherry producing areas has certain feasibility by utilizing the combined parameters of multiple elements and stable isotopes.
In order to further understand the distinguishing situation of each index on the production area of the cherries, the 11 parameters are subjected to stepwise distinguishing analysis (distinguishing analysis: stepwise distinguishing method and Fisher linear distinguishing analysis), effective variables for distinguishing the production areas of the cherries are screened out, and a distinguishing model is established.
At a 0.05 significance level, Mn, Co, Rb, Sr, Ba and delta2H, introducing the 6 index parameters into a model, and establishing a discriminant function as follows:
y (australia) ═ 185.044+4.177 Mn-5.773 Co +0.248Rb +5.761Sr +3.935 Ba-4.586 δ2H；
Y (USA) ═ 469.467+6.053 Mn-20.203 Co-1.736 Rb +12.548 Sr-17.171 Ba-7.743 delta2H；
Y (New Zealand) ═ 277.490+4.829 Mn-7.884 Co-1.164 Rb +18.949 Sr-17.518 Ba-5.679 delta2H；
Y (chile) ═ 246.049+4.451 Mn-13.999 Co +1.653Rb +11.061 Sr-21.685 Ba-5.415 δ2H；
Y (China) ═ 328.687+7.972 Mn-2.906 Co-1.312 Rb +9.055 Sr-30.025 Ba-6.153 delta2H；
When blind samples from different regions are identified and classified, parameter measurement results of the samples can be respectively substituted into the discrimination model, Y values of the regions are compared, and the blind samples belong to the regions with the maximum Y values.
The discrimination function is used for carrying out discrimination analysis on 38 cherry samples, the accuracy of cross check is 92.1%, wherein one sample in Australia, New Zealand and Chilean respectively has misjudgment.
It can be seen that, in this example, the content of 37 mineral elements and four stable isotopes in 38 cherry sample kernels from australia, usa, new zealand, chile and china was determined by inductively coupled plasma mass spectrometry (ICP-MS) and elemental analysis-stable isotope mass spectrometry (EA-IRMS), and variance analysis, principal component analysis, cluster analysis and discriminant analysis were performed by multivariate statistical analysis of each parameter index; the analysis results show that Mn, Co, Rb, Sr, Ba and delta2The overall rate of discrimination of the H combination index on the cherries is 92.1 percent; the method can provide necessary information for the origin tracing analysis of the cherry producing areas, and can be used as an effective reference index for the origin tracing of the cherry producing areas.
In summary, the embodiment provides a cherry producing area tracing method based on a kernel, which includes the following steps:
(1) taking kernels of a sample to be detected, drying and grinding to obtain a powder sample;
(2) determination of Mn, Co, Rb, Sr, Ba and delta in powder samples2The content of H element;
wherein, after the powder sample is pretreated, an inductively coupled plasma mass spectrometer (ICP-MS) is adopted to measure the contents of Mn, Co, Rb, Sr and Ba elements;
simultaneously adopting elemental analysis-isotope ratio mass spectrometry (EA-IRMS) to measure delta in powder sample2Testing the content of H;
(3) substituting the measured values into the following discriminant functions, calculating a Y value and judging the origin of the sample to be measured according to the Y value:
Y(Australia)＝–185.044+4.177Mn–5.773Co+0.248Rb+5.761Sr+3.935Ba–4.586δ2H；
Y(USA)＝–469.467+6.053Mn–20.203Co–1.736Rb+12.548Sr–17.171Ba–7.743δ2H；
Y(New Zealand)＝–277.490+4.829Mn–7.884Co–1.164Rb+18.949Sr–17.518Ba–5.679δ2H；
Y(Chilean)＝–246.049+4.451Mn–13.999Co+1.653Rb+11.061Sr–21.685Ba–5.415δ2H；
Y(China)＝–328.687+7.972Mn–2.906Co–1.312Rb+9.055Sr–30.025Ba–6.153δ2H；
The origin of the sample to be tested is the area corresponding to the maximum Y value.
Comparative example 1
The method of tracing the source of the producing area of the comparative example 1 is basically the same as that of the example 1, except that:
in the step (2), only detecting the contents of five mineral elements of Mn, Co, Rb, Sr and Ba;
in step (3), the discriminant function used is:
Y(Australia)＝–27.555+0.945Mn+10.680Co+1.467Rb+0.820Sr+29.398Ba；
Y(USA)＝–20.650+0.598Mn+7.572Co+0.323Rb+4.206Sr+24.126Ba；
Y(New Zealand)＝–36.070+0.827Mn+12.487Co+0.346Rb+12.832Sr+12.770Ba；
Y(Chilean)＝–26.544+0.636Mn+5.425Co+3.092Rb+5.228Sr+7.195Ba；
Y(China)＝–45.201+3.636Mn+19.168Co+0.324Rb+2.426Sr+2.795Ba；
Through calculation, the discrimination rate of the discrimination function on 38 cherry samples is 89.4%; wherein, China, Chile and Australia are well distinguished, but the United states (2) and New Zealand (2) have misjudgment.
Comparative example 2
The method of tracing the source of the origin of comparative example 2 is substantially the same as that of example 1, except that:
in the step (2), only detecting the contents of four mineral elements of Mn, Co, Sr and Ba;
in step (3), the discriminant function used is:
Y(Australia)＝–24.555+0.911Mn+11.140Co+2.095Sr+28.195Ba；
Y(USA)＝–20.505+0.590Mn+7.673Co+4.487Sr+24.081Ba；
Y(New Zealand)＝–35.903+0.819Mn+12.595Co+13.133Sr+12.722Ba；
Y(Chilean)＝–13.214+0.563Mn+6.395Co+7.916Sr+6.768Ba；
Y(China)＝–45.055+3.629Mn+19.270Co+2.707Sr+2.750Ba；
Through calculation, the discrimination rate of the discrimination function on 38 cherry samples is 84.2%; wherein, China and Chile are better distinguished, but Australia (2), America (3) and New Zealand (1) have misjudgment.
Claims (10)
1. A cherry producing area tracing method based on fruit stones is characterized by comprising the following steps:
(1) taking kernels of a sample to be detected, drying and grinding to obtain a powder sample;
(2) determination of Mn, Co, Rb, Sr, Ba and delta in powder samples2The content of H element;
(3) substituting the measured values into the following discriminant functions, calculating a Y value and judging the origin of the sample to be measured according to the Y value:
Y(Australia)＝–185.044+4.177Mn–5.773Co+0.248Rb+5.761Sr+3.935Ba–4.586δ2H；
Y(USA)＝–469.467+6.053Mn–20.203Co–1.736Rb+12.548Sr–17.171Ba–7.743δ2H；
Y(New Zealand)＝–277.490+4.829Mn–7.884Co–1.164Rb+18.949Sr–17.518Ba–5.679δ2H；
Y(Chilean)＝–246.049+4.451Mn–13.999Co+1.653Rb+11.061Sr–21.685Ba–5.415δ2H；
Y(China)＝–328.687+7.972Mn–2.906Co–1.312Rb+9.055Sr–30.025Ba–6.153δ2H。
2. The method for tracing the production place of cherries based on kernels of claim 1, wherein in step (1), the kernels do not contain nuts.
3. The method for tracing the production place of cherries based on the fruit stones, as claimed in claim 1, wherein in the step (2), the powder sample is pretreated and then the content of Mn, Co, Rb, Sr and Ba elements is measured by an inductively coupled plasma mass spectrometer.
4. The method for tracing the production place of cherries based on the fruit stones as claimed in claim 3, wherein the pre-treatment comprises the following steps in sequence:
1) powder samples were mixed with acid according to 1 g: (30-35) uniformly mixing the materials in a mass-volume ratio of mL, and performing microwave digestion;
2) after complete digestion, acid is removed to 1-2 mL;
3) adding water to desired volume, and mixing.
5. The method for tracing the production place of cherries based on the fruit pits of claim 3, wherein the detection parameters of the inductively coupled plasma mass spectrometer are as follows:
radio frequency power: 1600 w;
sampling depth: 10.0 mm;
argon carrier gas flow: 0.7L/min;
temperature of the atomization chamber: 2 ℃;
flow rate of argon diluted gas: 0.35L/min;
extraction of lens 1 voltage: 0V;
extraction of lens 2 voltage: -195V;
omega deflection voltage: -80V;
omega lens voltage: 8.6V;
collision cell entrance voltage: -40V;
collision cell exit voltage: -60V;
collision cell gas: helium gas;
collision cell gas flow rate: 5.0 mL/min;
eight-pole deflection voltage: -18.0V;
sample introduction time: 30 s;
peristaltic pump speed: 0.30 rps;
the stabilizing time is as follows: 35 s.
6. The method for tracing the production place of cherries based on the fruit stones of claim 1, wherein in the step (2), the elemental analysis-isotope ratio mass spectrometry is adopted to measure the delta in the powder sample2The content of H was tested.
7. The method for tracing the production place of cherries based on the fruit pits as claimed in claim 6, wherein the powder sample and the standard substance are put into the chamber to be tested together for balancing for at least 48h, and then the elemental analyzer and the stable isotope mass spectrometer are sequentially used for testing.
8. The method for tracing the production place of cherries based on the fruit stones of claim 7, wherein the detection parameters of the element analyzer are as follows:
H2pressure: 25 psi;
CO pressure: 13 psi;
he pressure: 2 psi.
9. The method for tracing the production place of cherries based on the cherries of claim 7, wherein the detection parameters of the stable isotope mass spectrometer are as follows:
H2acceleration voltage: 4184.67V;
CO acceleration voltage: 4492.39V;
an ionization mode: an EI ion source;
ion source voltage: 3.06 kV;
vacuum degree: 1.2X 10-6mBar；
Current: 1.5 mA;
δ2accuracy of continuous measurement of H<0.4 per mill, and the detection frequency n is 10.
10. The method for tracing the production place of cherries based on any one of claims 3-9, wherein if the number of samples to be tested is greater than 5, the equipment is checked by using a standard sample after at least every 5 samples to be tested are tested.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202210082187.9A CN114414653B (en) | 2022-01-24 | 2022-01-24 | Cherry origin tracing method based on kernel |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202210082187.9A CN114414653B (en) | 2022-01-24 | 2022-01-24 | Cherry origin tracing method based on kernel |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN114414653A true CN114414653A (en) | 2022-04-29 |
| CN114414653B CN114414653B (en) | 2023-12-01 |
Family
ID=81278289
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202210082187.9A Active CN114414653B (en) | 2022-01-24 | 2022-01-24 | Cherry origin tracing method based on kernel |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN114414653B (en) |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109752441A (en) * | 2018-11-12 | 2019-05-14 | 广东出入境检验检疫局检验检疫技术中心 | A kind of car li based on multielement/cherry place of production source tracing method |
| CN110068605A (en) * | 2019-04-11 | 2019-07-30 | 浙江省农业科学院 | A method of based on the PCA-LDA analysis and identification rice place of production |
| CN110687189A (en) * | 2019-10-08 | 2020-01-14 | 浙江省检验检疫科学技术研究院 | Method for distinguishing milk producing area tracing and identifying based on multiple elements and multiple indexes |
| CN111474275A (en) * | 2020-04-29 | 2020-07-31 | 内蒙古自治区农牧业科学院 | Potato origin tracing method based on mineral elements and stable isotopes |
| CN112766980A (en) * | 2020-12-29 | 2021-05-07 | 深圳海关食品检验检疫技术中心 | Durian traceability model modeling method and device and durian traceability method |
-
2022
- 2022-01-24 CN CN202210082187.9A patent/CN114414653B/en active Active
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109752441A (en) * | 2018-11-12 | 2019-05-14 | 广东出入境检验检疫局检验检疫技术中心 | A kind of car li based on multielement/cherry place of production source tracing method |
| CN110068605A (en) * | 2019-04-11 | 2019-07-30 | 浙江省农业科学院 | A method of based on the PCA-LDA analysis and identification rice place of production |
| CN110687189A (en) * | 2019-10-08 | 2020-01-14 | 浙江省检验检疫科学技术研究院 | Method for distinguishing milk producing area tracing and identifying based on multiple elements and multiple indexes |
| CN111474275A (en) * | 2020-04-29 | 2020-07-31 | 内蒙古自治区农牧业科学院 | Potato origin tracing method based on mineral elements and stable isotopes |
| CN112766980A (en) * | 2020-12-29 | 2021-05-07 | 深圳海关食品检验检疫技术中心 | Durian traceability model modeling method and device and durian traceability method |
Non-Patent Citations (1)
| Title |
|---|
| 周秀雯 等: ""基于矿物元素指纹差异的榴莲产地甄别"", 《食品科学》 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN114414653B (en) | 2023-12-01 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN108982692B (en) | Method for distinguishing milk powder producing area by element analysis-stable isotope mass spectrum | |
| Vaclavik et al. | Liquid chromatography–mass spectrometry-based metabolomics for authenticity assessment of fruit juices | |
| Paulsen et al. | Inductively coupled plasma mass spectrometric analysis of ultrapure acids | |
| CN111721857A (en) | Method for identifying litchi varieties by using extensive targeted metabonomics technology | |
| Husáková et al. | Analytical capabilities of inductively coupled plasma orthogonal acceleration time-of-flight mass spectrometry (ICP-oa-TOF-MS) for multi-element analysis of food and beverages | |
| CN110687189A (en) | Method for distinguishing milk producing area tracing and identifying based on multiple elements and multiple indexes | |
| Jarošová et al. | Elemental analysis of coffee: a comparison of ICP-MS and AAS methods. | |
| Faberi et al. | Fatty acid composition and δ13C of bulk and individual fatty acids as marker for authenticating Italian PDO/PGI extra virgin olive oils by means of isotopic ratio mass spectrometry | |
| CN113406247B (en) | Soybean origin tracing identification method based on combination of IRMS, LC-Q-TOF-MS and multi-element analysis | |
| CN110702832A (en) | Application of high-content 4-hydroxyquinoline as characteristic marker of jujube honey | |
| CN112162054A (en) | Authenticity evaluation method of sophora moorcroftianain honey | |
| Wang et al. | Simultaneous high sensitivity trace-element and isotopic analysis of gemstones using laser ablation inductively coupled plasma time-of-flight mass spectrometry | |
| Zhou et al. | Coupling neutral desorption sampling to dielectric barrier discharge ionization mass spectrometry for direct oil analysis | |
| CN114002368A (en) | Method for determining illegal added components in health food by ultra-high performance liquid chromatography-quadrupole-time-of-flight high resolution mass spectrometry | |
| Du et al. | Discrimination of the species and authenticity of Rhizoma Coptidis based on stable isotope and multielement fingerprinting and multivariate statistical analysis | |
| CN114414653A (en) | Cherry producing area tracing method based on fruit stones | |
| JP5758244B2 (en) | Origin determination method of Australian wheat using trace elements and isotope ratio | |
| CN106908527A (en) | A kind of method for differentiating the honey of lychee flowers place of production | |
| CN114420222B (en) | Distributed flow processing-based method for rapidly confirming fragment ion compound structure | |
| Weber et al. | Improvement of the chemometric variety characterization of wines by improving the detection limit for aroma compounds | |
| CN114924002A (en) | Method for identifying NFC and FC orange juice by combining metabonomics technology and machine learning | |
| Heumann | Trace determination and isotopic analysis of the elements in life sciences by mass spectrometry | |
| CN111175368A (en) | Method for identifying synthetic acetic acid added in brewed vinegar | |
| CN112014493B (en) | Method for identifying true and false fruit juice | |
| CN107219321B (en) | A kind of mixing mass spectrum screens out method |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |