CN109557014A - 一种快速检测发酵乳中乳酸菌数的方法 - Google Patents
一种快速检测发酵乳中乳酸菌数的方法 Download PDFInfo
- Publication number
- CN109557014A CN109557014A CN201910069041.9A CN201910069041A CN109557014A CN 109557014 A CN109557014 A CN 109557014A CN 201910069041 A CN201910069041 A CN 201910069041A CN 109557014 A CN109557014 A CN 109557014A
- Authority
- CN
- China
- Prior art keywords
- lactic acid
- acid bacteria
- sample
- acidified milk
- calibration
- 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
- JVTAAEKCZFNVCJ-UHFFFAOYSA-N lactic acid Chemical compound CC(O)C(O)=O JVTAAEKCZFNVCJ-UHFFFAOYSA-N 0.000 title claims abstract description 272
- 241000894006 Bacteria Species 0.000 title claims abstract description 141
- 239000004310 lactic acid Substances 0.000 title claims abstract description 136
- 235000014655 lactic acid Nutrition 0.000 title claims abstract description 136
- 238000000034 method Methods 0.000 title claims abstract description 119
- 235000020167 acidified milk Nutrition 0.000 title claims abstract description 95
- 238000001514 detection method Methods 0.000 title claims abstract description 46
- 238000002329 infrared spectrum Methods 0.000 claims abstract description 64
- 238000004458 analytical method Methods 0.000 claims abstract description 17
- 230000008676 import Effects 0.000 claims abstract description 5
- 241000186660 Lactobacillus Species 0.000 claims description 23
- 229940039696 lactobacillus Drugs 0.000 claims description 23
- 241000186000 Bifidobacterium Species 0.000 claims description 18
- 230000005540 biological transmission Effects 0.000 claims description 17
- 238000002790 cross-validation Methods 0.000 claims description 17
- 238000004611 spectroscopical analysis Methods 0.000 claims description 16
- 238000012360 testing method Methods 0.000 claims description 16
- 238000005070 sampling Methods 0.000 claims description 15
- 241000194020 Streptococcus thermophilus Species 0.000 claims description 13
- 238000012795 verification Methods 0.000 claims description 11
- 239000004698 Polyethylene Substances 0.000 claims description 7
- 238000012937 correction Methods 0.000 claims description 7
- 238000007789 sealing Methods 0.000 claims description 6
- 230000002159 abnormal effect Effects 0.000 claims description 4
- 229910000530 Gallium indium arsenide Inorganic materials 0.000 claims description 3
- 238000007689 inspection Methods 0.000 abstract description 7
- 238000003912 environmental pollution Methods 0.000 abstract description 4
- 235000021001 fermented dairy product Nutrition 0.000 abstract description 2
- 238000001228 spectrum Methods 0.000 description 32
- 230000003595 spectral effect Effects 0.000 description 16
- 244000005700 microbiome Species 0.000 description 15
- 238000010790 dilution Methods 0.000 description 14
- 239000012895 dilution Substances 0.000 description 14
- 230000008569 process Effects 0.000 description 11
- 239000001963 growth medium Substances 0.000 description 10
- 238000011156 evaluation Methods 0.000 description 9
- 235000013305 food Nutrition 0.000 description 9
- 238000005457 optimization Methods 0.000 description 8
- 238000005259 measurement Methods 0.000 description 7
- 238000012628 principal component regression Methods 0.000 description 7
- 238000010521 absorption reaction Methods 0.000 description 6
- 239000007788 liquid Substances 0.000 description 6
- 239000000463 material Substances 0.000 description 6
- -1 polyethylene Polymers 0.000 description 6
- 229920000573 polyethylene Polymers 0.000 description 6
- 238000004445 quantitative analysis Methods 0.000 description 6
- 230000001580 bacterial effect Effects 0.000 description 5
- 230000008901 benefit Effects 0.000 description 5
- 239000003153 chemical reaction reagent Substances 0.000 description 5
- 238000007865 diluting Methods 0.000 description 5
- 230000000694 effects Effects 0.000 description 5
- 239000000284 extract Substances 0.000 description 5
- 238000001914 filtration Methods 0.000 description 5
- 235000013336 milk Nutrition 0.000 description 5
- 239000008267 milk Substances 0.000 description 5
- 210000004080 milk Anatomy 0.000 description 5
- 239000000203 mixture Substances 0.000 description 5
- 230000005856 abnormality Effects 0.000 description 4
- 238000005516 engineering process Methods 0.000 description 4
- 238000009499 grossing Methods 0.000 description 4
- 238000012545 processing Methods 0.000 description 4
- 235000018102 proteins Nutrition 0.000 description 4
- 102000004169 proteins and genes Human genes 0.000 description 4
- 108090000623 proteins and genes Proteins 0.000 description 4
- FAPWRFPIFSIZLT-UHFFFAOYSA-M Sodium chloride Chemical compound [Na+].[Cl-] FAPWRFPIFSIZLT-UHFFFAOYSA-M 0.000 description 3
- 238000006243 chemical reaction Methods 0.000 description 3
- 238000010586 diagram Methods 0.000 description 3
- 235000015141 kefir Nutrition 0.000 description 3
- 238000012417 linear regression Methods 0.000 description 3
- 230000002906 microbiologic effect Effects 0.000 description 3
- 238000012544 monitoring process Methods 0.000 description 3
- 238000002360 preparation method Methods 0.000 description 3
- 238000012216 screening Methods 0.000 description 3
- 238000013207 serial dilution Methods 0.000 description 3
- 238000004566 IR spectroscopy Methods 0.000 description 2
- 239000000654 additive Substances 0.000 description 2
- 230000000996 additive effect Effects 0.000 description 2
- 238000004364 calculation method Methods 0.000 description 2
- 150000001720 carbohydrates Chemical class 0.000 description 2
- 230000008859 change Effects 0.000 description 2
- 239000006071 cream Substances 0.000 description 2
- 235000013399 edible fruits Nutrition 0.000 description 2
- 238000002474 experimental method Methods 0.000 description 2
- 238000000855 fermentation Methods 0.000 description 2
- 230000004151 fermentation Effects 0.000 description 2
- 239000000796 flavoring agent Substances 0.000 description 2
- 235000019634 flavors Nutrition 0.000 description 2
- 239000007789 gas Substances 0.000 description 2
- 238000011534 incubation Methods 0.000 description 2
- 150000002632 lipids Chemical class 0.000 description 2
- 229920002521 macromolecule Polymers 0.000 description 2
- 239000011159 matrix material Substances 0.000 description 2
- GJPZDZHEZDANAG-UHFFFAOYSA-N methyl n-(1h-benzimidazol-2-yl)carbamate;propan-2-yl n-(3,4-diethoxyphenyl)carbamate Chemical compound C1=CC=C2NC(NC(=O)OC)=NC2=C1.CCOC1=CC=C(NC(=O)OC(C)C)C=C1OCC GJPZDZHEZDANAG-UHFFFAOYSA-N 0.000 description 2
- 150000007523 nucleic acids Chemical class 0.000 description 2
- 102000039446 nucleic acids Human genes 0.000 description 2
- 108020004707 nucleic acids Proteins 0.000 description 2
- 238000002203 pretreatment Methods 0.000 description 2
- 238000003908 quality control method Methods 0.000 description 2
- 230000010076 replication Effects 0.000 description 2
- 239000011780 sodium chloride Substances 0.000 description 2
- 239000007787 solid Substances 0.000 description 2
- 238000007711 solidification Methods 0.000 description 2
- 230000008023 solidification Effects 0.000 description 2
- 239000000758 substrate Substances 0.000 description 2
- 238000010200 validation analysis Methods 0.000 description 2
- QIJRTFXNRTXDIP-UHFFFAOYSA-N (1-carboxy-2-sulfanylethyl)azanium;chloride;hydrate Chemical compound O.Cl.SCC(N)C(O)=O QIJRTFXNRTXDIP-UHFFFAOYSA-N 0.000 description 1
- GUBGYTABKSRVRQ-XLOQQCSPSA-N Alpha-Lactose Chemical compound O[C@@H]1[C@@H](O)[C@@H](O)[C@@H](CO)O[C@H]1O[C@@H]1[C@@H](CO)O[C@H](O)[C@H](O)[C@H]1O GUBGYTABKSRVRQ-XLOQQCSPSA-N 0.000 description 1
- OYPRJOBELJOOCE-UHFFFAOYSA-N Calcium Chemical compound [Ca] OYPRJOBELJOOCE-UHFFFAOYSA-N 0.000 description 1
- 241001478240 Coccus Species 0.000 description 1
- 241001269238 Data Species 0.000 description 1
- 241001468191 Lactobacillus kefiri Species 0.000 description 1
- GUBGYTABKSRVRQ-QKKXKWKRSA-N Lactose Natural products OC[C@H]1O[C@@H](O[C@H]2[C@H](O)[C@@H](O)C(O)O[C@@H]2CO)[C@H](O)[C@@H](O)[C@H]1O GUBGYTABKSRVRQ-QKKXKWKRSA-N 0.000 description 1
- 238000010222 PCR analysis Methods 0.000 description 1
- 241001052560 Thallis Species 0.000 description 1
- 238000002679 ablation Methods 0.000 description 1
- 238000002835 absorbance Methods 0.000 description 1
- 238000000862 absorption spectrum Methods 0.000 description 1
- 238000000149 argon plasma sintering Methods 0.000 description 1
- QVGXLLKOCUKJST-UHFFFAOYSA-N atomic oxygen Chemical compound [O] QVGXLLKOCUKJST-UHFFFAOYSA-N 0.000 description 1
- 238000003705 background correction Methods 0.000 description 1
- 238000010009 beating Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 239000011575 calcium Substances 0.000 description 1
- 229910052791 calcium Inorganic materials 0.000 description 1
- 235000001465 calcium Nutrition 0.000 description 1
- 210000004027 cell Anatomy 0.000 description 1
- 210000002421 cell wall Anatomy 0.000 description 1
- 230000001413 cellular effect Effects 0.000 description 1
- 235000013339 cereals Nutrition 0.000 description 1
- 238000002512 chemotherapy Methods 0.000 description 1
- HVYWMOMLDIMFJA-DPAQBDIFSA-N cholesterol Chemical compound C1C=C2C[C@@H](O)CC[C@]2(C)[C@@H]2[C@@H]1[C@@H]1CC[C@H]([C@H](C)CCCC(C)C)[C@@]1(C)CC2 HVYWMOMLDIMFJA-DPAQBDIFSA-N 0.000 description 1
- 150000001875 compounds Chemical class 0.000 description 1
- 235000020247 cow milk Nutrition 0.000 description 1
- 235000014048 cultured milk product Nutrition 0.000 description 1
- 229960001305 cysteine hydrochloride Drugs 0.000 description 1
- 238000007405 data analysis Methods 0.000 description 1
- 230000004069 differentiation Effects 0.000 description 1
- 238000009792 diffusion process Methods 0.000 description 1
- 230000029087 digestion Effects 0.000 description 1
- 238000009826 distribution Methods 0.000 description 1
- 238000000605 extraction Methods 0.000 description 1
- 239000003925 fat Substances 0.000 description 1
- 235000013373 food additive Nutrition 0.000 description 1
- 239000002778 food additive Substances 0.000 description 1
- 230000037406 food intake Effects 0.000 description 1
- 235000012055 fruits and vegetables Nutrition 0.000 description 1
- 125000000524 functional group Chemical group 0.000 description 1
- 150000004676 glycans Chemical class 0.000 description 1
- 230000036541 health Effects 0.000 description 1
- 230000036039 immunity Effects 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 208000015181 infectious disease Diseases 0.000 description 1
- 238000002129 infrared reflectance spectroscopy Methods 0.000 description 1
- 238000002347 injection Methods 0.000 description 1
- 239000007924 injection Substances 0.000 description 1
- 238000009434 installation Methods 0.000 description 1
- 230000000968 intestinal effect Effects 0.000 description 1
- 238000002955 isolation Methods 0.000 description 1
- 239000008101 lactose Substances 0.000 description 1
- XFIKMPRBSORKQP-JATHGWPISA-M lithium;9-[(e)-4-[(2s,3r,4r,5s)-3,4-dihydroxy-5-[[(2s,3s)-3-[(2s,3s)-3-hydroxybutan-2-yl]oxiran-2-yl]methyl]oxan-2-yl]-3-methylbut-2-enoyl]oxynonanoate Chemical class [Li+].C[C@H](O)[C@H](C)[C@@H]1O[C@H]1C[C@@H]1[C@@H](O)[C@@H](O)[C@H](C\C(C)=C\C(=O)OCCCCCCCCC([O-])=O)OC1 XFIKMPRBSORKQP-JATHGWPISA-M 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 230000013011 mating Effects 0.000 description 1
- 239000012528 membrane Substances 0.000 description 1
- 238000000491 multivariate analysis Methods 0.000 description 1
- 238000001320 near-infrared absorption spectroscopy Methods 0.000 description 1
- 235000016709 nutrition Nutrition 0.000 description 1
- 230000035764 nutrition Effects 0.000 description 1
- 239000001301 oxygen Substances 0.000 description 1
- 229910052760 oxygen Inorganic materials 0.000 description 1
- 238000011056 performance test Methods 0.000 description 1
- 229920001282 polysaccharide Polymers 0.000 description 1
- 239000005017 polysaccharide Substances 0.000 description 1
- 230000008092 positive effect Effects 0.000 description 1
- 239000000843 powder Substances 0.000 description 1
- 238000007781 pre-processing Methods 0.000 description 1
- 238000000746 purification Methods 0.000 description 1
- 238000004451 qualitative analysis Methods 0.000 description 1
- 239000002994 raw material Substances 0.000 description 1
- 238000009738 saturating Methods 0.000 description 1
- 230000035945 sensitivity Effects 0.000 description 1
- 238000000926 separation method Methods 0.000 description 1
- 238000010183 spectrum analysis Methods 0.000 description 1
- 229910001220 stainless steel Inorganic materials 0.000 description 1
- 239000010935 stainless steel Substances 0.000 description 1
- 230000001954 sterilising effect Effects 0.000 description 1
- 239000000126 substance Substances 0.000 description 1
- 238000010998 test method Methods 0.000 description 1
- 238000012549 training Methods 0.000 description 1
- 238000002235 transmission spectroscopy Methods 0.000 description 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N15/00—Investigating characteristics of particles; Investigating permeability, pore-volume or surface-area of porous materials
- G01N15/10—Investigating individual particles
-
- C—CHEMISTRY; METALLURGY
- C12—BIOCHEMISTRY; BEER; SPIRITS; WINE; VINEGAR; MICROBIOLOGY; ENZYMOLOGY; MUTATION OR GENETIC ENGINEERING
- C12Q—MEASURING OR TESTING PROCESSES INVOLVING ENZYMES, NUCLEIC ACIDS OR MICROORGANISMS; COMPOSITIONS OR TEST PAPERS THEREFOR; PROCESSES OF PREPARING SUCH COMPOSITIONS; CONDITION-RESPONSIVE CONTROL IN MICROBIOLOGICAL OR ENZYMOLOGICAL PROCESSES
- C12Q1/00—Measuring or testing processes involving enzymes, nucleic acids or microorganisms; Compositions therefor; Processes of preparing such compositions
- C12Q1/02—Measuring or testing processes involving enzymes, nucleic acids or microorganisms; Compositions therefor; Processes of preparing such compositions involving viable microorganisms
- C12Q1/04—Determining presence or kind of microorganism; Use of selective media for testing antibiotics or bacteriocides; Compositions containing a chemical indicator therefor
- C12Q1/06—Quantitative determination
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N15/00—Investigating characteristics of particles; Investigating permeability, pore-volume or surface-area of porous materials
- G01N15/10—Investigating individual particles
- G01N2015/1006—Investigating individual particles for cytology
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N15/00—Investigating characteristics of particles; Investigating permeability, pore-volume or surface-area of porous materials
- G01N15/10—Investigating individual particles
- G01N2015/1022—Measurement of deformation of individual particles by non-optical means
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N15/00—Investigating characteristics of particles; Investigating permeability, pore-volume or surface-area of porous materials
- G01N15/10—Investigating individual particles
- G01N2015/1024—Counting particles by non-optical means
Landscapes
- Chemical & Material Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Health & Medical Sciences (AREA)
- Organic Chemistry (AREA)
- Biochemistry (AREA)
- Proteomics, Peptides & Aminoacids (AREA)
- Physics & Mathematics (AREA)
- General Health & Medical Sciences (AREA)
- Zoology (AREA)
- Immunology (AREA)
- Analytical Chemistry (AREA)
- Wood Science & Technology (AREA)
- Engineering & Computer Science (AREA)
- Pathology (AREA)
- Biophysics (AREA)
- Biotechnology (AREA)
- Microbiology (AREA)
- Molecular Biology (AREA)
- Dispersion Chemistry (AREA)
- Toxicology (AREA)
- General Physics & Mathematics (AREA)
- Bioinformatics & Cheminformatics (AREA)
- General Engineering & Computer Science (AREA)
- Genetics & Genomics (AREA)
- Investigating Or Analysing Materials By Optical Means (AREA)
Abstract
本发明公开了一种快速检测发酵乳中乳酸菌数的方法,含以下步骤:(1)采集发酵乳作为定标集样品,采用平板计数法对定标集样品进行分析检测,得乳酸菌数实测值,并转换成对数值,建立基础数据库;同步对定标集样品进行近红外光谱扫描,录制其近红外光谱数据,建立近红外光谱数据库;将基础数据库与近红外光谱数据库进行对应,建立乳酸菌数定标模型;(2)对定标模型进行验证;(3)将待测发酵乳样品进行近红外光谱扫描,得其近红外光谱数据,将数据导入定标模型,经换算得待测发酵乳样品的乳酸菌数预测值。该方法操作简捷,快速、无损、高效、准确、成本低,不污染环境,可多组分同时检测,能满足发酵乳产品在线质量检验的高效性和及时性要求。
Description
技术领域
本发明属于食品质量与安全检测,具体涉及一种快速检测发酵乳中乳酸菌数的方法,更涉及一种基于近红外光谱技术实现发酵乳中乳酸菌数快速检测的方法。
背景技术
乳酸菌是一类能在可利用的碳水化合物发酵过程中产生大量乳酸的细菌,作为益生菌,被广泛地应用于发酵制品和生物制品领域,具有调节肠道菌群的平衡、抑制肠道致病菌的感染、防治肿瘤、降低血清胆固醇水平、促进消化、提高机体免疫功能等健康功效。乳酸菌要发挥其保健功能,必须满足三个条件:其一,乳酸菌必须是活的;其二,摄入足够的量,人体必须每日摄入1~100亿或以上活性乳酸菌才会产生积极的功效;其三,乳酸菌的功效和益处必须是经过临床验证的。发酵乳是以新鲜牛乳为原料,经巴氏杀菌后添加乳酸菌发酵而成,比普通牛乳更营养,蛋白质和钙更易消化吸收。发酵乳被认为是乳酸菌传递到人体内的最理想载体。
发酵乳中乳酸菌数量是产品发挥益生功效的重要条件,国标GB 19302-2010《食品安全国家标准发酵乳》中规定,乳酸菌数不得低于106CFU/g(mL)。另外,在产品的安全性和功能性试验中也常以乳酸菌数量为评价指标。因此,乳酸菌数检测对企业质量控制和日常监督至关重要。
目前发酵乳中乳酸菌数的常规检测方法为平板计数法,其相关标准GB 4789.35-2016《食品安全国家标准食品微生物学检验乳酸菌检验》在2017年6月23日实施,新国标对乳酸菌总数的检测进行了较大改动。培养基由单一的MRS培养基增加MRS、MC、莫匹罗星锂盐和半胱氨酸盐酸盐改良MRS三种培养基组合,培养条件由单一的厌氧培养增加为厌氧和需氧组合,培养时间由48h增加为72h。
该方法存在以下缺点:
(1)操作繁琐:包含样品制备、均质、6~8次10倍系列稀释、培养、计数、计算等步骤,对检测人员技术要求高;
(2)效率低,耗时长:培养时间长达72h,容易造成监控数据滞后现象,无法满足企业在线检测要求;
(3)检测成本高:无菌耗材消耗大,仪器设备要求高,需配套洁净室、厌氧箱等设备。
同时,由于发酵乳保质期短,只有21-28d,因此研发更加快速、准确的检测技术,对于产品的质量控制和在线监控具有重要意义。
近红外光谱方法(Near infrared reflection spectroscopy,NIRS)利用蛋白质、核酸、脂质、碳水化合物等细胞大分子物质的结构和组成对近红外光产生特征吸收的原理,利用功能基团的差异作为鉴定微生物的基础,通过样品选择、化学计量学方法选择、光谱预处理、建模区间选择等步骤,建立微生物与光谱的定性或定量模型,从而实现用近红外光谱信息对微生物含量的快速检测。
因此,近红外光谱技术具有快速、无损、高效、成本低、可多组分同时检测等优点,利用近红外光谱技术建立快速检测乳酸菌数的方法具有重要意义和应用前景。但能否将近红外光谱方法用于检测发酵乳中乳酸菌数中并能达到在线检测的高效性和及时性要求尚未可知。
发明内容
本发明的目的在于提供一种快速检测发酵乳中乳酸菌数的方法,该方法操作简捷,快速、无损、高效、准确、成本低,不污染环境,且可多组分同时检测,能满足发酵乳产品在线质量检验的高效性和及时性要求。
本申请发明人经过研究发现,近红外光谱不仅能反映微生物细胞壁、细胞膜、细胞质甚至细胞核中的蛋白质、多糖、脂质、核酸及其大分子、水分等混合成分子的分子震动信息,而且能敏锐地探测分子基团及其周围环境的变化。通过测定微生物的近红外光谱图即可获得微生物及其生物大分子结构的信息,用于鉴别微生物的种类、微生物的状态和数量。找出不同微生物图谱间的细微差别,确定不同微生物的特征谱峰和图谱带,为微生物的判别、分类、鉴定和大范围筛选提供依据。
鉴于此,本申请发明人尝试使用相应的NIR化学计量学软件,建立含量-光谱对应的NIR定标模型。录制未知样品的NIR谱,调用NIR定标模型,定量发酵乳中的乳酸菌数。
进一步的,本发明提供的一种快速检测发酵乳中乳酸菌数的方法,包括以下步骤:
(1)建立发酵乳中乳酸菌数定标模型:
采集发酵乳作为定标集样品,采用平板计数法对定标集样品进行分析检测,得出定标集样品的乳酸菌数实测值,并转换成对数值,建立基础数据库;
同步对定标集样品进行近红外光谱扫描,录制定标集样品的近红外光谱数据,建立近红外光谱数据库;
将所述定标集样品的乳酸菌数基础数据库与所述近红外光谱数据库进行一一对应,建立乳酸菌数定标模型;
(2)采用交叉验证法对步骤(1)中所述的乳酸菌数定标模型进行内部验证;
(3)待测发酵乳样品的检测:将待测发酵乳样品进行近红外光谱扫描,获得待测发酵乳样品的近红外光谱数据,将数据导入步骤(1)中建立的乳酸菌数定标模型,计算结果经换算得待测发酵乳样品的乳酸菌数的预测值。
在上述快速检测发酵乳中乳酸菌数的方法中:
优选的,步骤(1)中所述乳酸菌数包括双歧杆菌和乳杆菌计数A、嗜热链球菌计数B以及乳酸菌总数C,其中乳酸菌总数C为双歧杆菌和乳杆菌计数A与嗜热链球菌计数B二者结果之和;步骤(1)中所述乳酸菌数定标模型包括双歧杆菌和乳杆菌计数A、嗜热链球菌计数B以及乳酸菌总数C的定标模型。
乳酸菌数实测值单位为CFU/g,定标集样品乳酸菌数范围在106-109CFU/g,原始数据太大,直接对应建模效果不理想,故把数据转化成log10CFU/g的模式,建立基础数据库,再将所述定标集样品的乳酸菌数基础数据库数据与所述近红外光谱数据库进行一一对应,建立乳酸菌数定标模型,再采用偏最小二乘法建立乳酸菌数的定标模型,并采用交叉验证法对该定标模型进行内部验证。
优选的,步骤(1)中对定标集样品进行近红外光谱扫描时,采用的近红外光谱仪为Thermo ANTARISⅡ,配制InGaAs检测器,透射采样模块、Result 3数据采集软件和TQAnalyst 8数据分析软件。
优选的,采用Result 3数据采集软件采集时,以积分球漫透射方式和空气为背景,设置分辨率8cm-1,扫描范围4000~10000cm-1,扫描次数为32。
由于发酵乳样品的粘度大,考虑到弃样的后处理方便,以及聚乙烯的NIR光谱吸收峰较少,干扰少,优选的,步骤(1)中对定标集样品进行近红外光谱扫描时,所述定标集样品采用聚乙烯密封袋采样,采样后密封并应在15min内采集近红外光谱数据,防止乳酸菌数在常温下发生变化。
由于发酵乳样品是半固体状态,具有流动性,与固体样品的采集数据方法不同,本发明在采集发酵乳近红外光谱信息时,需要对透射采样模块进行改进,优选的,本发明中所述的透射采样模块优选采用凹字形压块,所述定标集样品采用聚乙烯密封袋采样后密封,置于所述近红外光谱仪的积分球漫透射光孔上,然后采用凹字形压块压紧,其中凹字形压块的凹槽部正对聚乙烯密封袋,所述凹字形压块中的凹槽的深度为3mm,然后录制定标集样品的近红外光谱数据。
采用该透射采样模块采样,可以使发酵乳采集条件以及密封袋的厚度保持一致。
具体的,该透漫射采样模块可以是不锈钢材质,“[”型,长10cm×宽3cm×高3cm,反射凹面经镜面抛光处理,凹面深度3mm。采集方法可以使得样品厚度一致。
优选的,步骤(1)中对定标集样品进行近红外光谱扫描时,在温度为(23±2)℃、湿度为(50±5)%的恒温恒湿室中进行。
近红外光谱采集的发酵乳原始光谱信息中除包含与样品结构组成有关的信息外,还可能受测试条件、环境温度、仪器状态、样品状态等因素的影响。而且样品中不同成分也会互相干扰,导致谱线重叠,低含量成分光谱峰被高含量成分光谱峰掩盖等问题。因此,在建模时候,需要考虑多种因素。
建立乳酸菌数模型时,设定的建模条件通常包括:
(1)化学计量学方法:逐步多元线性回归法(Stepwise Multiple LinearRegression,SMLR)、偏最小二乘法(Partial Least Squares,PLS)、主成分回归法(Principal Component Regression,PCR)等;
(2)光程类型(Pathlength Type):恒定光程(Constant)、多元散射校正(Multiplicative Signal Correction,MSC)、标准正态变量变换(Standard NormalVariate,SNV)等;
(3)数据格式(Data Format):原始光谱(Spectrum)、一阶导数(Firstderivative,1st Der)、二阶导数(Second derivative,2nd Der);
(4)平滑类型(Smoothing):不光滑(No Smoothing,NS)、卷积平滑滤波(Savizky-Golay Filter数据点为7,3项式平滑滤波,S-G)、Norris导数平滑滤波(Norris derivativefilter,ND)等;
(5)建模波段范围:4000~10000cm-1。
本发明通过对光谱信息的采集方法进行改良、数据的筛选、对检测数据进行多种预处理以达到消除随机噪声、样品背景干扰和测试条件引起的差异,通过对光谱信息进行有效提取和优化,提高了分辨率和运算效率,最终成功建立高效快速准确的定标模型。
即本申请发明人经过进一步的实验,发现以下条件建立模型时,获得的发酵乳中乳酸菌数的预测值与真实值更加接近。
优选的,步骤(1)中建立乳酸菌数定标模型时,采用偏最小二乘法,经过Savitzky-golay filter数据点为7,3项式平滑滤波处理,并结合一阶导数和多元散射校正(MSC)对近红外光谱的数据进行处理。
优选的,步骤(1)中建立乳酸菌数定标模型时,选取波段为5569-5716cm-1、5724-6403cm-1、7197-7506cm-1范围内的光谱数据建立模型。
步骤(2)中采用交叉验证法对该定标模型进行内部验证,交叉验证法的原理为:假设定标集有n个样品,可以每次从中取出m(m=1,2.3,…)个样品,作为临时验证集,以其余的(n-m)个样品作为校正集进行建模,然后对这m个样品进行预测,如此循环,则分别得到n个样品的交叉预测值,再以交叉预测值与标准值作相关图,同样类似的可以计算校正集交互验证的根均方误差(RMSECV)和相关系数(Rcv)。
优选的,步骤(1)中采用主成分分析-马氏距离方法剔除发酵乳中异常数据,获得定标集样品,提高近红外光谱定量分析的可靠性;步骤(1)中所述发酵乳为含活乳酸菌的发酵乳,所述定标集样品包括校正集样品和验证集样品,步骤(2)中采用交叉验证法对所述的乳酸菌数定标模型进行内部验证时,从定标集样品中取出85%的样品作为校正集(Calibration)样品进行建模,以剩下的15%样品作为验证集(Validation)样品进行验证。
作为本发明的一种优选的实施方式,所述定标集样品为67批次,其中校正集样品为57批次,验证集样品为10批次。
选择具有代表性的样品作为定标集样品,如选择不同风味、黏度、颜色、添加物、乳酸菌数含量等发酵乳样品。
步骤(3)中所述乳酸菌数预测值包括双歧杆菌和乳杆菌A、嗜热链球菌计数B以及乳酸菌总数C的预测值,其中乳酸菌总数C为双歧杆菌和乳杆菌计数A与嗜热链球菌计数B二者结果之和。
因此,本发明研究了采用近红外光谱快速检测发酵乳中乳酸菌数的方法,该方法利用具有代表性的发酵乳样品构成定标集样品,采用平板计数法对定标集样品的乳酸菌数进行分析检测,获得其实测值,再转换成对数值,建立基础数据库;同时采用近红外光谱仪在设定的建模条件下采集其近红外光谱信息,获得光谱数据,建立近红外光谱数据库;将基础数据库和近红外光谱数据库进行一一对应,采用PLS法建立发酵乳中乳酸菌数的定标模型,并验证;取待测样品,对其进行近红外光谱分析,将光谱数据导入定标模型,经换算得出待测样品的乳酸菌数。该方法具有快速、高效、准确、成本低、不污染环境等优点,能满足发酵乳产品在线质量检验的高效性和及时性要求。
本申请将上述近红外光谱方法直接应用于发酵乳中乳酸菌数检测时,在建模过程中遇到以下的困难:
(1)目前基于近红外技术快速检测微生物的方法研究中,由于微生物在样品中的分布中存在分布不均、含量低(一般低于103CFU/g)等特点,需要对微生物进行富集,常见的富集方法有对目标菌株进行分离、纯化、培养,制成粉末状菌体,或者采用膜过滤的方法收集菌体,其样品形态都是菌体粉末,通过对标准菌株或典型菌株的纯培养物制成,样品的纯度高,含量高,一般达109CFU/g以上,基质相似,近红外光谱信息不受样品基质的干扰,但以上方法不能实现真正意义上的在线快速检测。
(2)本发明定标集样品的乳酸菌数含量为106~109CFU/g,且分布均匀,具有在线快速检测的潜力。但发酵乳成分复杂,包含蛋白质、脂肪、乳糖、乳酸、盐分、食品添加剂等物质,乳酸菌数在发酵乳中所占含量较低。不同品牌不同口味的发酵乳还会添加果蔬或谷物,导致发酵乳的颗粒大小、黏度、颜色存在很大差异。因此利用近红外光谱分析技术对发酵乳中的微生物含量进行定量分析,属于从复杂、重叠、变动的背景中提取弱信息。直接扫描发酵乳获得的近红外光谱易受样品基质的干扰,多成分物质造成光谱的背景复杂,增加了近红外光谱分析的复杂性。
本发明通过收集各种类型样品,样品类型包括不同风味、黏度、颜色、添加物、乳酸菌数含量等,选择具有代表性的样品作为定标集样品,尽量保证样品的乳酸菌数在建模范围内分布均匀,且采用主成分分析-马氏距离方法剔除次异常数据,选择有效样品,提高模型的预测准确度。同时采用多元分析,利用特定波长短提取有效信息、运用校正技术等方法,成功提取有效信息,建立快速准确的乳酸菌数定标模型,可实现真正意义上的在线快速分析。
(3)近红外光谱采集的发酵乳原始光谱信息中除包含与样品结构组成有关的信息外,还可能受测试条件、环境温度、仪器状态、样品状态等因素的影响。而且样品中不同成分也会互相干扰,导致谱线重叠,低含量成分光谱峰被高含量成分光谱峰掩盖等问题。本发明通过对光谱信息的采集方法进行改良、数据的筛选、对检测数据进行多种预处理以达到消除随机噪声、样品背景干扰和测试条件引起的差异,通过对光谱信息进行有效提取和优化,提高了分辨率和运算效率,最终成功建立高效快速准确的定标模型。
(4)由于发酵乳是浑浊黏稠的半固体,样品的颗粒度不均匀会引起散射,当样品组成发生变化,其吸收系数和散射系数也随之变化,产生较大的散射误差,导致建模的准确性不高。本发明通过积分球漫透射的方式,同时采用凹字形透射采样模块采样,可以使发酵乳采集条件以及样品的厚度保持一致,大大增强了信号强度,改良了近红外吸收光谱采集形式,成功建立定标模型。
与现有技术相比,本发明具有以下优点:
(1)操作简单便捷:本发明方法仪器的高度自动化降低了对操作者的技能要求,无需前处理和稀释,只需进行简单操作即可获得检测结果,易于推广;
(2)效率高,耗时短:近红外技术检测耗时少于一分钟,极大地缩短了检测时间,保障生产过程的连续性,且能用于多组分的同时检测,尤其适用于企业在线监测或快速检测;
(3)无需试剂,成本低,不污染环境:近红外光谱分析中只需要取得样品的光谱信号,不使用任何试剂,可以脱离试验试剂的影响,测试过程中不产生污染;
(4)因此,本发明方法结果表明使用近红外透射光谱快速测定发酵乳中的乳酸菌数值是可行的,为今后发酵乳中乳酸菌数的快速检测提供了理论依据。
附图说明
图1是本发明实施例2中发酵乳中乳酸菌数的NIR透射光谱图;
图2是本发明实施例2中发酵乳中嗜热链球菌数定标模型中预测值和真实值(测定值)的相关性关系;
图3是本发明实施例2中发酵乳中双歧杆菌和乳杆菌数定标模型中预测值和真实值(测定值)的相关性关系;
图4是本发明实施例2中发酵乳中乳酸菌总数定标模型中预测值和真实值(测定值)的相关性关系;
图5是本发明实施例2中发酵乳中嗜热链球菌数定标模型中预测值和真实值(测定值)的相对误差;
图6是本发明实施例2中发酵乳中双歧杆菌和乳杆菌数定标模型中预测值和真实值(测定值)的相对误差;
图7是本发明实施例2中发酵乳中乳酸菌总数定标模型中预测值和真实值(测定值)的相对误差;
图8是本发明实施例2中发酵乳中嗜热链球菌数定标模型的十字交叉验证图(Cross validation);
图9是本发明实施例2中发酵乳中双歧杆菌和乳杆菌数定标模型的十字交叉验证图(Cross validation);
图10是本发明实施例2中发酵乳中乳酸菌总数定标模型的十字交叉验证图(Crossvalidation)。
具体实施方式
以下通过具体实施例来说明本发明中的发酵乳中乳酸菌数的检测方法:
实施例1
本实施例提供的快速检测发酵乳中乳酸菌数的方法,包括以下步骤:
1仪器和材料
1.1仪器与设备
高压灭菌锅:CL-40M ALP
生物安全柜:AC2-6S1 ESCO
隔水式恒温培养箱:GHP-9160上海一恒科技有限公司
拍击式均质器:Easymax AES Chemunex
电子稀释器:dilumat S AES Chemunex
厌氧培养罐:AnaeroPack日本三菱MGC
近红外光谱仪Thermo ANTARISⅡ:InGaAs检测器,安装Result 3数据采集软件和TQ Analyst 8数据分析软件。
1.2材料
材料:聚乙烯密实袋(7cm×5cm);无菌均质袋interscience;厌氧产气袋三菱瓦斯化学株式会社
样品:含活乳酸菌发酵乳样品;
试剂:氯化钠广州化学试剂厂
培养基:MRS(Man Rogosa Sharpe)培养基广东环凯微生物科技有限公司;MC(Modified Chalmers)培养基北京陆桥技术股份有限公司
2平板计数法检测发酵乳中乳酸菌数
采集具有代表性的发酵乳样品构成定标集样品,所有发酵乳样品均为含活乳酸菌的发酵乳,按照中国国家标准GB 4789.35-2016《食品安全国家标准食品微生物学检验乳酸菌检验》测定所述定标集样品的乳酸菌数,构成所述定标集样品的乳酸菌数平板计数法测定值,建立基础数据库。
其具体检验方法和计算方法如下:
2.1样品制备
样品的全部制备过程均应遵循无菌操作程序,以无菌生理盐水作为稀释液。称取25g样品,置于225mL稀释液的无菌均质袋中,用拍击式均质器拍打2min制成1:10的样品稀释液。
2.2 10倍系列稀释
吸取2.1中1:10样品稀释液1mL,注入装有9mL稀释液的试管中,使用旋涡混合器振摇试管1min,使其混合均匀,制成1:100的样品稀释液。按照上述操作步骤,按10倍系列依次制备样品稀释液。
2.3双歧杆菌和乳杆菌计数A
依据待检样品标签声称的乳酸菌类别和含量,估算待检样品中是否含有双歧杆菌和乳杆菌及其数量,选择2个~3个合适的连续稀释度,每个稀释度吸取1mL样品稀释液置于无菌培养皿内,每个稀释度做两个平行。然后注入约15mL冷却至48℃±1℃的MRS培养基,转动培养皿使之混合均匀。待培养基凝固后置于36℃±1℃培养箱中厌氧培养72h±2h,培养后计数与计算。从制备样品到倾注平板要求在15min内完成,以防双歧杆菌和乳杆菌数量发生变化。
2.4嗜热链球菌计数B
依据待检样品标签声称的乳酸菌类别和含量,估算待检样品中是否含有嗜热链球菌及其数量,选择2个~3个合适的连续稀释度,每个稀释度吸取1mL样品稀释液置于无菌培养皿内,每个稀释度做两个平行。然后注入约15mL冷却至48℃±1℃的MC培养基,转动培养皿使之混合均匀。待培养基凝固后置于36℃±1℃培养箱中需氧培养72h±2h,培养后计数与计算。从制备样品到倾注平板要求在15min内完成,以防嗜热链球菌数量发生变化。
2.5乳酸菌总数计数C
乳酸菌总数C为双歧杆菌和乳杆菌计数A和嗜热链球菌计数B二者结果之和,即C=A+B。双歧杆菌和乳杆菌计数A和嗜热链球菌计数B的计数范围均为30CFU~300CFU之间。
2.6结果计算
若只有一个稀释度平板上的菌落数在适宜计数范围内,计算两个平板菌落数的平均值,再将平均值乘以相应稀释倍数,作为每克中菌落数结果,按公式(1)计算。
式中:
N——样品中菌落数;
∑C——平板菌落数之和;
d——稀释因子。
若有两个连续稀释度的平板菌落数在适宜计数范围内时,按公式(2)计算。
式中:
N——样品中菌落数;;
∑C——平板(含适宜范围菌落数的平板)菌落数之和;
n1——第一稀释度(低稀释倍数)平板个数;
n2——第二稀释度(高稀释倍数)平板个数;
d——稀释因子(第一稀释度)。
2.7修约原则与报告
修约原则:菌落数小于100CFU时,按“四舍五入”原则修约,以整数报告。菌落数大于或等于100CFU时,第3位数字采用“四舍五入”原则修约后,用10的指数形式来表示,采用两位有效数字。
报告:根据菌落计数结果,报告单位以CFU/g表示。
3近红外光谱信息采集
采用近红外光谱仪在设定的建模条件下采集定标集样品的近红外光谱信息,构成所述定标集样品的近红外波段光谱信息,建立近红外光谱数据库。
3.1取样方法
使用聚乙烯密封袋(7cm×5cm)采样12g-15g,采样后密封并应在15min内完成检测,防止乳酸菌数在常温下发生变化而影响检测结果。
3.2采集方法
检测前开机预热30min,以空气为参比,扣除背景,性能测试后进行光谱扫描。测试过程环境温湿度控制在(23±2)℃,(50±5)%范围内。采用Result 3数据采集;采样方式:积分球漫透射;分辨率:8cm-1;波长范围:4000~10000cm-1;扫描次数:32次;信号增益:1X。
取均匀样品装进聚乙烯密实袋,密封后平铺于积分球漫透射光孔上,录制NIR谱。每个样品重复测定2次,取其平均光谱。采集的光谱数据用TQ Analyst 8数据分析软件进行处理和计算。采用积分球漫透射方法原因:减少散射误差。
4定标模型的评价指标
定标模型的评价指标直接影响模型的稳定性和准确性。本发明在模型建立和验证过程中使用的评价指标如下:
校正集根均方误差(Root Mean Square Error of Calibration,RMSEC)和预测集根均方误差(Root Mean Square Error of Prediction,RMSEP):
校正集交互验证的根均方误差(Root Mean Square Error of CrossValidation,RMSECV):
相关系数R:
式中:n表示样本数;yi表示第i个样品使用平板计数法获得的实测值;表示第i个样品的近红外预测值;表示样品使用平板计数法获得的实测值的平均值。RMSEC、RMSEP和RMSECV分别表示校正集、预测集样品和内部交叉验证时预测值和实测值偏离的大小,其值越小,表示所建模型的预测精度越大,且RMSEC、RMSEP和RMSECV之间的差异越小,模型的预测效果就越好。相关系数R表示样品参数指标的实测值和近红外预测值之间的相关程度,其值越接近于1,说明实测值与近红外预测值之间的相关程度越好。理想状态条件下,相关系数R的值为1,均方根误差为0。
5乳酸菌数实测值数据模式转换
按照中国国家标准GB 4789.35-2016《食品安全国家标准食品微生物学检验乳酸菌检验》测定所述定标集样品的乳酸菌数单位为CFU/g,定标集样品乳酸菌数范围在106-109CFU/g,原始数据太大,尝试直接建模效果不理想,故把数据转化成log10CFU/g的模式,再进行建模优化。利用TQ Analyst 8数据分析软件建立乳酸菌数的定标校正模型。
6建立模型的方法
6.1样品选择
在近红外光谱定量分析过程中,由于测量仪器、测试方法和环境等客观因素的影响以及技术人员本身主观因素的作用,再加上样本来源的多样性,用同一个模型分析不同来源的样品也可能存在模型不适应的问题。因此异常值的判别和处理是提高近红外分析质量的一个重要步骤。通过优化筛选,剔除偏离较大的样品,使其误判率为0,如采用主成分分析-马氏距离方法剔除异常数据,获得定标集样品,以提高近红外光谱定量分析的可靠性。
6.2化学计量学方法选择
在近红外分析中由于光谱吸收强度低,谱带较宽及光谱相互重叠非常严重,在进行定性及定量分析时需要化学计量学方法提取信息。在近红外光谱技术分析中经常采用的化学计量学方法有:逐步多元线性回归法(SMLR)、偏最小二乘法(PLS)、主成分回归法(PCR)等。本发明在建模过程中,以校正集根均方误差RMSEC及其相关系数Rc,预测集根均方误差RMSEP及其相关系数Rp为评价指标,对不同NIR分析方法比较,最终确定PLS为最佳处理方法。
6.3光谱数据预处理
近红外仪器自身的随机噪音,以及样品的不均匀、光散射等干扰,都会导致光谱偏移或漂移。光谱数据预处理可以降低样品表面不均匀和色差等因素的影响。光谱微分处理可以消除基线漂移、强化谱带特征、克服谱带重叠、提高分辨率和灵敏度。但与此同时导数的计算会引进噪音,降低信噪比,而平滑可以降低高频随机噪音。
发明以校正集根均方误差RMSEC及其相关系数Rc,预测集根均方误差RMSEP及其相关系数Rp为评价指标,比较了:(1)光程类型:恒定光程、多元散射校正(MSC)、标准正态变量变换(SNV);(2)数据格式:原始光谱、一阶导数(1st Der)、二阶导数(2nd Der);(3)平滑类型:不光滑(NS)、卷积平滑滤波(S-G)、Norris导数平滑滤波(ND)等三种光谱预处理方法,最终确定“MSC+1st Der+S-G”为最佳处理方法。
6.4建模波段选择
对反映样品参数信息突出的光谱区域进行优选,可以优化光谱范围,净化图谱信息,提高运算效率。获得近红外波段光谱信息中的波段范围为4000-10000cm-1。根据线性相关系数和NIR的吸光度,不断优化波段范围,最终选择定量光谱区间段为5569-5716cm-1、5724-6403cm-1、7197-7506cm-1。
7采用交叉验证法对定标模型进行内部验证
本实施例以定标集样品的85%作为校正集样品,定标集样品的15%作为验证集样品进行交叉验证。
8待测发酵乳样品的检测
取待测样品,对其进行近红外光谱分析,获得待测样品的近红外波段光谱信息,将光谱数据导入定标模型,结果经换算得待测样品的乳酸菌数预测值。
实施例2
本实施提供的快速检测发酵乳中乳酸菌数的方法,包括以下步骤:
(1)建立发酵乳中乳酸菌数定标模型:
采集发酵乳作为定标集样品,采用平板计数法对定标集样品进行分析检测,得出定标集样品的乳酸菌数实测值,并转换成对数值,建立基础数据库;
同步对定标集样品进行近红外光谱扫描,录制定标集样品的近红外光谱数据,建立近红外光谱数据库;
将所述定标集样品的乳酸菌数基础数据库与所述近红外光谱数据库进行一一对应,建立乳酸菌数定标模型;
(2)采用交叉验证法对步骤(1)中所述的乳酸菌数定标模型进行内部验证;
(3)待测发酵乳样品的检测:将待测发酵乳样品进行近红外光谱扫描,获得待测发酵乳样品的近红外光谱数据,将数据导入步骤(1)中建立的乳酸菌数定标模型,结果经换算得待测发酵乳样品的乳酸菌数的预测值。
具体实施过程如下:
1实验部分
1.1材料和仪器
含活乳酸菌发酵乳样品共70批次,由国家食品产品质量监督检验中心(广东)提供。
仪器同实施例1。
1.2平板计数法检测样品中乳酸菌数同实施例1。
1.3采集发酵乳样品近红外光谱信息
1.3.1取样方法
使用聚乙烯密封袋(7cm×5cm)采样12g-15g,采样后密封并应在15min内完成检测,防止乳酸菌数在常温下发生变化而影响检测结果。
采用聚乙烯密实袋装样原因:样品的粘度大、弃样的后处理方便,聚乙烯的NIR光谱吸收峰较少,干扰少。
1.3.2采集方法
ANTARISⅡ仪器采用Result 3操作系统采集光谱,智能透射方式,以积分球漫透射方式和空气为背景,设置分辨率8cm-1,扫描范围4000~10000cm-1,扫描次数为32。根据发酵乳的透明度和遮盖率,确定信号增益为1X。取均匀样品装进聚乙烯密实袋,密封后平铺于积分球漫透射光孔上,录制NIR谱。测试过程在温度(23±2)℃、湿度(50±5)%的恒温恒湿室中进行。每个样品重复测定2次,取其平均光谱。采集的光谱数据用TQ Analyst 8数据分析软件进行处理和计算。
1.4乳酸菌数实测值数据模式转换
把乳酸菌数实测值数据转化成log10CFU/g的模式,再进行建模优化。
1.5建立模型
将所述定标集样品的乳酸菌数基础数据库与所述近红外光谱数据库进行一一对应,以RMSEC、RMSEP、RMSECV及其相关系数Rc、Rp、Rcv为评价指标,通过剔除异常数据、选择化学计量学方法、光程类型、数据格式、平滑类型等光谱数据预处理、建模波段优化等步骤,建立乳酸菌数定标模型,并采用交叉验证法对定标模型进行内部验证。
2结果与讨论
2.1剔除异常数据
采用主成分分析-马氏距离方法剔除3批次异常数据,提高近红外光谱定量分析的可靠性,以剩余的67批次样品进行建模。NIR透射光谱叠加谱如图1所示。
2.2化学计量学方法选择
以RMSEC、RMSEP及其相关系数Rc、Rp为评价指标,对SMLR、PLS、PCR分析方法比较,结果如表1~2所示。由表1~2可知,PLS建模方法的RMSEC、RMSEP最小,Rc、Rp最接近1,故PLS为最佳处理方法。
表1 SMLR、PLS、PCR建模方法的RESEC和Rc
表2 SMLR、PLS、PCR建模方法的RESEP和Rp
2.3光谱数据预处理
以RMSEC、RMSEP及其相关系数Rc、Rp为评价指标,对恒定光程、MSC、SNV、一阶导数(S-G)、二阶导数(S-G)、一阶导数+MSC(S-G)、二阶导数+MSC(S-G)、一阶导数+SNV(S-G)、二阶导数+SNV(S-G)分析方法比较,结果如表3~4所示。综合考虑双歧杆菌和乳杆菌A、嗜热链球菌B和乳酸菌数C的RMSEC、RMSEP及其相关系数Rc、Rp,一阶导数+MSC(S-G)为最佳预处理方法。
表3近红外光谱数据预处理方法的RESEC和Rc
表4近红外光谱数据预处理方法的RESEP和Rp
2.4建模波段选择
发酵乳成分复杂,各种基质在近红外区域吸收丰富,吸收强度大,吸收峰宽,故相互干扰。通过选择高波数的定量波段5569-5716cm-1、5724-6403cm-1、7197-7506cm-1,不选择低波数干扰波段4000~5200cm-1和高波段干扰波段8900~9500cm-1,去除无关信息,提取有用信息,降低数据分析的计算量,达到优化光谱范围,提高运算效率的目的。
2.5采用交叉验证法对定标模型进行内部验证
以定标集样品中的57批次作为校正集样品,另外的10批次作为验证集样品,进行多变量建立乳酸菌数的定标模型,嗜热链球菌计数B、双歧杆菌和乳杆菌计数A、乳酸菌总数计数C的预测值和真实值的相关性关系如图2~4所示。
由图2~4可知校正集和预测集中实际测定值和模型预测值之间的相关系数都达到了0.91以上,模型的拟合度和预测精度都比较好。其预测值和真实值的相对误差如图5~7所示。
乳酸菌数定标模型的十字交叉验证图(Cross validation)如图8~10所示。由图8~10可知RMSEC、RMSEP、RMSECV的值都小于0.2,且相互之间差异很小,Rc、Rp、Rcv的值都大于0.89,说明模型的预测值和实测值相关性好。定标模型的RMSEC都小于RMSEP,说明样品代表性好,样品信息提取充分;RMSEP和RMSECV数值相近,说明建模样品与验证样品都具有代表性,模型信息拟合充分,模型预测性好。
3独立样品预测
从图2~10可知,经过剔除异常数据、化学计量学方法优选、光谱数据预处理、光谱波段优化等步骤建立的乳酸菌数定标模型,具有较好的预测能力,为进一步说明模型的可行性,本发明利用独立样品对近红外的准确性进行了考察。
准备10个含活乳酸菌发酵乳样品,按照1.2平板计数法检测样品中乳酸菌数的实测值,按1.3采集其近红外光谱信息,将光谱数据导入定标模型,得出样品的乳酸菌数预测值,同时将定标模型预测值和平板计数法实测值对比,分析其误差是否在允许范围内,要求两种方法的检测结果的绝对差值不超过0.45,具体结果见表5。
表5发酵乳样品的乳酸菌数定标模型预测值和平板计数法实测值比较(log10CFU/g)
由表5可知,10组数据的平板计数法实测值和模型预测值的绝对差值都小于0.45,均在允许范围内。
最后需要说明的是,虽然本发明以较佳实施例揭露如上,但并非用以限定本发明的范围。任何本领域的普通技术人员,在不脱离本发明的范围内,当可作些许的改进,即凡是依照本发明所做的同等改进,应为本发明的范围所涵盖。
Claims (10)
1.一种快速检测发酵乳中乳酸菌数的方法,包括以下步骤:
(1)建立发酵乳中乳酸菌数定标模型:
采集发酵乳作为定标集样品,采用平板计数法对定标集样品进行分析检测,得出定标集样品的乳酸菌数实测值,并转换成对数值,建立基础数据库;
同步对定标集样品进行近红外光谱扫描,录制定标集样品的近红外光谱数据,建立近红外光谱数据库;
将所述定标集样品的乳酸菌数基础数据库与所述近红外光谱数据库进行一一对应,建立乳酸菌数定标模型;
(2)采用交叉验证法对步骤(1)中所述的乳酸菌数定标模型进行内部验证;
(3)待测发酵乳样品的检测:将待测发酵乳样品进行近红外光谱扫描,获得待测发酵乳样品的近红外光谱数据,将数据导入步骤(1)中建立的乳酸菌数定标模型,结果经换算得待测发酵乳样品的乳酸菌数的预测值。
2.根据权利要求1所述的快速检测发酵乳中乳酸菌数的方法,其特征是:步骤(1)中所述乳酸菌数包括双歧杆菌和乳杆菌计数A、嗜热链球菌计数B以及乳酸菌总数C,其中乳酸菌总数C为双歧杆菌和乳杆菌计数A与嗜热链球菌计数B二者结果之和;步骤(1)中所述乳酸菌数定标模型包括双歧杆菌和乳杆菌计数A、嗜热链球菌计数B以及乳酸菌总数C的定标模型。
3.根据权利要求1所述的快速检测发酵乳中乳酸菌数的方法,其特征是:步骤(1)中对定标集样品进行近红外光谱扫描时,采用的近红外光谱仪为Thermo ANTARIS Ⅱ,配制InGaAs检测器,透射采样模块、Result 3数据采集软件和TQ Analyst 8数据分析软件。
4.根据权利要求3所述的快速检测发酵乳中乳酸菌数的方法,其特征是:采用Result 3数据采集软件采集时,以积分球漫透射方式和空气为背景,设置分辨率8cm-1,扫描范围4000~10000cm-1,扫描次数为32。
5.根据权利要求3所述的快速检测发酵乳中乳酸菌数的方法,其特征是:步骤(1)中对定标集样品进行近红外光谱扫描时,所述定标集样品采用聚乙烯密封袋采样,采样后密封并应在15min内采集近红外光谱数据,防止乳酸菌数在常温下发生变化而影响结果。
6.根据权利要求5所述的快速检测发酵乳中乳酸菌数的方法,其特征是:所述的透射采样模块为凹字形压块,所述定标集样品采用聚乙烯密封袋采样后密封,置于所述近红外光谱仪的积分球漫透射光孔上,然后采用凹字形压块压紧,其中凹字形压块的凹槽部正对聚乙烯密封袋,所述凹字形压块中的凹槽的深度为3mm,录制定标集样品的近红外光谱数据。
7.根据权利要求1所述的快速检测发酵乳中乳酸菌数的方法,其特征是:步骤(1)中对定标集样品进行近红外光谱扫描时,在温度为(23±2)℃、湿度为(50±5)%的恒温恒湿室中进行。
8.根据权利要求1所述的快速检测发酵乳中乳酸菌数的方法,其特征是:步骤(1)中建立乳酸菌数定标模型时,采用偏最小二乘法,经过savitzky-golay filter数据点为7,3项式平滑滤波处理,并结合一阶导数和多元散射校正对近红外光谱的数据进行处理。
9.根据权利要求1所述的快速检测发酵乳中乳酸菌数的方法,其特征是:步骤(1)中建立乳酸菌数定标模型时,选取波段为5569-5716cm-1、5724-6403cm-1、7197-7506cm-1范围内的光谱数据建立模型。
10.根据权利要求1所述的快速检测发酵乳中乳酸菌数的方法,其特征是:步骤(1)中采用主成分分析-马氏距离方法剔除发酵乳中异常数据,获得定标集样品,提高近红外光谱定量分析的可靠性;步骤(1)中所述发酵乳为含活乳酸菌的发酵乳,所述定标集样品包括校正集样品和验证集样品,步骤(2)中采用交叉验证法对所述的乳酸菌数定标模型进行内部验证时,从定标集样品中取出85%的样品作为校正集样品进行建模,以剩下的15%样品作为验证集样品进行验证。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910069041.9A CN109557014B (zh) | 2019-01-24 | 2019-01-24 | 一种快速检测发酵乳中乳酸菌数的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910069041.9A CN109557014B (zh) | 2019-01-24 | 2019-01-24 | 一种快速检测发酵乳中乳酸菌数的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109557014A true CN109557014A (zh) | 2019-04-02 |
CN109557014B CN109557014B (zh) | 2021-08-27 |
Family
ID=65873578
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910069041.9A Active CN109557014B (zh) | 2019-01-24 | 2019-01-24 | 一种快速检测发酵乳中乳酸菌数的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109557014B (zh) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110057777A (zh) * | 2019-05-24 | 2019-07-26 | 南京农业大学 | 一种面粉中呕吐毒素的定量检测方法 |
CN111259970A (zh) * | 2020-01-20 | 2020-06-09 | 江苏大学 | 一种馒头加工过程中面团发酵状态的智能化监测方法 |
CN111665216A (zh) * | 2020-06-02 | 2020-09-15 | 中南民族大学 | 一种判定速冻面米制品中大肠杆菌和金黄色葡萄球菌污染程度的方法 |
CN114384105A (zh) * | 2020-10-16 | 2022-04-22 | 仙乐健康科技股份有限公司 | 益生菌片剂稳定性试验预测模型的构建方法及应用方法 |
Citations (15)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5234904A (en) * | 1989-06-26 | 1993-08-10 | Kabushiki Kaisha Yakult Honsha | Antihypertensives containing polysaccharide peptidoglycan complexes extracted from the cell wall of lactic acid bacteria |
JPH07194394A (ja) * | 1993-12-28 | 1995-08-01 | Yakult Honsha Co Ltd | 乳酸菌含有物中の酵母類の検査用培地及び検査方法 |
US6054262A (en) * | 1995-07-31 | 2000-04-25 | Kabushiki Kaisha Yakult Honsha | Method for controlling culture of lactic bacteria |
KR20040090758A (ko) * | 2003-04-18 | 2004-10-27 | 박정호 | 유효 미생물 세라믹스 분말로된 원적외선 방사체 |
JP2006333838A (ja) * | 2005-06-06 | 2006-12-14 | Univ Of Tsukuba | 減圧噴霧乾燥法による発酵乳粉末の製造方法 |
JP2011036190A (ja) * | 2009-08-12 | 2011-02-24 | Enseki Aojiru Kk | 顆粒状経口食品の製造方法及びその顆粒状経口食品 |
CN103592257A (zh) * | 2013-12-02 | 2014-02-19 | 广东产品质量监督检验研究院 | 一种用近红外光谱快速测定木器涂料稀释剂中甲苯、乙苯和二甲苯含量的方法 |
CN103604771A (zh) * | 2013-12-02 | 2014-02-26 | 广东产品质量监督检验研究院 | 一种利用近红外光谱的主成分分析-马氏距离分类法鉴定水性墙面涂料常用乳液类型的方法 |
CN103616335A (zh) * | 2013-12-02 | 2014-03-05 | 江苏大学 | 发酵食品产酸菌产酸能力的快速鉴别方法 |
CN104031976A (zh) * | 2014-06-16 | 2014-09-10 | 贝因美婴童食品股份有限公司 | 一种检测发酵乳中乳酸菌数的方法 |
CN104824156A (zh) * | 2015-05-19 | 2015-08-12 | 江西阳光乳业股份有限公司 | 一种益生菌发酵乳 |
CN104950046A (zh) * | 2014-03-26 | 2015-09-30 | 浙江科技学院 | 一种鉴别发酵型和调配型酸乳饮品的游离氨基酸分析方法 |
CN106706555A (zh) * | 2016-11-21 | 2017-05-24 | 无锡迅杰光远科技有限公司 | 基于近红外光谱技术的奶粉测定方法及系统 |
CN207450476U (zh) * | 2017-09-26 | 2018-06-05 | 东莞市米目米信息技术有限公司 | 一种乳酸菌饮料生产用的贴标机 |
CN108562622A (zh) * | 2018-02-05 | 2018-09-21 | 西北农林科技大学 | 一种基于介电特性技术的生鲜羊乳菌落总数快速检测方法 |
-
2019
- 2019-01-24 CN CN201910069041.9A patent/CN109557014B/zh active Active
Patent Citations (15)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5234904A (en) * | 1989-06-26 | 1993-08-10 | Kabushiki Kaisha Yakult Honsha | Antihypertensives containing polysaccharide peptidoglycan complexes extracted from the cell wall of lactic acid bacteria |
JPH07194394A (ja) * | 1993-12-28 | 1995-08-01 | Yakult Honsha Co Ltd | 乳酸菌含有物中の酵母類の検査用培地及び検査方法 |
US6054262A (en) * | 1995-07-31 | 2000-04-25 | Kabushiki Kaisha Yakult Honsha | Method for controlling culture of lactic bacteria |
KR20040090758A (ko) * | 2003-04-18 | 2004-10-27 | 박정호 | 유효 미생물 세라믹스 분말로된 원적외선 방사체 |
JP2006333838A (ja) * | 2005-06-06 | 2006-12-14 | Univ Of Tsukuba | 減圧噴霧乾燥法による発酵乳粉末の製造方法 |
JP2011036190A (ja) * | 2009-08-12 | 2011-02-24 | Enseki Aojiru Kk | 顆粒状経口食品の製造方法及びその顆粒状経口食品 |
CN103592257A (zh) * | 2013-12-02 | 2014-02-19 | 广东产品质量监督检验研究院 | 一种用近红外光谱快速测定木器涂料稀释剂中甲苯、乙苯和二甲苯含量的方法 |
CN103604771A (zh) * | 2013-12-02 | 2014-02-26 | 广东产品质量监督检验研究院 | 一种利用近红外光谱的主成分分析-马氏距离分类法鉴定水性墙面涂料常用乳液类型的方法 |
CN103616335A (zh) * | 2013-12-02 | 2014-03-05 | 江苏大学 | 发酵食品产酸菌产酸能力的快速鉴别方法 |
CN104950046A (zh) * | 2014-03-26 | 2015-09-30 | 浙江科技学院 | 一种鉴别发酵型和调配型酸乳饮品的游离氨基酸分析方法 |
CN104031976A (zh) * | 2014-06-16 | 2014-09-10 | 贝因美婴童食品股份有限公司 | 一种检测发酵乳中乳酸菌数的方法 |
CN104824156A (zh) * | 2015-05-19 | 2015-08-12 | 江西阳光乳业股份有限公司 | 一种益生菌发酵乳 |
CN106706555A (zh) * | 2016-11-21 | 2017-05-24 | 无锡迅杰光远科技有限公司 | 基于近红外光谱技术的奶粉测定方法及系统 |
CN207450476U (zh) * | 2017-09-26 | 2018-06-05 | 东莞市米目米信息技术有限公司 | 一种乳酸菌饮料生产用的贴标机 |
CN108562622A (zh) * | 2018-02-05 | 2018-09-21 | 西北农林科技大学 | 一种基于介电特性技术的生鲜羊乳菌落总数快速检测方法 |
Non-Patent Citations (4)
Title |
---|
R. MUELAS ET AL.: ""Evaluation of individual lactic acid bacteria for the fermentation of goat milk: Quality parameters"", 《LWT - FOOD SCIENCE AND TECHNOLOGY》 * |
李宝志等: ""用塑料袋作为无菌采样容器的实验报告"", 《黑龙江医学》 * |
黄宝莹 等: ""发酵乳中乳酸菌数测定的不确定度评估"", 《中国乳品工业》 * |
黄常毅: ""近红外光谱技术在微生物发酵过程中的应用研究"", 《中国优秀硕士学位论文全文数据库 工程科技Ⅰ辑》 * |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110057777A (zh) * | 2019-05-24 | 2019-07-26 | 南京农业大学 | 一种面粉中呕吐毒素的定量检测方法 |
CN111259970A (zh) * | 2020-01-20 | 2020-06-09 | 江苏大学 | 一种馒头加工过程中面团发酵状态的智能化监测方法 |
CN111259970B (zh) * | 2020-01-20 | 2024-02-13 | 江苏大学 | 一种馒头加工过程中面团发酵状态的智能化监测方法 |
CN111665216A (zh) * | 2020-06-02 | 2020-09-15 | 中南民族大学 | 一种判定速冻面米制品中大肠杆菌和金黄色葡萄球菌污染程度的方法 |
CN114384105A (zh) * | 2020-10-16 | 2022-04-22 | 仙乐健康科技股份有限公司 | 益生菌片剂稳定性试验预测模型的构建方法及应用方法 |
Also Published As
Publication number | Publication date |
---|---|
CN109557014B (zh) | 2021-08-27 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109557014A (zh) | 一种快速检测发酵乳中乳酸菌数的方法 | |
CN109540838A (zh) | 一种快速检测发酵乳中酸度的方法 | |
JP6339244B2 (ja) | Ft−irスペクトルデータの多変量統計分析を用いた果実の糖度及び酸度の予測方法 | |
CN102621192B (zh) | 一种利用电子鼻检测芒果新鲜度的方法 | |
CN106950192A (zh) | 一种基于近红外光谱分析技术的植物蛋白饮料中主要成分含量快速检测的方法 | |
CN105136737A (zh) | 基于近红外光谱的馒头中马铃薯全粉含量快速测定的方法 | |
CN103411906A (zh) | 珍珠粉和贝壳粉的近红外光谱定性鉴别方法 | |
CN105044024A (zh) | 一种基于近红外光谱技术对葡萄果实进行无损检测的方法 | |
CN110487746A (zh) | 一种基于近红外光谱检测娃娃菜品质的方法 | |
CN104316489A (zh) | 一种近红外光谱检测灵芝提取物掺假的方法 | |
CN103018181A (zh) | 基于相关性分析和elm神经网络的软测量方法 | |
CN110006886A (zh) | 一种纳米化色敏传感器及其判别小麦霉变程度的方法 | |
US8679739B2 (en) | Method and apparatus for diagnostic analyses | |
CN101936894A (zh) | 一种基于近红外光谱及显微脂肪细胞数据融合的猪肉新鲜度无损检测技术 | |
Akiyama et al. | Monitoring the progress of lactic acid fermentation in yogurt manufacturing using terahertz time-domain–attenuated total-reflection spectroscopy | |
CN105675538A (zh) | 一种胡麻饼养分的检测方法 | |
Giovenzana et al. | Comparison of two immersion probes coupled with visible/near infrared spectroscopy to assess the must infection at the grape receiving area | |
CN108562622B (zh) | 一种基于介电特性技术的生鲜羊乳菌落总数快速检测方法 | |
CN110308226A (zh) | 一种发酵虫草菌粉生产全过程链质量控制快速检测方法 | |
Lénárt et al. | Application of FT-NIR spectroscopy on the determination of the fat and protein contents of lyophilized cheeses | |
CN108796005A (zh) | 一种实时监测谷氨酸棒杆菌发酵过程的方法 | |
CN114062299A (zh) | 一种乳果糖的定量检测方法 | |
CN108489934A (zh) | 一种检测花生油品质的方法 | |
CN113324941A (zh) | 原料牛奶保存时间的快速鉴定方法 | |
CN110296956A (zh) | 一种近红外光谱测定水稻秸秆发酵中有机质含量的方法 |
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 |