本发明属于兽用生物制品技术领域，具体涉及一种猪病毒性腹泻病原的多重液相芯片检测方法及其构建。
背景技术
由病毒引起的猪腹泻疫病流行范围广，仔猪发病率高，死亡率可高达100％。而且多病原感染、协同感染较为普遍。临床很难确定主要诱发病因，给诊断和控制带来困难。
引起猪腹泻的主要病毒有猪流行性腹泻病毒(porcineepidemicdiarrheavirus，pedv)、猪德尔塔冠状病毒(porcinedeltacoronavirus，pdcov)、猪传染性胃肠炎病毒(porcinetransmissiblegastroenteritisvirus，tgev)、猪捷申病毒(porcineteschvirus，ptv)、牛病毒毒性腹泻病毒(bovineviraldiarrheavirus，bvdv)、猪轮状病毒(porcinerotavirus，porv)、猪萨佩罗病毒(porcinesapelovirus，psv)、猪嵴病毒(porcinekobuvirus，pkv)、猪星状病毒(porcineastrovirus，pastv)、猪环曲病毒(porcinetorovirus，ptov)、猪札幌病毒(porcinesapovirus，psav)等，这些病毒引起不成程度腹泻，在多个病原共同感染情况下，可导致严重腹泻及猪只死亡。
目前，关于猪病毒性腹泻病原的检测研究最多的是多重pcr方法，但大多是2重或3重pcr方法。焦洋等(猪传染性胃肠炎病毒、猪流行性腹泻病毒和猪博卡病毒多重pcr检测方法的建立，动物医学进展，2013，34(08)：71-75)建立了针对pedv、tgev和pbov的3重pcr检测方法。
陈小芬等(pedv/tgev/pdcov快速检测方法的建立及pedvcoe基因和orf3基因遗传变异研究，华南农业大学，2016.)针对pedv、tgev和pdcov，建立了三重rt-pcr检测方法。
刘高鹏(多重pcr检测和鉴定五种猪腹泻病毒研究，浙江理工大学，2017)建立了针对pedv、tgev、猪轮状病毒a型(porcinegrouparotaviruses,gar)、猪轮状病毒c型(porcinegroupcrotaviruses，gcr)、猪圆环病毒2型(porcinecircovirus2，pcv2)的5重实时荧光定量pcr方法，检测172份临床样品，发现gcr、gar、tgev、pedv和pcv2阳性率分别为5.81％、2.32％、8.72％、20.92％和46.51％，总混合感染率为22.09％。
但是，这些现有技术中，常规单重rt-pcr能够检测低至0.1pg的病毒rna，而多重rt-pcr随着检测病原数量的增加，灵敏度下降十分明显。鉴于这种灵敏度的局限性，其检测病原范围也很有限。
近年来，luminex检测得到了广泛应用，根据各病原序列，设计合适的检测引物，并根据检测引物扩增的目标序列，设计特异性探针；进一步在微球上偶联探针序列，通过探针与目标序列之间的杂交，最后在检测mfi值。该方法中，探针序列需要自行设计，随着病原数量的增加，探针设计难度也增加，导致探针特异性不强；而且微球偶联过程中，又增加了人为的不确定因素，导致试验误差增加。
技术实现要素：
本发明的目的在于提供一种猪病毒性腹泻病原的多重液相芯片检测方法及其构建，能够对猪的1-11种腹泻病毒进行定性、定量检测，特异性好，灵敏度高，结果准确，为临床猪腹泻样品的大规模筛查提供有效方法。
为了达到上述目的，本发明提供以下技术方案：
一种猪病毒性腹泻病原的多重液相芯片检测方法，包括以下步骤：
1)根据猪病毒性腹泻病原的基因序列，设计扩增猪病毒性腹泻病原的pcr引物，所述pcr引物包括上游引物和下游引物；
2)根据猪病毒性腹泻病原的基因序列，选择合适的tag序列，将tag序列偶联于对应的微球上；将tag互补序列标记于pcr上游引物5’端，tag互补序列和pcr上游引物之间有12个碳原子的臂距；在pcr下游引物5’端标记生物素，获得修饰引物；
3)提取样品rna，反转录获得cdna；
4)以步骤3)的cdna为模板，利用修饰引物进行多重pcr扩增，退火温度50～58℃，获得多重pcr扩增产物；
5)步骤4)获得的扩增产物与偶联在微球上的tag序列进行杂交，40～45℃反应20～30min，设阴性对照，检测杂交反应液，分别获得样品及阴性对照的mfi值，根据mfi值，进而确定cmfi值；当cmfi≥3时，结果判定为阳性，否则为阴性；
其中，cmfi＝[mfi(样品)-mfi(阴性对照)]/mfi(阴性对照)。
进一步，还包括定量检测步骤：构建猪病毒性腹泻病毒的阳性质粒标准品，制作猪病毒性腹泻病毒的mfi值与浓度的标准曲线，利用该标准曲线和步骤5)测得的样品mfi值，对样品进行定量检测。
优选地，所述的猪病毒性腹泻病原中有猪的1-11种病毒性腹泻病原。
进一步，所述猪的1-11种病毒性腹泻病原选自ptv、pkov、pdcov、psv、psav、pastv、ptov、porv、pedv、bvdv和tgev。
进一步，所述猪的1-11种病毒性腹泻病原的pcr引物序列如下，从5’到3’：
ptov上游引物：aattgcttattggtggcttc；
下游引物：agcdatttgrgcdgcattc；
pdcov上游引物：cgcagttttcattgtgtcca；
下游引物：cctgtggcggatttctaact；
pastv上游引物：ccttwccccactgatgaaga；
下游引物：cctgtccatctgcctttctgt；
pkov上游引物：cgccgttcactctttgtcc；
下游引物：accaagcagcatcctaccag；
psv上游引物：ctggactgggcctatact；
下游引物：tyaggtacacacgggctc；
ptv上游引物：tgcagctctttcacatttrga；
下游引物：acttgtatgaggcccatcg；
psav上游引物：gtggcaacgtacaaygcrtgg；
下游引物：gcctccatcacgaacact；
bvdv上游引物：agccatgccchtagtaggac；
下游引物：cctctgcwrcaccctatc；
pedv上游引物：gtcaagatggctattctatgg；
下游引物：accaagaatgtgtcctgcg；
tgev上游引物：caaaggaataggtaacaggg；
下游引物：atctgcatgaggtccagt；
porv上游引物：ggattacttggyactacttt；
下游引物：aatgatgctgaytatggaagt。
进一步，所述猪的1-11种病毒性腹泻病原的tag序列如下，从5’到3’：
ptov-tag：ttgaaagagagagataagaataaa；
pdcov-tag：gatagtaaatagagaaagagttaa；
pastv-tag：gttaaatgtaaagtgaaagaatag；
pkov-tag：agtagtgaaagaaatgaaatgtaa；
psv-tag：ttaaagaagagaaagaaagtgtta；
ptv-tag：tgtgagtaaattgtgataaagtaa；
psav-tag：tgtttatagattgatgatagtgtt；
bvdv-tag：tatgaaatgtttgtttattgtgtg；
pedv-tag：gaaagagtatgaaagttgattaaa；
tgev-tag：aattgttgaatatgtttgtgtttg；
porv-tag：gtattagttaaagttgaaagatga。
又，所述猪的1-11种病毒性腹泻病原的tag互补序列如下，从5’到3’：
ptov-tag：aactttctctctctattcttattt；
pdcov-tag：ctatcatttatctctttctcaatt；
pastv-tag：caatttacatttcactttcttatc；
pkov-tag：tcatcactttctttactttacatt；
psv-tag：aatttcttctctttctttcacaat；
ptv-tag：acactcatttaacactatttcatt；
psav-tag：acaaatatctaactactatcacaa；
bvdv-tag：atactttacaaacaaataacacac；
pedv-tag：ctttctcatactttcaactaattt；
tgev-tag：ttaacaacttatacaaacacaaac；
porv-tag：cataatcaatttcaactttctact。
本发明提供所述多重猪病毒性腹泻病原的液相芯片检测方法的构建方法，包括以下步骤：
a)提取样品rna，反转录获得cdna；
b)根据猪病毒性腹泻病原的基因序列，设计扩增病毒性腹泻病原的pcr引物，该pcr引物包括上游引物和下游引物；
c)以cdna为模板，利用扩增猪病毒性腹泻病原的pcr引物进行pcr扩增，回收pcr扩增产物，构建猪病毒性腹泻病毒阳性质粒；
d)根据猪病毒性腹泻病原的基因序列，选择合适的tag序列，将tag序列偶联于微球；将tag互补序列标记于pcr上游引物5’端，tag互补序列和pcr上游引物之间有12个碳原子的臂距，并在pcr下游引物的5’端标记生物素，得到修饰引物；
e)多重pcr扩增
以步骤c)构建的猪病毒性腹泻病毒阳性质粒为模板，利用修饰引物，进行多重pcr扩增，获得多重pcr扩增产物，确定合适的引物浓度和退火温度；
f)杂交反应
步骤e)获得的扩增产物与偶联在微球上的tag序列进行杂交，设阴性对照，在37～45℃反应20～45min，检测杂交反应液，获得mfi值，根据mfi值的高低，确定合适的杂交温度和时间；根据mfi值确定cmfi值，判断样品的阴/阳性。
进一步，步骤e)中，确定的引物浓度和退火温度分别为：修饰引物的上下游引物浓度各为2.5μm，退火温度55℃；步骤f)中，确定的杂交温度和时间分别为：杂交温度40～45℃，杂交时间20～30min。
本发明利用luminex液相芯片xtag技术，实现对1-11种猪腹泻病毒的进行定性、定量检测，其以荧光编码的微球作为反应载体，根据微球内荧光染料配比的不同，产生出不同的编码微球，通过把luminex液相芯片特有的tag短序列结合到微球上，每个编码微球都可对应不同的目标物，从而实现多分子同时检测。
本发明中，首先根据猪腹泻病原的基因序列，分别筛选出合适的tag序列并偶联于各微球；其次，将tag互补序列标记于扩增各病原的pcr上游引物5’端，并在pcr下游引物5’端标记生物素；然后利用修饰引物进行多重pcr扩增，将pcr扩增产物与偶联在微球上的tag序列进行杂交；最后，在luminex200机器上进行检测；根据mfi值，确定各孔样品的阴/阳性，实现对猪病毒性腹泻多重病原进行定性定量检测。
本发明的检测方法中，杂交温度不仅对杂交效率有直接影响，而且对tag与anti-tag的特异性杂交也存在显著影响；为了进一步提高杂交效率，以及减少非特异性的结合，本发明对tag与anti-tag的杂交时间进行了筛选。
与现有技术相比，本发明具有如下有益效果：
本发明中，针对猪腹泻病原的基因序列，遵循与特定病原基因序列完全互补，而与其他病原基因序列无非特异性反应的原则，筛选出合适的tag序列，然后与偶联该tag序列的商品化微球进行偶联，将tag互补序列和生物素制作修饰引物，对待测样品进行扩增，扩增产物与偶联在微球上的tag序列进行杂交，检测mfi值，操作简单，无需进行探针序列的设计，既增加了特异性，又避免了不必要的误差。
本发明中，具有双重结合特性，偶联微球的tag序列和pcr扩增产物识别并结合之后，才能进一步被链霉亲和素所识别，最后发出特异荧光信号，因此，该方法的特异性强于其他检测方法；并且，本发明可同时使用多达上百种微球，最多可以同时检测上百种病原，还具有高通量性。
本发明中，微球可以通过tag序列多方位捕捉pcr扩增片段，再结合生物素-亲和素级联放大反应系统，大大提高了该方法的检测限量，比普通pcr灵敏度至少提高10倍。
附图说明
图1为本发明实施例1中猪的11种病毒性腹泻病原pcr扩增的目的片段。
图2为本发明实施例1中7重pcr引物在三种浓度下的扩增结果。
图3为本发明实施例1中4重pcr引物在三种浓度下的扩增结果。
图4为本发明实施例1中7重pcr引物在三个温度下的扩增结果。
图5为本发明实施例1中4重pcr引物在三个温度下的扩增结果。
图6-8为本发明实施例1中11重pcr产物的杂交中在三个温度及时间下的扩增结果。
图9-10为本发明实施例1中猪病毒性腹泻病原的luminex灵敏度检测结果。
图11为本发明实施例1中猪的11种病毒性腹泻病原的特异性检测结果。
具体实施方式
以下结合具体实施例对本发明作进一步说明。
magplex-tag微球、校验试剂盒和鞘液购luminex公司；链霉亲和素-藻红蛋白购自默克史密博公司；dnamarkerdl500、pmd18-t载体试剂盒和trizol购自takara公司；质粒提取试剂盒、dna胶回收试剂盒购自axygen公司；depc水购自sigma公司；goldview染料购自赛百胜公司；蛋白胨、酵母提取物购自oxoidltd公司；琼脂粉购自生工生物工程(上海)股份有限公司。
含有pedv、tgev、pdcov、porv、bvdv、psv、pkv、pastv、ptov、ptv以及psav病原的阳性样品由上海市农业科学院实验室保存。e.coli.dh5α感受态菌株由上海市农业科学院实验室制备保存.
实施例1猪病毒性腹泻病原的多重液相芯片检测方法的构建，包括以下步骤：
1.引物设计与修饰
首先，用oligo7软件分析猪的11种病毒性腹泻病原，分别为pedv、tgev、pdcov、porv、bvdv、psv、pkv、pastv、ptov、ptv以及psav，设计pcr引物，筛选合适的tag序列；
猪的11种病毒性腹泻病原的pcr引物的序列分别如下：
ptov上游引物：aattgcttattggtggcttc；
下游引物：agcdatttgrgcdgcattc；
pdcov上游引物：cgcagttttcattgtgtcca；
下游引物：cctgtggcggatttctaact；
pastv上游引物：ccttwccccactgatgaaga；
下游引物：cctgtccatctgcctttctgt；
pkov上游引物：cgccgttcactctttgtcc；
下游引物：accaagcagcatcctaccag；
psv上游引物：ctggactgggcctatact；
下游引物：tyaggtacacacgggctc；
ptv上游引物：tgcagctctttcacatttrga；
下游引物：acttgtatgaggcccatcg；
psav上游引物：gtggcaacgtacaaygcrtgg；
下游引物：gcctccatcacgaacact；
bvdv上游引物：agccatgccchtagtaggac；
下游引物：cctctgcwrcaccctatc；
pedv上游引物：gtcaagatggctattctatgg；
下游引物：accaagaatgtgtcctgcg；
tgev上游引物：caaaggaataggtaacaggg；
下游引物：atctgcatgaggtccagt；
porv上游引物：ggattacttggyactacttt；
下游引物：aatgatgctgaytatggaagt。
利用猪的11种腹泻病原的检测引物，分别进行单一pcr扩增，获得各阳性目的片段，参见图1，其中，m：dl500dna相对分子质量标准；1：psav；2：阴性对照；3：ptv；4：阴性对照；5：psv；6：阴性对照：7：pkv；8：阴性对照；9：pastv；10：阴性对照；11：pdcov；12：阴性对照；13：ptov；14：阴性对照；15：tgev；16：阴性对照；17：pedv；18：阴性对照；19：bvdv；20：阴性对照；21：porv；22：阴性对照；图1表明，各片段pcr扩增特异，大小均正确。
所述猪的11种病毒性腹泻病原的tag序列如下，从5’到3’：
ptov-tag：ttgaaagagagagataagaataaa；
pdcov-tag：gatagtaaatagagaaagagttaa；
pastv-tag：gttaaatgtaaagtgaaagaatag；
pkov-tag：agtagtgaaagaaatgaaatgtaa；
psv-tag：ttaaagaagagaaagaaagtgtta；
ptv-tag：tgtgagtaaattgtgataaagtaa；
psav-tag：tgtttatagattgatgatagtgtt；
bvdv-tag：tatgaaatgtttgtttattgtgtg；
pedv-tag：gaaagagtatgaaagttgattaaa；
tgev-tag：aattgttgaatatgtttgtgtttg；
porv-tag：gtattagttaaagttgaaagatga。
其次，将各病原的pcr引物进行修饰，在各上游引物5’端以tag互补序列进行修饰，两者之间增加了12个碳原子的臂距；下游引物5’端进行生物素修饰，11种腹泻病毒的修饰引物扩增长度以及所对应的微球编号见表1，引物的修饰由生工生物工程(上海)股份有限公司完成。
2.提取rna、反转录
吸取250μl猪粪便样品，加入750μltrizol裂解液，充分震荡混匀，静置10min；加入200μl三氯甲烷，充分震荡，4℃离心15min；小心吸取500μl上清液，加入等体积异丙醇，-20℃沉淀10min，然后4℃离心15min；弃上清液，加入1ml75％乙醇洗涤，4℃离心15min；弃去上清，倒扣纸上，风干，用20μl无rnasedepc水重悬，-80℃保存备用。
以提取的rna为模板，利用随机引物进行反转录，获得cdna，反应程序：42℃60min；75℃15min；4℃，∞。反转录产物-20℃保存，反转录反应体系参见表2。
表2
3.构建病毒的阳性质粒
以上述cdna为模板，利用猪的11种病毒性腹泻病原的pcr引物分别进行pcr扩增，具体反应体系为：2μlcdna，pcrmix12.5μl，上、下游引物(10μm)各0.5μl，用水补至25μl；反应程序：95℃，4min；(95℃，30s；55℃，30s；72℃，35s)，35个循环；72℃，7min；4℃，∞。
利用试剂盒回收pcr扩增产物，将回收的目的片段克隆入pmd18-t载体，并转化dh5α感受态，挑取单个菌落进行菌液pcr鉴定，将疑似阳性质粒送公司测序；在ncbi数据库内进行序列比对，将测序正确的阳性质粒于-20℃保存备用。
4.多重pcr反应体系的建立及条件筛选
将获得的猪的11种病毒阳性质粒进行稀释，以稀释好的猪的11种病毒阳性质粒为模板，用表1的修饰引物进行多重pcr；
反应体系为：各取2μl质粒、2.5μl10×pcrbuffer、2μldntp、0.5μlrtaq、11对病毒上、下游引物(10μm)各0.5μl，用水补至25μl。
反应程序为：95℃，4min；(95℃，30s；55℃，30s；72℃，30s)，35个循环；72℃，7min；4℃，∞。
进一步，设定不同引物浓度(2.5μm、5μm、10μm)和退火温度(50℃、55℃、58℃)进行pcr，根据mfi值，筛选引物浓度，结果参见图2-3，结果表明，3个引物浓度(2.5μm、5μm、10μm)均能扩增出特异性条带。
又，采用3个温度梯度，分别为50℃、55℃、58℃，进行多重pcr反应条件的优化，结果如图4-5，由图4-5可见，这3个退火温度均能扩增出目的条带，55℃为最佳退火温度。
5.多重pcr产物的杂交反应及条件筛选
将上述多重pcr反应的扩增产物与相对应的微球进行杂交，设立阴性对照，步骤如下：
(1)将配制好的微球工作液恢复至室温，加入八连管中，20μl/孔(每种微球个数为2500个)；
(2)在上述孔内加入5μlpcr扩增产物；
(3)用1×tmhybridization稀释sape，使其终浓度为10μg/ml；每孔加入75μlsape报告缓冲液，轻轻混匀；
(4)放置pcr仪上进行杂交反应，37～45℃反应20～45min；
(5)取100μl反应液在luminex200仪器上进行检测。
进一步，筛选杂交温度及时间，参照luminex推荐杂交温度37～45℃，选取40℃、42℃、45℃三个杂交温度，根据所得到的mfi值，确定最佳杂交温度；参照luminex推荐杂交时间20～45min，选取20min、25min、30min三个时间点进行杂交，根据所得到的mfi值，确定最佳杂交时间为30min。
结果参见图6-8，可见，42℃杂交30min可获得的mfi值最高，因此，判定luminex检测的最佳杂交条件为42℃杂交30min。
6.luminex检测
luminex200液相系统预热30min，采用xpoint3.1软件，对每孔样品进行读数。当读取的微球数≥100，说明读值有效；记录每孔的mfi值(平均荧光强度)。
选取48个阴性对照样品，分别进行luminex检测，获得各自的mfi值，从而来确定cmfi值。运用“信噪比”：cmfi＝[mfi(样品)-mfi(阴性对照)]/mfi(阴性对照)。当cut-offcmfi≥3时，结果判定为阳性，否则为阴性。
利用depc水将各病原的阳性质粒进行10倍倍比稀释，然后采用筛选的多重pcr反应条件，分别进行7重pcr和4重pcr反应进行扩增进行扩增，检测该方法的灵敏性，结果参见图9-10。
由图9-10可见，当质粒稀释至10-6倍时，luminex200仪器检测到的mfi值低于阳性判定值，可判定这两种多重pcr方法可检测ptv、pkov、pdcov、psv、psav、pastv、ptov、porv、pedv、bvdv、tgev的质粒最低检测限量分别为3.12×103copies/μl、2.92×103copies/μl、2.79×103copies/μl、3.37×103copies/μl、2.7×103copies/μl、3.02×103copies/μl、2.65×103copies/μl、2.57×103copies/μl、1.74×103copies/μl、2.41×103copies/μl、2.75×103copies/μl。
以单一病毒的cdna作为模板，分别进行7重pcr和4重pcr反应，进行扩增，检测该方法的特异性，结果如图11，可见，各反应中，只有相对应的病原能扩增出目的条带，其余病原扩增均为阴性，说明该方法的特异性良好。
7.绘制定量检测标准曲线
将11种猪病毒性腹泻病毒的标准阳性质粒进行10倍倍比稀释，稀释成100-10-7这8个浓度梯度。利用表1的修饰引物，将每个稀释度的混合质粒标准品分别进行pcr扩增反应，体系中同时加入11种病毒的修饰引物。将各pcr扩增产物与相对应的微球进行杂交，然后在luminex机器上进行检测，系统会自动生成11条标准曲线，11种病原的标准曲线的r2均大于0.99，说明所有标准曲线的线性回归分析均具有高的相关系数。
实施例2临床猪病毒性腹泻样品的的液相芯片检测与分析
利用实施例1中的液相芯片检测方法，对173份猪腹泻粪便样品检测，检测设阴性对照，以复孔方式进行，获得平行实验的两组数据，取两组数据的平均值，对每个样品进行数据分析，以其中一个样品为例，检测获得的两组mfi值检测值及阴性对照值参见表3。
1.定性分析过程如下：
12#xtag微球：样品cmfi＝[mfi(样品)-mfi(阴性对照)]/mfi(阴性对照)＝0.99＜3，为porv阴性；
14#xtag微球：样品cmfi＝[mfi(样品)-mfi(阴性对照)]/mfi(阴性对照)＝0.15＜3，为psv阴性；
19#xtag微球：样品cmfi＝[mfi(样品)-mfi(阴性对照)]/mfi(阴性对照)＝8.56＞3，为bvdv阳性；
20#xtag微球：样品cmfi＝[mfi(样品)-mfi(阴性对照)]/mfi(阴性对照)＝0.64＜3，为pedv阴性；
39#xtag微球：样品cmfi＝[mfi(样品)-mfi(阴性对照)]/mfi(阴性对照)＝2.25＜3，为psav阴性；
43#xtag微球：样品cmfi＝[mfi(样品)-mfi(阴性对照)]/mfi(阴性对照)＝6.67＞3，为ptov阳性；
44#xtag微球：样品cmfi＝[mfi(样品)-mfi(阴性对照)]/mfi(阴性对照)＝7.58＞3，为pkov阳性；
51#xtag微球：样品cmfi＝[mfi(样品)-mfi(阴性对照)]/mfi(阴性对照)＝0.52＜3，为pastv阴性；
53#xtag微球：样品cmfi＝[mfi(样品)-mfi(阴性对照)]/mfi(阴性对照)＝0.48＜3，为tgev阴性；
72#xtag微球：样品cmfi＝[mfi(样品)-mfi(阴性对照)]/mfi(阴性对照)＝0.89＜3，为pdcov阴性；
74#xtag微球：样品cmfi＝[mfi(样品)-mfi(阴性对照)]/mfi(阴性对照)＝5.56＞3，为ptv阳性。
可见，该检测样品中同时感染病毒bvdv、ptov、pkov和ptv四种病毒。
173份样品的检测分析结果参见表4-5。
由表4可见，porv的阳性检出率为14.5％、psv的阳性检出率为6.9％、bvdv的阳性检出率为6.4％、pedv的阳性检出率为32.4％、pkov的阳性检出率为26.0％、psav的阳性检出率为11.6％、ptv的阳性检出率为1.7％、pastv的阳性检出率为26.0％、tgev的阳性检出率为3.5％、pdcov的阳性检出率为0.6％、ptv的阳性检出率为1.7％。
pedv检出率最高，其次是pkov和pastv，接着是porv和psav。tgev、pdcov和ptv检出率较低。
由表5可见，在173份临床样品中有16.8％(29/173)的样品中存在2种病毒感染，主要是pkov和pastv；10.9％(19/173)的样品中存在3种病毒感染，主要是pkov、pastv和ptov；3.4％(6/173)的样品中存在4种病毒感染，主要是pkov、pastv、ptov与pedv；1.7％(3/173)的样品中存在5种病毒感染，主要是pkov、pastv、psav、ptov与ptv混合感染，3份样品中存在pkov、pedv、porv、pastv、psav、ptov与psv混合感染。
2.定量分析
根据测得的mfi值，以定性判断标准为参照，结合各病毒的标准曲线，进一步进行定量判定，部分样品的定量检测结果参见表6。
由表6可知，porv、psv、bvdv、pedv、pkov、psav、pastv、tgev、pdcov和ptv的病毒含量分别≥100.29pg/ml、2.23pg/ml、156.53pg/ml、10.97pg/ml、34.69pg/ml、24.73pg/ml、367.62pg/ml、54.66pg/ml、10.64pg/ml、14.68pg/ml时，可判定为阳性。
进一步，可利用聚类分析原理，将各样品中的病原含量以颜色表示出来，可以很直观的估算各样品中病原含量。