Calibration of contact parameters of ‘Luli’ apple for simulation based on discrete element method
-
摘要:目的
构建基于离散元法的苹果颗粒模型，并精确标定其仿真接触参数，为振动收获中接近承接或采后处理装置的设计优化提供理论依据。
方法采用球形颗粒组成的方式构建‘鲁丽’苹果离散元模型，通过对比分析确定最合适的球形颗粒半径。采用试验与仿真相结合的方法，确定恢复系数和摩擦系数等接触参数。通过方差分析评估苹果接触部位、果实质量和材料类型对接触参数的影响。调整仿真参数以复现不同试验条件，结合数据拟合得到参数方程，并利用试验数据进行验证。
结果苹果的接触部位对其恢复系数的影响不显著。确定了颗粒半径为2 mm的苹果离散元模型，并以此为准，标定了苹果与特高密度和高密度泡棉的恢复系数、静摩擦系数、滚动摩擦系数分别为0.61和0.47、0.46和0.61、
0.0166 和0.0288 ，而苹果相互之间的对应参数分别为0.65、0.42和0.032 0。结论成功构建了苹果离散元模型，通过无底圆筒提升试验验证了标定参数的有效性。
Abstract:ObjectiveThe aim of this study was to construct a discrete element model of apples and precisely calibrate the corresponding contact parameters, providing a theoretical basis for the design and optimization of close-range catching mechanisms in apple vibration harvesting or post-harvest processing devices.
MethodThe discrete element model of ‘Luli’ apples was constructed using a spherical particle bonding method, and the optimal particle radius was identified through comparative analysis. A combined method of experimental testing and simulation was adopted to determine the contact parameters, such as the restitution coefficient and friction coefficients. Analysis of variance (ANOVA) was employed to evaluate the effects of collision zone, fruit mass, and foam type on the contact parameters. Various test conditions were simulated to obtain data, which were then used to fit parameter equations. The parameter equations were subsequently verified using experimental data.
ResultThe collision zone had no significant effect on the coefficient of restitution. A discrete element model witha 2 mm particle radius was established. The calibrated restitution coefficients, static and rolling friction coefficients between the apple particle models and the super-high and high-density foams were 0.61 and 0.47, 0.46 and 0.61,
0.0166 and0.0288 , respectively. The corresponding values for interactions between apple particle models were 0.65, 0.42 and 0.032 0.ConclusionThe discrete element model of apple is successfully constructed, and the effectiveness of the calibrated parameters is verified through a bottomless cylinder lifting experiment.
-
Keywords:
- Mechanics /
- Discrete element method /
- Apple particle model /
- Contact parameter /
- Virtual calibration
-
桑树是一种药食同源的多年生木本植物，广泛应用于医药、食品、畜牧养殖等各个行业中[1-3]。桑树病害的发生一直制约着蚕桑产业的健康稳定发展，其中，桑树青枯病作为一种毁灭性的细菌性检疫病害，严重影响了桑树的产量和质量[4]。1969年广东省顺德市首次报道了桑树青枯病的发生[5]，随后传播至广东大部分的桑树种植区。至今，桑青枯病仍在广东、广西等地的蚕桑主产区流行发生，并在全国范围内的桑园种植地多有发生报道[6-7]。桑树青枯病的病原菌为茄科雷尔氏菌Ralstonia solanacearum，茄科雷尔氏菌的寄主范围广，能感染桑树在内的450余种植物的根或茎，根据其对不同寄主致病性差异的分类报道，侵染桑树的茄科雷尔氏菌多为5号生理小种[8-9]。由于茄科雷尔氏菌种类多，存在不同的致病变种，且致病机理复杂又危害严重，被世界各国的检疫部门列入植物生产和出口贸易的重点检疫对象[10-12]。
建立高效快速的桑树青枯病病原茄科雷尔氏菌检测技术，不仅对桑树苗木的流通和检疫提供技术支持、为桑树产业的健康稳定发展保驾护航，还对各种作物青枯病的防控有重要的参考价值。目前，有关青枯病的病原菌茄科雷尔氏菌的分子检测技术主要有各种的聚合酶链式反应(PCR)法[13-14]和环介导等温扩增 (Loop-mediated isothermal amplification，LAMP)技术[15]。然而这2种检测技术在实际应用中均存在一定的局限性，如PCR法检测时间相对较长，灵敏度低，在桑树青枯病的检测中效果不明显[16]；而LAMP技术虽有检测灵敏度高、反应时间短等优点，但环境的气溶胶往往影响检测结果的可靠性或引起假阳性[17-18]。等温多自配引发扩增 (Isothermal multiple self-matching-initiated amplification，IMSA)技术是2015年马学军团队发明的一种新型核酸扩增技术[19-20]，该技术主要针对靶基因设计6条引物，并能特异性地识别7个基因位点，与LAMP和荧光定量PCR相比具有更强的特异性及更高的灵敏度，在食品安全、动物健康养殖、人畜共患病病原的检测应用中证实是可靠的[21-24]，但有关IMSA检测植物病原菌特别是桑树青枯病病原的报道较少。
本研究拟以茄科雷尔氏菌果胶裂解酶基因(Pectate lyase gene)为靶标，尝试融合IMSA引物设计策略和LAMP扩增反应的技术体系，建立一种快速高效的茄科雷尔氏菌IMSA-LAMP检测技术，并对其反应体系和实用性等进行优化，以期为桑树细菌病检测、桑树青枯病检疫和桑树种苗的流通提供新的检测技术和手段。
1. 材料与方法
1.1 供试材料
2020―2021年从广东、广西两地的桑园采集了疑似桑青枯病症状的茎及根共24份，如表1所示。
表 1 桑树青枯病田间病样Table 1. Field samples of mulberry bacterial wilt地区 District 坐标 Coordinate 品种 Cultivar 数量 Quantity 广东阳山 Yangshan, Guangdong E112.641°，N24.465° 伦教40 Lunjiao 40 2 广东罗定 Luoding, Guangdong E111.570°，N22.768° 伦教40 Lunjiao 40 2 广东广州 Guangzhou, Guangdong E113.359°，N23.170° 伦教40 Lunjiao 40 2 广东英德 Yingde, Guangdong E113.454°，N24.179° 伦教40 Lunjiao 40 3 广西都安 Du’an, Guangxi E108.105°，N23.931° 粤桑11 Yuesang 11 3 广西来宾 Laibin, Guangxi E109.234°，N23.727° 粤桑11 Yuesang 11 4 广西环江 Huanjiang, Guangxi E108.258°，N24.826° 农桑14 Yuesang 11 2 广西忻城 Xincheng, Guangxi E108.665°，N24.066° 粤桑11 Yuesang 11 2 广西宜州 Yizhou, Guangxi E108.636°，N24.485° 粤桑11 Yuesang 11 4 桑青枯病病原菌茄科雷尔氏菌菌株YZqk1、YLqk1(GDMCC No.：1.1615)、XZqk2(GDMCC No.：1.1619)、XCqk4(GDMCC No.：1.1616)、LZqk5(GDMCC No.：1.1617)、YDqk6(GDMCC No.：1.1620)和其他桑源细菌菌株：阴沟肠杆菌XCYG-001(GDMCC No.：1.1600)、克雷伯氏菌LCKL-001(GDMCC No.：1.1602)、菠萝泛菌LCFJ-001(GDMCC No.：1.1601)、铜绿假单胞菌YD-001(GDMCC No.：60613)、丁香假单胞菌ZJDX-003(GDMCC No.：61844)、多黏类芽孢杆菌YD-002(GDMCC No.：61095)、短小杆菌YDDX-001、芽孢杆菌YD-003，均由华南农业大学亚太地区蚕桑培训中心实验室分离鉴定，含GDMCC No.的菌种由广东省菌种保藏中心保藏。
1.2 供试菌株培养
30 ℃条件下，将茄科雷尔氏菌菌株置于TTC琼脂培养基[25]培养5 d，将非茄科雷尔氏菌菌株接种于LB琼脂培养基[25]培养2 d。
1.3 细菌、疑似桑树青枯病病样及健康桑树DNA提取
供试细菌根据Ezup柱式细菌基因组DNA抽提试剂盒[生工生物工程(上海)股份有限公司，以下简称为“生工”]说明书提取总DNA。供试病样及健康桑树样本根据Ezup柱式植物总DNA抽提试剂盒(生工)说明书提取总DNA。DNA提取后，利用超微量分光光度计(NDlife)进行浓度测定，用TE缓冲液稀释到10 ng/μL左右，–20 ℃保存备用，作为后续检测的模板。
1.4 IMSA-LAMP引物设计与合成
IMSA-LAMP引物设计，如图1A所示，在IMSA-LAMP检测中使用了6种引物，包括2种茎引物(SteF和SteR)和2对嵌套杂交引物(2个外引物DsF和DsR和2个内引物FIT和RIT)；引物特异性识别靶基因的7个不同区域，从5′端标记为F3、F2、F1、T、R1c、R2c和R3c；DsF和DsR引物由F3和R3以及F1c和R1c序列组成，FIT和RIT引物由F2和R2以及Tc和T序列组成，SteF和SteR引物分别是R1c和F1c序列[26]。IMSA-LAMP扩增原理如图1B所示，为便于说明，DNA合成从DsF开始，FIT被设置为起始过程(DNA合成以类似方式在DsR和RIT中进行)；带箭头的水平直线代表底物延长方向，带箭头的斜线代表与目标位点退火的引物，带箭头的弧线表示2个区域的自动匹配。在该步骤中，生成了4种不同长度的基本自匹配结构(SMS-1~SMS-4)[26]。根据NCBI报道的茄科雷尔氏菌果胶裂解酶基因(WP_197360060.1)序列，利用引物在线设计软件( http://primerexplorer.jp/elamp5.0.0/index.html)设计之后，根据引物可能存在的二聚体和相近的退火温度，调整并设计出6套IMSA-LAMP检测的引物组。所有引物均由擎科生物有限公司合成，引物均用HPLC法纯化，引物序列见表2。
![]()
表 2 用于扩增茄科雷尔氏菌果胶裂解酶基因的引物序列Table 2. Primer sequences used for amplifying pectate lyase genes of Ralstonia solanacearum引物组编号
Primer mix number引物名称
Primer name引物序列(5′→3′)
Primer sequence (5′→3′)长度/bp
Length1 DsF-1 GTGTCTCAGTCCCAGTGTGGCCTAGTTGGTGGGGTAAAGGC 41 DsR-1 GACAATGGGGGCAACCCTGATCCCAGGTATTACCCAGTGCGA 42 FIT-1 CTGCCTCCCGTAGGAGCCTGCCAAGGCGAACGATCAGTAG 40 RIT-1 CAGGCTCCTACGGGAGGCAGTCTTTCCGGACAAAAGTGCT 40 SteF-1 GACAATGGGGGCAACCCTGATC 22 SteR-1 GTGTCTCAGTCCCAGTGTGGC 21 2 DsF-2 GTACTGCGCGCTGGAGCTGGGAACAACCATTTCCAGACGC 40 DsR-2 GGAACATGCTGAACCAGCCGGGGCCGTTTGCACCGAACC 39 FIT-2 AGCTTGCTCGCTTTGAGCTGCCCTCGACGGGGAATCAC 38 RIT-2 GCAGCTCAAAGCGAGCAAGCTGCCGCCATGATGACCG 37 SteF-2 GGAACATGCTGAACCAGCCGGG 22 SteR-2 GTACTGCGCGCTGGAGCTGG 20 3 DsF-3 GTACTGCGCGCTGGAGCTGGACCATTTCCAGACGCC 36 DsR-3 GGAACATGCTGAACCAGCCGGGGCCGTTTGCACCGAACC 39 FIT-3 GAACAGCTTGCTCGCTTTGAGCTCTCGACGGGGAATCACGA 41 RIT-3 AGCTCAAAGCGAGCAAGCTGTTCGCCGCCATGATGACCG 39 SteF-3 GGAACATGCTGAACCAGCCGGG 22 SteR-3 GTACTGCGCGCTGGAGCTGG 20 4 DsF-4 AACATGCTGAACCAGCCGGGCGCCGTTTGCACCGAACC 38 DsR-4 GTACTGCGCGCTGGAGCTGGACAACCATTTCCAGACGCC 39 FIT-4 GAACAGCTTGCTCGCTTTGAGCTTCGACGGGGAATCACGA 40 RIT-4 AGCTCAAAGCGAGCAAGCTGTTCCGAGCCGCCATGATGAC 40 SteF-4 AACATGCTGAACCAGCCGGGC 21 SteR-4 GTACTGCGCGCTGGAGCTGG 20 5 DsF-5 AGCAGTTGTACTGCGCGCTGGCCATTTCCAGACGCCCTC 39 DsR-5 CGACAACGCCGGCCAGAAGGCCGTTTGCACCGAAC 35 FIT-5 CGAACAGCTTGCTCGCTTTGAGGAATCACGATGCGGGTTCC 41 RIT-5 TCAAAGCGAGCAAGCTGTTCGTCGAGCCGCCATGATGAC 39 SteF-5 CGACAACGCCGGCCAGAA 18 SteR-5 AGCAGTTGTACTGCGCGCTGG 21 6 DsF-6 AGCAGTTGTACTGCGCGCTGCCATTTCCAGACGCCCTC 38 DsR-6 CGACAACGCCGGCCAGAAGGCCGTTTGCACCGAAC 35 FIT-6 CGAACAGCTTGCTCGCTTTGAGAATCACGATGCGGGTTCC 40 RIT-6 TCAAAGCGAGCAAGCTGTTCGGGAATCGAGCCGCCATGATGAC 43 SteF-6 CGACAACGCCGGCCAGAA 18 SteR-6 AGCAGTTGTACTGCGCGCTG 20 1.5 IMSA-LAMP反应体系与反应条件
以提取的供试细菌、疑似桑树青枯病病样及健康桑树样本DNA为模板，对茄科雷尔氏菌果胶裂解酶基因片段进行荧光定量PCR特异性扩增。扩增体系：2× LAMP Master Mix 12.5 μL(生工)，DNA Polymerase 0.5 μL(生工)，Sterilized ddH2O 8.0 μL(生工)，5×103 mmol SYTOTM 9核酸染料1 μL(Thermo Fisher Scientific)，引物Mix 1.0 μL(外引物、茎引物、内引物体积比为1∶4∶8)，DNA模板2.0 μL，共25.0 μL。反应程序：63 ℃ 1 min，80 ℃ 1 min；将63 ℃ 1 min作为1个循环，设计45个循环。
1.6 IMSA-LAMP引物特异性
为了验证“1.4”筛选的引物检测桑树青枯病病菌的特异性，以桑源茄科雷尔氏菌及非茄科雷尔氏菌菌株进行IMSA-LAMP引物特异性试验，以健康桑树DNA作为阴性对照、无菌水为空白对照、茄科雷尔氏菌DNA为阳性对照，反应条件同“1.5”。
1.7 IMSA-LAMP不同引物浓度比例优化
为保证茄科雷尔氏菌DNA能在IMSA-LAMP体系中达到最佳扩增效果，在上述“1.5”反应扩增体系基础上，对引物比例进行筛选。以10 ng/μL的茄科雷尔氏菌YZqk1 DNA为模板，以“1.6”最佳引物组为扩增引物，设置4组不同的引物浓度比例，即外引物(100 μmol/L)、茎引物(100 μmol/L)、内引物(100 μmol/L)体积比依次为1∶4∶8、1∶2∶4、1∶1∶2、1∶8∶4，见表3。使用荧光定量PCR仪进行IMSA-LAMP扩增，根据起始扩增时间及荧光强度等选择最佳引物浓度比例。
表 3 外、茎、内引物体积配比Table 3. Volume ratio of outer, stem and inner primersμL 外、茎、内引物体积比
Volume ratio of primers
(Outer∶Stem∶Inner)外引物
Outer primer茎引物
Stem primer内引物
Inner primer无菌水
Sterilized ddH2O总和
TotalDsF DsR SteF SteR FIT RIT 1∶4∶8 2 2 8 8 16 16 48 100 1∶2∶4 2 2 4 4 8 8 72 100 1∶1∶2 2 2 2 2 4 4 84 100 1∶8∶4 2 2 16 16 8 8 48 100 1.8 IMSA-LAMP反应温度优化
以10 ng/μL的茄科雷尔氏菌YZqk1 DNA为模板，在61.5、62.0、62.5、63.0、63.5、64.0、64.5和65.0 ℃ 8个反应温度下，通过荧光定量PCR仪进行IMSA-LAMP扩增，反应体系见“1.5”。通过比较不同反应温度的荧光扩增曲线，获得最佳反应温度。
1.9 IMSA-LAMP检测法与PCR法的灵敏度对比
将茄科雷尔氏菌YZqk1纯培养后置于25 mL LB肉汤中扩增培养，于28 ℃条件下培养48 h，用稀释涂布计数方法将菌液浓度稀释成1×106 CFU/mL，取1 mL于离心机中12 000 r/min离心4 min，取沉淀，根据Ezup柱式细菌基因组DNA抽提试剂盒说明书提取总DNA。利用超微量分光光度计进行浓度测定，按10倍稀释法用TE缓冲液将DNA梯度稀释7次，具体计算公式：稀释后DNA浓度=稀释前DNA浓度÷10稀释次数，将每个梯度作为模板进行IMSA-LAMP检测，反应体系见“1.5”。选择优化后的引物、引物浓度比例及反应温度，根据荧光扩增曲线判定检测结果；同时，采用Primer Premier 5.0软件对茄科雷尔氏菌果胶裂解酶基因设计1对PCR检测引物：PL-F(AGCTTGCTCGCTTTGAGCTGCTGCAAGACC)/ PL-R(TGATCGCCCGAACAACCATTTCCAGACGCC)，拟扩增的目的片段长度为135 bp，并利用该引物对模板进行PCR检测，反应体系为2× Taq PCR Master Mix 12.5 μL(生工)，10 μmol/L引物PL-F、PL-R各1 μL，DNA模板2 μL，Sterilized ddH2O补足至25 μL；扩增条件为94 ℃预变性5 min；94 ℃变性30 s、55 ℃退火30 s、72 ℃延伸30 s，30个循环；72 ℃最终延伸8 min，于4 ℃保存。PCR扩增产物取5 μL经12 g/L琼脂糖凝胶电泳检测。
1.10 IMSA-LAMP检测方法的验证
将广东、广西地区收集的24份疑似桑树青枯病病样(表1)提取DNA后，进行IMSA-LAMP检测，并取1份健康桑树样本为阴性对照、茄科雷尔氏菌YZqk1 DNA为阳性对照。反应体系见“1.5”，其中引物、引物浓度比例及反应温度均选择优化后的结果，并通过荧光扩增曲线进行判读；同时利用“1.9”的PCR反应体系和参数进行同步平行检测。
2. 结果与分析
2.1 IMSA-LAMP引物筛选
将本研究设计的6组引物(表2)进行IMSA-LAMP检测，检测结果如图2所示。引物组1~6均出现了扩增曲线，引物组6于反应20 min时最先出现上升趋势，35 min后反应进入平缓期；而其他引物组均在反应25 min后才出现扩增。因此，选择引物组6作为IMSA-LAMP法最佳引物组；最佳引物组6的6条引物结合位置及7个特异识别位点，如图3所示。
![]()
图 2 IMSA-LAMP 6组引物对茄科雷尔氏菌的检测结果Figure 2. Detection results of Ralstonia solanacearum under six primer groups of IMSA-LAMP2.2 IMSA-LAMP引物特异性试验
利用IMSA-LAMP对6株茄科雷尔氏菌及8株非茄科雷尔氏菌DNA进行检测，结果如图4所示。6株茄科雷尔氏菌(YZqk1、YLqk1、XZqk2、XCqk4、LZqk5、YDqk6)均在45 min内出现了“S”扩增曲线(图4：曲线1~6)。其他8株非茄科雷尔氏菌、无菌水及健康桑树DNA均未出现扩增曲线(图4：曲线7~17)。表明，IMSA-LAMP检测方法筛选的引物能有效地检测茄科雷尔氏菌。
![]()
图 3 IMSA-LAMP引物组6的6条检测引物的位置及7个特异识别位点F1：CAGCGCGCAGTACAACTGCT，F2：GAATCACGATGCGGGTTCC， F3：CCATTTCCAGACGCCCTC，T：TCAAAGCGAGCAAGCTGTTCGGGAA，R1：CGACAACGCCGGCCAGAA，R2：GTCATCATGGCGGCTCGA，R3：GTTCGGTGCAAACGGCCFigure 3. Locations of six primers of primer group 6 and their seven specific recognition sites in IMSA-LAMP![]()
图 4 茄科雷尔氏菌IMSA-LAMP特异性检测结果1~6：茄科雷尔氏菌(YZqk1、YLqk1、XZqk2、XCqk4、LZqk5、YDqk6)；7：阴沟肠杆菌(XCYG-001)；8：克雷伯氏菌(LCKL-001)；9：菠萝泛菌(LCFJ-001)；10：铜绿假单胞菌(YD-001)；11：丁香假单胞菌(ZJDX-003)；12：多黏类芽孢杆菌(YD-002)；13：短小杆菌(YDDX-001)；14：芽孢杆菌(YD-003)；15：无菌水；16、17：健康桑树DNAFigure 4. Specific detection results of IMSA-LAMP for Ralstonia solanacearum1−6: Ralstonia solanacearum (YZqk1, YLqk1, XZqk2, XCqk4, LZqk5, YDqk6); 7: Enterobacter cloacae (XCYG-001); 8: Klebsiella sp. (LCKL-001); 9: Pantoea ananas (LCFJ-001); 10: Pseudomonas aeruginosa (YD-001); 11: Pseudomonas syringae (ZJDX-003); 12: Paenibacillus polymyxa (YD-002); 13: Bacillus pumilus (YDDX-001); 14: Bacillus sp. (YD-003); 15: Sterile water; 16−17: DNA of healthy mulberry2.3 IMSA-LAMP 引物浓度比例优化
设置外、茎、内引物浓度比例为1∶4∶8、1∶2∶4、1∶1∶2、1∶8∶4，以10 ng/μL的茄科雷尔氏菌YZqk1 DNA为模板进行IMSA-LAMP扩增，扩增温度为63 ℃，反应时间为45 min，结果如图5所示。4种不同引物浓度比例均出现了“S”曲线，当外、茎、内引物浓度比例为1∶4∶8时，起始扩增时间最短，为12 min，25 min到达扩增的平缓期；当外、茎、内引物浓度比例为1∶8∶4时，起始扩增时间为14 min，30 min达平缓期；当外、茎、内引物浓度比例为1∶2∶4时，起始扩增时间为23 min，35 min达到平缓期；当外、茎、内引物浓度比例为1∶1∶2时，起始扩增时间为25 min，45 min内未出现平缓期。最终可确定最佳外、茎、内引物比例为1∶4∶8。
![]()
图 5 不同外、茎、内引物浓度比例下茄科雷尔氏菌IMSA-LAMP检测结果Figure 5. IMSA-LAMP detection results ofRalstonia solanacearum under different concentration ratios of outer, stem and inner primers2.4 IMSA-LAMP反应温度优化
根据最佳引物比例进行IMSA-LAMP反应温度筛选，设置温度梯度为61.5、62.0、62.5、63.0、63.5、64.0、64.5及65.0 ℃，结果如图6所示。在61.5~65.0 ℃时，均有扩增曲线。在64.5 ℃时，起始扩增时间最短，为13 min，于24 min到达扩增的平缓期；在65.0 ℃时，在14 min曲线开始扩增，于25 min达到扩增平缓期；在61.5~64.0 ℃时，均在反应15 min后开始扩增，于30 min后达到扩增平缓期。故选择64.5 ℃作为最佳反应温度。
![]()
图 6 不同反应温度茄科雷尔氏菌IMSA-LAMP检测结果1：61.5 ℃，2：62.0 ℃，3：62.5 ℃，4：63.0 ℃，5：63.5 ℃，6：64.0 ℃，7：64.5 ℃，8：65.0 ℃Figure 6. Detection results of IMSA-LAMP for Ralstonia solanacearum at different reaction temperatures2.5 IMSA-LAMP法的灵敏度
茄科雷尔氏菌YZqk1纯培养后，取1 mL 1×106 CFU/mL的菌悬液，根据Ezup柱式细菌基因组DNA抽提试剂盒说明书提取总DNA，并利用超微量分光光度计进行浓度测定，其DNA质量浓度为2 ng/μL，10倍浓度梯度稀释7次，DNA质量浓度为2 ng/μL~200 ag/μL，对应菌为1×106~1×10−2 CFU/mL，将每个梯度作为模板进行IMSA-LAMP检测，并利用PCR进行平行检测作为对照。IMSA-LAMP检测结果如图7所示，茄科雷尔氏菌YZqk1 DNA质量浓度为2 ng/μL~200 fg/μL(对应菌为1×106~1×102 CFU/mL)，均出现了扩增曲线；而DNA质量浓度在20 fg/μL~200 ag/μL(对应菌为1×10~1×10−2 CFU/mL)未出现扩增曲线。因此，IMSA-LAMP对茄科雷尔氏菌最低DNA检测质量浓度为200 fg/μL，对应菌为1×102 CFU/mL。另外，PCR检测到的桑树茄科雷尔氏菌YZqk1最低DNA质量浓度为20 pg/μL，对应菌为1×104 CFU/mL(图8)。IMSA-LAMP检测桑树茄科雷尔氏菌具有良好的灵敏度，是PCR的100倍。
![]()
图 7 IMSA-LAMP对茄科雷尔氏菌灵敏度的检测结果1：2 ng/μL (1×106 CFU/mL)；2：200 pg/μL (1×105 CFU/mL)；3：20 pg/μL (1×104 CFU/mL)；4：2 pg/μL (1×103 CFU/mL)；5：200 fg/μL (1×102 CFU/mL)；6：20 fg/μL (1×10 CFU/mL)；7：2 fg/μL (1×10−1 CFU/mL)；8：200 ag/μL (1×10−2 CFU/mL)；9：无菌水；10：健康桑树DNAFigure 7. Sensitivity test results of IMSA-LAMP for Ralstonia solanacearum1: 2 ng/μL (1×106 CFU/mL); 2: 200 pg/μL (1×105 CFU/mL); 3: 20 pg/μL (1×104 CFU/mL); 4: 2 pg/μL (1×103 CFU/mL); 5: 200 fg/μL (1×102 CFU/mL); 6: 20 fg/μL (1×10 CFU/mL); 7: 2 fg/μL (1×10−1 CFU/mL); 8: 200 ag/μL (1×10−2 CFU/mL); 9: Sterilized water; 10: DNA of healthy mulberry![]()
图 8 PCR对茄科雷尔氏菌灵敏度检测结果1：2 ng/μL (1×106 CFU/mL)；2：200 pg/μL (1×105 CFU/mL)；3：20 pg/μL (1×104 CFU/mL)；4：2 pg/μL (1×103 CFU/mL)；5：200 fg/μL (1×102 CFU/mL)；6：20 fg/μL (1×10 CFU/mL)；7：2 fg/μL (1×10−1 CFU/mL)；8：200 ag/μL (1×10−2 CFU/mL)；9：无菌水；10：健康桑树DNAFigure 8. Sensitivity test results of PCR for Ralstonia solanacearum1: 2 ng/μL (1×106 CFU/mL); 2: 200 pg/μL (1×105 CFU/mL); 3: 20 pg/μL (1×104 CFU/mL); 4: 2 pg/μL (1×103 CFU/mL); 5: 200 fg/μL (1×102 CFU/mL); 6: 20 fg/μL (1×10 CFU/mL); 7: 2 fg/μL (1×10−1 CFU/mL); 8: 200 ag/μL (1×10−2 CFU/mL); 9: Sterilized water; 10: DNA of healthy mulberry2.6 IMSA-LAMP法对生产样本的检测
将收集的24份疑似桑树青枯病病样提取DNA后，利用IMSA-LAMP法与PCR法进行试验，评估2种检验方法的检测结果。IMSA-LAMP检测结果如图9所示，收集的疑似桑树青枯菌病样中，有21份样本出现了“S”扩增曲线(曲线2~22)，广西壮族自治区环江市及忻城市3份疑似桑树青枯病病样未出现扩增现象(曲线23~25)，检出率为87.5%。同步进行的PCR法，其中9份样品为阳性，15份样品为阴性(图10)，检出率为37.5%，9份PCR检测阳性的样本在IMSA-LAMP法检测结果中也为阳性，说明2种方法有一定的重复性，但IMSA-LAMP法的检出率明显高于PCR。
![]()
图 9 IMSA-LAMP对不同地区疑似桑树青枯病病样的检测结果1：茄科雷尔氏菌YZqk1；2~5：广西来宾；6~9：广西宜州；10~12：广东英德；13~14：广东阳山；15~16：广东罗定；17~18：广东广州；19~21：广西都安；22~23：广西环江；24~25：广西忻城；26：健康桑树DNAFigure 9. Detection results of suspected mulberry bacterial wilt samples from different districts by IMSA-LAMP1: Ralstonia solanacearum YZqk1; 2–5: Laibin of Guangxi; 6–9: Yizhou of Guangxi; 10–12: Yingde of Guangdong; 13–14: Yangshan of Guangdong; 15–16: Luoding of Guangdong; 17–18: Guangzhou of Guangdong; 19–21: Du'an of Guangxi; 22–23: Huanjiang of Guangxi; 24–25: Xincheng of Guangxi; 26: Healthy mulberry DNA![]()
图 10 PCR对不同地区疑似桑树青枯病病样的检测结果1：茄科雷尔氏菌YZqk1；2~5：广西来宾；6~9：广西宜州；10~12：广东英德；13~14：广东阳山；15~16：广东罗定；17~18：广东广州；19~21：广西都安；22~23：广西环江；24~25：广西忻城；26：健康桑树DNAFigure 10. Detection results of suspected mulberry bacterial wilt samples from different districts by PCR1: Ralstonia solanacearum YZqk1; 2–5: Laibin of Guangxi; 6–9: Yizhou of Guangxi; 10–12: Yingde of Guangdong; 13–14: Yangshan of Guangdong; 15–16: Luoding of Guangdong; 17–18: Guangzhou of Guangdong; 19–21: Du'an of Guangxi; 22–23: Huanjiang of Guangxi; 24–25: Xincheng of Guangxi; 26: Healthy mulberry DNA3. 讨论与结论
桑树青枯病作为一种毁灭性的细菌性病害，被国家列为重点检疫对象[27]。建立一种桑树茄科雷尔氏菌快速有效的检测方法，对预防桑树青枯病的发生和扩散具有重要的意义。桑树青枯病病原检测方法分为传统分离技术及分子检测技术。传统分离技术主要是通过病样组织分离病原菌的方式分离出桑树茄科雷尔氏菌，但操作繁琐，分离时间长，结果受人为影响大；随着分子检测技术的快速发展，核酸检测技术被广泛应用于医疗诊断、动植物安全生产、食品安全及动植物检疫等检测中，大大提高检测效率[28-30]。其中，LAMP是常见的核酸检测技术，本研究前期将LAMP技术应用于蚕、桑病害的检测研究，能有效检测出蚕桑病毒、细菌及真菌等病原物[31-33]。而IMSA技术的优点在于引物结合位点多，相较于LAMP技术拥有更高的灵敏度及特异性。Chen等[34]应用IMSA法对牛奶样品中金黄色葡萄球菌Staphylococcus aureus进行检测，灵敏度比LAMP方法高10倍；Liu等[35]对人畜共患病病原肠毒素大肠埃希菌(Enterotoxigenic Escherichia coli，ETEC)进行了IMSA的检测应用研究，相比于荧光定量PCR、LAMP、交叉引物等温扩增技术具有更高的灵敏度。同样地，本研究将LAMP反应体系与IMSA引物设计原理相结合应用于桑树青枯病的茄科雷尔氏菌检测，试验结果显示，IMSA-LAMP检测技术明显优于PCR的检测结果，本方法节约时间，且结果准确可靠。
果胶酶是一类分解果胶物质的酶的总称，在细菌生长和繁殖代谢过程中都能合成，而果胶裂解酶基因在青枯菌侵染植物过程中发挥了重要作用[36-37]。本研究在以NCBI登录的茄科雷尔氏菌基因组分析的基础上，以茄科雷尔氏菌果胶裂解酶基因(WP_197360060.1)序列为靶基因序列，进行了引物设计与筛选，经试验结果比较确定引物组6为最终的检测引物组。通过对茄科雷尔氏菌及桑源其他细菌的进行IMSA-LAMP检测，发现仅有茄科雷尔氏菌DNA样本出现扩增曲线，表明本研究建立的IMSA-LAMP检测技术特异性明显。对反应温度及引物比例进行优化，结果表明外、茎、内引物浓度比例为1∶4∶8，恒温64.5 ℃的荧光定量PCR仪在45 min内可完成检测，该方法比常规PCR方法快60 min以上，反应结果可通过荧光扩增曲线出现的先后进行快速判定，本研究建立的IMSA-LAMP法检测茄科雷尔氏菌DNA最低质量浓度达200 fg/μL，对应的菌为1×102 CFU/mL，比Okiro等[38]与李信申等[39]报道的青枯病LAMP检测方法的灵敏度高10倍左右，是常规PCR检测法灵敏度的100倍，同样的结果在Gou等[40]使用IMSA检测猪圆环病毒PCV3的应用中得到佐证。
在对采自田间的24份疑似桑树青枯病样本检测中，IMSA-LAMP检测出21份样本，3份样本未检测出桑茄科雷尔氏菌；同时利用PCR法对24份病样进行平行检测，PCR仅检测出9份，这9份样本也同时印证了IMSA-LAMP的检测结果。对于其他IMSA-LAMP检测出而PCR未能检测出的阳性样本，我们分析有2种原因，一是说明IMSA-LAMP检测方法更加灵敏，这个推测已经在本研究纯茄科雷尔氏菌模板检测中得到验证；二是说明导致桑树青枯和枯萎症状的病原菌种类很多，本研究课题组已经相继在桑树枯萎的病样中分离出阴沟肠杆菌[25, 41]、菠萝泛菌[42]及克雷伯氏菌[43]。而阴沟肠杆菌、菠萝泛菌及克雷伯氏菌作为机会性植物致病病原菌，在生产中导致植物枯萎的机理仍有待深入研究。综上，本研究建立的茄科雷尔氏菌的IMSA-LAMP检测方法准确可靠，可用于大批量桑树青枯病的快速分子诊断，为桑树细菌病检测和防疫提供了新技术支持。
致谢：试验过程中得到华南农业大学动物科学学院硕士研究生孟芳、王亚琴等的帮助，LAMP检测试剂得到广州双螺旋基因技术有限公司张璜博士与林梦蝶实验员的大力支持。
-
![]()
图 1 苹果离散元颗粒模型的建立过程
Figure 1. Establishment process of apple discrete element particle model
![]()
图 2 不同颗粒半径下的苹果离散元模型
Figure 2. Apple discrete element models with different particle radii
![]()
图 4 苹果跌落在泡棉的恢复系数的实测试验(a)和仿真试验(b)
1：吸盘；2：苹果；3：泡棉；4：直尺。
Figure 4. Actual measurement (a) and simulation test (b) of restitution coefficient of apple dropping onto foam
1: Suction cup; 2: Apple; 3: Foam; 4: Ruler.
![]()
图 5 不同试验条件下苹果与泡棉间的碰撞恢复系数
Figure 5. Restitution coefficients of collisions between apples and foams under different test conditions
![]()
图 6 仿真试验中苹果与不同密度泡棉间碰撞恢复系数与反弹高度的关系
Figure 6. Relationship between restitution coefficient and rebound height of collisions between apples and foams with different densities in simulation test
![]()
图 7 苹果间碰撞恢复系数的实测试验(a)和仿真试验(b)
1：棉线；2：苹果A；3：苹果B；4：模拟管道。
Figure 7. Actual measurement (a) and simulation test (b) of the restitution coefficient for collision between apples
1: Cotton thread; 2: Apple A 3: Apple B; 4: Simulated pipeline.
![]()
图 8 摆升高度与苹果间碰撞恢复系数的关系
Figure 8. Relationship between swing height and restitution coefficient of collision between apples
![]()
图 9 仿真试验中苹果间的碰撞恢复系数与苹果最大摆升高度的关系
Figure 9. Relationship between the restitution coefficient of collision between apples and the maximum swing height of apples in simulation test
![]()
图 10 仿真试验中苹果与泡棉间的静摩擦系数与其临界倾角的关系
Figure 10. Relationship between static friction coefficient and critical tilting angle between apples and foams in simulation test
![]()
图 11 苹果与泡棉的滚动摩擦系数测定试验
1：卷尺；2：苹果；3：泡棉。
Figure 11. Determination test of rolling friction coefficient between apple and foam
1: Tape measure; 2: Apple; 3: Foam.
![]()
图 12 仿真试验中苹果与泡棉间滚动摩擦系数与苹果水平滚动距离的关系
Figure 12. Relationship between rolling friction coefficient between apple and foam and horizontal rolling distance of apples in simulation test
![]()
图 13 苹果间静摩擦系数的实测试验(a)和仿真试验(b)
1：苹果；2：数显倾角仪；3：泡棉。
Figure 13. Actual measurement (a) and simulation test (b) of static friction coefficient between apples
1: Apple; 2: Digital inclinometer; 3: Foam.
![]()
图 14 仿真试验中苹果间静摩擦系数与临界倾角的关系
Figure 14. Relationship between static friction coefficient and critical tilt angle of apples in simulation test
![]()
图 15 实测试验(a)和仿真试验(b)下苹果的堆积情景
1：挡板；2：苹果；3：泡棉。
Figure 15. Apple stacking scenarios in actual measurement (a) and simulation test (b)
1: Baffle; 2:Apple; 3: Foam.
![]()
图 16 基于图像处理的苹果堆积角测定过程
Figure 16. Determination of apple repose angle based on image processing
![]()
图 17 苹果堆积角的测定结果及其与抽提速度的关系
Figure 17. The results of the apple stacking angle and its relationship with the lifting speed
![]()
图 18 仿真试验中苹果间滚动摩擦系数与堆积角的关系
Figure 18. Relationship between rolling friction coefficient and repose angle of apples in simulation test
![]()
图 19 圆筒抬升实际验证试验(a)与仿真试验(b)下苹果的堆积场景
1：圆筒，2：苹果，3：泡棉
Figure 19. Apple stacking scenarios in actual verification test (a) and simulation test (b) using cylinder lifting
1: Cylinder, 2:Apple, 3: Foam
![]()
图 20 圆筒抬升速度与苹果堆积角的关系
Figure 20. Relationship between cylinder lifting speed and repose angle of apples
表 1 不同质量的苹果在不同泡棉上的水平滚动距离及滚动摩擦系数
Table 1 Horizontal rolling distance and rolling friction coefficient of apples with different weights on different foams
果实质量/g
Fruit weight泡棉类型
Foam type水平滚动距离/mm
Horizonal rolling distance滚动摩擦系数
Rolling friction coefficient<170 高密度 533.32±30.77 0.0289 ±0.0010 特高密度 1071.80 ±35.600.0167 ±0.0040 179~190 高密度 524.54±37.24 0.0290 ±0.0010 特高密度 1071.80 ±40.800.0167 ±0.0005 >190 高密度 533.35±42.00 0.0289 ±0.0010 特高密度 1043.90 ±22.900.0170 ±0.0003 
下载: 导出CSV
表 2 不同质量苹果间的静摩擦系数
Table 2 Static friction coefficients between apples with different weights
苹果质量/g
Apple weight临界倾角
Critical tilt angle静摩擦系数
Static friction coefficient< 170 22.48±0.8 0.41±0.02 179~190 21.66±0.4 0.40±0.01 >190 22.5±0.3 0.41±0.01
下载: 导出CSV
-
[1] ZHANG Z, IGATHINATHANE C, LI J, et al. Technology progress in mechanical harvest of fresh market apples[J]. Computers and Electronics in Agriculture, 2020, 175: 105606 doi: 10.1016/j.compag.2020.105606
[2] 文恩杨, 李玉华, 牛子孺, 等. 蒜种颗粒离散元模型参数标定[J]. 农机化研究, 2021, 43(5): 160-167. doi: 10.3969/j.issn.1003-188X.2021.05.028 [3] SCHEFFLER O C, COETZEE C J, OPARA U L, et al. A discrete element model (DEM) for predicting apple damage during handling[J]. Biosystems Engineering, 2018, 172: 29-48. doi: 10.1016/j.biosystemseng.2018.05.015
[4] ZHAO H, HUANG Y, LIU Z, et al. Applications of discrete element method in the research of agricultural machinery: A review[J]. Agriculture, 2021, 11(5): 425. doi: 10.3390/agriculture11050425
[5] 徐泳, 李艳洁, 李红艳, 等. 离散元法在农业机械化中应用评述[J]. 农机化研究, 2004(5): 26-30. doi: 10.3969/j.issn.1003-188X.2004.05.008 [6] KAFASHAN J, WIACEK J, RAMON H, et al. Modelling and simulation of fruit drop tests by discrete element method[J]. Biosystems Engineering, 2021, 212: 228-240. doi: 10.1016/j.biosystemseng.2021.08.007
[7] 侯占峰, 戴念祖, 陈智, 等. 冰草种子物性参数测定与离散元仿真参数标定[J]. 农业工程学报, 2020, 36(24): 46-54. doi: 10.11975/j.issn.1002-6819.2020.24.006 [8] FANG W, WANG X, HAN D, et al. Review of material parameter calibration method[J]. Agriculture, 2022, 12(5): 706. doi: 10.3390/agriculture12050706
[9] 张国忠, 陈立明, 刘浩蓬, 等. 荸荠离散元仿真参数标定与试验[J]. 农业工程学报, 2022, 38(11): 41-50. doi: 10.11975/j.issn.1002-6819.2022.11.005 [10] 牛智有, 孔宪锐, 沈柏胜, 等. 颗粒饲料破损离散元仿真参数标定[J]. 农业机械学报, 2022, 53(7): 132-140. doi: 10.6041/j.issn.1000-1298.2022.07.013 [11] 王宪良, 胡红, 王庆杰, 等. 基于离散元的土壤模型参数标定方法[J]. 农业机械学报, 2017, 48(12): 78-85. doi: 10.6041/j.issn.1000-1298.2017.12.009 [12] 贾富国, 韩燕龙, 刘扬, 等. 稻谷颗粒物料堆积角模拟预测方法[J]. 农业工程学报, 2014, 30(11): 254-260. doi: 10.3969/j.issn.1002-6819.2014.11.031 [13] 王云霞, 梁志杰, 张东兴, 等. 基于离散元的玉米种子颗粒模型种间接触参数标定[J]. 农业工程学报, 2016, 32(22): 36-42. doi: 10.11975/j.issn.1002-6819.2016.22.005 [14] YAN D, YU J, WANG Y, et al. A review of the application of discrete element method in agricultural engineering: A case study of soybean[J]. Processes, 2022, 10(7): 1305. doi: 10.3390/pr10071305
[15] 任甲辉, 武涛, 刘庆庭, 等. 蔗段离散元仿真建模方法与参数标定[J]. 华南农业大学学报, 2022, 43(3): 124-132. doi: 10.7671/j.issn.1001-411X.202108015 [16] LI C, LI Z, WANG T, et al. Parameter optimization of column-comb harvesting of litchi based on the EDEM[J]. Scientia Horticulturae, 2023, 321: 112216. doi: 10.1016/j.scienta.2023.112216
[17] FAN G, WANG S, SHI W, et al. Simulation parameter calibration and test of typical pear varieties based on discrete element method[J]. Agronomy, 2022, 12(7): 1720. doi: 10.3390/agronomy12071720
[18] 马云海. 农业物料学[M]. 北京: 化学工业出版社, 2015. [19] FU H, YANG J, DU W, et al. Determination of coefficient of restitution of fresh market apples caused by fruit-to-fruit collisions with a sliding method[J]. Biosystems Engineering, 2022, 224: 183-196. doi: 10.1016/j.biosystemseng.2022.10.010
[20] 哈尔滨工业大学理论力学教研室. 理论力学 II[M]. 北京: 高等教育出版社, 2016. [21] FU H, HE L, MA S, et al. ‘Jazz’ apple impact bruise responses to different cushioning materials[J]. Transactions of the ASABE, 2017, 60(2): 327-336. doi: 10.13031/trans.11946
[22] PANG W, STUDMAN C J, WARD G T. Bruising damage in apple-to-apple impact[J]. Journal of Agricultural Engineering Research, 1992, 52: 229-240. doi: 10.1016/0021-8634(92)80063-X
[23] 刘文政, 何进, 李洪文, 等. 基于离散元的微型马铃薯仿真参数标定[J]. 农业机械学报, 2018, 49(5): 125-135. doi: 10.6041/j.issn.1000-1298.2018.05.014 [24] 吴孟宸, 丛锦玲, 闫琴, 等. 花生种子颗粒离散元仿真参数标定与试验[J]. 农业工程学报, 2020, 36(23): 30-38. doi: 10.11975/j.issn.1002-6819.2020.23.004 [25] 李勤良. 颗粒堆积性质和散状物料转载过程的DEM仿真研究[D]. 沈阳: 东北大学, 2010. [26] SCHEFFLER O C, COETZEE C J, OPARA U L. A discrete element model (DEM) for predicting apple damage during handling[J]. Biosystems Engineering, 2018, 172: 29-48. doi: 10.1016/j.biosystemseng.2018.05.015
-
期刊类型引用(1)
1. 陈文静，董章勇，宋汉达，罗梅. 环介导等温扩增技术在植物病原物检测中的应用. 仲恺农业工程学院学报. 2024(04): 48-54 . 
百度学术
其他类型引用(2)
计量
- 文章访问数: 45
- HTML全文浏览量: 7
- PDF下载量: 42
- 被引次数: 3
