Molecular detection of PRRSV and variation analysis of ORF5 gene in Fujian province in 2017
-
摘要:目的
研究2017年福建地区中小规模猪场猪繁殖与呼吸综合征病毒(PRRSV)毒株的分子流行情况。
方法从福建地区中小规模猪场收集83份疑似蓝耳阳性病料，对编码GP5蛋白的ORF5基因进行了RT-PCR检测和遗传进化分析。
结果采集的83份疑似蓝耳病料中阳性率为73%，测得的序列之间ORF5的核苷酸相似性达80.3%~100.0%，其中所测的各个分支毒株与其所属分支的代表毒株NADC30、GM2和JXA1的核苷酸相似性分别为92.6%~93.9%、89.8%~93.0%和93.4%~99.3%，检测到新分支中新出现的毒株与JXA1的核苷酸相似性为91.0~93.0%。GP5氨基酸比对结果显示不同亚群GP5信号肽区、诱导表位、主要中和表位及高变区等功能区发生了显著变异。
结论福建地区出现了PRRSV新分支，使福建PRRSV变异种类更加多样化，建议加强PRRSV流行及变异的监测。
-
关键词:
- 猪繁殖与呼吸道综合征 /
- 蓝耳病 /
- 分子流行病学 /
- ORF5基因 /
- 遗传变异
Abstract:ObjectiveTo study the molecular prevalence of the porcine peproductive and respiratory syndrome virus (PRRSV) strain in small and medium-sized pig farms of Fujian province in 2017.
MethodWe collected 83 suspected blue ear positive disease materials from small and medium-sized pig farms in Fujian province. RT-PCR detection and genetic evolution analysis of ORF5 gene which encoded GP5 protein was performed.
ResultThe positive rate of 83 suspected blue ear disease samples was 73%, and the nucleotide similarity among the measured ORF5sequences was 80.3%−100.0%. The nucleotide similarities between tested strains of different branches and the representative strains of those branches, NADC30, GM2 and JXA1, were 92.6%−93.9%, 89.8%−93.0% and 93.4%−99.3%, respectively. The new identified strain in the new branch had nucleotide similarity of 91.0%−93.0% to JXA1. The GP5 amino acid alignment showed significant variation in functional regions such as GP5 signal peptide region, induction epitope, principal neutralizing epitope and hypervariable region in different subgroups.
ConclusionA new branch of PRRSV has emerged in Fujian area, making the PRRSV variants in Fujian more diverse. It is necessary to strengthen the monitoring of PRRSV epidemics and mutations.
-
在现代畜牧业中，精确测量动物体尺对评估其生长和健康状况至关重要。传统测量方法主要靠人工，不仅耗时，还因测量者技能和经验不同而导致误差。计算机技术的进步使得数字化畜牧业成为推进行业发展的关键[1]，其中视觉技术得到了广泛应用。众多研究利用计算机视觉技术来检测动物行为[2-3]、计数牲畜[4-6]以及评估肉质[7-8]，展现了其在畜牧业中的潜力，为动物体尺的自动测量提供了重要支持。
传统视觉技术测量动物体尺的主要方法包括基于三维重建和基于关键点提取的技术。三维重建技术通过从多个角度捕获动物的图像，以克服动物姿态多样性的问题[9-11]。然而，在实际的猪场环境中，这种方法面临多种挑战。例如，猪的活动可能导致重建的模型扭曲；表面纹理和光照的变化可能影响重建的精度；选择和配置硬件设备，以及满足大量的计算资源需求，也是这种方法所面临的挑战。另一方面，基于关键点提取的技术需要图像尺度在旋转和光照变化时保持稳定[12-13]，其准确性主要依赖于关键点提取模型的鲁棒性；其次，猪的多种姿势(如站立、行走或躺卧)会影响图像中提取到的关键点信息的分析；此外，图像的畸变性也可能进一步造成测量结果不准确。无论是基于三维重建的技术，还是基于关键点提取的技术，都需要进一步的研究和改进，以提高在实际环境中的应用效果。
随着农业智能化迅速发展，YOLO系列的视觉算法得到了广泛应用[14-16]。殷献博等[17]提出了采用多注意力机制改进的YOLOX-Nano模型，在柑橘梢期监测任务中达到了88.07%的平均精度，实现了自动化监测；另外，杨秋妹等[18]应用YOLO技术的高精度特点，开发了一种自动化猪只盘点系统，即使在复杂环境下也能准确地进行盘点，减少了对人工操作的依赖。此外，ByteTrack作为一种多目标跟踪算法[19]，与YOLO结合在农业领域广泛应用。Zheng等[20]利用Yolov7和ByteTrack实现了自然场景下奶牛的自动监控，平均分析速度达到47帧/s；同时，Qian等[21]提出了基于Yolov5的鱼检测算法，结合ByteTrack算法实时跟踪，解决了遮挡问题。在这些研究的启发下，本研究提出了一种基于Yolov8-Pose和ByteTrack追踪算法的猪只体尺估计方法，旨在实现实际环境中猪场猪只体尺的自动测量。
1. 图像采集与数据集制作
1.1 图像来源
本研究采用的猪只图像数据源自双河基地实验猪场，通过场内安装的室内监控系统进行采集。监控系统采用的是海康威视公司生产的型号为DS-2CD3146FWD-IS的半球形摄像头，该摄像头以顶置方式安装在猪栏上方，距离地面高度为2.0~2.5 m，以获得最佳的俯拍视角。在数据收集阶段，随机选取猪场中保育舍、后备舍以及营养调控舍的场景进行视频裁剪，共裁剪238段视频片段，平均每个片段为30 s，每间隔6帧进行抽取，共抽取7 324张原始图片，所有图像的分辨率均设定为2 560像素×1 440像素。经人工筛选后，选出1 822张合适的图片进行标注，共包含6 452头猪只，用于本研究的数据分析。
1.2 数据集制作
本研究采用了一种详细的关键点标注策略，为图像中的每头猪标注6个主要关键点(图1展示了4个场景猪只的关键点标注)：头部、右肩、右臀、尾巴、左臀和左肩，通过测量这些关键点间的距离，估算出猪的臀宽、体长和肩宽等重要的体尺数据。为了保证数据集的独立性和模型的泛化性，按照约8∶1∶1的比例划分训练集(含1 622张图像)、验证集(200张)和测试集(200张)，以支持模型训练和评估的需要。
![图 1 4个不同场景猪只6个关键点的标注]()
图 1 4个不同场景猪只6个关键点的标注Figure 1. Annotation of six key points of pigs in four different scenarios为了进一步提高模型的检测能力，在训练过程中引入了Mosaic、色域变换和左右翻转3种数据增强策略。其中，Mosaic增强是指抽取训练样本集中的4张图片进行随机裁剪并拼接形成一个新的训练样本，图2展示了经过Mosaic增强后的样本图片。左右翻转的概率设为0.5，色域变换HSV分别设为0.015、0.700和0.400。为避免过拟合，更加关注原始数据，训练的最后10个epochs中不再采用任何数据增强策略。
2. 方法与算法设计
2.1 研究方法
本研究针对猪自然静止与运动状态下姿态的多样性，提出了一种结合关键点检测算法和追踪算法的体尺估计方案。通过分析视频流中连续帧的图像，实现体尺估计并完成相应记录。引入感兴趣区域(Region of interest，ROI)减少相机畸变的影响并提高处理速度，采用基于帧关联性的异常检测算法优化预测结果，提高体尺估计的精度和稳定性。图3为本研究工作流程。
2.2 体尺算法设计
2.2.1 摄像机标定
在精确测量动物体尺的过程中，需将图像上的像素距离转换为实际世界中的距离。这种转换涉及计算一个缩放因子(实际距离与图像中测量到的像素距离的比值)，使得从图像数据中直接测量物理尺寸成为可能。为实现这一转换，本研究首先通过固定摄像头并使用标定板从多个角度和位置拍摄图像，来获得摄像机的内外参数。内参数包括焦距、主点坐标以及畸变系数，这些参数反映了摄像机如何将三维世界映射到二维图像平面。外参数则描述了摄像机相对于世界坐标系的位置和方向，为确定物体在现实世界中的位置提供必要信息。通过MATLAB标定工具箱处理标定图像，可以得到这些参数。随后，利用内参数对图像进行畸变校正，并使用外参数将摄像机坐标系下的点转换到世界坐标系下。图4展示了经过相机矫正后的图像，其中的坐标系对应真实世界中的坐标系，进而通过对比棋盘格的已知实际尺寸(每个正方形小格的边长为70 mm)与其在矫正图像中的像素尺寸，计算出缩放因子。
2.2.2 关键点检测
本研究选用2023年新推出的Yolov8-Pose作为关键点检测模型。该模型具有以下优势：通过端到端的训练框架直接从图像中检测目标框和关键点，速度更快，利用得到的目标框更易于进行追踪任务；深度卷积网络和多尺度特征提取能力使其能够有效地识别和定位猪的关键点，精度高，无需对算法网络进行调整或优化；集成了注意力机制和上下文信息融合，增强了模型在复杂背景下尤其是在遮挡和变化光照条件下检测关键点的能力；模型轻量化、易部署，适用于各种环境。
2.2.3 结合ByteTrack追踪算法的体尺估计
ByteTrack算法是一种基于目标检测的多目标跟踪技术，区别于依赖重识别(ReID)的方法，它仅使用目标检测步骤中获得的目标框信息来执行跟踪。相对于其他追踪算法，该算法保留了低置信度的检测框，采用卡尔曼滤波器来预测高低分目标框的未来状态，并运用匈牙利算法实现当前目标与已存在轨迹之间的匹配。由于其精度和速度的优势，本研究选择结合ByteTrack追踪算法和关键点检测方法，使用视频流作为输入，通过自动识别并追踪视频中的每一只猪，为其分配一个唯一ID，随后连续捕捉每只猪只的动态变化，并测量每一时刻的关键点信息，利用像素转化技术将关键点的像素距离转换为实际的体长、肩宽和臀宽，并将这些数据记录到相应的ID中。通过对视频流中猪只的连续观察和分析，能够更多地获取其体态信息，减少单次测量可能引入的误差，从而确保数据的稳定性和可靠性；同时，相对人工选取特定姿态照片的处理方式，减少了人工依赖和人工误差。以下是体尺计算公式，具体参考图5。
![图 5 猪只6个关键点和体尺测量示意图]()
图 5 猪只6个关键点和体尺测量示意图Figure 5. Schematic diagram of six key points and body measurements for pigs$$ {L_{{\mathrm{body}}}} = \sqrt {{{\left[ {{\theta _1}({x_2} - {x_3})} \right]}^2} + {{\left[ {{\theta _2}({y_2} - {y_3})} \right]}^2}} \text{，} $$ (1) $$ {L_{{\mathrm{shoulder}}}} = \sqrt {{{\left[ {{\theta _1}({x_2} - {x_6})} \right]}^2} + {{\left[ {{\theta _2}({y_2} - {y_6})} \right]}^2}} \text{，} $$ (2) $$ {L_{{\mathrm{hip}}}} = \sqrt {{{\left[ {{\theta _1}({x_3} - {x_5})} \right]}^2} + {{\left[ {{\theta _2}({y_3} - {y_5})} \right]}^2}} \text{，} $$ (3) 式中， Lbody、Lshoulder、Lhip分别为猪只的体长、肩宽、臀宽，(xi, yi)为第i个关键点的坐标，按顺序i=1~6分别对应头、右肩、右臀、尾、左臀、左肩；
$ {\theta _1} $ 、$ {\theta _2} $ 分别为横、纵坐标方向的缩放因子，以实现像素距离到实际尺寸的转换(缩放因子受摄像头配置参数和安装位置影响)。2.3 ROI筛选与姿态和异常检测过滤
2.3.1 ROI筛选
在实际生产中猪群倾向于聚集在摄像范围边缘区域，图片畸变性导致同一猪只在不同位置的体尺估计存在显著偏差，在远离中心位置时呈现缩小的趋势。同时，猪只的聚集和拥挤情况也会导致猪的姿态变形，进一步影响体尺估计的准确性。此外，猪只数量增加会显著降低处理速度。为提高估计准确度和处理速度，本研究设置ROI，仅对中央部分的猪只进行评估，既减少了边缘畸变和拥挤姿态变形的影响，也保证了足够的检测范围以便一段时间内能够统计到全部猪只。在本试验设置中，选定的ROI是图像中心的红色矩形框，其大小和形状可根据具体场景灵活调整，如图6所示。
2.3.2 姿态和异常检测过滤
在测量位于ROI区域内猪只体尺时，分2个阶段进行。第1阶段是猪只姿态的检测过滤，在这一阶段中需要考虑猪只是否处于站立状态以及站立时是否存在歪曲。当猪只处于站立状态时，肩宽与体长的比值以及臀宽与体长的比值均大于设定阈值(阈值的大小仅与猪舍场景有关，与猪只的实际大小无关，不同场景设定不同阈值)，并且头部关键点与双肩形成的三角形以及头部关键点与左右臀部关键点形成的三角形均为锐角三角形时，表示猪只姿态端正，测得数据有效，进入第2阶段，反之数据视为无效，丢弃。图7为需要被过滤掉的异常猪只姿态的示例。第2阶段是异常检测过滤，进一步排除由于运动模糊等因素可能导致的同一猪只在连续帧中体尺测量结果有较大差异的情况，本研究利用视频流中连续帧的关键点位置关联性，采用异常检测处理测量误差：对当前有效数据的前10帧内的数据进行分析比较，若当前测得数据与各数据最小差距小于4 cm，则保存到猪只字典中，反之则丢弃。经过一段时间运行，每头猪体尺数据为当前该猪只字典所有数据的平均值。
3. 模型训练与结果
3.1 模型训练
3.1.1 试验平台与模型参数设置
在Linux操作系统下，通过Python 3.8.1和PyTorch 1.8.0搭建的深度学习框架进行网络模型的训练与测试。试验配置包括Intel Xeon Silver 4110 CPU和2张NVIDIA GeForce RTX 2080 Ti GPU，利用CUDA 11.1并行计算框架以及CuDNN 11.7深度神经网络加速库进行性能优化。Yolov8-Pose模型训练采用640像素×640像素的输入图像分辨率，使用Adam优化器，设置训练批次大小为32，总共进行300次迭代训练。训练阶段，初始学习率设置为0.001，通过余弦退火策略进行学习率更新，针对单一目标类别，并关注6个关键点的检测。
3.1.2 评价指标
在执行关键点检测任务时，通常通过4个核心指标来评估模型的性能：精确度(Precision，P)、召回率(Recall，R)、平均精确度(Average precision，AP)和平均精确度均值(mean Average precision，mAP)。这些指标共同提供了一个全面的框架，量化模型在定位和识别关键点方面的效能，计算公式如下。
$$ P = \dfrac{{{\mathrm{TP}}}}{{{\mathrm{TP}} + {\mathrm{FN}}}} \times 100{\text{%}} \text{，} $$ (4) $$ R = \dfrac{{{\mathrm{TP}}}}{{{\mathrm{TP}} + {\mathrm{FN}}}} \times 100{\text{%}} \text{，} $$ (5) $$ {\mathrm{AP}} = \displaystyle\int\limits_0^1 {P(R){\mathrm{d}}R \times } 100{\text{%}} \text{，} $$ (6) $$ {\mathrm{mAP}} = \dfrac{1}{M}{\displaystyle\sum}_{k } {{\mathrm{AP}}_k \times } 100{\text{%}} \text{，} $$ (7) 式中，TP表示猪只位置被正确预测出来的数量，FP表示环境被错误预测为猪只位置的数量，FN为猪只位置未被预测出来的数量，M为检测目标的类别总数，APk为第k类AP。模型的性能通过P-R曲线(图8)进行可视化评估，曲线图标明了模型在精确度和召回率方面的表现，其下方的面积即为模型的平均精确度。图8A是目标框的P-R曲线，图8B是关键点的P-R曲线。最终，选出表现最优的模型权重，并在测试集上进行效果验证，图9展示了模型在验证集中随机抽取的4张图像中的推理效果。
![图 8 目标框(A)和关键点(B)的P-R曲线]()
图 8 目标框(A)和关键点(B)的P-R曲线Figure 8. P-R curves for bounding boxes (A) and key points (B)3.2 结果分析
3.2.1 ByteTrack追踪效果
在评价追踪任务性能时，主要关注3个指标：每秒处理帧数(Frames per second，FPS)；IDF1分数(Identity F1 score，IDF1)；多目标跟踪准确率(Multiple object tracking accuracy，MOTA)。FPS衡量系统每秒能处理的图像帧数，直接反映系统的处理速度；IDF1反映跟踪过程中正确维持目标身份的比例，表明在多目标跟踪中目标身份连续性的保持情况；MOTA是评估多目标跟踪性能的综合指标，通过综合考量正确追踪到的目标数量、漏检目标数量、误检目标数量及目标定位误差等多个方面，提供对追踪系统整体性能的全面评价。IDF1和MOTA的计算公式如下。
$$ {\rm{I D F}} 1=\dfrac{2 {\rm{I D T P}}}{2 {\rm{I D T P}}+{\rm{I D F P}}+{\rm{I D F N}}}\text{，} $$ (8) $$ {\rm{MOTA}} = 1 - \dfrac{{\displaystyle\sum\nolimits_{{t}} {\left( {{\rm{F}}{{\rm{N}}_t} + {\rm{F}}{{\rm{P}}_t} + {\rm{IDS}}{{\rm{W}}_t}} \right)} }}{{\displaystyle\sum\nolimits_t {{\rm{G}}{{\rm{T}}_t}} }} \text{，} $$ (9) 式中，IDTP表示真实目标与跟踪结果正确匹配的数量，IDFP表示真实目标与跟踪结果错误匹配的数量，IDFN表示未匹配的真实目标数量，GTt表示第t帧中真实目标的个数，FNt表示第t帧中漏检的个数，FPt表示第t帧中虚检的个数，IDSWt表示第t帧中轨迹发生转变的个数。
本研究设计了4个测试场景评估3种多目标跟踪算法(FairMOT、ByteTrack、BoT-SORT)的表现。场景1展示了5头成年大猪在无遮挡的环境中自由活动；场景2聚焦于5头小猪，同样在无遮挡条件下；场景3增加了复杂性，引入遮挡物，同时仍有5头成年大猪；最后，场景4呈现了一个挑战性的大规模场景，超过20头猪只共同出现且环境拥挤，没有遮挡物。通过表1可知，在无遮挡和适量猪只情况下，选择ByteTrack算法能够保证猪只跟踪的准确性和稳定性，且能保持较高的处理速度。
表 1 3种多目标追踪算法对猪只追踪的验证结果1)Table 1. Validation results of three multiple target tracking algorithms for pig tracking场景
Scene猪数量
Number
of pigs有无遮挡
Occlusion
situationFairMOT ByteTrack BoT-SORT MOTA/% IDF1/% FPS/(帧·s−1) MOTA/% IDF1/% FPS/(帧·s−1) MOTA/% IDF1/% FPS/(帧·s−1) 1 5 无 86.1 87.7 43 91.3 92.9 142 96.7 90.7 59 2 5 无 84.6 85.4 42 94.7 97.3 142 95.3 97.6 58 3 5 有 73.7 72.3 43 67.6 82.0 144 82.8 88.3 59 4 >20 无 76.3 74.9 36 79.8 86.9 129 84.5 91.4 52 1) MOTA：多目标跟踪准确率，IDF1：IDF1分数，FPS：每秒处理帧数
1) MOTA: Multiple object tracking accuracy, IDF1: Identity F1 score, FPS: Frames per second3.2.2 体尺估计评价指标
本研究采用决定系数(R2)、平均绝对误差(Mean absolute error，MAE)、平均绝对百分比误差(Mean absolute percentage error，MAPE)作为评估指标，衡量模型预测值与真实值之间的偏差程度；决定系数越大，说明模型预测值与真实值拟合程度越高，平均绝对误差、平均绝对百分比误差较小表明误差水平较低，反映模型的高准确性。
3.2.3 猪只体尺估计结果分析
随机抽取5个包含多头猪只的圈栏进行体尺估计(单位为cm)，每个圈栏分3个组别评估不同数据处理策略的效果。第1组采用ROI筛选，仅记录摄像头中心范围内猪只的体长数据，防止摄像头畸变、猪只拥挤导致的误差。第2组在ROI筛选的基础上引入姿态和异常检测过滤，识别并过滤由运动模糊、姿态歪曲等因素导致的误差数据。第3组未采用以上机制。分别对各圈栏3个组别猪只的体长、肩宽、臀宽数据求决定系数、平均绝对误差、平均绝对百分比误差，如表2所示。
表 2 猪只体长、肩宽、臀宽数据偏差对比1)Table 2. Variance comparison of pig body length, shoulder width and hip width data性状
Trait圈栏号(猪只数量)
Circle ID
(number of pigs)第1组 Group 1 第2组 Group 2 第3组 Group 3 R2 MAE/cm MAPE/% R2 MAE/cm MAPE/% R2 MAE/cm MAPE/% 体长
Body
length1(7) 0.83 4.71 5.99 0.89 4.00 4.95 0.50 8.14 9.97 2(6) 0.92 2.75 3.84 0.93 2.50 3.48 0.45 6.50 8.65 3(5) 0.67 2.00 2.26 0.82 1.40 1.50 −1.65 5.80 6.34 4(2) 0.20 12.00 13.98 0.81 6.50 7.10 0.92 4.00 4.15 5(4) 0.00 2.50 3.66 0.89 1.00 1.47 −2.13 5.25 7.65 总计 Total(24) 0.85 3.96 5.08 0.94 2.88 3.61 0.75 6.13 7.74 肩宽
Shoulder
width1(7) 0.81 1.57 6.15 0.83 1.71 7.50 −0.54 4.71 18.82 2(6) 0.61 1.00 5.39 0.81 0.83 4.57 0.04 1.50 7.89 3(5) 0.71 1.00 4.54 0.77 1.00 4.83 0.58 1.40 6.93 4(2) −0.06 4.00 14.59 0.19 3.00 10.24 −1.31 6.00 22.16 5(4) −1.13 2.25 13.77 0.28 3.50 12.00 −5.05 3.75 21.87 总计 Total(24) 0.77 1.63 7.60 0.86 1.25 5.65 0.13 3.17 14.40 臀宽
Hip
width1(7) 0.71 1.86 8.76 0.74 1.86 8.79 −0.33 4.14 18.50 2(6) 0.72 0.67 3.76 0.81 0.67 4.30 0.67 0.83 5.40 3(5) 0.71 1.60 7.45 0.84 1.40 6.65 0.27 2.80 13.99 4(2) 0.27 3.00 11.99 0.59 2.00 7.45 −1.49 5.50 22.60 5(4) 0.25 1.25 8.06 0.80 0.75 5.12 −4.36 4.25 27.14 总计 Total(24) 0.81 1.54 7.80 0.88 1.29 6.50 0.25 3.13 15.63 1)R2：决定系数，MAE：平均绝对误差，MAPE：平均绝对百分比误差
1) R2: Determination coefficient, MAE: Mean absolute error, MAPE: Mean absolute percentage error试验结果表明，针对5个栏24头猪只，通过将跟踪算法和关键点检测算法结合，能够确保猪只体尺测量的平均绝对误差小于7 cm，平均绝对百分比误差保持在16%以内。在引入ROI筛选、姿态和异常检测过滤后，24头猪只体尺的平均绝对误差小于3 cm，体长的平均绝对百分比误差保持在4%以内，肩宽保持在6%以内，臀宽保持在7%以内，同时，体长、肩宽和臀宽的决定系数分别为0.94、0.86和0.88，展示了模型良好的拟合效果。此外，ROI筛选显著提高了数据处理速度，从原先的3.0帧/s提升至19.3帧/s。
4. 结论
本研究针对猪场摄像头俯拍场景(无遮挡环境)，提出了基于关键点检测和多目标跟踪的猪只体尺估计方法，在Yolov8-Pose的基础上融合了ByteTrack追踪算法、ROI筛选、姿态和异常检测过滤机制，主要优点如下。
1)相较于其他特征提取或轮廓测量方法，该方法无需人工筛选图片，同时能够捕捉更多猪只的姿态变化信息，提供更精确的结果，且适用于猪场实际生产环境，在无需人工干预的情况下可以对多头猪只进行自动测量。与三维重建测量方法相比，这种方法无需大量计算资源，无需建立特定测量场景，同时避免了猪只剧烈活动时进行三维重建的困难。
2)引入ROI机制可以降低猪群拥挤因素、边缘畸变性对体尺估计的影响，并提高数据处理速度。
3)加入姿态和异常检测过滤算法，避免了运动模糊、姿态不正等干扰因素导致的误差增加。
综上所述，本研究提出的体尺估计方法为猪场生产场景提供了一个轻量化、易部署、高精度的自动体尺测量解决方案，能够在养殖场边缘和终端设备上顺畅运行，助力养殖管理的智能化。
-
![]()
图 1 PRRSV GP5氨基酸序列比对分析结果
PS：信号肽；D：诱导表位；PNE:主要中和表位；HVR1、HVR2：高变区；T1、T2：T细胞表位；B：B细胞表位；TM1、TM2、TM3：跨膜区域；黄色高亮部分表示糖基化位点
Figure 1. Comparison analysis of amino acid locus variation in PRRSV GP5
PS: Peptide signal; D: Decoy epitope; PNE: Principal neutralizing epitope; HVR1, HVR2: Hypervariable region; T1, T2: T cell epitope; B: B cell epitope; TM1, TM2, TM3: Transmembrane regions; Yellow highlighted region refers to glycosylation site
![]()
图 2 ORF5基因遗传进化分析
采用NJ法建树，自举值检验为1 000
Figure 2. Genetic evolution analysis of the ORF5 gene
The tree is built using the NJ method, and the bootstrap value test is 1 000
表 1 引物和扩增片段长度
Table 1 Primer and amplified fragment length
引物名称
Primer name引物序列(5'→3')
Primer sequence长度/bp
LengthORF5-1-F TGAGACCATGAGGTGGGC 726 ORF5-1-R GAAAACGCCAAAAGCACC ORF5-2-F TGCTCCATTTCATGACAC 974 ORF5-2-R GCATCTGGAGGTGATGAAT 
下载: 导出CSV
-
[1] SNIJDER E J, KIKKERT M, FANG Y. Arterivirus molecular biology and pathogenesis[J]. J Gen Virol, 2013, 94(10): 2141-2163.
[2] FENG Y, ZHAO T, NGUYEN T, et al. Porcine respiratory and reproductive syndrome virus variants, Vietnam and China, 2007[J]. Emerg Infect Dis, 2008, 14(11): 1774-1776. doi: 10.3201/eid1411.071676
[3] TIAN K, YU X, ZHAO T, et al. Emergence of fatal PRRSV variants: Unparalleled outbreaks of atypical PRRS in China and molecular dissection of the unique hallmark[J]. PLoS One, 2007, 2(6): e526. doi: 10.1371/journal.pone.0000526
[4] TONG G Z, ZHOU Y J, HAO X F, et al. Highly pathogenic porcine reproductive and respiratory syndrome, China[J]. Emerg Infect Dis, 2007, 13(9): 1434-1436. doi: 10.3201/eid1309.070399
[5] ZHOU Y J, HAO X F, TIAN Z J, et al. Highly virulent porcine reproductive and respiratory syndrome virus emerged in China[J]. Transbound Emerg Dis, 2008, 55(3/4): 152-164.
[6] 宁宜宝, 郑杰, 张纯萍, 等. 我国南方猪高热病的研究: Ⅱ: 猪繁殖与呼吸障碍综合征病毒的分离、鉴定和致病性测定[J]. 中国兽药杂志, 2007(1): 14-18. doi: 10.3969/j.issn.1002-1280.2007.01.005 [7] LUNNEY J K, FANG Y, LADINIG A, et al. Porcine reproductive and respiratory syndrome virus (PRRSV): Pathogenesis and interaction with the immune system[J]. Annu Rev Anim Biosci, 2016, 4: 129-154. doi: 10.1146/annurev-animal-022114-111025
[8] 邓蓉, 郑文堂, 肖红波. 中国生猪生产区域布局与经营规模分析[J]. 现代化农业, 2016(10): 44-46. doi: 10.3969/j.issn.1001-0254.2016.10.024 [9] ZHAO K, YE C, CHANG X, et al. Importation and recombination are responsible for the latest emergence of highly pathogenic porcine reproductive and respiratory syndrome virus in China[J]. J Virol, 2015, 89(20): 10712-10716. doi: 10.1128/JVI.01446-15
[10] ZHANG Q, JIANG P, SONG Z, et al. Pathogenicity and antigenicity of a novel NADC30-like strain of porcine reproductive and respiratory syndrome virus emerged in China[J]. Vet Microbiol, 2016, 197: 93-101. doi: 10.1016/j.vetmic.2016.11.010
[11] LIU J K, WEI C H, DAI A L, et al. Complete genomic characterization of two European-genotype porcine reproductive and respiratory syndrome virus isolates in Fujian province of China[J]. Arch Virol, 2017, 162(3): 823-833. doi: 10.1007/s00705-016-3136-9
[12] BAI X, WANG Y, XU X, et al. Commercial vaccines provide limited protection to NADC30-like PRRSV infection[J]. Vaccine, 2016, 34(46): 5540-5545. doi: 10.1016/j.vaccine.2016.09.048
[13] 张金辉. 生猪禁养政策持续加码, 养猪企业何去何从?[J]. 猪业科学, 2016(11): 142-143. doi: 10.3969/j.issn.1673-5358.2016.11.051 [14] 王亚辉, 彭华, 李娟. 我国21个地区公布畜禽禁养时间表2016年环保政策将更严[J]. 中国猪业, 2016(5): 13-18. doi: 10.3969/j.issn.1673-4645.2016.05.002 [15] 刘岳. 猪繁殖与呼吸综合征病毒GP5蛋白分段表达及鉴定[D]. 泰安: 山东农业大学, 2010. [16] XIE J, CUI T, CUI J, et al. Epidemiological and evolutionary characteristics of the PRRSV in Southern China from 2010 to 2013[J]. Microb Pathog, 2014, 75: 7-15. doi: 10.1016/j.micpath.2014.08.001
[17] SHI M, LAM T T, HON C C, et al. Molecular epidemiology of PRRSV: A phylogenetic perspective[J]. Virus Res, 2010, 154(1/2): 7-17. doi: 10.1016/j.virusres.2010.08.014
[18] LU W H, TUN H M, SUN B L, et al. Re-emerging of porcine respiratory and reproductive syndrome virus (lineage 3) and increased pathogenicity after genomic recombination with vaccine variant[J]. Vet Microbiol, 2015, 175(2/3/4): 332-340. doi: 10.1016/j.vetmic.2014.11.016
[19] 申思. HP-PRRSV囊膜糖蛋白5糖基化位点突变及功能初探[D]. 泰安: 山东农业大学, 2016. -
期刊类型引用(3)
1. 刘继芳，周向阳，李敏，韩书庆，郭雷风，迟亮，杨璐，吴建寨. 人工智能驱动畜牧新质生产力高质量发展：制约因素、生成逻辑与推进路径. 智慧农业(中英文). 2025(01): 165-177 . 
百度学术
2. 王亚彬，徐爱俊，周素茵，叶俊华. 基于Byte的生猪多目标跟踪算法. 农业工程学报. 2025(07): 145-155 . 百度学术
3. 李振波，孙浩翔，郭倩男，张涵钰，刘皓南. 机器视觉技术在大体型家畜无接触体尺测量中的研究进展. 农业工程学报. 2025(07): 1-12 . 百度学术
其他类型引用(1)
