An Indirect ELISA Assay for Detecting Specific SIgA Antibody of Porcine Epidemic Diarrhea Virus
-
摘要:目的 快速了解猪群感染猪流行性腹泻病毒或免疫猪流行性腹泻病毒相关疫苗后猪群特异性SIgA抗体水平。方法 对猪流行性腹泻病毒（Porcine epidemic diarrhea virus，PEDV）ZJ08株进行纯化，采用纯化的病毒作为包被抗原，以辣根过氧化物酶标记的SC（Secretory component，SC）多克隆抗体为标记抗体，通过优化检测条件，建立乳汁中猪流行性腹泻病毒特异性SIgA抗体间接ELISA检测方法。结果 成功建立了猪流行性腹泻病毒SIgA抗体检测方法。PEDV抗原的最佳包倍浓度为0.6 μg·孔−1，样本最佳稀释比例为1∶5，作用时间1 h；SC标记抗体最适稀释比例为1∶100，作用时间1.5 h。建立的ELISA方法能够特异性检测样本中抗PEDV SIgA抗体，与PoRV（Porcine rotavirus，PoRV）、TGEV（Transmissible gastroenteritis virus，TGEV）的阳性乳汁均无交叉反应；批内、批间重复试验的变异系数小于10%，具有良好的特异性和重复性。通过临床采集的乳汁样品测定发现，样本中PEDV IgA抗体阳性率达到98.1%，而采用针对特异性SIgA的抗体方法检测，其阳性率仅为69.2%，样本中IgA和SIgA抗体符合率仅为70.6%，提示仅仅依据PEDV IgA抗体水平并不能很好地反映猪群免疫该疫苗后特异性黏膜免疫的真实水平。结论 建立的PEDV特异性SIgA抗体检测方法可为分析、评估猪群免疫PEDV疫苗后乳汁中特异SIgA抗体产生及猪群的保护情况提供可靠依据。Abstract:Objective A method for detecting SIgA antibodies of porcine epidemic diarrhea virus (PEDV) was established.Methods Using a purified PEDV ZJ08 strain as coating antigen and the antibody for secretory component(SC) labeled by horseradish peroxidase (HRP) as secondary antibody, an indirect ELISA assay with optimized conditions was developed to detect the specific SIgA antibody of PEDV.Result The optimized assay applied a working antigen concentration of 0.6 μg per well, a 1∶5 sample dilution, a coating at 37 ℃ for 1 h, and an 100× dilution of the labeled antibody for 1.5 h incubation at 37 ℃. Other than the target detection, the assay did not react with porcine rotavirus or transmissible gastroenteritis virus. The inter- or intra-batch variations of the assay on test results were less than 10%. On the milk from vaccinated pigs, a 98.1% positive detection on IgA, which was significantly higher than 69.2% on SIgA and merely 70.6% on combined IgA and SIgA, was achieved indicating the deficiency of a single IgA test for an accurate determination on the mucosal immunity of the animal.Conclusion The established assay appeared to appropriately reveal the specific SIgA antibody produced in pigs for protection of the livestock from PEDV by the vaccination.
-
Keywords:
- Porcine epidemic diarrhea virus /
- indirect ELISA /
- SIgA antibody
-
0. 引言
【研究意义】 蝴蝶兰（Phalaenopsis equestris）花型独特、花色淡雅、花期持久，具有较高的观赏价值和市场经济价值，是我国销售量最大的盆花花卉，是目前市场上非常有竞争力和发展潜力的一种价值较高的花卉。但是病毒性病害严重影响和危害了蝴蝶兰的发展前景，已经成为蝴蝶兰产业持续、健康和快速发展的一个巨大绊脚石。福建漳州蝴蝶兰花卉产业在全国具有举足轻重的重要地位。福建漳州作为全国最大的兰花培育集散中心，年产出几千万盆兰花，其中蝴蝶兰约500万～800万盆，投放市场的成品兰占全国兰花市场份额的10%，同时也已成为向美国、日本、韩国和欧洲等国出口蝴蝶兰的重要基地[1]。植物病毒病被称“植物癌症”，一旦发病，尚无有效方法挽救，防治非常困难。近几年，由于全球气候变化、生态环境恶化、农业种植结构、耕种模式改变等原因，蓟马和蚜虫等虫媒传播的花卉病毒病日益严重，对我国蝴蝶兰等重要花卉的安全生产造成了严重的威胁。因此，亟需加强病毒-蝴蝶兰的互作研究，从而为制定防治策略提供理论依据。【前人研究进展】目前，研究表明感染蝴蝶兰植株的病毒种类较多，主要包括建兰花叶病毒（cymbidium mosaic virus, CymMV）[2]、齿兰环斑病毒（odontoglossum ringspot virus, ORSV）[3]、黄瓜嵌纹病毒（cucumber mosaic virus, CMV）[4]、番茄斑萎病毒（tomato spotted wilt virus, TSWV）[5]、建兰环斑病毒（cymbidium ringspot virus, CymRV）、辣椒褪绿病毒（capsicum chlorosis virus, CaCV）[6]、凤仙花坏死斑病毒（impetiens necrotic spot virus, INSV）[7]、兰花斑点病毒（ovrchid fleck rhabdovirus, OFV）[8]、蝴蝶兰黄化斑点病毒 （phalaenopsis chlorosisspot virus, PhCSV）[9]、康乃馨斑驳病毒（carnation mottle virus, CaMV）[10]和石斛兰叶脉坏疽病毒（dendrobium veinnecrosis closterovirus, DVNV） [11]等。其中CymMV和ORSV是蝴蝶兰生产中的主要病毒，且两者经常复合侵染蝴蝶兰植株。CymMV和ORSV分别属于甲型线形病毒科（Alphaflexiviridae） 和帚状病毒科（Virgaviridae）的正单链RNA病毒[12]。病毒感染蝴蝶兰会影响植株的生长发育，具体表现为生长速度减缓，叶片出现规则或者不规则的褪绿黄化、炭疽坏死斑块，花器出现畸形、色泽不均的斑块、提早凋萎等不同症状[13−14]，对蝴蝶兰影响巨大。RNA干扰（RNA interference, RNAi）是植物、昆虫和真菌等生物保守的抵御病毒等病原物侵染的免疫防御机制[15−16]。研究表明，病毒侵染模式植物拟南芥后，植物的Dicer-like (DCL)蛋白识别并切割长的病毒来源的双链基因组为21、22、24 nt等不同长度的小干扰RNA（small interference RNA, siRNA），这些病毒来源的siRNA（vsiRNAs）与Argonaute结合形成RNA诱导的沉默复合体，随后复合物中vsiRNAs通过碱基匹配进行基因定位，Argonaute蛋白通过核酸切割活性沉默基因[17]，从而沉默外源病毒等病原物基因的表达，由此可知，RNAi免疫途径被激活或者病毒与寄主植物发生互作的一个典型特征是病毒侵染寄主植物细胞后能够诱导产生大量病毒来源的vsiRNAs。由于病原系统的差异，病毒vsiRNAs在丰度、长度分布、5′端碱基偏好性、正负义来源、基因组高低频切割位点等方面呈现广泛的多样性[15]。不同种的病毒侵染不同寄主细胞后产生的 vsiRNAs 丰度存在较大差异，丰度一般在0.1%～20%[15]。vsiRNAs的长度分布主要集中在21 nt和22 nt[15]。vsiRNAs 5′端碱基偏好性根据vsiRNAs结合的AGO蛋白而有所不同[15]。vsiRNAs近乎等比例来源于正负义基因组或者主要来源于正义或负义基因组[15]。vsiRNAs在病毒基因组上呈现连续但不均匀的分布，即形成高低频切割位点[15]。 【本研究切入点】运用深度测序技术，已经在众多不同的植物细胞中鉴定到大量病毒来源的vsiRNAs，vsiRNAs的鉴定已成为研究病毒和寄主互作的有效途径或方法[18−19]。然而关于蝴蝶兰与病毒病原之间的互作尚有待深入探讨。【拟解决的关键问题】通过电镜和PCR等技术鉴定CymMV和ORSV感染蝴蝶兰，并通过深度测序技术鉴定分析病毒来源的vsiRNAs特征来分析CymMV、ORSV和蝴蝶兰寄主之间的相互作用，以期为防治蝴蝶兰病毒病提供依据。
1. 材料与方法
1.1 试验材料
病毒侵染和健康的蝴蝶兰样品均由福建漳州台湾农民创业园绿乐园艺有限公司提供。多糖多酚植物总RNA提取试剂盒（RNAprep pure plant plus kit）、一步法RT-PCR试剂盒（fastking one step RT-PCR kit）和DNA Marker分子量标准（MD102）均购自天根生化科技（北京）有限公司。
1.2 RT-PCR
根据CymMV和ORSV外壳蛋白基因序列（Accession登录号分别为: AB693982和AB693988）设计特异扩增引物CymMV-F：ATGGGAGAGCCCACTCCAACT，CymMV-R：TATCCTCCTGGAAACCAGCCT；ORSV-F：ATGTTTTACACTATTACAGA，ORSV-R：TGAGACTTGATTGTACATACCA，两个基因片段大小均约400 bp。用多糖多酚植物总RNA提取试剂盒（RNAprep pure plant plus kit）提取病毒侵染和健康蝴蝶兰叶片总RNA后，根据一步法RT-PCR试剂盒（fastking one step RT-PCR kit）说明进行RT-PCR检测。一步RT-PCR体系（50 μL）为：2×fastking one step RT-PCR mastermix 25 μL，25×RT-PCR Enzyme Mix 2 μL，上下游引物（10 μmol·L−1）各1.25 μL，RNA模版0.01～1.00 μg，ddH2O补至50 μL。RT-PCR程序：42 ℃，30 min；95 ℃，3 min；94 ℃，30 s，55～65 ℃，30 s，72 ℃，30 s，35个循环；72 ℃，5 min。RT-PCR扩增产物用琼脂糖凝胶电泳检测。
1.3 电镜观察
电镜负染观察主要步骤：获取病毒侵染和健康蝴蝶兰叶片的汁液，将汁液滴加在覆有Formvar膜的铜网上，2 min后用滤纸吸去铜网膜上多余的液体，用1%磷钨酸（pH 7.2）进行负染，用透射电子显微镜H-
7650 Hitachi（日立）观察病毒粒体的形态和大小等特征。电镜超薄切片观察主要步骤：病毒侵染和健康蝴蝶兰叶片切成0.5 cm的小块后，依次用2.5%戊二醛进行固定，用酒精梯度脱水，用Spurr试剂进行渗透，用环氧树脂进行包埋；随后在聚合器中70 ℃进行聚合；聚合后的样品在超薄切片机Leica UC7（德国莱卡）上进行切片制备；超薄切片转移至铜网进行染色；最后将铜网置于透射电镜H-
7650 Hitachi （日立）观察病毒粒体的形态和大小等特征。1.4 siRNAs深度测序和特征分析
参照Lan的方法[18]进行siRNAs的深度测序和鉴定分析。主要步骤为：用多糖多酚植物总RNA提取试剂盒（RNA preppure plant plus kit）提取病毒侵染和健康蝴蝶兰叶片总RNA，10 μg总RNA经8%聚丙烯酰胺凝胶电泳后，回收纯化18～32 nt的siRNAs，用T4 RNA酶分别连接3′接头和5′接头后，SuperScript Ⅱ酶反转录后进行PCR反应，通过聚丙烯酰胺凝胶电泳回收纯化PCR产物，Illumina HiSeq 2000测序仪（Illumina，美国）进行深度测序。用软件Bowtie将siRNAs与CymMV和ORSV基因组进行匹配，再进行病毒来源vsiRNAs丰度、长度、碱基偏好性和正负义链分布等特征分析。
2. 结果与分析
2.1 CymMV和ORSV在蝴蝶兰中的侵染
采集蝴蝶兰感病样品（图1A），以蝴蝶兰叶片总RNA为模板，CymMV-F/R和ORSV-F/R为引物，RT-PCR扩增两种病毒的外壳蛋白基因部分片段，经琼脂糖凝胶电泳检测。结果显示，RT-PCR扩增获得了两条特异的目的片段，片段大小约400 bp（图1B），与引物设计预期结果一致，初步表明CymMV和ORSV两种病毒侵染蝴蝶兰。
![]()
图 1 CymMV和ORSV在蝴蝶兰叶片的共侵染鉴定A：病毒侵染的蝴蝶兰样品；B：CymMV和ORSV外壳蛋白基因的RT-PCR扩增产物的琼脂糖凝胶电泳， M为DNA marker MD102；1和2分别为 CymMV和ORSV外壳蛋白基因产物； C：CymMV和ORSV病毒电镜观察的负染图；D：CymMV和ORSV病毒电镜观察的超薄切片图。Figure 1. Identification of CymMV/ORSV co-infections in P. equestrisA: P. equestris plant infected by viruses; B: agarose gel electrophoresis of RT-PCR products of CP genes of CymMV and ORSV; M: DNA marker MD102; 1 and 2: RT-PCR products of CP genes of CymMV and ORSV, respectively; C: electron microscopic negative staining images of CymMV and ORSV; D: electron microscopic images on thin sections of CymMV and ORSV.为了进一步验证RT-PCR结果的准确性和可靠性，采集蝴蝶兰感病叶片进行负染和超薄切片制备，然后用电镜观察结果。电镜结果表明，负染和超薄切片中都能够观察到长约300 nm线状的病毒粒子和长约500 nm杆状的病毒粒子（图1C、D），300 nm线状的病毒粒子和500 nm杆状的病毒粒子分别代表CymMV和ORSV两种病毒，在健康蝴蝶兰叶片样品中并未观察到该种大小和形态特征的病毒，因此电镜结果进一步证实CymMV和ORSV两种病毒侵染了蝴蝶兰。
2.2 CymMV和ORSV来源的vsiRNAs特征分析
2.2.1 siRNAs深度测序总体概况
病毒侵染证实后，提取感病和健康植株的总RNA进行siRNAs的深度测序，鉴定和分析病毒来源的vsiRNAs特征。siRNAs深度测序结果显示，健康和病毒侵染的蝴蝶兰植株分别获得7 563 892和6 133 689个读数的原始种类的小RNA（total siRNAs）；将这些total siRNAs匹配到CymMV和ORSV两种病毒基因组后，结果表明，健康和病毒侵染的蝴蝶兰植株中能匹配到病毒基因组的total vsiRNAs读数分别为906和826 682；在病毒侵染的蝴蝶兰样品中，匹配到CymMV和ORSV病毒基因组的total vsiRNAs读数分别为369 524（约44.7%）和457 104（55.3%）。
2.2.2 vsiRNAs的长度特征
获得siRNAs读数后进行特征分析，长度特征分析表明，在健康蝴蝶兰样品中，24、21、22 、23 nt长度的total siRNAs丰度依次较高（图2A），然而在病毒侵染的蝴蝶兰样品中，21、24、22、23 nt长度的total siRNAs丰度依次较高（图2A），total siRNAs最大丰度的差异表明病毒侵染改变了total siRNAs的长度分布形式或特征。另外，在病毒侵染的蝴蝶兰样品中，CymMV和ORSV来源的vsiRNAs长度丰度均主要为21 nt和22 nt（图2B），其中21 nt vsiRNAs丰度均为最高，22 nt vsiRNAs次之，该结果与Pai等的研究结果一致[19]。21 nt vsiRNAs高丰度不仅表明21 nt vsiRNAs在RNAi免疫过程中起主要的抵御病毒作用，还暗示在蝴蝶兰细胞中DCL4和DCL2蛋白可能是vsiRNAs分子形成过程中的主要参与者。vsiRNAs是双链小RNA分子，vsiRNAs 5′或3′端突出的碱基数量称为悬挂（overhang）。对于21 nt的vsiRNAs分子，悬挂的碱基数量=21−碱基匹配数量。21 nt vsiRNAs的overhang分析表明，CymMV和ORSV来源的21 nt的vsiRNAs中，具有2个碱基overhang的vsiRNAs丰度都最高（图2C），暗示在蝴蝶兰细胞中，CymMV和ORSV来源的21 nt的vsiRNAs中具有2个碱基overhang的vsiRNAs可能是RNAi免疫反应最有效的激发者。
![]()
图 2 蝴蝶兰细胞中CymMV和ORSV来源的vsiRNAs特征A：健康蝴蝶兰和病毒侵染蝴蝶兰中siRNAs长度分布；B： 病毒侵染蝴蝶兰中vsiRNAs长度分布； C：病毒侵染的蝴蝶兰中5'端距离为1~21 nt vsiRNAs读数；D： Total vsiRNAs（左）和21 nt vsiRNAs（右） 5'端第一个碱基偏好性。Figure 2. Characterization of CymMV- and ORSV-derived vsiRNAs in P. equestrisA: size distribution of total small RNAs in virus-infected and viruses-free P. equestris plants; B: size distribution of vsiRNAs matching CymMV and ORSV genomes in virus-infected P. equestris plants; C: reads of 21-nt vsiRNAs with 1-21 nt distance between 5′ ends from CymMV and ORSV genome in P. equestris plants; D: relative frequency of 5′-terminal nucleotide of total vsiRNAs (left) and 21 nt vsiRNAs (right).2.2.3 vsiRNAs的碱基偏好性
vsiRNAs 5′端的第一个碱基种类会影响vsiRNAs结合到哪种Argonaute蛋白。Total vsiRNAs 5′端碱基偏好性分析表明，CymMV和ORSV来源的vsiRNAs 5′端碱基偏好性都是依次偏好于U、A、C和G（图2D，左Panel），该结果表明CymMV和ORSV来源的vsiRNAs因第一个碱基的差异其可能结合各自不同的Argonaute蛋白，进而介导RNAi免疫效应。长度丰度分析表明，21 nt vsiRNAs丰度最高，其可能是蝴蝶兰细胞RNAi免疫途径的主要激发者，因此，进一步分析了21 nt vsiRNAs的5′端第一个碱基偏好性分析，分析结果表明，21 nt vsiRNAs的5′端第一个碱基偏好性也是依次偏好于U、A、C和G（图2D，右Panel），该结果进一步暗示vsiRNAs 5′端第一个碱基的类型与其结合的Argonaute蛋白之间具有密切的联系。
2.2.4 vsiRNAs的正负义链来源
total vsiRNAs正负义链来源（极性）分布表明，CymMV来源的vsiRNAs中，分别有
326662 个读数（88.4%）和42 862个读数（11.6%）的vsiRNAs来自于正义链和负义链病毒基因组（图3A）；与此类似，ORSV来源的vsiRNAs中，分别有425 548个读数（93.1%）和31 520个读数（6.9%）的vsiRNAs来自于正义链和负义链病毒基因组（图3A），该结果表明蝴蝶兰细胞中CymMV和ORSV vsiRNAs主要来源都是病毒基因组的正义链。![]()
图 3 CymMV和ORSV vsiRNAs的生成前体分析A：CymMV和ORSV vsiRNAs的正负义链分布；B：vsiRNAs在CymMV基因组上高低频切割位点；C：vsiRNAs在ORSV基因组上高低频切割位点； D：CymMV正链基因组的二级结构预测；E：ORSV正链基因组的二级结构预测。方框I为图3B和图3D中高频切割位点和基因组茎环结构的对应关系。方框II为图3C和图3E中高频切割位点和基因组茎环结构的对应关系。Figure 3. Biogenesis precursors of CymMV and ORSV vsiRNAsA: polarity distribution of vsiRNAs matching CymMV and ORSV genomes from co-infected P. equestris plants; B and C: vsiRNAs hotspots and cold spots distribution along CymMV and ORSV genomes, respectively; D and E: secondary structures of CymMV and ORSV RNAs predicted with RNAfold server. Box I shows the corresponding relationship between the high frequency cutting sites in Figures 3B and the stem ring structure in Figure 3D. Box II shows the corresponding relationship between the high frequency cutting sites in Figures 3C and the stem ring structure in Figure 3E.2.2.5 高低频切割位点及其与基因组结构的关系分析
将total vsiRNAs定位到CymMV和ORSV病毒基因组具体位点后，高低频切割位点分析表明，vsiRNAs沿着CymMV和ORSV病毒基因组呈现连续但不均匀的分布特征（图3B、C）；大部分（95%）total vsiRNAs产生位点位于CymMV和ORSV病毒基因组内部，而不是基因组的5′和3′两端（图3B、3C）；在病毒基因组不同位点形成了高频切割点（hotspot）和低频切割点（coldspot）（图3B、C）。
用RNAfold server在线平台预测CymMV和ORSV病毒基因二级结构，结构预测结果表明在CymMV和ORSV病毒基因内部形成了大量的茎环二级结构（图3D、E），且茎环二级结构的形成与高频切割位点是相对应的，表明CymMV和ORSV病毒基因内部茎环二级结构是RNAi免疫途径关键因子DCL蛋白的切割底物，从而形成不同大小的vsiRNAs，即茎环结构是vsiRNAs的主要生成前体。
3. 讨论
蝴蝶兰是具有重要观赏和经济价值的花卉之一。但是病毒性病害已经成为制约蝴蝶兰产业持续健康发展的较大因素。福建漳州蝴蝶兰花卉产业在全国具有举足轻重的重要地位，但是CymMV和ORSV两种病毒对蝴蝶兰危害严重[12−14]。因此，亟需加强蝴蝶兰病毒性病害尤其是CymMV和ORSV两种病毒的检测调查。本研究采集了漳州漳浦农民产业园区的健康和疑似病毒侵染的蝴蝶兰植株，依次采用RT-PCR和传统的电镜技术明确了漳州蝴蝶兰植株感染了CymMV和ORSV两种主要病毒。
植物病毒素有“植物癌症”之称，目前并无有效的防治办法。为了探索有效的防治措施，务必加强病毒和寄主之间相互作用的研究。RNAi是植物、昆虫和真菌等生物保守的抵御病毒等病原物侵染的免疫防御机制[15−16]。RNAi免疫途径被激活的一个典型特征是病毒侵染寄主植物细胞后能够诱导产生大量病毒来源的vsiRNAs。因此，病毒来源的vsiRNAs的鉴定和特征分析已经成为加深理解病毒和寄主相互作用的重要途径和方法。运用高通量深度技术已经在众多的病毒-植物病原系统中鉴定了大量的vsiRNAs[18]。在本研究中，利用siRNAs深度测序技术在CymMV和ORSV复合侵染的蝴蝶兰叶片中获得了大量（约13.48%）病毒来源的vsiRNAs，该结果表明CymMV和ORSV复合侵染蝴蝶兰叶片细胞后激活了寄主植物RNAi免疫途径，病毒与寄主确实存在相互作用。
病毒与寄主互作的具体机制可以通过分析vsiRNAs特征加以理解。研究表明24 nt siRNA主要在miRNA途径中调控寄主本身的生长发育等生命活动中发挥作用，而21 nt的siRNAs主要在siRNA介导的抗病毒免疫途径中发挥作用[15]。本研究中，高丰度值21 nt siRNAs的出现表明CymMV和ORSV的侵染激活了蝴蝶兰细胞中siRNAs介导的抗病毒免疫途径。在拟南芥等模式植物细胞中，21、22、24 nt的vsiRNAs分别由DCL4、DCL2和DCL1蛋白切割底物产生[17]，因此推测在蝴蝶兰细胞中该长度的vsiRNAs可能也有其同源蛋白切割底物产生。研究表明许多病原系统产生的21 nt的vsiRNAs中，其5′或3′端具有2个碱基的突出（overhang）siRNAs丰度最高[20−22]，与此相似，本研究的蝴蝶兰-病毒系统中，21 nt的vsiRNAs中，5′或3′端具有2个碱基的突出（overhang）vsiRNAs丰度也最高，该结果暗示CymMV和ORSV来源的21 nt的vsiRNAs中具有2个碱基overhang的vsiRNAs可能是RNAi免疫反应最有效的激发者。
研究表明siRNAs 5′端的碱基在介导siRNAs与Argonaute蛋白的结合密切相关[23]。比如，在拟南芥模式植物中，5′端为碱基A的siRNAs 更易与Argonaute2和Argonaute4结合，5′端为碱基U的siRNAs 更易与Argonaute1结合，5′端为碱基C的siRNAs 更易与Argonaute5结合，5′端为碱基G的siRNAs 更易与其他Argonaute结合。本研究中vsiRNAs 5′端碱基偏好性分析表明，CymMV和ORSV来源的vsiRNAs中，5’端碱基偏好性都是依次偏好于U、A、C和G；另外21 nt vsiRNAs 5′端碱基偏好性同样是依次偏好于U、A、C和G，该碱基偏好性结果与其他病原系统的vsiRNAs碱基偏好性一致[24−25]；该结果表明CymMV和ORSV来源的vsiRNAs可能结合各自不同的Argonaute蛋白，进而介导RNAi免疫效应，但是由于目前蝴蝶兰基因组无法获得，也就无法获得蝴蝶兰相关的Argonaute蛋白信息，因此后续获得蝴蝶兰基因组中编码Argonaute的基因序列是一项非常重要的工作。
在不同病毒激活寄主RNAi免疫途径中，vsiRNAs产生的前体各不相同[15]。vsiRNAs的形成前体主要通过分析其在病毒基因组正负义链分布、核苷酸产生的位点以及频率分布来推测。在本研究发现，CymMV和ORSV复合侵染的蝴蝶兰细胞中，CymMV和ORSV来源的vsiRNAs主要来自于病毒基因组的正义链，该结果与其他一些正单链RNA病毒在寄主细胞中形成的vsRNAs的正负义偏好性相反，这些病原系统中vsiRNAs大约等比例来源于病毒基因组正负义链，因此他们推测病毒基因组复制过程中形成的双链为vsiRNAs的形成前体[26−29]。进一步将 vsiRNAs定位到CymMV和ORSV病毒基因组具体位点后，高低频切割位点分析表明，vsiRNAs沿着CymMV和ORSV病毒基因组呈现持续但不均匀的分布特征；大部分（95%）total vsiRNAs产生位点位于CymMV和ORSV病毒基因组内部，而不是基因组的5′和3′两端；在病毒基因组不同位点形成了高频切割点和低频切割点。随后用RNAfold server在线平台预测CymMV和ORSV病毒基因二级结构，结构预测结果表明在CymMV和ORSV病毒基因内部形成了大量的茎环二级结构，且茎环二级结构的形成与高频切割位点相对应，因此推测CymMV和ORSV病毒基因组正义链通过碱基匹配自身折叠产生的 pre-miRNA（Precursor miRNA）的茎环结构区域为DCL蛋白的切割底物，即vsiRNAs产生的前体，该预测结果与Flynt和Molnár 等 [30−31] 的研究结果一致。
4. 结论
本研究首先通过电镜和PCR技术鉴定了CymMV和ORSV两种病毒在蝴蝶兰植株中的复合侵染；然后通过深度测序技术鉴定分析了来源于CymMV和ORSV两种病毒基因组的vsiRNAs的丰度、长度、悬挂、碱基偏好性、正负义链分布、Hotspot和Coldspot等特征，从而加深了蝴蝶兰抗病毒RNAi免疫途径防御CymMV和ORSV两种病毒的机制的理解，对研究蝴蝶兰病毒性病害的防治措施具有重要的理论意义。
-
![]()
图 1 不同浓度蔗糖层PEDV RT-PCR检测结果
M：1000 DNA Marker；1～3：40%～60%层带未稀释、100倍稀释、1000倍稀释；4～6为30%层带未稀释、100倍稀释、1000倍稀释；7～9为20%层带未稀释、100倍稀释、1000倍稀释。
Figure 1. RT-PCR detections of PEDV on sucrose layers
M：DNA ladder DL1000; 1–3: on not diluted 40%–60% sucrose layer and on those diluted 100× and 1 000×, respectively; 4–6: on not diluted 30% sucrose layer and on those diluted 100× and 1 000×, respectively; 7–9 on not diluted 20% sucrose layer and on those diluted 100× and 1 000×, respectively.
![]()
图 2 Western Blot分析PET32a-SC重组蛋白的表达
M：蛋白 Marker；1：SC 蛋白通过 His 标签和羊抗鼠二抗的反应检测。
Figure 2. PET32a-SC protein expression by western blot
M: Protein marker; 1: SC protein tested by His tag and goat anti-mouse secondary antibody.
![]()
图 3 PET32a-SC蛋白SDS-PAGE分析
M：蛋白 marker；1、2：纯化后 PET32a-SC。
Figure 3. SDS-PAGE on expressed PET32a-SC protein
M: Protein marker; 2-3: purified protein of PET32a-SC.
表 1 抗原最佳包被浓度及乳汁最佳稀释度的确定（OD450 nm）
Table 1 Optimal antigen concentration and milk dilution determined according to OD450 nm
样本稀释度
Sample dilution抗原稀释度
Antigen dilution/(μg·孔−1)P N P/N 1∶5 1.2 0.821 0.130 6.32 0.6 0.663 0.083 7.99 0.3 0.627 0.096 6.53 0.15 0.584 0.075 7.79 1∶10 1.2 0.725 0.131 5.53 0.6 0.572 0.096 5.96 0.3 0.494 0.093 5.31 0.15 0.455 0.077 5.91 1∶20 1.2 0.69 0.132 5.23 0.6 0.563 0.105 5.36 0.3 0.472 0.115 4.10 0.15 0.436 0.099 4.40 1∶30 1.2 0.519 0.128 4.05 0.6 0.412 0.107 3.85 0.3 0.324 0.117 2.77 0.15 0.271 0.111 2.44 
下载: 导出CSV
表 2 最佳包被条件的确定
Table 2 Optimum conditions for coating
指标 Index 包被条件 Coating condition 4 ℃包被过夜
4 ℃ wrap overnight37 ℃ 1 h 37 ℃ 2 h P 0.657 0.572 0.620 N 0.130 0.082 0.081 P/N 5.05 6.98 7.65
下载: 导出CSV
表 3 最佳封闭条件的确定
Table 3 Optimum conditions for blocking
指标 Index 封闭浓度
Closure condition封闭条件（37 ℃）
Closure condition1% 2% 5% 1 h 2 h 3 h P 1.180 1.121 1.034 1.110 1.070 1.077 N 0.241 0.213 0.220 0.282 0.226 0.255 P/N 4.89 5.26 4.68 3.94 4.73 4.22
下载: 导出CSV
表 4 最佳乳汁样本工作条件的确定
Table 4 Optimum working conditions for milk specimens
指标
Index孵育时间 Incubation time 37 ℃ 45 min 37 ℃ 60 min 37 ℃ 90 min P 0.638 0.639 0.654 N 0.096 0.085 0.087 P/N 6.65 7.52 7.52
下载: 导出CSV
表 5 最佳标记抗体工作条件的确定
Table 5 Optimum working conditions for HRP-labeled secondary antibody
指标
Index稀释比例
dilution作用时间
Reaction conditions1∶50 1∶100 1∶200 1∶300 37 ℃ 45 min 37 ℃ 60 min 37 ℃ 90 min P 0.819 0.891 0.709 0.706 0.746 0.796 0.867 N 0.168 0.162 0.172 0.159 0.146 0.150 0.160 P/N 4.88 5.5 4.12 4.44 5.11 5.31 5.42
下载: 导出CSV
表 6 底物最佳反应时间的确定
Table 6 Optimum reaction time for TMB
指标 Index 孵育时间 Incubation time 15 min 20 min 25 min 30 min P 0.869 0.934 1.019 1.089 N 0.136 0.136 0.215 0.187 P/N 6.39 6.87 4.93 5.82
下载: 导出CSV
表 7 乳汁样本的重复性检验结果
Table 7 Repeatability of whey samples
编号
Number组内重复
Intra group repetition组间重复
Intergroup duplication¯X±SD 变异系数CV/% ¯X±SD 变异系数CV/% 1 0.167±0.009 0.3 0.168±0.040 1 2 0.225±0.005 0.2 0.223±0.052 2 3 0.249±0.006 0.2 0.252±0.009 1 4 0.352±0.017 0.6 0.342±0.032 2 5 0.304±0.006 0.2 0.298±0.028 1 6 0.340±0.009 0.3 0.335±0.025 2 7 0.159±0.040 1.3 0.162±0.059 3 8 0.641±0.002 0.1 0.623±0.066 6 9 0.174±0.021 0.8 0.165±0.032 5 10 0.161±0.013 0.5 0.157±0.035 1
下载: 导出CSV
-
[1] 杨朋霞, 闫若潜, 王华俊, 等. 基于猪流行性腹泻病毒N蛋白间接ELISA抗体检测方法的建立及初步应用 [J]. 中国预防兽医学报, 2019, 41(10):1021−1026. YANG P X, YAN R Q, WANG H J, et al. Establishment and preliminary application of an indirect ELISA for detecting antibody based on N protein of porcine epidemic diarrhea variant virulent strain [J]. Chinese Journal of Preventive Veterinary Medicine, 2019, 41(10): 1021−1026.(in Chinese)
[2] 陈建飞. 猪流行性腹泻病毒全基因组序列分析及感染性cDNA克隆的构建[D]. 北京: 中国农业科学院, 2012: 3. CHEN J F. Analysis of complete genome sequences and construction of an infectious cDNA clone of porcine epidemic diaririea virus[D]. Beijing: Chinese Academy of Agricultural Sciences, 2012: 3. (in Chinese)
[3] LI Y C, WU Q X, HUANG L L, et al. An alternative pathway of enteric PEDV dissemination from nasal cavity to intestinal mucosa in swine [J]. Nature Communications, 2018, 9: 3811. DOI: 10.1038/s41467-018-06056-w
[4] LI W T, LI H, LIU Y B, et al. New variants of porcine epidemic diarrhea virus, China, 2011 [J]. Emerging Infectious Diseases, 2012, 18(8): 1350−1353. DOI: 10.3201/eid1803.120002
[5] CHEN Q, LI G W, STASKO J, et al. Isolation and characterization of porcine epidemic diarrhea viruses associated with the 2013 disease outbreak among swine in the United States [J]. Journal of Clinical Microbiology, 2014, 52(1): 234−243. DOI: 10.1128/JCM.02820-13
[6] SONG Q Y, STONE S, DREBES D, et al. Characterization of anti-porcine epidemic diarrhea virus neutralizing activity in mammary secretions [J]. Virus Research, 2016, 226: 85−92. DOI: 10.1016/j.virusres.2016.06.002
[7] JUNG K, SAIF L J, WANG Q H. Porcine epidemic diarrhea virus (PEDV): An update on etiology, transmission, pathogenesis, and prevention and control [J]. Virus Research, 2020, 286: 198045. DOI: 10.1016/j.virusres.2020.198045
[8] PARK J E, KANG K J, RYU J H, et al. Porcine epidemic diarrhea vaccine evaluation using a newly isolated strain from Korea [J]. Veterinary Microbiology, 2018, 221: 19−26. DOI: 10.1016/j.vetmic.2018.05.012
[9] POONSUK K, GIMÉNEZ-LIROLA L G, ZHANG J Q, et al. Does circulating antibody play a role in the protection of piglets against porcine epidemic diarrhea virus? [J]. PLoS One, 2016, 11(4): e0153041. DOI: 10.1371/journal.pone.0153041
[10] 陈冰清, 葛艳, 司伏生, 等. 猪流行性腹泻病毒诱导的猪肠道-乳腺-sIgA轴以及免疫机制探讨 [J]. 微生物学通报, 2021, 48(11):4240−4249. DOI: 10.13344/j.microbiol.china.210041 CHEN B Q, GE Y, SI F S, et al. Gut-mammary-sIgA axis and immunologic mechanisms induced by porcine epidemic diarrhea virus [J]. Microbiology China, 2021, 48(11): 4240−4249.(in Chinese) DOI: 10.13344/j.microbiol.china.210041
[11] 孙立旦, 陈立功, 周双海, 等. 猪流行性腹泻病毒血清IgA抗体间接ELISA方法的建立 [J]. 中国兽医学报, 2019, 39(3):381−386. DOI: 10.16303/j.cnki.1005-4545.2019.03.01 SUN L D, CHEN L G, ZHOU S H, et al. Indirect ELISA for detecting specific IgA antibody in serum to porcine epidemic diarrhea virus [J]. Chinese Journal of Veterinary Science, 2019, 39(3): 381−386.(in Chinese) DOI: 10.16303/j.cnki.1005-4545.2019.03.01
[12] 谷长维, 谷长乐, 胡博. 猪流行性腹泻病毒血清中特异性IgG抗体间接ELISA检测方法的建立 [J]. 中国动物传染病学报, 2022, 30(3):106−112. DOI: 10.3969/j.issn.1674-6422.2022.3.zgsyjscb202203015 GU C W, GU C L, HU B. Development of indirect ELISA assay for detection of porcine epidemic diarrhea virus antibodies in sera of naturally infected pigs [J]. Chinese Journal of Animal Infectious Diseases, 2022, 30(3): 106−112.(in Chinese) DOI: 10.3969/j.issn.1674-6422.2022.3.zgsyjscb202203015
[13] DIEL D G, LAWSON S, OKDA F, et al. Porcine epidemic diarrhea virus: An overview of current virological and serological diagnostic methods [J]. Virus Research, 2016, 226: 60−70. DOI: 10.1016/j.virusres.2016.05.013
[14] 高冉, 刘建波, 陈建飞, 等. 抗猪SIgA SC片段单克隆抗体的制备及鉴定 [J]. 中国预防兽医学报, 2018, 40(5):455−458. GAO R, LIU J B, CHEN J F, et al. Preparation and identification of monoclonal antibodies against porcine SIgA SC fragment [J]. Chinese Journal of Preventive Veterinary Medicine, 2018, 40(5): 455−458.(in Chinese)
[15] 丛广义, 刘建波, 时洪艳, 等. 猪流行性腹泻病毒特异性SIgA抗体间接ELISA检测方法的建立及应用 [J]. 中国预防兽医学报, 2019, 41(8):812−817. CONG G Y, LIU J B, SHI H Y, et al. Establishment and application of indirect ELISA for detection of porcine epidemic diarrhea virus specific SIgA [J]. Chinese Journal of Preventive Veterinary Medicine, 2019, 41(8): 812−817.(in Chinese)
[16] 李春喜, 姜丽娜, 邵云. 生物统计学[M]. 5版. 北京: 科学出版社, 2016. [17] SUN R Q, CAI R J, CHEN Y Q, et al. Outbreak of porcine epidemic diarrhea in suckling piglets, China [J]. Emerging Infectious Diseases, 2012, 18(1): 161−163. DOI: 10.3201/eid1801.111259
[18] 王永恒, 周孟云, 李昱辰, 等. 猪流行性腹泻病毒黏液SIgA抗体ELISA检测方法的建立 [J]. 中国兽医科学, 2020, 50(1):10−19. DOI: 10.16656/j.issn.1673-4696.2020.0004 WANG Y H, ZHOU M Y, LI Y C, et al. Establishment of ELISA for detecting porcine epidemic diarrhea virus specific SIgA in mucus [J]. Chinese Veterinary Science, 2020, 50(1): 10−19.(in Chinese) DOI: 10.16656/j.issn.1673-4696.2020.0004
[19] 朱利塞, 李菲, 王娟, 等. 母猪初乳中猪流行性腹泻病毒SIgA抗体ELISA检测方法的建立及应用 [J]. 中国动物传染病学报, 2023, 31(1):55−60. DOI: 10.19958/j.cnki.cn31-2031/s.20211103.006 ZHU L S, LI F, WANG J, et al. Development of an ELISA method for detection of SIgA antibodies against porcine epidemic diarrhea virus in sow colostrum [J]. Chinese Journal of Animal Infectious Diseases, 2023, 31(1): 55−60.(in Chinese) DOI: 10.19958/j.cnki.cn31-2031/s.20211103.006
[20] 丁振江, 库旭钢, 闫贵伟, 等. 猪流行性腹泻IgA抗体间接ELISA检测方法的建立 [J]. 畜牧与兽医, 2014, 46(11):1−5. DING Z J, KU X G, YAN G W, et al. Development of an indirect ELISA method for detection of IgA antibody against porcine epidemic diarrhea virus [J]. Animal Husbandry & Veterinary Medicine, 2014, 46(11): 1−5.(in Chinese)
[21] 王治平, 张云静, 孙玉洁, 等. 基于单克隆抗体捕获的猪血清和乳汁中PEDVIgA抗体间接ELISA方法的建立和评估 [J]. 中国畜牧兽医, 2020, 47(7):2248−2255. WANG Z P, ZHANG Y J, SUN Y J, et al. Development and evaluation of an indirect ELISA for detection of porcine epidemic diarrhea virus IgA in pig serum and milk captured by monoclonal antibody [J]. China Animal Husbandry & Veterinary Medicine, 2020, 47(7): 2248−2255.(in Chinese)
[22] 柏家果, 杜琛, 王若木, 等. 猪流行性腹泻病毒特异性SIgA 间接ELISA检测方法的建立及初步应用 [J]. 中国兽医学报, 2023, 43(3):455−461. BAI J G, DU C, WANG R M, et al. Establishment and preliminary application ofindirect ELISA method for detection of specific SIgA against porcine epidemic diarrhea virus [J]. Chinese Journal of Veterinary Science, 2023, 43(3): 455−461.(in Chinese)
-
期刊类型引用(0)
其他类型引用(1)
计量
- 文章访问数: 478
- HTML全文浏览量: 220
- PDF下载量: 14
- 被引次数: 1
