CN115436531A - 一种基于大曲非挥发性物质鉴别大曲质量的方法 - Google Patents
一种基于大曲非挥发性物质鉴别大曲质量的方法 Download PDFInfo
- Publication number
- CN115436531A CN115436531A CN202211287619.6A CN202211287619A CN115436531A CN 115436531 A CN115436531 A CN 115436531A CN 202211287619 A CN202211287619 A CN 202211287619A CN 115436531 A CN115436531 A CN 115436531A
- Authority
- CN
- China
- Prior art keywords
- daqu
- quality
- substances
- nonvolatile
- identifying
- 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
- 239000000126 substance Substances 0.000 title claims abstract description 62
- 238000000034 method Methods 0.000 title claims abstract description 46
- 240000004808 Saccharomyces cerevisiae Species 0.000 claims abstract description 33
- 238000012216 screening Methods 0.000 claims abstract description 20
- 238000005516 engineering process Methods 0.000 claims abstract description 13
- 238000007637 random forest analysis Methods 0.000 claims abstract description 12
- 230000000694 effects Effects 0.000 claims abstract description 6
- 239000000203 mixture Substances 0.000 claims abstract description 6
- 238000004422 calculation algorithm Methods 0.000 claims abstract description 5
- 238000013375 chromatographic separation Methods 0.000 claims abstract description 5
- 238000005259 measurement Methods 0.000 claims abstract description 5
- 238000004895 liquid chromatography mass spectrometry Methods 0.000 claims description 8
- 238000000926 separation method Methods 0.000 claims description 7
- 239000000463 material Substances 0.000 claims description 6
- 238000005457 optimization Methods 0.000 claims description 5
- 239000012071 phase Substances 0.000 description 14
- 239000000796 flavoring agent Substances 0.000 description 13
- OKKJLVBELUTLKV-UHFFFAOYSA-N Methanol Chemical compound OC OKKJLVBELUTLKV-UHFFFAOYSA-N 0.000 description 9
- 239000002253 acid Substances 0.000 description 9
- 235000019634 flavors Nutrition 0.000 description 8
- 150000002500 ions Chemical class 0.000 description 7
- 238000004458 analytical method Methods 0.000 description 5
- 238000000855 fermentation Methods 0.000 description 5
- 230000004151 fermentation Effects 0.000 description 5
- 241001480003 Chaetothyriales Species 0.000 description 4
- 238000007599 discharging Methods 0.000 description 4
- 238000004519 manufacturing process Methods 0.000 description 4
- 238000001819 mass spectrum Methods 0.000 description 4
- BDAGIHXWWSANSR-UHFFFAOYSA-N methanoic acid Natural products OC=O BDAGIHXWWSANSR-UHFFFAOYSA-N 0.000 description 4
- 244000005700 microbiome Species 0.000 description 4
- 238000012360 testing method Methods 0.000 description 4
- WEVYAHXRMPXWCK-UHFFFAOYSA-N Acetonitrile Chemical compound CC#N WEVYAHXRMPXWCK-UHFFFAOYSA-N 0.000 description 3
- 239000007864 aqueous solution Substances 0.000 description 3
- 150000001875 compounds Chemical class 0.000 description 3
- 238000010276 construction Methods 0.000 description 3
- 230000001419 dependent effect Effects 0.000 description 3
- 230000006870 function Effects 0.000 description 3
- 238000011160 research Methods 0.000 description 3
- 238000003860 storage Methods 0.000 description 3
- OSWFIVFLDKOXQC-UHFFFAOYSA-N 4-(3-methoxyphenyl)aniline Chemical compound COC1=CC=CC(C=2C=CC(N)=CC=2)=C1 OSWFIVFLDKOXQC-UHFFFAOYSA-N 0.000 description 2
- 241000894006 Bacteria Species 0.000 description 2
- PVNIIMVLHYAWGP-UHFFFAOYSA-N Niacin Chemical compound OC(=O)C1=CC=CN=C1 PVNIIMVLHYAWGP-UHFFFAOYSA-N 0.000 description 2
- 240000007594 Oryza sativa Species 0.000 description 2
- 235000007164 Oryza sativa Nutrition 0.000 description 2
- 235000012538 ammonium bicarbonate Nutrition 0.000 description 2
- 238000002790 cross-validation Methods 0.000 description 2
- 238000001514 detection method Methods 0.000 description 2
- 238000010494 dissociation reaction Methods 0.000 description 2
- 230000005593 dissociations Effects 0.000 description 2
- 238000000605 extraction Methods 0.000 description 2
- 235000019253 formic acid Nutrition 0.000 description 2
- 239000007788 liquid Substances 0.000 description 2
- 238000004949 mass spectrometry Methods 0.000 description 2
- 239000002207 metabolite Substances 0.000 description 2
- 238000002360 preparation method Methods 0.000 description 2
- 235000009566 rice Nutrition 0.000 description 2
- 238000002470 solid-phase micro-extraction Methods 0.000 description 2
- 238000012549 training Methods 0.000 description 2
- 238000004704 ultra performance liquid chromatography Methods 0.000 description 2
- 239000003643 water by type Substances 0.000 description 2
- YPZUZOLGGMJZJO-XRGAULLZSA-N (3as,5as,9as,9br)-3a,6,6,9a-tetramethyl-2,4,5,5a,7,8,9,9b-octahydro-1h-benzo[e][1]benzofuran Chemical compound CC([C@@H]1CC2)(C)CCC[C@]1(C)[C@@H]1[C@@]2(C)OCC1 YPZUZOLGGMJZJO-XRGAULLZSA-N 0.000 description 1
- TUSDEZXZIZRFGC-UHFFFAOYSA-N 1-O-galloyl-3,6-(R)-HHDP-beta-D-glucose Natural products OC1C(O2)COC(=O)C3=CC(O)=C(O)C(O)=C3C3=C(O)C(O)=C(O)C=C3C(=O)OC1C(O)C2OC(=O)C1=CC(O)=C(O)C(O)=C1 TUSDEZXZIZRFGC-UHFFFAOYSA-N 0.000 description 1
- 241000228197 Aspergillus flavus Species 0.000 description 1
- 229920002786 Corilagin Polymers 0.000 description 1
- 102100024452 DNA-directed RNA polymerase III subunit RPC1 Human genes 0.000 description 1
- 206010013911 Dysgeusia Diseases 0.000 description 1
- 102000004190 Enzymes Human genes 0.000 description 1
- 108090000790 Enzymes Proteins 0.000 description 1
- 101000689002 Homo sapiens DNA-directed RNA polymerase III subunit RPC1 Proteins 0.000 description 1
- FBOZXECLQNJBKD-ZDUSSCGKSA-N L-methotrexate Chemical compound C=1N=C2N=C(N)N=C(N)C2=NC=1CN(C)C1=CC=C(C(=O)N[C@@H](CCC(O)=O)C(O)=O)C=C1 FBOZXECLQNJBKD-ZDUSSCGKSA-N 0.000 description 1
- OUYCCCASQSFEME-QMMMGPOBSA-N L-tyrosine Chemical compound OC(=O)[C@@H](N)CC1=CC=C(O)C=C1 OUYCCCASQSFEME-QMMMGPOBSA-N 0.000 description 1
- RFMMMVDNIPUKGG-YFKPBYRVSA-N N-acetyl-L-glutamic acid Chemical compound CC(=O)N[C@H](C(O)=O)CCC(O)=O RFMMMVDNIPUKGG-YFKPBYRVSA-N 0.000 description 1
- 241000209140 Triticum Species 0.000 description 1
- 235000021307 Triticum Nutrition 0.000 description 1
- 150000001413 amino acids Chemical class 0.000 description 1
- 239000011324 bead Substances 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000009835 boiling Methods 0.000 description 1
- 238000013124 brewing process Methods 0.000 description 1
- 238000012790 confirmation Methods 0.000 description 1
- 238000007796 conventional method Methods 0.000 description 1
- TUSDEZXZIZRFGC-XIGLUPEJSA-N corilagin Chemical compound O([C@H]1[C@H](O)[C@H]2OC(=O)C3=CC(O)=C(O)C(O)=C3C3=C(O)C(O)=C(O)C=C3C(=O)OC[C@@H](O1)[C@H]2O)C(=O)C1=CC(O)=C(O)C(O)=C1 TUSDEZXZIZRFGC-XIGLUPEJSA-N 0.000 description 1
- CPWYQGWOJMNXGJ-UHFFFAOYSA-N corilagin Natural products OC1C2COC(=O)c3c(O)c(O)c(O)c(O)c3c4c(O)c(O)c(O)c(O)c4C(=O)OC1C(O)C(OC(=O)c5cc(O)c(O)c(O)c5)O2 CPWYQGWOJMNXGJ-UHFFFAOYSA-N 0.000 description 1
- 238000007405 data analysis Methods 0.000 description 1
- 230000003247 decreasing effect Effects 0.000 description 1
- 238000000113 differential scanning calorimetry Methods 0.000 description 1
- 201000010099 disease Diseases 0.000 description 1
- 208000037265 diseases, disorders, signs and symptoms Diseases 0.000 description 1
- 229940079593 drug Drugs 0.000 description 1
- 239000003814 drug Substances 0.000 description 1
- 238000013210 evaluation model Methods 0.000 description 1
- 230000007717 exclusion Effects 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 238000013467 fragmentation Methods 0.000 description 1
- 238000006062 fragmentation reaction Methods 0.000 description 1
- 239000003205 fragrance Substances 0.000 description 1
- 125000003630 glycyl group Chemical group [H]N([H])C([H])([H])C(*)=O 0.000 description 1
- 238000000227 grinding Methods 0.000 description 1
- 238000012165 high-throughput sequencing Methods 0.000 description 1
- 230000002045 lasting effect Effects 0.000 description 1
- 238000004811 liquid chromatography Methods 0.000 description 1
- 230000014759 maintenance of location Effects 0.000 description 1
- 235000020068 maotai Nutrition 0.000 description 1
- 239000003550 marker Substances 0.000 description 1
- 229960003512 nicotinic acid Drugs 0.000 description 1
- 235000001968 nicotinic acid Nutrition 0.000 description 1
- 239000011664 nicotinic acid Substances 0.000 description 1
- 235000015097 nutrients Nutrition 0.000 description 1
- 150000007524 organic acids Chemical class 0.000 description 1
- 235000005985 organic acids Nutrition 0.000 description 1
- WLJVXDMOQOGPHL-UHFFFAOYSA-N phenylacetic acid Chemical compound OC(=O)CC1=CC=CC=C1 WLJVXDMOQOGPHL-UHFFFAOYSA-N 0.000 description 1
- 239000000843 powder Substances 0.000 description 1
- 239000000047 product Substances 0.000 description 1
- 238000013441 quality evaluation Methods 0.000 description 1
- 238000011002 quantification Methods 0.000 description 1
- 239000002994 raw material Substances 0.000 description 1
- 230000035945 sensitivity Effects 0.000 description 1
- 238000012163 sequencing technique Methods 0.000 description 1
- 239000000243 solution Substances 0.000 description 1
- 239000007921 spray Substances 0.000 description 1
- 238000007619 statistical method Methods 0.000 description 1
- 239000000758 substrate Substances 0.000 description 1
- 239000006228 supernatant Substances 0.000 description 1
- OUYCCCASQSFEME-UHFFFAOYSA-N tyrosine Natural products OC(=O)C(N)CC1=CC=C(O)C=C1 OUYCCCASQSFEME-UHFFFAOYSA-N 0.000 description 1
- 238000010200 validation analysis Methods 0.000 description 1
- 238000012795 verification Methods 0.000 description 1
- 235000013343 vitamin Nutrition 0.000 description 1
- 239000011782 vitamin Substances 0.000 description 1
- 229940088594 vitamin Drugs 0.000 description 1
- 229930003231 vitamin Natural products 0.000 description 1
- 150000003722 vitamin derivatives Chemical class 0.000 description 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N30/00—Investigating or analysing materials by separation into components using adsorption, absorption or similar phenomena or using ion-exchange, e.g. chromatography or field flow fractionation
- G01N30/02—Column chromatography
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N30/00—Investigating or analysing materials by separation into components using adsorption, absorption or similar phenomena or using ion-exchange, e.g. chromatography or field flow fractionation
- G01N30/02—Column chromatography
- G01N30/62—Detectors specially adapted therefor
- G01N30/72—Mass spectrometers
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N30/00—Investigating or analysing materials by separation into components using adsorption, absorption or similar phenomena or using ion-exchange, e.g. chromatography or field flow fractionation
- G01N30/02—Column chromatography
- G01N30/86—Signal analysis
- G01N30/8658—Optimising operation parameters
Landscapes
- Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Chemical & Material Sciences (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- General Physics & Mathematics (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Measuring Or Testing Involving Enzymes Or Micro-Organisms (AREA)
Abstract
本发明公开了一种基于大曲非挥发性物质鉴别大曲质量的方法，采用色谱分离技术获取大曲样品中的非挥发性物质，对大曲非挥发性物质变量指标进行筛选与优化，采用随机森林算法中变量重要性度量对大曲非挥发性物质变量进行优化筛选，获取对大曲质量分类贡献大的物质组成变量，构建大曲质量等级鉴别的判别模型，确定判别准则，对建模样品进行回判，以验证判别效果，将待测大曲样品按照获取的大曲非挥发性物质信息并输入上述构建的判别模型中，输出待测大曲样品的质量分类。本发明通过提取大曲中的非挥发性物质后，通过对非挥发性物质的分析直接判断大曲类型，实现了大曲的标准化鉴定，具有准确性好、稳定高效等优势。
Description
技术领域
本发明属于大曲生产技术领域，特别涉及一种基于大曲非挥发性物质鉴别大曲质量的方法。
背景技术
酱香型白酒是我国主要香型白酒之—，具有酒体醇厚、酱香突出、空杯留香持久和回味悠长等特点，这种独特的风味与其特殊且充满古风遗韵的酿造工艺息息相关。其中，以小麦为原料制成的酱香型髙温大曲含有多种菌种和酶，在酱香型白酒生产过程中起着菌种剂、糖化、发酵、生香以及提倛营养底物的功能，直接关系着酒的产量、风味和质量，堪称白酒的灵魂。在酱香型大曲的制备过程中，由于大曲存储仓库不同位置温度、湿度及环境中气体成分的差异，导致不同发酵程度大曲的产生。在实际生产过程中通常通过颜色将不同发酵程度的大曲分为白曲、黄曲和黑曲。然而，传统的通过目视判断大曲颜色的方法过于依赖生产人员的经验，亟需开发一种数字化的标准对不同发酵程度的大曲进行区分与鉴定。
现有的大曲鉴别方法有根据大曲优势微生物、挥发性风味物质和热焓值等指标，并结合判别分析方法进行数据分析实现大曲质量的判别。如在专利号为201310201060.5、专利名称为“一种构建判别模型鉴别大曲质量的方法”的中国发明专利提供了一种应用固相微萃取和判别分析鉴别大曲质量，该方法是通过大曲风味物质对大曲质量进行评价，但是这种鉴别方法时间长，成本高，且在储存过程大曲中的风味物质易挥发，导致无法准确鉴别。在专利号为201410419220.8、专利名称为“一种构建判别模型鉴别大曲质量的方法”的中国发明专利提供了一种应用454高通量测序技术和多元统计分析鉴别大曲质量，但这种鉴别方法测序时间长，且无法对大曲微生物中的活菌和死菌进行区分，实验操作繁琐；并且采用这种方法不能反映微生物的绝对定量情况。在专利号为CN202010518502.9、专利名称为“一种鉴别大曲质量的方法”的中国发明专利提供了一种差示扫描量热分析技术用于大曲样品分析，该方法是通过热焓值来对大曲质量进行判别，而热焓值是表征大曲系统能量的一个状态函数，是用于大曲性能的间接表征值，因此，难以及时、快速反映整个大曲的质量水平，对于大曲鉴别的判准率难以有效提高。
非挥发性物质是大曲发酵过程中形成的重要产物，对白酒发酵的风味有重要影响，例如白酒中常见的氨基酸类和有机酸类，因此在一定程度上能代表大曲的质量水平。这类化合物在大曲发酵过程由于其不易挥发的化学特性，在大曲发酵和储存过程中能稳定存在于大曲中，因此对于大曲质量水平的鉴别更具有代表性和稳定性。液相色谱-质谱联用仪(Liquid Chromatograph Mass Spectrometer)，简称LC-MS，结合了液相色谱仪有效分离热不稳性及高沸点化合物的分离能力与质谱仪很强的组分鉴定能力，是一种分离分析复杂有机混合物的有效手段，具有灵敏度高、适用范围广等特点，对于难挥发、强极性、热不稳定的代谢物分析具有极大优势。近年来，LC-MS作为代谢组学重要的研究手段，广泛应用于药物研究、临床疾病标志物研究等领域中。
发明内容
本发明意在提供一种基于大曲非挥发性物质鉴别大曲质量的方法，通过提取大曲中的非挥发性物质后，通过对非挥发性物质的分析直接判断大曲类型，实现了大曲的标准化鉴定，具有准确性好、稳定高效等优势。
本发明的方案:一种基于大曲非挥发性物质鉴别大曲质量的方法，包括以下步骤：
步骤(1)，大曲非挥发性物质的提取：采用色谱分离技术获取大曲样品中的非挥发性物质；
步骤(2)，变量筛选与优化：基于步骤(1)所得的大曲非挥发性物质与样品质量等级，对大曲非挥发性物质变量指标进行筛选与优化，采用随机森林算法中变量重要性度量对大曲非挥发性物质变量进行优化筛选，获取对大曲质量分类贡献大的物质组成变量；
步骤(3)，建立判别模型：基于步骤(2)所优化的变量信息与样品质量等级，构建大曲质量等级鉴别的判别模型，确定判别准则；
步骤(4)，验证模型：对建模样品进行回判，以验证判别效果；
步骤(5)，鉴别：将待测大曲样品按照步骤(1)获取的大曲非挥发性物质信息并输入上述构建的判别模型中，输出待测大曲样品的质量分类。
进一步，所述步骤(1)中所述分离技术为LC-MS分离技术。
进一步，所述步骤(1)中大曲样品为出仓黄曲、出仓白曲和出仓黑曲的任一种。
进一步，所述步骤(2)中对大曲非挥发性物质变量指标进行筛选与优化，采用变量重要性度量Mean Decrease Accuracy(平均精度下降)或Mean Decrease Gini(平均基尼下降)。
进一步，所述步骤(2)中优化筛选的原则为下降Mean Decrease Gini变量排序前20种非挥发性物质。
进一步，所述步骤(2)中变量筛选采用R语言构建。
进一步，所述步骤(3)中判别模型的构建及判别模型的构建采用随机森林法。
进一步，所述步骤(3)中判别模型采用R语言构建。
本发明的有益效果为：(1)利用LC-MS分离技术，解决了大曲挥发性风味物质不能全面准确反映大曲质量的优劣。(2)通过随机森林对大曲非挥发性物质进行筛选与优化，同时结合随机森林判别分析构建大曲质量鉴别评价模型并应用于大曲质量等级鉴别，构建了一套适用于出仓黄曲、出仓白曲和出仓黑曲等级鉴别的判别方法，该判别模型一经构建完成，以后在用该模型判别大曲质量时，可以随时调用，无需再进行重复构建，具有操作使用简单、判别准确率高等优点。
附图说明
图1为本发明实施例2重要性排序排名前20的大曲非挥发性物质变量图示。
具体实施方式
下面通过具体实施方式进一步详细说明：
实施例1：一种基于大曲非挥发性物质鉴别大曲质量的方法，包括以下步骤：
步骤(1)，大曲非挥发性物质的提取：采用色谱分离技术获取大曲样品中的非挥发性物质，具体的，色谱分离技术为LC-MS分离技术；大曲样品为拆仓曲，即出仓黄曲、出仓白曲和出仓黑曲的任一种。
步骤(2)，变量筛选与优化：基于步骤(1)所得的大曲非挥发性物质与样品质量等级，对大曲非挥发性物质变量指标进行筛选与优化，采用随机森林算法中变量重要性度量对大曲非挥发性物质变量进行优化筛选，获取对大曲质量分类贡献大的物质组成变量；本实施例中，对大曲非挥发性物质变量指标进行筛选与优化，采用变量重要性度量MeanDecrease Accuracy(平均精度下降)或Mean Decrease Gini(平均基尼下降)，优化筛选的原则为下降Mean Decrease Gini变量排序前20种非挥发性物质，具体为；步骤(2)中，优化筛选得到的大曲非挥发性物质组成变量为Indolelactic.acid(C169)、N.Acetylglutamic.acid(C184)、L.3.Phenyllactic.acid(C209)、Cantharidin(C185)、Sebacic.acid(C145)、Corilagin(C183)、Nicotinic.acid(C161)、Octadecanedioic.acid(C93)、Glycyl.L.leucine(C186)、Isoferulic.acid(C181)、o.Tyrosine(C150)、Ambroxane(C166)、Pentadecanoic.acid(C194)、X2.Furoic.acid(C192)、X3.Indolepropionic.acid(C170)、Vitamin.B5(C86)、Clotrimazole(C126)、Aspartic.acid(C177)、Glutamic.acid(C202)和Isoliensinine(C160)。
步骤(3)，建立判别模型：基于步骤(2)所优化的变量信息与样品质量等级，采用R语言通过随机森林法构建大曲质量等级鉴别的判别模型，确定判别准则；
步骤(4)，验证模型：对建模样品进行回判，以验证判别效果；
步骤(5)，鉴别：将待测大曲样品按照步骤(1)获取的大曲非挥发性物质信息并输入上述构建的判别模型中，输出待测大曲样品的质量分类。
实施例2：实验过程
1.样品信息
本实施例中的样品来源于贵州茅台酒股份有限公司提供的90个高温大曲，其中样品包括：
建模样品：出仓黄曲30个、出仓白曲30个、出仓黑曲30个；
测试样品：出仓黄曲15个、出仓白曲15个、出仓黑曲15个。
2.大曲非挥发性物质的提取
称取20mg大曲细粉，加入800μL甲醇提取剂珠打研磨2分钟，取300μL上清夜冻干后，用100μL 30％甲醇/水复溶。
LC-MS分析：在正电离模式下的色谱柱为Waters ACQUITY UPLC BEH C8(1.7μm,2.1×50mm)，柱温为60℃，流动相流速为0.4mL/min。流动相A为含有0.1％甲酸(体积比)的水溶液，流动相B为含有0.1％甲酸(体积比)的乙腈溶液。流动相初始比例是5％ B相，维持0.5分钟。随后在1.5分钟内线性至40％B相，在6分钟内线性升至100％B相。在0.1min内降至5％ B相，并维持1.9分钟。升总梯度时间为12分钟。
在负电离模式下的色谱柱为Waters ACQUITY UPLC HSS T3(1.8μm,2.1×50mm)，柱温为55℃，流动相流速为0.35mL/min。流动相A为含有6.5mM NH4HCO3的水溶液，流动相B为含有6.5mM NH4HCO3的95％甲醇/水溶液(体积比)。流动相初始比例是2％ B相，维持0.5分钟。随后在1.5分钟内线性升至40％B相，在6分钟内线性升至100％B相。在0.1min内降至2％ B相，并维持1.9分钟。总梯度时间为12分钟。
质谱参数为离子源gas 1，离子源gas 2，Curtain gas分别为50psi，50psi，35psi。离子源温度为500℃，CAD gas为7，Spray voltage在正负离子模式下分别为5500V和-4500V。一级质谱扫描范围是80–1250Da，去簇电压(DP)在正负离子模式下分别为80V和-80V，DP spread为0V。一级质谱中，碰撞能(CE)在正负离子模式下分别为10V和-10V，CEspread为0V。富集时间0.1s，开启动态背景排除。二级质谱采集模式选择数据依赖的采集方式(IDA)，解离模式选择诱导碰撞解离(CID)。Q1分辨率为单位分辨率，二级质谱扫描范围为50–1250Da。去簇电压(DP)在正负离子模式下分别为80V和-80V，DP spread为0V。碰撞能设置低，中，高三种，分别对应于15V(-15V)，30V(-30V)，45V(-45V)。富集时间为0.04s。
基于MetExDB数据库，利用MetEx定性方法对所得特征峰进行分析，经精确质量，色谱保留规律和质谱碎裂规律共同确证，从不同拆仓曲中鉴定出221种代谢物(物质编号为C1，C2，C3，…，C221)与相对含量。
3.变量筛选与优化
将变量与样品信息输入，利用R语言随机森林算法，采用变量重要性度量MeanDecrease Gini对获得的大曲非挥发性物质变量进行筛选与优化，获得重要性排序前20种非挥发性物质，如图1所示。
4.建立判别模型
通过R语言，采用随机森林分析法，以上述20类大曲非挥发性物质为自变量，以3种高温大曲的判断为因变量，进行随机森林判别分析，构建大曲质量等级的判别模型。
5.验证模型
首先利用交叉验证对建模样品进行回判，以验证判别效果，结果见表1。
表1建模样品交叉验证结果
注：Y：出仓黄曲；W：出仓白曲；B：出仓黑曲
从表1可以看出，用于建模的样品共90个：出仓黄曲30个、出仓白曲30个、出仓黑曲30个。该判别模型对出仓黄曲的判别准确率为29/30＝96.7％，对出仓白曲的判别准确率为30/30＝100％，对出仓黑曲的判别准确率为30/30＝100％，综合判别准确率为(29+30+30)/(30+30+30)＝98.9％，表明所建立的判别模型对大曲质量鉴别的判别较为有效。
6.鉴别
采用判别模型对45个外部的测试样品(出仓黄曲15个、出仓白曲15个、出仓黑曲15个)进行判别，结果见表2。
表2外部测试样品验证结果
注：Y：出仓黄曲；W：出仓白曲；B：出仓黑曲
从表2可以看出，该二次判别对外部测试样品中出仓黄曲的判别准确率为14/15＝93.3％，对出仓白曲的判别准确率为15/15＝100％，对出仓黑曲的判别准确率为15/15＝100.0％，综合判别准确率为(14+15+15)/(15+15+15)＝97.8％，该结果进一步表明所建立的判别模型对大曲质量鉴别的判别效果较好。
显然，相对于固相微萃取的方式，本申请提供的方法所用时间极少，2小时之内即可完成，而且成本极低，除步骤(1)中LC-MS检测需要用到少许耗材外，整个过程基本无需其他耗材，而该方法的可靠性和准确性方面也远远超过通过微生物和挥发性风味物质建立的质量评价方法。更为重要的是R语言中提供了很多用于构建判别模型的程序包或函数，因此采用R语言构建判别模型可保证非专业人士通过较少的培训也可轻易操作，极大的避免了操作人员不够专业所带来的误判，也避免了企业为培养专业的操作人员所付出的高昂成本。
以上所述的仅是本发明的实施例，方案中公知的特性等常识在此未作过多描述。本发明要求的保护范围应当以其权利要求的内容为准，说明书中的具体实施方式等记载可以用于解释权利要求的内容。
Claims (8)
1.一种基于大曲非挥发性物质鉴别大曲质量的方法，其特征在于：包括以下步骤：
步骤(1)，大曲非挥发性物质的提取：采用色谱分离技术获取大曲样品中的非挥发性物质；
步骤(2)，变量筛选与优化：基于步骤(1)所得的大曲非挥发性物质与样品质量等级，对大曲非挥发性物质变量指标进行筛选与优化，采用随机森林算法中变量重要性度量对大曲非挥发性物质变量进行优化筛选，获取对大曲质量分类贡献大的物质组成变量；
步骤(3)，建立判别模型：基于步骤(2)所优化的变量信息与样品质量等级，构建大曲质量等级鉴别的判别模型，确定判别准则；
步骤(4)，验证模型：对建模样品进行回判，以验证判别效果；
步骤(5)，鉴别：将待测大曲样品按照步骤(1)获取的大曲非挥发性物质信息并输入上述构建的判别模型中，输出待测大曲样品的质量分类。
2.根据权利要求1所述的一种基于大曲非挥发性物质鉴别大曲质量的方法，其特征在于：所述步骤(1)中的分离技术为LC-MS分离技术。
3.根据权利要求2所述的一种基于大曲非挥发性物质鉴别大曲质量的方法，其特征在于：所述步骤(1)中大曲样品为出仓黄曲、出仓白曲和出仓黑曲的任一种。
4.根据权利要求1所述的一种基于大曲非挥发性物质鉴别大曲质量的方法，其特征在于：所述步骤(2)中对大曲非挥发性物质变量指标进行筛选与优化，采用变量重要性度量Mean Decrease Accuracy或Mean Decrease Gini。
5.根据权利要求4所述的一种基于大曲非挥发性物质鉴别大曲质量的方法，其特征在于：所述步骤(2)中优化筛选的原则为下降Mean Decrease Gini变量排序前20种非挥发性物质。
6.根据权利要求5所述的一种基于大曲非挥发性物质鉴别大曲质量的方法，其特征在于：所述步骤(2)中变量筛选采用R语言构建。
7.根据权利要求1所述的一种基于大曲非挥发性物质鉴别大曲质量的方法，其特征在于：所述步骤(3)中判别模型的构建及判别模型的构建采用随机森林法。
8.根据权利要求7所述的一种基于大曲非挥发性物质鉴别大曲质量的方法，其特征在于：所述步骤(3)中判别模型采用R语言构建。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202211287619.6A CN115436531B (zh) | 2022-10-20 | 2022-10-20 | 一种基于大曲非挥发性物质鉴别大曲质量的方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202211287619.6A CN115436531B (zh) | 2022-10-20 | 2022-10-20 | 一种基于大曲非挥发性物质鉴别大曲质量的方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN115436531A true CN115436531A (zh) | 2022-12-06 |
| CN115436531B CN115436531B (zh) | 2024-06-25 |
Family
ID=84252782
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202211287619.6A Active CN115436531B (zh) | 2022-10-20 | 2022-10-20 | 一种基于大曲非挥发性物质鉴别大曲质量的方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN115436531B (zh) |
Citations (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20030139885A1 (en) * | 2001-11-05 | 2003-07-24 | Irm, Llc | Methods and devices for proteomics data complexity reduction |
| CN103571728A (zh) * | 2013-08-16 | 2014-02-12 | 江南大学 | 脂肽类化合物用于改善白酒香气的用途及其测定方法 |
| CN104372075A (zh) * | 2014-08-22 | 2015-02-25 | 贵州茅台酒股份有限公司 | 一种构建判别模型鉴别大曲质量的方法 |
| CN105861346A (zh) * | 2016-06-14 | 2016-08-17 | 江南大学 | 一株低产尿素、产风味的异常威克汉姆酵母菌株及其在食品生产中的应用 |
| CN109828044A (zh) * | 2019-02-21 | 2019-05-31 | 安徽古井贡酒股份有限公司 | 一种超高效合相色谱串接QDa同时快速检测酒类产品中8种酚酸类物质的方法 |
| CN109949863A (zh) * | 2019-02-18 | 2019-06-28 | 贵州茅台酒股份有限公司 | 一种基于随机森林模型鉴别大曲质量的方法 |
| CN113219101A (zh) * | 2021-05-20 | 2021-08-06 | 贵州茅台酒股份有限公司 | 一种葡萄酒中外源性风险物质的检测与筛查方法及其在检测中的用途 |
| CN113533549A (zh) * | 2021-01-28 | 2021-10-22 | 岛津企业管理（中国）有限公司 | 白酒口味物质鉴定分析系统 |
-
2022
- 2022-10-20 CN CN202211287619.6A patent/CN115436531B/zh active Active
Patent Citations (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20030139885A1 (en) * | 2001-11-05 | 2003-07-24 | Irm, Llc | Methods and devices for proteomics data complexity reduction |
| CN103571728A (zh) * | 2013-08-16 | 2014-02-12 | 江南大学 | 脂肽类化合物用于改善白酒香气的用途及其测定方法 |
| CN104372075A (zh) * | 2014-08-22 | 2015-02-25 | 贵州茅台酒股份有限公司 | 一种构建判别模型鉴别大曲质量的方法 |
| CN105861346A (zh) * | 2016-06-14 | 2016-08-17 | 江南大学 | 一株低产尿素、产风味的异常威克汉姆酵母菌株及其在食品生产中的应用 |
| CN109949863A (zh) * | 2019-02-18 | 2019-06-28 | 贵州茅台酒股份有限公司 | 一种基于随机森林模型鉴别大曲质量的方法 |
| CN109828044A (zh) * | 2019-02-21 | 2019-05-31 | 安徽古井贡酒股份有限公司 | 一种超高效合相色谱串接QDa同时快速检测酒类产品中8种酚酸类物质的方法 |
| CN113533549A (zh) * | 2021-01-28 | 2021-10-22 | 岛津企业管理（中国）有限公司 | 白酒口味物质鉴定分析系统 |
| CN113219101A (zh) * | 2021-05-20 | 2021-08-06 | 贵州茅台酒股份有限公司 | 一种葡萄酒中外源性风险物质的检测与筛查方法及其在检测中的用途 |
Non-Patent Citations (2)
| Title |
|---|
| 刘露等: "白酒分析检测研究进展", 酿酒, vol. 45, no. 6, pages 28 - 32 * |
| 吴江，吴一平: "基于ggplot的政经数据可视化", 31 July 2020, 华夏出版社, pages: 219 - 222 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN115436531B (zh) | 2024-06-25 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| KR101965293B1 (ko) | 한국산과 중국산 백미 판별용 Biomarker 및 이를 이용한 판별 모델 | |
| Martínez-Jarquín et al. | Automated chemical fingerprinting of Mexican spirits derived from Agave (tequila and mezcal) using direct-injection electrospray ionisation (DIESI) and low-temperature plasma (LTP) mass spectrometry | |
| CN111272931A (zh) | 一种茶叶的原产地溯源方法 | |
| CN112578039B (zh) | 一种燕窝中掺假物的分析鉴定方法 | |
| CN108535375B (zh) | 一种基于液质联用代谢组学分析玛咖标志性代谢产物的方法 | |
| He et al. | Comprehensive two-dimensional gas chromatography–time of flight mass spectrometry (GC× GC-TOFMS) in conventional and reversed column configuration for the investigation of Baijiu aroma types and regional origin | |
| CN111044638A (zh) | 一种对不同品种烤烟烟叶分类鉴别的方法 | |
| CN111060642A (zh) | 一种对同品种不同产地烟叶分类鉴别的方法 | |
| Chen et al. | Classification of vinegar with different marked ages using olfactory sensors and gustatory sensors | |
| Mutz et al. | Feasibility of NIR spectroscopy coupled with chemometrics for classification of Brazilian specialty coffee | |
| CN111208251B (zh) | 一种白茶年份判别方法 | |
| JP2020511671A (ja) | 改善された精度、同定および定量化のためのｉｒｏａメタボロミクスワークフロー | |
| CN109214423B (zh) | 一种基于动静态数据融合的食品质量判别分析方法 | |
| Jin et al. | Tracing the origin of Taiping Houkui green tea using 1H NMR and HS-SPME-GC–MS chemical fingerprints, data fusion and chemometrics | |
| Lias et al. | A study on the application of electronic nose coupled with DFA and statistical analysis for evaluating the relationship between sample volumes versus sensor intensity of agarwood essential oils blending ratio | |
| CN111337605B (zh) | 一种评价荷花蜂花粉真实性的方法 | |
| CN113390980A (zh) | 一种煎饼加工中风味物质变化的评价方法 | |
| CN106442589A (zh) | 一种采用能量色散x荧光光谱指纹识别技术鉴别化妆品品质的方法 | |
| CN115436531A (zh) | 一种基于大曲非挥发性物质鉴别大曲质量的方法 | |
| CN113049722A (zh) | 一种挥发性香气物质指纹图谱的构建方法及其应用 | |
| CN108287200A (zh) | 质谱参照数据库的建立方法及基于其的物质分析方法 | |
| CN110221010B (zh) | 一种测定酱香型白酒大曲中吡喃酮类化合物的方法 | |
| CN109100461A (zh) | 一种利用蛋白组学技术区分有机大米和非有机大米的方法 | |
| CN113533587A (zh) | 基于气相离子迁移谱鉴别辣椒粉品种的方法 | |
| Bambina et al. | 1H-NMR Spectroscopy Coupled with Chemometrics to Classify Wines According to Different Grape Varieties and Different Terroirs |
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 |