Discrete element simulation modeling method and parameter calibration of sugarcane segment
-
摘要:目的
蔗段作为大长径比秆状物料,其离散元模型的构建方法与仿真参数的设定尚不清楚,仿真模型精度对颗粒间的动力学响应特性有较大影响,需通过参数标定提高仿真参数的准确性。
方法以蔗段物理堆积角为响应值,采用仿真试验方法优化标定离散元参数。首先,采用物理试验测定蔗段的基本物性参数,并基于多球聚合模型和XML的方法构建蔗段仿真模型;然后,应用Plackett-Burman试验对蔗段离散元仿真中的8个初始参数进行显著性筛选,并对显著性参数进行最陡爬坡试验,确定最优参数区间;最后,基于Box-Behnken试验建立显著性参数与堆积角的二阶回归方程,以物理试验堆积角42.70°为目标值,对回归方程进行优化求解。
结果显著性筛选试验得出蔗段泊松比、蔗段−蔗段静摩擦系数、蔗段−蔗段滚动摩擦系数对仿真堆积角影响显著;最优参数组合为:蔗段泊松比0.35、蔗段−蔗段静摩擦系数0.53、蔗段−蔗段滚动摩擦系数0.04。最优参数组合的仿真试验结果表明,仿真堆积角与物理试验堆积角无显著性差异,两者相对误差为0.99%,进一步验证了蔗段离散元标定参数的可靠性。
结论蔗段离散元模型与最优仿真参数可用于蔗段离散元仿真试验,并可为大长径比秆状农业物料的离散元参数标定提供参考。
Abstract:ObjectiveThe construction method of discrete element model and the setting of simulation parameters of sugarcane segment as stalk material with large length and diameter ratio are not clear, the accuracy of the model has a great influence on the dynamic response characteristics between particles, and the accuracies of the simulation parameters need to be improved by parameter calibration.
MethodTaking the physical repose angle of sugarcane segment as the response value, the discrete element parameters were optimized and calibrated by simulation test method. Firstly, the basic physical parameters of sugarcane segment were measured by physical test, and the sugarcane segment simulation model was constructed based on the multi-sphere polymerization model and XML method. Then, the Plackett-Burman test was used to screen the significance of eight initial parameters in the discrete element simulation of sugarcane segment, and the optimal value ranges of significant parameters were determined by the steepest ascent search test. Further, the second-order regression equation between the significant parameters and repose angle was established based on the Box-Behnken test, and taking the physical repose angle of 42.70° as the target value, the regression equation was optimized.
ResultThe significance screening test showed that the Poisson’s ratio of sugarcane segment, the static and dynamic friction coefficients between two sugarcane segments had significant influence on the simulated repose angle. The optimal parameter combination was as follows: The Poisson’s ratio of sugarcane segment was 0.35, the static and dynamic friction coefficients between sugarcane segments were 0.53 and 0.04 respectively. The simulation test results of the optimal parameter combination showed that there was no significant difference between the simulated repose angle and the physical repose angle, and the relative error between them was 0.99%, which further verifies the reliability of sugarcane segment parameters calibrated by discrete element method.
ConclusionThe sugarcane segment discrete element model and the optimal simulation parameters can be used in the sugarcane segment discrete element simulation test, and can provide a reference for calibration of the discrete element parameter of stalk agricultural materials with large length and diameter ratio.
-
Keywords:
- Discrete element /
- Calibration /
- Optimization /
- Simulation parameter /
- Repose angle /
- Sugarcane segment
-
蓝舌病病毒(Bluetongue virus,BTV)、流行性出血病病毒(Epizootic hemorrhagic disease virus,EHDV)和帕利亚姆血清群病毒(Palyam serogroup virus,PALV)均为通过库蠓Culicoides传播的呼肠孤病毒科Reoviridae环状病毒属Orbivirus成员,广泛分布于热带、亚热带与温带地区,牛是3种病毒的重要宿主[1-3]。BTV、EHDV和PALV在全球引起的动物疫病的暴发和流行给牛羊养殖业造成了严重经济损失,导致牛羊进出口贸易受限,影响了国际畜产品正常贸易。据估计,2009年全球畜牧业仅因BTV引发蓝舌病而导致的经济损失就高达30亿美元[4];近年来EHDV在中东多国特别是在阿尔及利亚、突尼斯、摩洛哥和以色列等国家的牛群中多次引发流行性出血病的暴发,严重影响了当地畜牧业的健康发展[5];PALV感染妊娠期母牛能引起新生犊牛积水性无脑和小脑发育不全,曾在日本和韩国暴发,给两国的肉牛养殖业造成了严重的经济损失[6]。目前,蓝舌病和流行性出血病已经被世界动物卫生组织(Office international des epizooties,OIE)列为法定报告的跨境动物疫病。
BTV、EHDV和PALV血清型众多,目前世界范围内已发现27种BTV血清型[7]、9种EHDV血清型[8]和10种PALV血清型[9]。2013年在公益性行业(农业)科研专项的资助下,在全国范围内开展了虫媒病毒的流行病学调查,发现我国主要流行12种BTV血清型(BTV-1、-2、-3、-4、-5、-7、-9、-12、-15、-16、-21和-24)[10],5种EHDV血清型(EHDV-1、-5、-6、-7和-10)[11]和3种PALV血清型[Chuzan disease virus (CHUV)、D’Aguilar virus (DAV)与Bunyip Creek virus (BCV)][12]等,表明流行于我国的BTV、EHDV和PALV血清型具有丰富多样性。对3种病毒的血清学调查发现,我国至少12个省区的牛羊群中同时存在BTV、EHDV和PALV抗体阳性,且阳性率由北向南逐渐升高,特别是在南方地区,如云南、广西和广东等地的牛群中,3种病毒的血清抗体阳性率均高于30%;说明BTV、EHDV和PALV在我国已经广泛存在,严重威胁我国畜牧产业的健康发展[13-15]。
快速诊断病原是动物疫病防控的重要一环,BTV、EHDV和PALV不仅具有共同的传播媒介与易感动物,在流行区域上也存在重叠;因此,建立一种能同时诊断这3种病原的检测方法具有重要的应用价值。多重RT-qPCR可在1个反应管中同时检测多种目的基因,克服了普通荧光PCR只能检测单一基因的缺点,具有特异性强、敏感性高和短时间内检测多种病原的优点,被广泛应用于动物传染病的诊断[16]。本研究旨在运用TaqMan探针技术建立一种三重RT-qPCR方法,提高病毒检出率,简化操作流程,降低检测成本,为临床快速鉴别检测BTV、EHDV和PALV提供有效的技术手段。
1. 材料与方法
1.1 毒株、血液和核酸样品
24个BTV血清型(BTV-1~BTV-24)南非参考毒株由OIE南非蓝舌病参考实验室(Onderste-poort Veterinary Institute)提供,6个EHDV血清型(EHDV-1、-2、-5、-6、-7、-8)参考毒株由澳大利亚麦克阿瑟·伊丽莎白农业研究所(Elizabeth Macarthur Agricultural Institute,EMAI)提供,中国分离的12种BTV血清型、5种EHDV血清型和3种PALV血清型代表毒株信息见表1。阿卡斑病毒(Akabane virus, AKAV)核酸由云南省畜牧兽医科学院热带亚热带动物病毒病重点实验室保存,小反刍兽疫病毒(Peste des petits ruminants virus, PPRV)核酸由中国农业科学院哈尔滨兽医研究所惠赠,口蹄疫病毒(Foot-and-mouth disease virus, FMDV)、牛流行热病毒(Bovine ephemeral fever virus, BEFV)核酸由中国农业科学院兰州兽医研究所惠赠。BTV Seg-10、EHDV Seg-5和PALV Seg-7基因节段重组质粒pLB-BTV-NS3-T7、pLB-EHDV-NS1-T7和pLB-CHUV-VP7-T7由云南省畜牧兽医科学院热带亚热带动物病毒病重点实验室前期制备并保存。
1.2 主要试剂与仪器
病毒RNA抽提试剂盒MagMAX-96 Viral RNA Isolation Kit购自Thermo公司,DNA胶回收试剂盒、荧光定量PCR试剂盒购自宝生物工程(大连)有限公司,T7启动子RNA体外转录试剂盒购自Promega公司,RNA纯化试剂盒、质粒DNA提取试剂盒、pLB平末端克隆载体购自天根生化科技有限公司,7500实时荧光定量PCR系统购自ABI公司。
1.3 引物和探针的设计与合成
从NCBI检索并下载BTV、EHDV、PALV病毒的全基因组序列,并根据云南省畜牧兽医科学院热带亚热带动物病毒病重点实验室分离的3种病毒的全基因组序列,利用MEGA-X进行序列对比分析,确定并选择BTV NS3、EHDV NS1和PALV VP7基因保守区序列,利用Primer Express 6.0软件设计相应的特异性引物,分别用FAM、VIC和TAMRA荧光基团标记TaqMan探针,使用BHQ1作为淬灭基团(表2)。引物和探针均由宝生物工程(大连)有限公司合成,并用RNase-free 水稀释至20 μmol/L。
表 2 蓝舌病病毒、流行性出血病病毒和帕利亚姆血清群病毒三重RT-qPCR的引物与探针Table 2. Primers and probes used for detecting bluetongue virus, epizootic hemorrhagic disease virus and Palyam serogroup virus by triplex RT-qPCR引物/探针
Primer/Probe扩增或标记位置
Amplification or
label location扩增或标记长度/bp
Amplification or
label length引物或探针序列(5′→3′)
Primer or probe sequenceBTV-NS3-F Seg-10 189~207 TGGAYAAAGCGATGTCAAA BTV-NS3-R Seg-10 266~285 ACATCATCACGAAACGCTTC BTV-NS3-Probe Seg-10 243~264 FAM-ARGCTGCATTCGCATCGTACGC-BHQ1 EHDV-NS1-F Seg-5 12~32 TCTTCGTCGACTGCCATCGAG EHDV-NS1-R Seg-5 66~94 AACATTTTGACATGATTKGCRTARTAACT EHDV-NS1-Probe Seg-5 34~63 TAMRA-TGGAGCGCTTTTTGAGAAAATACAACATGA-BHQ1 PALV-VP7-F Seg-7 275~295 CATCAATGGCAACAATCGGTG PALV-VP7-R Seg-7 346~369 ATTCAGCATACCTGTAATTCGTAC PALV-VP7-Probe Seg-7 314~341 VIC-TTCCATATACAACGTCGGCAATGACAAG-BHQ1 1.4 标准品的制备
提取pLB-BTV-NS3-T7、pLB-EHDV-NS1-T7和pLB-CHUV-VP7-T7基因序列重组质粒,使用XbaⅠ内切酶进行质粒的酶切线性化。取1 μg纯化后的DNA作为模板,使用T7体外转录试剂盒进行BTV Seg-10 (822 bp)、EHDV Seg-5 (533 bp)和PALV Seg-7 (950 bp) ssRNA的体外转录,反应体系为:DNA模板5 μL、RiboMAXTM Express T7 2× Buffer 10 μL、T7 Express Enzyme Mix 2 μL,加入RNase-free 水补齐至20 μL;37 ℃反应30 min。使用RNA纯化试剂盒纯化转录产物,Nanodrop ND-2000测定核酸浓度后,根据相对分子质量计算纯化后ssRNA的拷贝数。
1.5 三重RT-qPCR反应条件的优化
以体外转录的BTV Seg-10、EHDV Seg-5和PALV Seg-7的ssRNA为模板,采用方阵法分别对三重RT-qPCR方法的引物浓度(0.1~0.5 μmol·L−1)、探针浓度(0.1~0.8 μmol·L−1)、退火温度(50~60 ℃)和模板体积(2~8 μL)进行筛选,确定最佳反应条件。
1.6 特异性试验
以制备好的BTV Seg-10、EHDV Seg-5和PALV Seg-7的ssRNA模板作为阳性对照,与AKAV、PPRV、FMDV和BEFV核酸一同进行三重RT-qPCR检测,设立阴性牛羊血液总RNA及DEPC水分别作为阴性和空白对照,重复3次。
1.7 敏感性试验和标准曲线的绘制
分别将3种病毒的体外转录ssRNA以10倍梯度稀释获得拷贝数为1×1010~1×10−1 µL−1的12个梯度,分别以阴性牛羊血液总RNA及DEPC水作为阴性和空白对照。所有样品及对照均设2个复孔,同时进行单重和三重RT-qPCR检测,并以RNA浓度的对数值为横坐标、三重RT-qPCR结果Ct值为纵坐标绘制标准曲线。
1.8 重复性试验
以BTV、EHDV和PALV体外转录的ssRNA的3个稀释度(拷贝数分别为1×103、1×105和1×107 μL−1)为模板,每个稀释度各取10份样本分别在不同的时间重复检测3次,同时设立阴性牛羊血液总RNA及DEPC水分别作为阴性和空白对照,并通过计算Ct值的标准差和变异系数,验证三重RT-qPCR的批内和批间重复性。
1.9 混合病毒核酸样品的检测试验
将BTV、EHDV和PALV 3种病毒高、中、低浓度(拷贝数分别为1×106、1×104和1×102 µL−1)的ssRNA分别按体积比1∶1混合,同时进行单重和三重RT-qPCR检测,并设立阴性牛羊血液总RNA及DEPC水分别作为阴性和空白对照。在3种病毒核酸混合的情况下分析高浓度的核酸模板是否会影响三重RT-qPCR对低浓度核酸模板的检测结果。
1.10 毒株与临床血液样本的检测
提取2012—2016年在我国分离的BTV、EHDV和PALV不同血清型毒株及已分离毒株对应的阳性血液样本核酸,与24个BTV血清型及6个EHDV血清型参考毒株的核酸同时进行单重和三重RT-qPCR检测。比较试验结果,分析和评价该方法的临床实用性。
表 1 中国分离的不同蓝舌病病毒、流行性出血病病毒和帕利亚姆血清群病毒血清型代表毒株信息Table 1. Representative strains of different serotypes of bluetongue virus, epizootic hemorrhagic disease virus and Palyam serogroup virus isolated in China血清型 Serotype 毒株 Strain 分离年份 Isolation year 分离地点 Isolation place 来源动物 Animal source BTV-1 V147 2012 云南师宗 Shizong, Yunnan 绵羊 Sheep BTV-2 V146 2013 云南普洱 Puer, Yunnan 黄牛 Cattle BTV-3 V025 2013 云南德宏 Dehong, Yunnan 黄牛 Cattle BTV-4 V030 2013 云南普洱 Puer, Yunnan 黄牛 Cattle BTV-5 V084 2012 云南师宗 Shizong, Yunnan 黄牛 Cattle BTV-7 V089 2014 广东汕头 Shantou, Guangdong 黄牛 Cattle BTV-9 V008 2013 云南德宏 Dehong, Yunnan 黄牛 Cattle BTV-12 V007 2012 广西南宁 Nanning, Guangxi 黄牛 Cattle BTV-15 V061 2014 云南普洱 Puer, Yunnan 黄牛 Cattle BTV-16 V109 2012 云南师宗 Shizong, Yunnan 绵羊 Sheep BTV-21 V057 2013 江苏须臾 Xuyu, Jiangsu 黄牛 Cattle BTV-24 V015 2013 云南德宏 Dehong, Yunnan 黄牛 Cattle EHDV-1 V128 2013 广东汕头 Shantou, Guangdong 奶牛 Cow EHDV-5 V127 2015 广西融水 Rongshui, Guangxi 水牛 Buffalo EHDV-6 V274 2014 云南普洱 Puer, Yunnan 黄牛 Cattle EHDV-7 V296 2013 云南师宗 Shizong, Yunnan 水牛 Buffalo EHDV-10 V280 2017 云南师宗 Shizong, Yunnan 黄牛 Cattle CHUV V144 2012 云南师宗 Shizong, Yunnan 黄牛 Cattle DAV V157 2015 云南红河 Honghe, Yunnan 水牛 Buffalo BCV V252 2016 广东汕头 Shantou, Guangdong 奶牛 Cow 2. 结果与分析
2.1 标准品的制备
以构建好的3种重组质粒为模板进行体外转录,获得BTV Seg-10 (822 bp)、EHDV Seg-5 (533 bp)和PALV Seg-7 (950 bp)的ssRNA,经纯化后测定质量浓度分别为656.5、408.6和785.7 ng·μL−1,计算获得拷贝数分别为1.41×1012、1.36×1012和1.46×1012 μL−1。将获得的ssRNA样品作为三重RT-qPCR灵敏性评价的标准品。
2.2 三重RT-qPCR方法的建立与优化
通过引物和探针浓度的不同配比试验及反应条件的优化试验,最终确定反应体系。总反应体积为20 μL,其中2× One Step RT-PCR Buffer Ⅲ 10 μL,Ex Taq HS (5 U·μL−1) 0.4 μL,PrimeScript RT Enzyme Mix Ⅱ 0.4 μL,Rox Reference Dye Ⅱ (50×) 0.4 μL,BTV上、下游引物各0.2 μL及探针0.4 μL,EHDV的上、下游引物各0.3 μL及探针0.6 μL,PALV上、下游引物各0.2 μL及探针0.4 μL,变性模板RNA 4 μL,双蒸水补充至20 μL。反应条件:42 ℃ 5 min;92 ℃ 10 s;92 ℃ 5 s,60 ℃ 34 s,40个循环。
2.3 特异性试验
将3种病毒体外转录ssRNA作为阳性对照,同时对AKAV、PPRV、FMDV和BEFV的核酸及阴性对照进行三重RT-qPCR检测。结果(图1)显示,本研究建立的三重RT-qPCR具有高度特异性,仅与对应的BTV、EHDV和PALV核酸产生扩增曲线,其他病原(AKAV、PPRV、FMDV、BEFV)和阴性对照均无扩增曲线,判定为阴性。
2.4 敏感性试验和标准曲线的绘制
以拷贝数为1×1010~1×10−1 µL−1的12个梯度ssRNA作为模板,根据已优化的条件同时进行三重和单重RT-qPCR反应。结果显示,所有阴性对照及空白对照均无阳性扩增(图2a、2b、2c);在核酸拷贝数为1×1010~1×101 µL−1的范围内,标准曲线均具有良好的线性,相关系数(R)均在0.998以上,引物扩增效率(E)为91.99%~103.42%(图2d、2e、2f)。各浓度体外转录ssRNA单重与三重RT-qPCR检测的Ct均值(表3)表明本研究建立的三重 RT-qPCR最低可检测到拷贝数为1×101 µL−1的BTV、EHDV和PALV核酸,与单重检测结果一致。
图 2 蓝舌病病毒、流行性出血病病毒和帕利亚姆血清群病毒体外转录ssRNA对三重RT-qPCR灵敏性评价曲线上数字0~10分别代表ssRNA拷贝数为1×100~1×1010 μL−1Figure 2. Sensitivity evaluation of bluetongue virus, epizootic hemorrhagic disease virus and Palyam serogroup virus in vitro transcribed ssRNA to triplex RT-qPCRThe numbers from 0 to 10 on the curves indicate ssRNA copy numbers of 1×100−1×1010 μL−1 respectively表 3 单重和三重RT-qPCR方法Ct值比较1)Table 3. Ct value comparison between monoplex and triplex RT-qPCR核酸拷贝数/μL−1
Nucleic acid copy number蓝舌病病毒
Bluetongue virus流行性出血病病毒
Epizootic hemorrhagic disease virus帕利亚姆血清群病毒
Palyam serogroup virus单重 Monoplex 三重 Triplex 单重 Monoplex 三重 Triplex 单重 Monoplex 三重 Triplex 1×1010 8.57 9.00 9.05 9.15 6.47 6.80 1×109 11.63 12.17 12.23 12.26 9.22 9.20 1×108 15.43 15.72 15.71 15.76 13.75 13.45 1×107 18.83 19.11 19.35 19.04 16.83 17.05 1×106 22.30 22.34 22.59 22.48 20.92 21.01 1×105 25.40 25.75 26.15 26.12 24.59 24.55 1×104 28.50 28.87 29.47 29.59 28.09 28.13 1×103 31.75 32.10 32.65 33.17 31.03 30.96 1×102 34.38 35.04 35.51 36.03 34.28 34.61 1×101 37.78 38.11 38.28 38.64 37.48 37.66 1×100 — — — — — — 1×10−1 — — — — — — 1) “—”表示未检测出
1) “—”indicates not detectable2.5 混合病毒核酸样品的检测
将BTV、EHDV和PALV 3种病毒的ssRNA分别按高、中、低3个浓度等体积混合,同时进行单重和三重RT-qPCR检测。结果(表4)显示,一种病毒核酸模板浓度较高(拷贝数为106 μL−1)而另一种病毒核酸模板浓度较低(拷贝数为102 μL−1)时,三重RT-qPCR依然可以同时检测到3种病毒,并且与单重RT-qPCR结果无明显差异。
表 4 混合病毒核酸样品单重和三重RT-qPCR检测结果Table 4. Results of mixed virus nucleic acid samples detected by monoplex and triplex RT-qPCR样品
Sample病毒1)
Virus核酸拷贝数/ μL−1
Nucleic acid copy numberCt 单重 Monoplex 三重 Triplex 1 BTV 1×106 23.20 23.24 EHDV 1×104 30.45 30.57 PALV 1×102 35.33 35.66 2 BTV 1×104 29.48 29.85 EHDV 1×102 37.01 37.08 PALV 1×106 21.82 21.91 3 BTV 1×102 35.43 36.09 EHDV 1×106 23.49 23.38 PALV 1×104 29.07 29.11 1) BTV:蓝舌病病毒;EHDV:流行性出血病病毒;PALV:帕利亚姆血清群病毒
1) BTV: Bluetongue virus; EHDV: Epizootic hemorrhagic disease virus; PALV: Palyam serogroup virus2.6 重复性试验
以高、中、低3个稀释度的BTV、EHDV和PALV的ssRNA作为模板,分3个批次进行检测,每批次每个稀释度各10份样本,通过计算批内、批间变异系数评价多重体系的稳定性。结果(表5)显示,BTV、EHDV和PALV病毒核酸高中低3个浓度所对应的批内和批间Ct值的标准差均在0.5以内,变异系数均在2%以下,可见三重体系对于3种病毒核酸的检测均具有良好的稳定性。
表 5 三重RT-qPCR重复性试验中Ct值与变异系数Table 5. Means and coefficients of variation in repeatability tests of triplex RT-qPCR病毒
Virus核酸拷贝数/μL−1
Nucleic acid copy numberCt1) 变异系数/% Coefficient of variation 批内 Intra-assay 批间 Inter-assay 批内 Intra-assay 批间 Inter-assay BTV 1×107 18.61±0.19 18.89±0.29 1.05 1.54 1×105 25.39±0.26 25.32±0.40 1.03 1.58 1×103 31.37±0.43 32.54±0.46 1.37 1.40 EHDV 1×107 19.32±0.25 19.12±0.30 1.30 1.59 1×105 26.28±0.32 26.66±0.48 1.22 1.82 1×103 32.62±0.35 33.17±0.40 1.08 1.19 PALV 1×107 16.85±0.20 17.15±0.24 1.16 1.38 1×105 24.58±0.27 24.43±0.40 1.10 1.64 1×103 30.91±0.32 30.55±0.50 1.03 1.44 1)表中数据为平均值±标准差
1)The data are mean ± standard deviation2.7 毒株与临床血液样本的检测
对2012—2016年在我国分离的BTV、EHDV和PALV不同血清型毒株及分离出病毒的阳性血液样本和BTV、EHDV参考毒株核酸同时进行单重和三重RT-qPCR检测。结果(表6和表7)显示,本研究建立的三重RT-qPCR方法对参考毒株和在我国分离培养的毒株核酸检测Ct值为14.67~27.62,对血液样本检测Ct值为27.69~35.24;单重RT-qPCR检测病毒毒株Ct值为15.08~26.04,对血液样品检测Ct值为28.03~34.67。对2种方法的检测结果进行单因素方差分析,P>0.05,无显著差异。本研究建立的三重RT-qPCR方法可准确检测出国内外BTV、EHDV和PALV不同血清型毒株及临床血液样本中对应的不同血清型病毒核酸,具有较好的群特异性,与单重检测结果无明显差异。
表 6 单重和三重RT-qPCR检测中国分离的蓝舌病病毒、流行性出血病病毒和帕利亚姆血清群病毒毒株及对应阳性血液样本的验证试验Table 6. Demonstration tests of bluetongue virus, epizootic hemorrhagic disease and Palyam serogroup virus strains isolated in China and corresponding blood samples by monoplex and triplex RT-qPCR血清型
Serotype毒株
StrainCt(病毒样品 Virus sample) Ct(血液样品 Blood sample) 单重 Monoplex 三重 Triplex 单重 Monoplex 三重 Triplex BTV-1 V147 16.21 16.79 32.01 32.63 BTV-2 V146 22.12 22.60 30.89 30.67 BTV-3 V025 19.03 18.87 30.02 29.82 BTV-4 V030 23.13 23.45 31.52 31.76 BTV-5 V084 20.57 21.67 31.87 32.21 BTV-7 V089 25.23 25.12 31.13 30.53 BTV-9 V008 17.55 17.63 30.21 29.98 BTV-12 V007 21.06 20.15 28.03 27.69 BTV-15 V061 18.73 18.67 32.89 33.13 BTV-16 V109 20.12 19.32 33.34 34.27 BTV-21 V057 24.32 25.07 32.08 31.45 BTV-24 V015 17.55 17.34 34.67 35.24 EHDV-1 V128 18.67 18.69 30.34 29.68 EHDV-5 V127 20.23 20.12 29.76 30.20 EHDV-6 V274 22.95 23.32 34.32 34.46 EHDV-7 V296 15.15 14.67 33.58 32.14 EHDV-10 V280 16.07 16.12 30.31 29.07 CHUV V144 21.22 20.69 31.68 30.57 DAV V157 16.32 15.32 30.57 31.75 BCV V252 23.46 24.54 31.76 32.62 表 7 单重和三重RT-qPCR检测蓝舌病病毒和流行性出血病病毒参考毒株的验证试验Table 7. Demonstration test of reference strains of bluetongue virus and epizootic hemorrhagic disease virus by monoplex and triplex RT-qPCR血清型
SerotypeCt 血清型
SerotypeCt 单重 Monoplex 三重 Triplex 单重 Monoplex 三重 Triplex BTV-1 17.09 16.72 BTV-16 19.66 20.07 BTV-2 20.62 21.45 BTV-17 15.08 14.86 BTV-3 21.09 22.63 BTV-18 18.60 17.32 BTV-4 18.32 17.82 BTV-19 22.23 21.12 BTV-5 20.44 19.31 BTV-20 21.15 20.67 BTV-6 21.16 20.08 BTV-21 26.04 27.62 BTV-7 22.67 23.43 BTV-22 22.62 23.51 BTV-8 15.96 15.67 BTV-23 25.34 26.02 BTV-9 22.26 21.30 BTV-24 20.26 19.09 BTV-10 20.68 21.67 EHDV-1 19.68 18.93 BTV-11 21.75 22.52 EHDV-2 21.79 22.22 BTV-12 21.63 20.03 EHDV-5 18.60 17.45 BTV-13 22.78 23.12 EHDV-6 19.73 18.96 BTV-14 24.06 23.56 EHDV-7 20.62 19.65 BTV-15 23.12 22.67 EHDV-8 19.63 20.05 3. 讨论与结论
作为虫媒病毒的BTV、EHDV和PALV血清型多样,传播途径相同,流行地域重叠,并且在动物感染早期存在部分相似的临床症状。国外研究曾报道3种病毒交叉感染的情况[17]。2013—2016年在公益性行业(农业)科研专项的资助下,云南省畜牧兽医科学院热带亚热带动物病毒病重点实验室在我国多个虫媒病毒监控点的不同哨兵牛上同时分离出了BTV、EHDV和PALV,对这3种病原的早期诊断是一项重大挑战,往往仅能检测出其中1种病毒核酸阳性忽略另外2种病毒,增加病毒检测的漏检率,不利于动物疫病的防控。基于此,本试验在原有RT-qPCR对虫媒病毒的检测基础上,建立了针对BTV、EHDV和PALV的三重RT-qPCR检测方法。
多重RT-qPCR方法易受引物、试剂和反应条件等诸多因素影响[18]。因此本研究首先根据BTV NS3、EHDV NS1和PALV VP7基因的高度保守性和较好群特异性等特点[19-21],选择作为引物和探针设计的靶基因,同时适当引入简并碱基,并进行BLAST序列对比分析确保引物和探针的保守性和特异性。随后对反应体系的退火温度,模板体积、引物和探针浓度进行了一系列优化,确定了最优反应条件,将复合扩增中引物间、模板间和引物与模板间可能存的竞争抑制和相互干扰降到最低,从而实现从单种病毒检测到同时进行3种病毒的检测,节省样品用量、试验资源和操作时间。
本研究利用BTV、EHDV和PALV的体外转录ssRNA对三重RT-qPCR灵敏度进行测试并绘制标准曲线,结果显示,三重反应对病毒体外转录ssRNA拷贝数最低检测限在10 μL−1级别,与单重反应相比,灵敏度无明显差异。龚雪蕊等[22]建立的寨卡病毒、登革病毒和基孔肯雅病毒三重RT-PCR检测方法的最低拷贝数检出限为10 μL−1,Nunes等[23]报道的用于南美洲汉坦4种流行病毒的RT-qPCR检测方法的拷贝数灵敏度为10 μL−1,二者均与本研究相近。本研究构建的标准曲线的相关系数均大于等于0.998,扩增效率基本相同且均大于90%,符合RT-qPCR扩增效率在90%~110%的合理范围。
前期调查发现,BTV、EHDV和PALV在临床上出现混合感染时,3种病毒在动物血液中的载量会存在较大差异;因此,本试验将3种病毒不同浓度的核酸模板进行混合,分析高浓度核酸模板是否对低浓度核酸模板存在干扰。结果显示,一种病毒核酸模板浓度较高而另一种病毒模板浓度较低时,所建立的三重RT-qPCR方法依然可以同时检测到3种病毒,并且与单重法的结果无明显差异,表明高浓度模板对低浓度模板的扩增检测干扰不明显。重复性试验中,对3种病毒高、中、低3个浓度的模板ssRNA进行检测,变异系数均小于2%,表明本研究建立的三重RT-qPCR检测方法具有较好的稳定性和重复性,有良好的检测效果。
2017年刘佳佳等[24]建立了BTV和EHDV双重RT-qPCR检测方法,但仅进行了1种BTV和EHDV血清型的特异性验证。本研究以单重RT-qPCR检测方法作为参照,使用建立的三重RT-qPCR同时对不同BTV和EHDV血清型参考毒株和我国主要流行的不同BTV、EHDV和PALV血清型毒株及其阳性血液样品进行检测。结果显示,本研究建立的三重RT-qPCR方法具有较好的群特异性,可以检测出不同血清型的参考毒株和在我国主要流行的毒株及其血液样品中对应病毒的核酸,与单重法的检检测结果无明显差异。
综上所述,本研究将在我国已经存在、主要流行于南方热带亚热带地区的3种虫媒病毒作为一个组合,首次建立了一套三重RT-qPCR检测方法,填补了同时检测BTV、EHDV和PALV 3种病毒核酸方法的空缺。该方法可应用于疫区或出入境动物带毒状况监测,对早期临床疫病鉴定和防疫指导具有参考价值。此外,多重RT-qPCR检测方法试验步骤简化、样本用量少、试验成本低,可为疫情防控争取宝贵时间。
-
表 1 Plackett-Burman试验参数列表
Table 1 List of Plackett-Burman test parameters
参数水平
Parameter level蔗段泊
松比(X1)
Poisson’s ratio of
sugarcane segment蔗段剪切
模量(X2)
Shear modulus of
sugarcane segment蔗段−蔗段
恢复系数(X3)
Sugarcane segment -
sugarcane segment
restitution coefficient蔗段−蔗段静
摩擦系数(X4)
Sugarcane segment -
sugarcane segment
static friction coefficient蔗段−蔗段滚
动摩擦系数(X5)
Sugarcane segment -
sugarcane segment
rolling friction coefficient蔗段−钢板
恢复系数(X6)
Sugarcane segment -
steel restitution
coefficient蔗段−钢板
静摩擦系数(X7)
Sugarcane segment -
steel static
friction coefficient蔗段−钢板滚
动摩擦系数(X8)
Sugarcane segment-
steel rolling
friction coefficient−1 0.30 3.0 0.30 0.30 0.01 0.20 0.10 0 0 0.34 10.8 0.45 0.55 0.35 0.45 0.35 0.03 1 0.38 18.6 0.60 0.80 0.06 0.70 0.60 0.05 表 2 Plackett-Burman试验方案及结果1)
Table 2 Plackett-Burman test scheme and results
序号
No.X1 X2 X3 X4 X5 X6 X7 X8 堆积角/(°)
Repose angle1 1 1 −1 1 1 1 −1 −1 49.39 2 −1 1 1 −1 1 1 1 −1 34.99 3 1 −1 1 1 −1 1 1 1 40.19 4 −1 1 −1 1 1 −1 1 1 48.03 5 −1 −1 1 −1 1 1 −1 1 31.80 6 −1 −1 −1 1 −1 1 1 −1 37.07 7 1 −1 −1 −1 1 −1 1 1 36.45 8 1 1 −1 −1 −1 1 −1 1 28.39 9 1 1 1 −1 −1 −1 1 −1 31.59 10 −1 1 1 1 −1 −1 −1 1 35.83 11 1 −1 1 1 1 −1 −1 −1 48.14 12 −1 −1 −1 −1 −1 −1 −1 −1 26.28 13 0 0 0 0 0 0 0 0 40.02 1) X1、X2分别表示蔗段泊松比和剪切模量;X3~X5分别表示蔗段−蔗段恢复系数、静摩擦系数和滚动摩擦系数;X6~X8分别表示蔗段−钢板恢复系数、静摩擦系数和滚动摩擦系数
1) X1 and X2 are Poisson’s ratio and shear modulus of sugarcane segment; X3−X5 are restitution coefficient, static friction coefficient and rolling friction coefficient of sugarcane segment-sugarcane segment, respectively; X6−X8 are restitution coefficient, static friction coefficient and rolling friction coefficient of sugarcane segment-steel, respectively表 3 Plackett-Burman试验参数显著性分析
Table 3 Analysis of significance of parameters in Plackett-Burman test
参数1)
Parameter平方和
Sum of square自由度
Degree of freedomF P2) X1 33.835 1 20.372 0.020* X2 5.727 1 3.448 0.160 X3 0.785 1 0.473 0.541 X4 398.477 1 239.917 0.001** X5 203.775 1 122.690 0.002** X6 1.680 1 1.012 0.389 X7 6.007 1 3.617 0.153 X8 3.819 1 2.300 0.227 1) X1、X2分别表示蔗段泊松比和剪切模量;X3~X5分别表示蔗段−蔗段恢复系数、静摩擦系数和滚动摩擦系数;X6~X8分别表示蔗段−钢板恢复系数、静摩擦系数和滚动摩擦系数;2) “*”“**”分别表示在0.05和0.01水平影响显著
1) X1 and X2 are Poisson’s ratio and shear modulus of sugarcane segment, respectively; X3−X5 are restitution coefficient, static friction coefficient and rolling friction coefficient of sugarcane segment-sugarcane segment, respectively; X6−X8 are restitution coefficient, static friction coefficient and rolling friction coefficient of sugarcane segment-steel, respectively. 2) “*” and “**” indicate significant effects at 0.05 and 0.01 levels respectively表 4 最陡爬坡试验方案设计及结果
Table 4 Design scheme and results of the steepest ascent search test
序号
No.蔗段泊松比(X1)
Poisson’s ratio of
sugarcane segment蔗段−蔗段静摩擦系数(X4)
Sugarcane segment-sugarcane
segment static friction coefficient蔗段−蔗段滚动摩擦系数(X5)
Sugarcane segment-sugarcane
segment rolling friction coefficient堆积角/(°)
Repose angle相对误差/%
Relative error1 0.300 0.30 0.01 32.09 24.85 2 0.316 0.40 0.02 34.98 18.08 3 0.332 0.50 0.03 40.36 5.49 4 0.348 0.60 0.04 42.97 0.62 5 0.364 0.70 0.05 43.80 2.56 6 0.380 0.80 0.06 48.32 13.16 表 5 Box-Behnken试验设计方案及结果1)
Table 5 Box-Behnken test design scheme and results
序号
No.X1 X4 X5 堆积角/(°)
Repose angle1 −1 0 −1 45.41 2 −1 0 1 46.64 3 1 0 −1 41.95 4 1 0 1 43.72 5 −1 −1 0 40.99 6 1 −1 0 39.20 7 −1 1 0 43.90 8 1 1 0 41.33 9 0 −1 −1 41.24 10 0 −1 1 41.91 11 0 1 −1 44.82 12 0 1 1 45.46 13 0 0 0 43.01 14 0 0 0 42.21 15 0 0 0 42.98 1) X1:蔗段泊松比;X4、X5分别表示蔗段−蔗段静摩擦系数和滚动摩擦系数
1) X1: Poisson’s ratio of sugarcane segment; X4 and X5 are static friction coefficient and rolling friction coefficient of sugarcane segment-sugarcane segment, respectively表 6 Box-Behnken试验回归模型方差分析
Table 6 Analysis of variance of the Box-Behnken test regression model
方差来源
Source of variance平方和
Sum of square自由度
Degree of freedom均方
Mean squareP1) 模型 Model 55.20 9 6.13 0.0031** X1 14.42 1 14.42 0.0015** X4 18.51 1 18.51 0.0008** X5 2.32 1 2.32 0.0530 X1X4 0.1521 1 0.1521 0.5470 X1X5 0.0729 1 0.0729 0.6736 X4X5 0.0002 1 0.0002 0.9811 X12 0.0863 1 0.0863 0.6472 X42 5.54 1 5.54 0.0114* X52 12.63 1 12.63 0.0020** 残差 Residual 1.82 5 0.3649 失拟项 Lack of fit 1.41 3 0.4710 0.3183 纯误差 Pure error 0.4113 2 0.2056 总和 Sum 57.03 14 1) “*”“**”分别表示在0.05和0.01水平影响显著
1) “*” and “**” indicate significant effects at 0.05 and 0.01 levels,respectively -
[1] 樊秋菊, 黄清玲, 吴合槟, 等. 国内外甘蔗收获机械化发展概况与前景[J]. 甘蔗糖业, 2020, 49(6): 1-11. doi: 10.3969/j.issn.1005-9695.2020.06.001 [2] 刘庆庭, 刘晓雪, 武涛, 等. 我国甘蔗机械化“十四五”规划与甘蔗产业可持续发展[J]. 现代农业装备, 2020, 41(6): 2-9. doi: 10.3969/j.issn.1673-2154.2020.06.001 [3] 邢浩男, 马少春, 王风磊, 等. 切段式甘蔗收割机排杂风机结构优化与试验[J]. 农业工程学报, 2020, 36(20): 67-75. doi: 10.11975/j.issn.1002-6819.2020.20.009 [4] HORABIK J, MOLENDA M. Parameters and contact models for DEM simulations of agricultural granular materials: A review[J]. Biosystems Engineering, 2016, 147(7): 206-225.
[5] OINONEN A, MARQUIS G. A constitutive model for interface problems with frictional contact and cohesion[J]. European Journal of Mechanics - A/Solids, 2015, 49(1): 205-213.
[6] 曾智伟, 马旭, 曹秀龙, 等. 离散元法在农业工程研究中的应用现状和展望[J]. 农业机械学报, 2021, 52(4): 1-20. doi: 10.6041/j.issn.1000-1298.2021.04.001 [7] 刘凡一, 张舰, 李博, 等. 基于堆积试验的小麦离散元参数分析及标定[J]. 农业工程学报, 2016, 32(12): 247-253. doi: 10.11975/j.issn.1002-6819.2016.12.035 [8] 侯占峰, 戴念祖, 陈智, 等. 冰草种子物性参数测定与离散元仿真参数标定[J]. 农业工程学报, 2020, 36(24): 46-54. doi: 10.11975/j.issn.1002-6819.2020.24.006 [9] 吴孟宸, 丛锦玲, 闫琴, 等. 花生种子颗粒离散元仿真参数标定与试验[J]. 农业工程学报, 2020, 36(23): 30-38. doi: 10.11975/j.issn.1002-6819.2020.23.004 [10] 武涛, 黄伟凤, 陈学深, 等. 考虑颗粒间黏结力的黏性土壤离散元模型参数标定[J]. 华南农业大学学报, 2017, 38(3): 93-98. doi: 10.7671/j.issn.1001-411X.2017.03.015 [11] 王黎明, 范盛远, 程红胜, 等. 基于EDEM的猪粪接触参数标定[J]. 农业工程学报, 2020, 36(15): 95-102. doi: 10.11975/j.issn.1002-6819.2020.15.012 [12] 刘文政, 何进, 李洪文, 等. 基于离散元的微型马铃薯仿真参数标定[J]. 农业机械学报, 2018, 49(5): 125-135. doi: 10.6041/j.issn.1000-1298.2018.05.014 [13] 张荣芳, 焦伟, 周纪磊, 等. 不同填充颗粒半径水稻种子离散元模型参数标定[J]. 农业机械学报, 2020, 51(S1): 227-235. doi: 10.6041/j.issn.1000-1298.2020.S1.026 [14] 鹿芳媛, 马旭, 谭穗妍, 等. 水稻芽种离散元主要接触参数仿真标定与试验[J]. 农业机械学报, 2018, 49(2): 93-99. doi: 10.6041/j.issn.1000-1298.2018.02.012 [15] 张涛, 刘飞, 赵满全, 等. 大豆种子与排种器接触物理参数的测定与离散元仿真标定[J]. 中国农业大学学报, 2017, 22(9): 86-92. doi: 10.11841/j.issn.1007-4333.2017.09.11 [16] 王云霞, 梁志杰, 张东兴, 等. 基于离散元的玉米种子颗粒模型种间接触参数标定[J]. 农业工程学报, 2016, 32(22): 36-42. doi: 10.11975/j.issn.1002-6819.2016.22.005 [17] 王美美, 王万章, 杨立权, 等. 基于响应面法的玉米籽粒离散元参数标定[J]. 华南农业大学学报, 2018, 39(3): 111-117. doi: 10.7671/j.issn.1001-411X.2018.03.017 [18] 马帅, 徐丽明, 袁全春, 等. 葡萄藤防寒土与清土部件相互作用的离散元仿真参数标定[J]. 农业工程学报, 2020, 36(1): 40-49. doi: 10.11975/j.issn.1002-6819.2020.01.005 [19] VAN ZEEBROECK M, TIJSKENS E, DINTWA E, et al. The discrete element method (DEM) to simulate fruit impact damage during transport and handling: Model building and validation of DEM to predict bruise damage of apples[J]. Postharvest Biology and Technology, 2006, 41(1): 85-91. doi: 10.1016/j.postharvbio.2006.02.007
[20] 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
[21] 于庆旭, 刘燕, 陈小兵, 等. 基于离散元的三七种子仿真参数标定与试验[J]. 农业机械学报, 2020, 51(2): 123-132. doi: 10.6041/j.issn.1000-1298.2020.02.014 [22] FAVIER J F, ABBASPOUR-FARD M H, KREMMER M, et al. Shape representation of axi-symmetrical, non-spherical particles in discrete element simulation using multi-element model particles[J]. Engineering Computations, 1999, 16(4): 467-480. doi: 10.1108/02644409910271894
[23] ABBASPOUR-FARD M H. Theoretical validation of a multi-sphere discrete element model suitable for biomaterials handling simulation[J]. Biosystems Engineering, 2004, 88(2): 153-161. doi: 10.1016/j.biosystemseng.2004.03.010
[24] KRUGGEL-EMDEN H, RICKELT S, WIRTZ S, et al. A study on the validity of the multi-sphere discrete element method[J]. Powder Technology, 2008, 188(2): 153-165. doi: 10.1016/j.powtec.2008.04.037
[25] SALOT C, GOTTELAND P, VILLARD P. Influence of relative density on granular materials behavior: DEM simulations of triaxial tests[J]. Granular Matter, 2009, 11(4): 221-236. doi: 10.1007/s10035-009-0138-2
[26] KOZHAR S, DOSTA M, ANTONYUK S, et al. DEM simulations of amorphous irregular shaped micrometer-sized titania agglomerates at compression[J]. Advanced Powder Technology, 2015, 26(3): 767-777. doi: 10.1016/j.apt.2015.05.005
[27] 韩丹丹, 张东兴, 杨丽, 等. 内充气吹式玉米排种器工作性能EDEM-CFD模拟与试验[J]. 农业工程学报, 2017, 33(13): 23-31. doi: 10.11975/j.issn.1002-6819.2017.13.004 [28] 原建博, 李骅, 吴崇友, 等. 基于离散单元法的水稻籽粒快速颗粒建模研究[J]. 南京农业大学学报, 2018, 41(6): 1151-1158. doi: 10.7685/jnau.201801004 [29] 王美美, 王万章, 杨立权, 等. 基于EDEM的玉米子粒建模方法的研究[J]. 河南农业大学学报, 2018, 52(1): 80-84. [30] 张涛, 刘飞, 赵满全, 等. 玉米秸秆接触物理参数测定与离散元仿真标定[J]. 中国农业大学学报, 2018, 23(4): 120-127. doi: 10.11841/j.issn.1007-4333.2018.04.15 [31] 夏蕊, 杨兆建, 李博, 等. 基于离散元法的煤散料堆积角试验研究[J]. 中国粉体技术, 2018, 24(6): 36-42. [32] 韩燕龙, 贾富国, 唐玉荣, 等. 颗粒滚动摩擦系数对堆积特性的影响[J]. 物理学报, 2014, 63(17): 173-179. [33] 罗菊川, 文晟, 李涵光, 等. 甘蔗尾茎泊松比的试验分析[J]. 华南农业大学学报, 2017, 38(6): 118-124. doi: 10.7671/j.issn.1001-411X.2017.06.018 [34] 罗菊川, 区颖刚, 刘庆庭. 甘蔗尾茎力学特性试验[J]. 农机化研究, 2016, 38(7): 220-223. doi: 10.3969/j.issn.1003-188X.2016.07.045 [35] 刘庆庭, 区颖刚, 卿上乐, 等. 甘蔗茎秆在扭转、压缩、拉伸荷载下的破坏试验[J]. 农业工程学报, 2006, 22(6): 201-204. doi: 10.3321/j.issn:1002-6819.2006.06.043 [36] 刘庆庭, 区颖刚, 袁纳新. 甘蔗茎在弯曲荷载下的破坏[J]. 农业工程学报, 2004, 20(3): 6-9. doi: 10.3321/j.issn:1002-6819.2004.03.002 [37] 王美美. 单芽段甘蔗种植机排种机理与种植均匀性研究[D]. 广州: 华南农业大学, 2014. [38] 温翔, 杨望, 郭无极, 等. 切段式甘蔗收割机排杂离散元仿真参数标定及验证[J]. 中国农机化学报, 2020, 41(1): 12-18. [39] 黄積佰. 切段式甘蔗收获机除杂系统的流固耦合仿真与优化[D]. 广州: 华南农业大学, 2020. -
期刊类型引用(7)
1. 杨冬,黄亚萍. 收割机脱粒装置模块化设计——基于经验预测和德尔菲法. 农机化研究. 2024(07): 106-110 . 百度学术
2. 旷娜,苗雪雪,邹丹,梁玉刚,唐启源,方宝华. 同等控温下头季稻和再生季稻稻米品质差异分析. 中国稻米. 2024(02): 43-46 . 百度学术
3. 吴泽栋,刘浩蓬,张国忠,周勇,宿衍强,刘冉. 振动挤压式荸荠收获机的设计与试验. 华中农业大学学报. 2024(02): 205-214 . 百度学术
4. 王芳,张伟,申珂,刘中峰. 基于改进鱼群算法的履带式收割机路径跟踪. 山西农业大学学报(自然科学版). 2024(04): 109-117 . 百度学术
5. 任学志,蔡孟媛. 脱粒分离装置的研究现状及发展趋势. 农机使用与维修. 2023(02): 35-37 . 百度学术
6. 杨德生,黄见良,彭少兵. 机收再生稻高产优质栽培技术研究进展. 中国稻米. 2023(05): 1-8 . 百度学术
7. 刘婉茹,周勇,徐红梅,付建伟,张妮,谢干,张国忠. 水稻鼓形纵轴流脱粒滚筒杆齿优化与试验. 农业工程学报. 2023(15): 34-45 . 百度学术
其他类型引用(5)