Simulation and experiment of air-blowing precision seed-metering device for Panax notoginseng
-
摘要:目的
为了满足三七Panax notoginseng的机械化种植需求，减少机械式清种过程中对种子造成的损伤，设计了一种适用于播种三七种子的气吹式精密排种器。
方法确定了排种器的主要结构参数，并建立了清种过程中的力学模型。通过建立排种器内部流场模型，运用Fluent软件对不同清种风压条件下流场进行仿真分析，验证了清种风压范围。为了检验仿真确定的风压范围的可行性并寻求最佳工作参数组合，选取合格指数、漏播指数和重播指数作为试验指标，作业速度、种层高度、清种风压作为试验影响因素，采用正交试验方法，对排种器进行了台架试验研究。
结果最优参数组合：作业速度为0.6 m/s、种层高度为90 mm、清种风压为0.5 kPa，此时试验合格指数为90.48，漏播指数为4.24，重播指数为5.28。
结论该气吹式排种器能够满足三七的播种要求，为进行排种器的田间试验提供了参考。通过试验结果与之前仿真分析的过程对比，清种风压变化对于排种器充填性能的影响一致，验证了利用Fluent模拟确定清种风压的可行性。
Abstract:ObjectiveA novel kind of air-blowing precision seed-metering device was designed in order to fulfill the mechanized planting needs of Panax notoginseng and reduce seed damage in the process of mechanical seeding.
MethodThe main structural parameters of the seed-metering device were determined. The mechanical model in seed clearing process was established. The internal flow field model of the seed-metering device was established, Fluent software was used to simulate the flow field under different air-blowing pressure, and the range of air-blowing pressure for seed clearing was verified. We used orthogonal design and performed bench test to further examine the feasibility of the range of air-blowing pressure and find out the optimal combination of working parameters. Eligible index, leakage sowing index and repeat sowing index were selected as experimental indexes. Operating speed, seed layer height and air-blowing pressure were selected as the influencing factors in the test.
ResultThe optimal combination of parameters was the operating speed of 0.6 m/s, seed layer height of 90 mm and air-blowing pressure of 0.5 kPa. Under these conditions, the qualified index was 90.48, the leakage sowing index was 4.24, and the repeat sowing index was 5.28.
ConclusionThe air-blowing seed-metering device can meet the requirements of Panax notoginseng seeding. This study provides a theoretical basis for the field experiment of the seed-metering device. Comparison between the test results and the simulation analysis indicates that the influence of air-blowing pressure changes on the filling performance of the seed-metering device is consistent, and Fluent simulation is feasible for determining the air-blowing pressure for seed clearing.
-
代谢组学利用现代分析技术对生物体内小分子化合物进行定性定量研究[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 ℃，几乎所有物质都已经被气化分离。
本文对金龟子绿僵菌菌液的前处理技术及气相色谱−质谱联用分析方法进行研究，对前处理过程中的淬灭、提取和衍生时间进行了优化。所建方法对微生物主要代谢物如糖、氨基酸、醇、有机酸、烷烃、酯等均可检出。此方法简单、方便、高效，且稳定性好，重复性良好，该方法的建立有助于开展更深入的代谢机制相关研究，并可为农、林领域开展病原微生物代谢组相关研究提供参考。
-
![]()
图 1 三七气吹式排种器结构示意图
1：排种轴，2：排种轮，3：壳体，4：型孔，5：清种气嘴，6：风管接口，7：进种管，8：种层高度调节板，9：充种室，10：投种口
Figure 1. Structure diagram of air-blowing seed-metering device for Panax notoginseng
1：Seeding shaft，2：Seeding wheel，3：Shell，4：Type hole，5：Seed clearing nozzle，6：Duct interface，7：Seed tube，8：Seed height adjustment plate，9：Seed filling room，10：Seeding port
![]()
图 2 排种轮结构参数示意图
D：排种轮直径，ω：排种轮角速度，ξ：充种角，τ：清种角，ψ：护种角，φ：投种角
Figure 2. Schematic of structural parameters of seeding disc
D：Diameter of the seeding wheel，ω：Angular velocity of the seeding wheel，ξ：Seed filling angle，τ：Seed clearing angle，ψ：Seed protecting angle，φ：Seed throwing angle
![]()
图 3 型孔结构示意图
d1：型孔上端直径，d2：型孔下端直径，h：型孔深度
Figure 3. Structure diagram of type hole
d1：Diameter of the upper end of type hole, d2：Diameter of the lower end of type hole, h：Depth of type hole
![]()
图 4 清种区种子受力简化模型
Fq：种子受到的清种压力，Fl：惯性离心力，FN：型孔侧面对种子的支持力，Ff：种子与型孔侧面的摩擦力，α：压力Fq与x轴的夹角，β：y轴与水平面夹角，δ：摩擦力Ff与y轴的夹角，ε：型孔锥角，ω：排种轮角速度，G：种子的重力，O：排种轮转动中心，1、2：种子
Figure 4. Simplified model of force on the seed in seed clearing zone
Fq: Air-blowing force on the seed, Fl: Inertial centrifugal force, FN: Supporting force on the side of type hole, Ff: Friction force between the seed and the side of type hole, α: Angle between Fq and x-axis, β: Angle between y-axis and the horizontal plane, δ: Angle between Ff and y-axis, ε: Cone angle of the type hole, ω: Angular velocity of the seeding wheel, G: Seed gravity, O: Rotation center of the seeding wheel, 1 and 2: Seed
![]()
图 5 不同清种风压下排种器内部流场模拟
Figure 5. Simulation of internal flow field in the seed-metering device under different air-blowing pressure
![]()
图 6 试验装置组成图
1：气力管道，2：排种器，3：种床带，4：三七种子
Figure 6. Composition diagram of experimental device
1：Pneumatic pipeline，2：Seed-metering device，3：Seed bed belt，4：Panax notoginseng seed
表 1 排种性能试验因素与水平
Table 1 Factors and levels of seeding performance test
水平
Level作业速度/(m·s–1)
Operating
speed种层高度/mm
Seed layer
height清种风压/kPa
Air-blowing pressure1 0.4 50 0.4 2 0.6 70 0.5 3 0.8 90 0.6 4 1.0 110 0.7 
下载: 导出CSV
表 2 排种性能试验设计与结果
Table 2 Design and result of seeding performance test
试验序号
Test No.因素 Factor 指标 Index 作业速度
Operating speed (A)种层高度
Seed layer height (B)清种风压
Air-blowing pressure (C)合格指数
Eligible index漏播指数
Leakage sowing index重播指数
Repeat sowing index1 1 1 1 57.23 37.03 5.74 2 1 2 2 76.49 18.52 4.99 3 1 3 3 87.91 3.42 8.67 4 1 4 4 68.63 12.16 19.21 5 2 1 2 79.37 13.56 7.07 6 2 2 1 79.43 14.74 5.83 7 2 3 4 89.34 4.24 6.42 8 2 4 3 82.52 3.53 13.95 9 3 1 3 66.71 28.16 5.13 10 3 2 4 70.45 27.32 2.23 11 3 3 1 76.37 7.95 15.68 12 3 4 2 69.72 3.34 26.94 13 4 1 4 49.67 36.83 13.5 14 4 2 3 67.52 23.95 8.53 15 4 3 2 79.47 3.39 17.14 16 4 4 1 55.45 8.12 36.43
下载: 导出CSV
表 3 排种性能试验结果的极差分析1)
Table 3 Range analysis of seeding performance test result
指标
Index项目
ItemA B C 合格指数
Eligible
indexk1 72.565 63.245 67.120 k2 82.665 73.472 76.263 k3 70.813 83.273 76.165 k4 63.028 69.080 69.523 R 19.637 20.028 9.143 漏播指数
Leakage
sowing indexk1 17.782 28.895 16.960 k2 9.018 21.133 9.703 k3 16.693 4.750 14.765 k4 18.073 6.787 20.137 R 9.055 24.145 10.434 重播指数
Repeat
sowing indexk1 9.652 7.860 15.920 k2 8.317 5.395 14.035 k3 12.495 11.977 9.070 k4 18.900 24.133 10.340 R 10.583 18.738 6.850 1) A、B 和 C 分别表示作业速度、种层高度和清种风压
1) A, B and C represent operating speed, seed layer height and air blowing pressure, respectively
下载: 导出CSV
表 4 排种性能试验结果的方差分析1)
Table 4 Variance analysis of seeding performance test result
指标
Index变异来源
Variation sourceSS DF F 合格指数
Eligible
indexA 782.77 3 61.151** B 856.51 3 66.912** C 260.73 3 20.369** R 25.60 6 漏播指数
Leakage
sowing indexA 220.90 3 5.959* B 1 610.30 3 43.437** C 230.97 3 6.230* R 74.14 6 重播指数
Repeat
sowing indexA 265.84 3 11.454** B 829.99 3 35.762** C 121.53 3 5.236* R 46.42 6 1) A、B、C 和 R 分别表示作业速度、种层高度、清种风压和残差；“*” 和 “**” 分别影响显著和极显著
1) A, B, C and R represent operating speed, seed layer height, air blowing pressure and residual, respectively; “*” and “**” represent significant effects at 0.05 and 0.01 levels, respectively
下载: 导出CSV
-
[1] 崔秀明. 依靠科技进步促进三七产业创新发展[J]. 中国现代中药, 2018, 20(3): 247-252. [2] 张维. 气吹式大豆精量播种机的设计[J]. 中国农机化学报, 2014, 35(4): 6-8. [3] 李洪昌, 高芳, 赵湛. 国内外精密排种器研究现状与发展趋势[J]. 中国农机化学报, 2014, 35(2): 12-16. [4] 史嵩. 气压组合孔式玉米精量排种器设计与试验研究[D]. 北京: 中国农业大学, 2015. [5] ZULIN Z, UPADHVAVA S K. Hydropneumatic seeder for primed seed[J]. Trans ASABE, 1998, 41(2): 307-314. doi: 10.13031/2013.17175
[6] 张泽平, 马成林. 精播排种器及排种理论研究进展[J]. 吉林工业大学学报, 1995, 25(4): 112-117. [7] 张守勤, 王成和. 气力轮式排种器型孔的流场及作用[J]. 农业机械学报, 1991, 22(4): 26-31. [8] 心男. 基于EDEM-FLUENT耦合的气吹式排种器工作过程仿真分析[D]. 长春: 吉林大学, 2013. [9] 马成林. 气吹排种器充填原理的研究[J]. 农业机械学报, 1981, 12(4): 1-13. [10] 胡树荣, 马成林, 李慧珍, 等. 气吹式排种器的排种频率和气流速度对排种质量影响的研究[J]. 吉林大学学报(工学版), 1981(4): 26-35. [11] 刘佳, 崔涛, 张东兴, 等. 气吹式精密排种器工作压力试验研究[J]. 农业工程学报, 2011, 27(12): 18-22. doi: 10.3969/j.issn.1002-6819.2011.12.004 [12] 崔涛, 韩丹丹, 殷小伟,等. 内充气吹式玉米精量排种器设计与试验[J]. 农业工程学报, 2017, 33(1): 8-16. [13] 孙裕晶, 马成林, 张勇智. 基于均匀设计的精密排种器结构优化方法[J]. 吉林大学学报(工学版), 2004, 34(4): 569-572. [14] 廖庆喜, 杨波, 李旭, 等. 内充气吹式油菜精量排种器气室流场仿真与试验[J]. 农业机械学报, 2012, 43(4): 51-54. doi: 10.6041/j.issn.1000-1298.2012.04.011 [15] 雷小龙, 廖宜涛, 张闻宇, 等. 油麦兼用气送式集排器输种管道气固两相流仿真与试验[J]. 农业机械学报, 2017, 48(3): 57-68. [16] 李成华, 高玉芝. 倾斜圆盘气吹式播种器设计及试验研究[J]. 沈阳农业大学学报, 2007, 38(6): 830-836. doi: 10.3969/j.issn.1000-1700.2007.06.014 [17] 胡靖明. 气吹式小粒种子精量排种器气室仿真[J]. 林业机械与木工设备, 2016, 44(7): 12-14. [18] 张波屏. 播种机械设计原理[M]. 北京: 机械工业出版社, 1982. [19] ORTH R J, MARION S R, GRANGER S, et al. Evaluation of a mechanical seed planter for transplanting Zostera marina (eelgrass) seeds[J]. Aquat Bot, 2009, 90: 204-208. doi: 10.1016/j.aquabot.2008.07.004
[20] 张德文, 李林, 王惠民. 精密播种机械[M]. 北京: 农业出版社, 1982. [21] 胡树荣, 马成林, 李慧珍, 等. 气吹式排种器锥孔的结构参数对排种质量影响的研究[J]. 农业机械学报, 1981, 12(3): 21-31. [22] 任闯, 高筱钧, 苏微, 等. 三七种子的物理机械特性试验[J]. 湖南农业大学学报(自然科学版), 2015, 41(1): 109-112. [23] 左春柽, 马成林, 张守勤, 等. 新型气力精密排种器的空气动力学原理[J]. 农业工程学报, 1997, 13(3): 110-114. doi: 10.3321/j.issn:1002-6819.1997.03.023 [24] 欧阳洁, 李静海, 崔俊芝. 颗粒轨道模型中相间耦合关系及曳力计算的研究[J]. 动力工程, 2004, 24(6): 857-862. doi: 10.3321/j.issn:1000-6761.2004.06.021 [25] 刘素. 气吹式三七精密排种器设计与试验研究[D]. 昆明: 昆明理工大学, 2018. [26] 全国农业机械标准化技术委员会. GB/T6973-2005单粒(精密)播种机试验方法[S]. 北京: 中国标准出版社, 2006. [27] 袁志发, 周静芋. 试验设计与分析[M]. 北京: 高等教育出版社, 2000. -
期刊类型引用(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)
计量
- 文章访问数: 1519
- HTML全文浏览量: 1
- PDF下载量: 1653
- 被引次数: 27
