Optical path simulation and experiment for spectral nondestructive detection of pomelo fruit internal quality based on ZEMAX
-
摘要:目的
解决柚果Citrus maxima在无损检测中采集的光谱信息较弱、柚果内部品质检测不准确的问题。
方法采用ZEMAX软件对光谱无损检测光路进行光学仿真，并搭建实物装置进行测试与验证。为进一步验证该光路对柚果内部品质的无损检测效果，采集400~1 100 nm波段132个柚果光谱数据，理化测定其可溶性固形物含量并进行训练建模和测试验证。
结果最优光路参数为：灯珠数量为7、柚果与光源间距为2 cm，探测器与柚果间距为0.5 cm。通过Savitzky-Golay (SG)平滑和标准正态变换(Standard normal variate，SNV)预处理，竞争自适应重加权采样(Competitive adaptive reweighted sampling，CARS)特征选择，建立柚果糖度偏最小二乘回归(Partial least squares regression，PLSR)预测模型，训练集R2和RMSE分别为0.81和0.85，测试集R2和RMSE分别为0.81和0.89。
结论本研究得到的柚果光谱无损检测光路参数可用于柚果内部品质的无损检测，检测效果较好，也为其他大尺寸水果内部无损检测提供了参考。
Abstract:ObjectiveThis study is aimed to address the issue that the spectral information collected in the non-destructive testing of pomelo fruit (Citrus maxima) is weak, and the detected internal quality of pomelo fruit is not accurate.
MethodThe optical simulation of the optical path of spectral nondestructive testing was carried out by using ZEMAX software, and tested and verified by the physical device. In order to further verify the nondestructive detection effect of the optical path on the internal quality of pomelo fruit, the spectral data of 132 pomelo fruits in 400−1 100 nm band were collected, and their soluble solid contents were determined physically and chemically for training modeling and testing verification.
ResultThe optimal parameters of the optical path were as following: The number of lamp beads was 7, the distance between the fruit and the light source was 2 cm, and the distance between the detector and the fruit was 0.5 cm. Through Savitzky-Golay (SG) smoothing and standard normal variate (SNV) preprocessing, competitive adaptive reweighted sampling (CARS) feature selection, a partial least squares regression (PLSR) prediction model was established. The R2 and RMSE of the training set were 0.81 and 0.85, and those of the test set were 0.81 and 0.89, respectively.
ConclusionThe optical path parameters of pomelo fruit obtained in this study are applicable for non-destructive detection of pomelo fruit internal quality, yielding satisfactory detection results. The results provide a reference for the internal non-destructive testing of other large fruits.
-
代谢组学利用现代分析技术对生物体内小分子化合物进行定性定量研究[1]，结合现代生物信息学可实现对代谢物全面地、无偏倚地分析[2]。研究微生物代谢物的变化动态，有助于我们更加了解微生物内部的生理变化。目前代谢组分析技术主要有气相色谱−质谱联用(GC-MS)[3]、液相色谱−质谱联用(LC-MS)[4]和核磁共振(NMR)[5]等。基于GC-MS技术的样品前处理步骤主要包括：淬灭、提取、衍生化和检测。在代谢物淬灭过程中，应尽可能地保证微生物细胞的完整性，避免胞内代谢物外流[6]，可以通过迅速降低温度或者在菌体培养液中加入有机溶剂等手段抑制代谢酶活性。冷甲醇淬灭是普遍使用的方法[7]，它可以在亚秒时间内阻断代谢反应，将胞内外的代谢物区分开[8-9]。常用的提取方法有冷甲醇液氮反复冻融法[10]和75%乙醇煮沸法[11]等。常用的提取液有热乙醇、高氯酸或碱、甲醇和三氯甲烷等。在破碎细胞时，提取液将胞内的代谢物溶解，为了确保提取到稳定可靠的代谢物，可以加入4−羟乙基哌嗪乙磺酸(HEPES)和哌嗪乙磺酸(PIPES)等盐离子进行缓冲[12]。此外超临界流体萃取、微波辅助萃取[13]、固相萃取、固相微萃取等方法也可用于微生物代谢物提取。但至今，仍然没有一种适用于所有微生物的提取方法[14]，不同微生物、不同提取方法提取的代谢物倾向性和分析技术亦不同。因此在每一项新的代谢组学研究中，提取方法一定要有针对性的优化和测试。
近年来，以绿僵菌Metarhizium spp.、白僵菌Beauveria spp.为代表的真菌生物农药已成为生物防治领域的热门研究[15]。在菌种保存和应用过程中，继代培养造成菌株形态变异和毒力减退等退化现象普遍发生，极大地影响了真菌杀虫剂在生产上的应用。继代培养后分生孢子产量与毒力的下降常常伴随着次生代谢产物、代谢路径和代谢网络等层面的变化[16-19]。目前微生物代谢组学研究主要集中于大肠埃希菌Escherichia coli[20]、酿酒酵母Saccharomyces cerevisiae[21]及枯草芽孢杆菌Bacillus subtilis[22-23]等几种模式微生物，并且没有一种方法可以适用于所有微生物的代谢研究。为了获得真实、可靠和全面的代谢组学信息，需要一个准确、稳定的代谢组学检测方法，但目前关于绿僵菌代谢组仅有少量基于LC-MS检测的研究[24]，鲜见基于GC-MS的相关研究。本文以绿僵菌Ma09为研究对象，采用高分辨率的GC-MS技术，对绿僵菌代谢物前处理技术进行优化，以期建立适用于绿僵菌的代谢组检测方法，为获得和解析绿僵菌的“代谢物指纹图谱”提供技术支持，为绿僵菌的深入研究及应用提供参考。
1. 材料与方法
1.1 供试菌株
绿僵菌Ma09，由广东省林业科学研究院提供，绿僵菌斜面保存于4 ℃冰箱备用。转接于PPDA液体培养基，200 r/min、25 ℃条件下培养72 h。
1.2 方法
1.2.1 淬灭
参照Tsuchido等[25]的方法采用D260 nm和D280 nm回收率法评价淬灭剂对细胞膜损伤的程度。核酸和蛋白的特异吸光波长分别为260和280 nm，淬灭后菌液的D260 nm和D280 nm增加，意味着细胞内的核酸和蛋白已渗出到菌液中，说明细胞膜已产生损伤。取72 h菌体培养液10 mL，4 ℃、6000 r/min离心10 min，离心所得沉淀再用9 g/L NaCl溶液润洗1次去除残留培养基。分别用10 mL不同淬灭剂(0 ℃)将菌体重悬，以9 g/L NaCl溶液作为阴性对照、热9 g/L NaCl (70 ℃, 30 min)溶液作为阳性对照，4 ℃、6000 r/min离心10 min，得到上清液，试验重复3次。测定上清液的D260 nm和D280 nm，根据如下公式计算核酸的回收率(η260 nm)和蛋白的回收率(η280 nm)。
$$ \begin{array}{l} {\eta }_{260{\rm{\;nm}}} = \dfrac{{D}_{{260}{\rm{\;nm}}}{(}{{\text{处理}}}{)}-{D}_{{260}{\rm{\;nm}}}{({\text{阴性对照}})}}{{D}_{{260}{\rm{\;nm}}}{({\text{阳性对照}})}-{D}_{{260}{\rm{\;nm}}}{({\text{阴性对照}})}}{ \times 100{\text{%}}}{,} \\ {\eta }_{{280}{\rm{\;nm}}} = \dfrac{{D}_{{280}{\rm{\;nm}}}{(}{{\text{处理}}}{)} - {D}_{{280}{\rm{\;nm}}}{({\text{阴性对照}})}}{{D}_{{280}{\rm{\;nm}}}{({\text{阳性对照}})} - {D}_{{280}{\rm{\;nm}}}{({\text{阴性对照}})}}{ \times 100{\text{%}}}{\text{。}} \end{array}$$ 淬灭剂参考前人的研究[7, 26]，选择甲醇、乙醇和甘油溶液，体积分数分别为20%、40%和60%。添加的缓冲盐为9 g/L的NaCl、NaCN和(NH4)2CO3，10 mmol/L EDTA及10 mmol/L HEPES溶液。
1.2.2 绿僵菌代谢物气相色谱质谱检测
将绿僵菌Ma09菌液以40%(φ)乙醇溶液淬灭，−40 ℃冷甲醇提取后N2吹干，加入20 mg/mL的甲氧基胺盐酸盐吡啶溶液在37 ℃反应60 min，加入含1%(φ)三甲基氯硅烷(TMCS)的N, O−双(三甲基硅烷基)三氟乙酰胺(φ为96%的BSTFA )衍生1 h。采用Agilent 8890+5977B GC/MS系统，毛细管色谱柱为HP-5ms(30 m×0.25 mm×0.25μm) (φ为5%的苯基甲基硅氧烷)非极性弹性石英毛细管柱。质谱条件：电子轰击源(EI)，离子源温度为230 ℃，四极杆温度为150 ℃，接口温度为280 ℃，电子能量为70 eV，调谐方式为标准调谐，质谱扫描方式为全扫描，扫描范围为30~500 aum。使用NIST谱库进行图谱检索，试验重复3次。
气相色谱以纯度大于99.99%(φ)的高纯氦气(Ne)为载气，载气流速1 mL/min。采用不分流进样方式，进样量l μL。进样口温度280 ℃。气相色谱条件如下：
条件A：起始温度50 ℃，保持3 min；以10 ℃/min升至150 ℃，保持5 min；然后以5 ℃/min升至200 ℃，保持5 min；以10 ℃/min升至280 ℃，保持5 min。
条件B：参考明明[27]的GC条件。起始温度70 ℃，保持4 min；以3 ℃/min升至200 ℃；以10 ℃/min升至280 ℃，保持5 min。
条件C：参考王洪彬等[23]的GC条件。初始温度70 ℃，恒温保持2 min；以3 ℃/min升到133 ℃；以2 ℃/min升至200 ℃；再以3 ℃/min升至220 ℃；以5 ℃/min升温至280 ℃。
条件D：参考胡志宏等[28]的GC条件，起始温度65 ℃，保持2 min；以5 ℃/min升至185 ℃；以1 ℃/min升至200 ℃；以15 ℃/min升至280 ℃，保持5 min。
条件E：参考郭刚等[29]的GC条件，起始温度70 ℃，保持4 min；以5 ℃/min升至280 ℃，保持4 min。
1.2.3 提取
热乙醇提取法参考Castrillo等[11]的方法略作修改，淬灭后的绿僵菌沉淀加入5 mL75%(φ)热乙醇溶液(含0.25 mol/L HEPES溶液)，pH为7.5，95 ℃水浴3 min；然后迅速在冰浴中冷却3 min，12000 r/min离心10 min，取上清液在90 ℃条件下旋转蒸发至1 mL。N2吹干后待衍生处理。
热甲醇提取法参照Maharjan等[12]的方法。淬灭后的绿僵菌沉淀中加入2 mL的甲醇溶液[V(甲醇)∶V(H2O)=2∶1]混合均匀后迅速放入70 ℃水中平衡30 min，然后以12000 r/min离心5 min，取上清液N2吹干待后衍生化处理。
冷甲醇提取法参考前人研究[30-32]并作相应修改，淬灭后的绿僵菌沉淀用2 mL预冷的甲醇(−20 ℃)溶解，迅速放入−20 ℃冰箱中保存20 min，冰浴冻融10 min后以12000 r/min离心5 min。上清液转移至另一离心管，取上清液1 mL将沉淀再次悬起，重复上述步骤，2次上清液合并，N2吹干后待衍生化处理。
甲醇氯仿提取法参考Maharjan等[12]的方法。淬灭后的绿僵菌沉淀用1 mL甲醇溶解后，加入1 mL冷三氯甲烷(−40 ℃)涡旋30 s后置于−20 ℃中保存60 min，此过程中每隔15 min涡旋1次，最后12000 r/min离心10 min，取上清液N2吹干待衍生化处理。
氯仿−甲醇−缓冲液提取法参考Maharjan等[12]的方法略作修改。淬灭后的绿僵菌沉淀用5 mL冷三氯甲烷(−40 ℃)溶解后，加入2 mL含有3 mmol/L HEPES和3 mmol/L EDTA的冷缓冲溶液(0 ℃)，混合均匀后12000 r/min离心5 min。收集上层液并在剩余的三氯甲烷相中加入2 mL预冷的甲醇(−40 ℃)和2 mL冷缓冲液，涡旋30 s后10000 r/min离心5 min，收集上层水相并与前一次提取液合并，N2吹干待衍生化处理。
1.2.4 衍生化
代谢物经淬灭和提取后必须进行衍生化处理后才能进行GC-MS检测分析。对衍生时间进行比较，取吹干后的样品，加入20 mg/mL的甲氧基胺盐酸盐吡啶溶液80 μL，剧烈振荡30 s后，在37 ℃条件下反应90 min，反应结束后冷却至室温，然后加入80 μL含1%(φ)TMCS的BSTFA衍生剂，在70 ℃条件下分别衍生0.5、1.0和1.5 h，冷却至室温后进行GC-MS检测，3次重复。
1.3 数据处理
利用Mass Hunter软件对结果进行分析，提取离子峰并用NIST MS数据库进行定性分析，将所有离子峰的峰面积数据归一化，采用SPSS 20.0统计软件进行数据处理及分析，采用Duncan’s法进行差异显著性分析。
2. 结果与分析
2.1 整体水平的细胞损伤及代谢物的渗出
D260 nm和D280 nm是评价细胞渗出物的指标。图1的结果显示，D260 nm最小的淬灭剂为40%(φ)乙醇溶液和20%(φ)甘油溶液，两者差异不显著，但明显低于其他淬灭剂处理的(P<0.05)。D280 nm最小的淬灭剂为60%(φ)甲醇溶液和20%(φ)甘油溶液，且两者间差异不显著，但明显低于其他淬灭剂处理的(P<0.05)。
![]()
图 1 不同淬灭剂处理后菌液的D260 nm和D280 nm1、2、3分别表示体积分数为60%, 40% 和20%的甘油溶液；4、5、6分别表示体积分数为60%, 40% 和20%的甲醇溶液；7、8分别表示体积分数为60%和40%的乙醇溶液；各图中的不同小写字母表示差异显著(P<0.05，Duncan’s法)Figure 1. D260 nm and D280 nm of bacterial fluid with different quenching solutions1, 2 and 3 represent glycerol solutions with volume fractions of 60%, 40% and 20%, respectively; 4, 5 and 6 represent methanol solutions with volume fractions of 60%, 40% and 20%, respectively; 7, 8 represent ethanol solutions with volume fractions of 60% and 40%, respectively; Different lowercase letters in each figure indicate significant differences (P<0.05，Duncan’s test)D260 nm和D280 nm越小说明淬灭过程中的代谢物渗出越少，淬灭效果越好。回收率计算结果(表1)显示，40%(φ)乙醇溶液和20%(φ)甘油溶液2种处理的核酸和蛋白回收率最低。
表 1 不同淬灭剂处理后核酸和蛋白的回收率Table 1. The recovery rates of nucleic acid and protein using different quenching solutions淬灭剂
Quenching solution回收率/% Recovery rate η260 nm η280 nm 40%(φ)乙醇溶液 Ethanol solution 9.63 11.61 20%(φ)甘油溶液 Glycerol solution 14.53 14.68 60%(φ)甲醇溶液 Methanol solution 23.37 23.24 60%(φ)乙醇溶液 Ethanol solution 28.71 35.78 40%(φ)甘油溶液 Glycerol solution 31.55 39.96 60%(φ)甘油溶液 Glycerol solution 52.00 61.23 40%(φ)甲醇溶液 Methanol solution 58.94 56.12 20%(φ)甲醇溶液 Methanol solution 71.03 59.10 2.2 不同离子盐对淬灭效果的影响
以40%(φ)乙醇溶液为空白对照，分别加入不同种类的离子盐，从图2可以看出缓冲盐的加入降低了菌液的D260 nm和D280 nm，减少了菌体细胞的损伤，但是多重比较并未发现几种离子盐之间差异显著[F(D260 nm)=0.893，F(D280 nm)=2.009]。操作中发现，20%(φ)甘油溶液由于黏稠度比较大，不利用后期的浓缩，因此本试验选择40%(φ)乙醇溶液且不加入任何离子盐作为最佳淬灭剂。
![]()
图 2 不同离子盐加入后菌液的D260 nm和D280 nm各图中不同小写字母表示差异显著(P<0.05，Duncan’s法)Figure 2. D260 nm and D280 nm of bacterial fluid added with different ionic saltDifferent lowercase letters in each figure indicate significant differences (P<0.05，Duncan’s test)2.3 绿僵菌代谢物气相色谱质谱条件的确定
图3为5种气相色谱条件的总离子流图。由于条件不同，各处理所用的时间也不同，条件C的总时最长为75 min，条件A仅需46 min。此外出峰时间和出峰数量存在明显差异。
![]()
图 3 不同气相色谱条件检测出的总离子流图A：起始温度50 ℃保持3 min，以10 ℃/min升至150 ℃保持5 min，5 ℃/min升至200 ℃保持5 min，10 ℃/min升至280 ℃保持5 min；B：起始温度70 ℃保持4 min，以3 ℃/min升至200 ℃，10 ℃/min升至280 ℃保持5 min；C：初始温度70 ℃保持2 min，以3 ℃/min至133 ℃，以2 ℃/min升至200 ℃，以3 ℃/min升至220 ℃，以5 ℃/min升至280 ℃；D：起始温度65 ℃保持2 min，以5 ℃/min升至185 ℃，以1 ℃/min升至200 ℃，以15 ℃/min升至280 ℃保持5 min；E：起始温度70 ℃保持4 min，以5 ℃/min升至280 ℃保持4 minFigure 3. Graph of total ion current (TIC) under different chromatographic conditionsA: The initial temperature was 50 ℃ for 3 min, then increased to 150 ℃ for 5 min at 10 ℃/min, then increased to 200 ℃ at 5 ℃/min and maintained for 5 min, and finally increased to 280 ℃ at 10 ℃/min and maintained for 5 min; B: The initial temperature was 70 ℃ for 4 min, then increased to 200 ℃ at 3 ℃ /min and then increased to 280 ℃ at 10 ℃ /min for 5 min; C: The initial temperature was 70 ℃ for 2 min, the temperature rose to 133 ℃ at 3 ℃/min, then rose to 200 ℃ at 2 ℃/min, then rose to 220 ℃ at 3 ℃/min, and finally rose to 280 ℃ at 5 ℃/min; D: The initial temperature was 65 ℃ for 2 min, then increased to 185 ℃ at 5 ℃/min, 200 ℃ at 1 ℃/min, and finally increased to 280 ℃ at 15 ℃/min for 5 min; E: The initial temperature was 70 ℃ for 4 min, and then increased to 280 ℃ at 5 ℃/min and maintained for 4 min根据检测出来的有效峰的数量可以判断出条件A的效果最好，条件A所获得的有效峰为97.33个，显著多于其他几种条件的(P<0.05)。
2.4 提取方法对绿僵菌代谢物提取效果的影响
选用了GC-MS对5种不同方法的提取效果进行分析，根据有效峰数量评判最优方法，各处理获得有效峰数量差异显著(P<0.05)。各处理样本TIC谱图见图4。
![]()
图 4 不同提取方法检测出的总离子流图Figure 4. Graph of total ion current (TIC) using different extraction methods根据有效峰的数量和强度判断，冷甲醇法效果最佳，可以检测出有效峰109.4个，其次是热乙醇法。
2.5 不同衍生物比例和时间的影响
不同衍生时间的总离子流图见图5，以图5中物质数量和峰强度作为判断标准。根据图6有效峰数量可知，衍生0.5 h时大多数物质浓度较低且有多种物质无法检测到，而衍生1.0和1.5 h时检测到的物质数量显著提高了(P<0.05)，衍生1.5 h后可以检测到90.33个有效峰。
![]()
图 5 不同衍生时间的总离子流图(TIC)Figure 5. Graph of total ion current (TIC) at different derivative time![]()
图 6 不同衍生时间的有效峰图中不同小写字母表示差异显著(P<0.05，Duncan’s法)Figure 6. Effective peak at different derivative timeDifferent lowercase letters indicate significant difference (P<0.05，Duncan’s test)2.6 样品前处理最佳方法的确定
取10 mL培养72 h的绿僵菌菌液，40%冷乙醇溶液淬灭后用2 mL冷甲醇提取，离心后将上清液用N2吹干，加入80 μL 20 mg/mL的甲氧基胺盐酸盐吡啶溶液，剧烈振荡30 s后，在37 ℃环境中反应90 min，反应结束后冷却至室温，然后加入80 μL BSTFA(φ为1% TMCS)衍生剂，在70 ℃条件下反应，衍生1.5 h。在气相色谱条件A下进行代谢物检测鉴定。采用该方法进行3次重复测定，根据其保留时间相对标准偏差进行重复性评价，公式如下：
$$ {\text{相对标准偏差}}=\frac{{\text{保留时间标准偏差}}}{{\text{保留时间平均值}}}\times 100{\text{%}} {\text{。}}$$ 从表2可以看出该方法检测保留时间相对标准偏差为0.000002~0.005789，均小于0.05，说明该方法稳定性和重复性较好，且精密度较高，可以满足代谢组学分析的要求。采用优化方法处理的样品经GC-MS分析后共鉴定到50种代谢物。
表 2 绿僵菌代谢物检测最适方法保留时间的重复性Table 2. Repeatability of the retention time for the optimal method of Metarhizium metabolites detection名称
Name分子式
Formula保留时间/s Retention time 相对标准偏差/%
Relative standard deviation重复1 Repeat 1 重复2 Repeat 2 重复3 Repeat 3 硼烷四氢吡啶 C5H8BN 3.1749 3.1769 3.1805 0.000893 N-Difuorophosphoxy-O-trimethylsilylhydroxylamine C3H10F2NO2PSi 3.2892 3.2755 3.2987 0.003547 N,N−二甲基乙醇胺 C4H11NO 3.5455 3.5479 3.5472 0.000348 2−二甲胺基乙硫醇 C4H11NS 4.0018 3.9991 3.9949 0.000869 氨基磺酸 H3NO3S 4.0513 4.0543 4.0528 0.000370 N−乙烯基吡啶溴化铵 C7H8BrN 4.0649 4.0792 4.0990 0.004195 甜菜碱 C5H11NO2 5.1147 5.1380 5.1264 0.002272 顺丁烯二酸二丁基锡 C12H20O4Sn 5.3026 5.2618 5.2822 0.003862 2−甲氨基乙醇 C4H11NO 6.4913 6.4917 6.5010 0.000845 1,2−丙二烯−1,3−二酮 C3O2 7.0261 6.9981 7.0121 0.001996 3,6,9−三噁十一烷二酸 C8H14O7 7.9223 7.9183 7.9203 0.000252 甘露糖胺 C6H13NO5 8.6969 8.6872 8.6921 0.000557 鸟嘌呤 C10H13N5O5 9.4114 9.3677 9.3682 0.002673 1−氧−3−氨基吡嗪 C4H5N3O 9.5307 9.4940 9.5124 0.001929 1,2,3−三唑−4−联苯甲醛 C3H3N3O 10.0717 10.0699 10.0708 0.000089 二乙基二甲基锡烷 C6H16Sn 10.4308 10.4281 10.4295 0.000129 L−高丝氨酸 C6H11NO2 10.4546 10.4630 10.4622 0.000443 1，3−二癸炔 C6H16Sn 10.6731 10.6717 10.6714 0.000085 1−哌嗪乙醇 C8H9BN2 10.9266 10.9298 10.9288 0.000149 甘露醇 C6H10O4 11.6747 11.6409 11.6275 0.002088 亮氨酸 C4H10Si 11.6378 11.6347 11.6393 0.000201 2,5−二甲基苯甲醛 C9H10O 11.6923 11.6955 11.6939 0.000136 异硫氰酸异丁酯 C5H9NS 14.0125 14.0328 14.0531 0.001446 2−甲基 −2−异氰酸丙烷 C5H9NO 14.1314 14.1048 14.2224 0.004357 硼酸三乙酯 C6H15BO3 14.2657 14.2609 14.2807 0.000723 N−二乙氨乙基−端粒酶−丁基−异丙基磷 C10H24NP 14.3573 14.3526 14.3351 0.000815 3−(二乙基硼氧基)− 1−丙硫醇 C7H17BOS 15.8069 15.7658 15.8480 0.002600 2,4−二−3−丁基酚 C14H22O 16.8852 16.8860 16.8856 0.000023 10,12−二十三碳二炔酸 C6H11N3O4 18.0846 18.0822 18.0834 0.000066 9,12−十八烯酸 C8H12N2O2 22.6812 22.8048 22.7430 0.002717 精氨酸 C6H14N4O2 23.5295 23.2586 23.3941 0.005789 L−脯氨酰基−L− 缬氨酸 C10H16N2O2 25.2445 25.2353 25.2399 0.000182 花生四烯酸 C20H32O2 27.6091 27.6176 27.6018 0.000286 5,8,11−二十碳三烯酸 C4H6N6O 27.8460 27.8163 27.9191 0.001898 棕榈酸甲酯 C17H34O2 27.9531 27.9876 27.9704 0.000616 D−α−阿拉伯呱喃糖 C17H42O5Si4 28.1760 28.1761 28.1760 0.000002 癸酸 C10H19AgO2 28.5201 28.5239 28.5220 0.000066 D−盐藻糖醇 C21H54O5Si5 28.5610 28.6237 28.5924 0.001096 肌醇 C24H60O6Si6 33.1030 33.1003 33.1017 0.000040 大麦芽碱 C10H15NO 36.3261 36.5231 36.3994 0.002734 (8, 11−十七二烯)4,5−二氰恶唑 C20H35NO 37.0850 37.0833 37.0842 0.000022 9−十八烯酰胺 C18H35NO 37.8143 37.8053 37.8251 0.000262 2−异氰−1 3−二甲基苯 C9H9N 39.2796 39.4483 39.3640 0.002142 氯化磷酸二乙酯 C6H14ClO3P 39.4663 39.4687 39.4495 0.000265 白桦脂醇 C30H50O2 40.3501 40.3509 40.3545 0.000058 十八烯酸单甘油酯 C21H40O4 41.0059 41.0365 41.0212 0.000372 环己烷 C3H12Si3 41.2616 41.2442 41.2529 0.000210 苯甲酰甘氨酸 C9H9NO3 41.2707 41.2523 41.2615 0.000222 十四烷酸乙酯 C16H32O2 41.3497 41.3536 41.3517 0.000047 十八碳二烯酸甲酯 C19H34O2 42.2349 42.2371 41.3508 0.012185 3. 讨论与结论
代谢物分析的关键要选择高敏感性的检测手段，要求最大程度、最大范围地提取细胞内代谢物，从而更加真实的反映菌体胞内代谢情况。虽然微生物代谢组分析已经成为一种广泛使用的方法，但只使用一个已建立起来的方法进行代谢组分析所能检测出的代谢物覆盖率很低[33]，目前为止还没有一种普遍适用的淬灭方法。另外丝状真菌的生理和形态与细菌、酵母菌有很大的不同，培养物通常高度黏稠且不均匀，液体培养基中有时会形成菌丝团或菌球，很难从发酵过程中获得具有代表性的样品[34]。因此不能采取与细菌、酵母菌相同的方法，需要对代谢物样品制备步骤进行系统优化。
理想的淬灭应具有以下几个条件：灭活新陈代谢的速度比样品中代谢物的转化速度更快；在淬灭步骤中保护样品的完整性；淬灭尽量不引起代谢物的化学和物理性质的变化；样品淬灭以后能够适合接下来的分析步骤；有可靠的重复性[35]。Wittmann等[36]认为微生物突然进入低于−60 ℃的环境中会产生一个冷休克的状态，但是冷甲醇会造成细胞内代谢物的大量泄漏。对于甲醇是否真的能够引起细胞代谢物泄漏的说法不一，Villas等[37]用冷甲醇处理酿酒酵母很少发生泄漏，而Maharjan[12]用冷甲醇处理谷氨酸棒杆菌胞内代谢物则易泄漏。目前大部分研究都选用60%冷甲醇溶液作为淬灭剂，也有选用甘油和乙醇作为淬灭剂的研究。Villas等[37]认为60%(φ)冷甘油−盐水溶液是淬灭酵母和细菌最好的淬灭剂。Spura等[7]发现40%(φ)的乙醇加8 g/L氯化钠溶液(−20 ℃)作为淬灭剂效果要远远好于60%甲醇溶液。本试验比较了甲醇、乙醇和甘油作为淬灭剂的效果，其中，效果最好的是40%(φ)乙醇溶液和20%(φ)甘油溶液，60%(φ)甲醇溶液效果次之。此外对于使用甘油作为淬灭剂也一直存在着争议，Villas-bôas等[26]认为在60%(φ)冷甘油加盐水溶液淬灭酵母和细菌可以很好的解决细胞泄露问题，本试验中发现，甘油虽然可以很好的保护细胞膜的完整性，但由于黏稠度较大，易被色谱柱吸附，影响检测质量，在色谱总离子图(TIC)中前端产生大片覆盖区域，不利于分离。缓冲盐对淬灭效果的影响也存在2种说法。Faijes等[9]认为添加了HEPES和 (NH4)2CO3可以维持菌体内外渗透压平衡，显著减少胞内代谢物的泄露。Canelas等[38]则认为加入缓冲盐会加速代谢物渗出。由此可见，微生物淬灭方法和淬灭剂的选择是一个艰难的过程，需要根据研究对象进行优化评价，建立适合相应微生物淬灭的方法。本试验中在淬灭剂加入几种离子盐后，D260 nm、D260 nm与对照相比较均有小幅度降低，但是多重比较结果表明与对照差异并不显著。这说明针对不同的菌种，是否需要加入缓冲盐也是应该考量的。本试验认为40%乙醇溶液不加入缓冲盐更适合作为绿僵菌代谢组的淬灭剂。
在不改变理化性质的基础上最大限度地将微生物代谢物从样品中提取出来；提取后的样品形式要与分析检测技术的需求一致；尽可能地去除干扰物质；便于浓缩[39]。提取剂的极性决定着最后获得代谢物的极性范围[40]，利用提取液与代谢物的相似相容原理，可以破坏菌体细胞结构，从而提取到代谢物质。Maharjan等[12]使用6种提取剂提取大肠埃希菌代谢物，包括酸、碱、冷甲醇、热甲醇、热乙醇和甲醇+三氯甲烷，结果表明冷甲醇提取到了更多的代谢物。此外，利用热乙醇提取氨基酸和途径中间产物，简单快速且无盐加入，易蒸发和浓缩，但高温对热不稳定的代谢物有破坏作用。Castrillo等[11]用4种提取方法分别提取了酵母菌的代谢组，不同方法产生了明显不同的代谢物。Meyer等[41]对比了6种方法提取金黄色葡萄球菌代谢物，即冷甲醇、冷乙醇、冷甲醇、热乙醇、冷氯仿+水+乙醇和热水，发现冷乙醇效果最好。Canelas等[10]认为甲醇反复冻融法和酸性甲腈乙醇不适用于酵母的代谢物提取。本试验比较了绿僵菌代谢组的5种提取方法，效果最好的还是冷甲醇法，这与大部分报道是一致的。
生物体内的低极性难挥发的代谢物质如氨基酸、脂肪、糖类和核苷酸类等需要GC-MS分析时，必须先进行衍生化。目前常用的衍生化方法为肟化硅烷化衍生法[42]。肟化硅烷化衍生法实际包括2个步骤：首先是甲氧基胺与醛酮类物质的加成反应。但对于种类不同的醛酮物质影响也是不同的。对于还原性糖类物质主要影响因素为反应时间[11]。而对于α−酮酸类物质，pH是主要因素[43]。其次，肟化后需要进行全面的硅烷化反应，常用的硅烷化试剂有N−(叔丁基二甲基硅烷基)−N−甲基三氟乙酰胺(MTBSTFA)、双(三甲基硅烷基)三氟乙酰胺(BSTFA)等，通常会在使用中加入1%(φ)的三甲基氯硅(TMCS)。Kanani等[44]认为BSTFA衍生化时间对回收的代谢物种类产生较大影响，本试验选取了甲氧基胺盐酸盐37 ℃条件下肟化90 min，加入BSTFA(φ为1% TMCS)70 ℃条件下进行硅烷化，衍生时间为1.5 h时效果最好。
气相色谱主要用于代谢物质的分离，质谱则是对代谢物的相对分子质量和分子结构进行鉴定，为了能够更有效地分离样品中的代谢物，对气相色谱升温程序做了具体的优化选择，起始温度过高不利于代谢物的检测，因此本试验中色谱条件的起始温度定为50 ℃，采取程序升温法。程序开始时柱温较低，低沸点物质得到分离，中等沸点物质缓慢移动，随着柱温升高，低沸点和高沸点物质得到分离，最后升温至280 ℃，几乎所有物质都已经被气化分离。
本文对金龟子绿僵菌菌液的前处理技术及气相色谱−质谱联用分析方法进行研究，对前处理过程中的淬灭、提取和衍生时间进行了优化。所建方法对微生物主要代谢物如糖、氨基酸、醇、有机酸、烷烃、酯等均可检出。此方法简单、方便、高效，且稳定性好，重复性良好，该方法的建立有助于开展更深入的代谢机制相关研究，并可为农、林领域开展病原微生物代谢组相关研究提供参考。
-
![]()
图 3 不同灯珠数量和不同柚果与光源间距光学仿真结果
Figure 3. Optical simulation results with different number of lamp beads and distance between pomelo fruit and light source
![]()
图 4 不同灯珠数量和不同柚果与光源间距试验结果
Figure 4. Test results with different number of lamp beads and distance between pomelo fruit and light source
![]()
图 5 柚果与光源间距试验下柚果损伤情况
Figure 5. Damage of pomelo fruit under the test of distance between pomelo fruit and light source
-
[1] HOU B, LI G. HFD: Hierarchical feature detector for stem end of pomelo with transformers[J]. Applied Sciences, 2023, 13(8): 4976. doi: 10.3390/app13084976
[2] 李小婷, 徐世荣, 潘东明, 等. 7种柚子果实品质分析与模糊综合评判[J]. 安徽农业科学, 2016, 44(27): 78-80. [3] FAN S, GUO Z, ZHANG B, et al. Using vis/NIR diffuse transmittance spectroscopy and multivariate analysis to predicate soluble solids content of apple[J]. Food Analytical Methods, 2016, 9(5): 1333-1343. doi: 10.1007/s12161-015-0313-5
[4] 孙潇鹏, 徐赛, 陆华忠. 基于光谱技术与光学仿真的柚果在线检测托盘设计与试验[J]. 食品与机械, 2019, 35(12): 56-62. [5] ITO H, FUKINO-ITO N, HORIE H, et al. Non-destructive detection of physiological disorders in melons using near infrared (NIR) spectroscopy[J]. Acta Horticulturae, 2004(654): 229-234.
[6] MAGWAZA L S, OPARA U L. Analytical methods for determination of sugars and sweetness of horticultural products: A review[J]. Scientia Horticulturae, 2015, 184: 179-192. doi: 10.1016/j.scienta.2015.01.001
[7] 吴爽. 基于高光谱漫透射和神经网络柚子内部品质的无损检测模型研究[D]. 福州: 福建农林大学, 2021. [8] 张雨辰. 蜜柚可食率无损检测方法与在线检测分级装备研究[D]. 杭州: 浙江大学, 2022. [9] 孙潇鹏, 刘灿灿, 陆华忠, 等. 基于可见−近红外透射光谱的蜜柚检测中影响因素分析[J]. 包装与食品机械, 2022, 40(4): 1-7. [10] 徐赛, 陆华忠, 丘广俊, 等. 水果品质无损检测研究进展及应用现状[J]. 广东农业科学, 2020, 47(12): 229-236. [11] 徐赛, 孙潇鹏, 张倩倩. 大型厚皮水果的无损检测技术研究[J]. 农产品质量与安全, 2019(5): 30-35. [12] 张海辉, 陈克涛, 苏东, 等. 基于特征光谱的苹果霉心病无损检测设备设计[J]. 农业工程学报, 2016, 32(18): 255-262. [13] FRASER D G, JORDAN R B, KÜNNEMEYER R, et al. Light distribution inside mandarin fruit during internal quality assessment by NIR spectroscopy[J]. Postharvest Biology and Technology, 2003, 27(2): 185-196. doi: 10.1016/S0925-5214(02)00058-3
[14] SUN C, AERNOUTS B, VAN BEERS R, et al. Simulation of light propagation in citrus fruit using Monte Carlo multi-layered (MCML) method[J]. Journal of Food Engineering, 2021, 291: 110225. doi: 10.1016/j.jfoodeng.2020.110225
[15] QIN J, LU R. Monte Carlo simulation for quantification of light transport features in apples[J]. Computers and Electronics in Agriculture, 2009, 68(1): 44-51. doi: 10.1016/j.compag.2009.04.002
[16] SUN J, KÜNNEMEYER R, MCGLONE A, et al. Investigations of optical geometry and sample positioning in NIRS transmittance for detecting vascular browning in apples[J]. Computers and Electronics in Agriculture, 2018, 155: 32-40. doi: 10.1016/j.compag.2018.09.041
[17] GOBRECHT A, BENDOULA R, ROGER J M, et al. Combining linear polarization spectroscopy and the representative layer theory to measure the beer-lambert law absorbance of highly scattering materials[j]. Analytica Chimica Acta, 2015, 853: 486-494.
[18] 田喜, 陈立平, 王庆艳, 等. 全透射近红外光谱的苹果整果糖度在线检测模型优化[J]. 光谱学与光谱分析, 2022, 42(6): 1907-1914. [19] 王风丽, 马晓英, 侯虎, 等. 基于Zemax软件的光栅光谱仪分辨能力仿真实验[J]. 实验室研究与探索, 2021, 40(11): 111-116. [20] 齐先豪, 韩森, 张凌华, 等. 抑制相干噪声的照明系统仿真与误差分析[J]. 应用光学, 2023, 44(3): 533-538. [21] 郑建锋. 激光导引头中基于四象限探测器的光学设计[J]. 光学技术, 2023, 49(1): 11-16. [22] 高升, 徐建华. 基于近红外光谱的红提内部品质无损检测研究[J]. 食品工业科技, 2022, 43(22): 7-14. [23] 张伏, 曹炜桦, 崔夏华, 等. 基于SG-CARS-IBP的圣女果可溶性固形物可见/近红外光谱无损检测[J]. 光谱学与光谱分析, 2023, 43(3): 737-743. [24] 吴晨凯, 张亮, 沈黄通, 等. 光在柑橘组织中的入射深度和分布情况初探[J]. 光谱学与光谱分析, 2014, 34(3): 601-604. [25] 孙晓林, 周童童, 孙志忠, 等. 光学特性在果蔬品质无损检测中的研究进展[J]. 食品研究与开发, 2022, 43(4): 209-218. [26] 王晨曦, 雷鸣, 雷志勇. 光通量变化对光电探测器信噪比影响研究[J]. 国外电子测量技术, 2019, 38(1): 102-105. [27] 田昊. 蜜柚光学特性与可溶性固形物在线无损检测方法研究[D]. 杭州: 浙江大学, 2022. [28] 葛晴. 近红外差分光谱法无创血糖检测中的光源−探测器距离优选[D]. 天津: 天津大学, 2020. [29] 江山. 单像素计算成像的信噪比及其应用研究[D]. 济南: 山东大学, 2019. [30] 郑少锋, 冯小军, 韩强. 被常见光源聚光灯照射物体温升的分析及在标准中应用[J]. 照明工程学报, 2012, 23(3): 102-105. [31] 第五鹏瑶, 卞希慧, 王姿方, 等. 光谱预处理方法选择研究[J]. 光谱学与光谱分析, 2019, 39(9): 2800-2806. [32] 阚相成, 李耀翔, 王立海, 等. 基于光谱预处理的低温水曲柳原木含水率检测[J]. 中南林业科技大学学报, 2022, 42(11): 154-163. [33] 刘燕德, 施宇, 蔡丽君, 等. 基于CARS算法的脐橙可溶性固形物近红外在线检测[J]. 农业机械学报, 2013, 44(9): 138-144. doi: 10.6041/j.issn.1000-1298.2013.09.025 -
期刊类型引用(12)
1. 谌生，文家燕，高远. 基于改进A*算法和Bezier曲线平滑的单AGV路径规划. 广西科技大学学报. 2025(02): 59-68 . 
百度学术
2. 蔡畅，陈建峰，闫青丽，刘芬. 考虑特殊区域的离散区域划分方法. 系统工程与电子技术. 2024(05): 1656-1664 . 百度学术
3. 衣杰，王锐霖，杜登云，王旭瑞，于复生. 自适应稳定平衡单足机器人设计. 制造业自动化. 2024(04): 158-161 . 百度学术
4. 孙月平，方正，袁必康，孙杰，孟祥汶，汪彦彤，赵德安. 基于FIA~*-APF算法的蟹塘投饵船动态路径规划. 农业工程学报. 2024(09): 137-145 . 百度学术
5. 李务太，宫金良，张彦斐. 基于改进人工势场-PSO算法的农业机器人全区域覆盖研究. 山东理工大学学报(自然科学版). 2024(05): 47-53 . 百度学术
6. 邢慧杰，张晓滨，张宏伟. APF与A*融合的多目标点路径规划算法. 计算机技术与发展. 2024(08): 116-121 . 百度学术
7. 侯宇翔，高焕兵，汪子健，杜传胜. 基于改进跳点搜索法的移动机器人路径规划. 组合机床与自动化加工技术. 2023(03): 54-58 . 百度学术
8. 周哲，欧阳勇. 融合神经网络GRU和RRT算法的机械臂路径规划. 组合机床与自动化加工技术. 2023(09): 7-10 . 百度学术
9. 魏鹏，闫晓静，徐军，杨代斌，袁会珠. 植物保护施药技术创新与装备智能化的研究现状与展望. 现代农药. 2023(05): 1-8+27 . 百度学术
10. 王春平，汪家琪，高金凤. 基于A~*算法和TEB算法的智能小车自主导航避障研究. 计算机与数字工程. 2023(12): 2836-2840+2917 . 百度学术
11. 南江萍. 改进蚁群算法在物流机器人路径规划中的应用实践研究. 科技经济导刊. 2021(19): 13-14 . 百度学术
12. 万远，胡同，贺彪，杨彪. 多条件约束三维动态行人疏散路径优化方法. 火灾科学. 2021(04): 192-201 . 百度学术
其他类型引用(15)
下载:
计量
- 文章访问数: 111
- HTML全文浏览量: 30
- PDF下载量: 19
- 被引次数: 27
