《Transboundary and Emerging Diseases》：Development and Application of a Rep-Targeted SYBR Green I qPCR Assay for Porcine Circovirus 4: Etiological Investigation in China (2022–2025) and Molecular Evolutionary Analysis
编辑推荐：
猪圆环病毒（PCV）4型（PCV4）是一种新出现的圆环病毒属病毒，其流行病学和进化动态尚未完全阐明。研究人员开发了一种靶向复制相关蛋白（Rep）基因的SYBR Green I实时荧光定量PCR检测方法，该方法对测试的常见猪病原体具有分析特异性，估计95%检出限
猪圆环病毒（PCV）4型（PCV4）是一种新出现的圆环病毒属病毒，其流行病学和进化动态尚未完全阐明。研究人员开发了一种靶向复制相关蛋白（Rep）基因的SYBR Green I实时荧光定量PCR检测方法，该方法对测试的常见猪病原体具有分析特异性，估计95%检出限（LOD95%）为8.22 copies/μL。在2022年至2025年间从中国29个省份采集的1920份诊断样本中，96份检测为PCV4阳性（5.00%；95%置信区间[CI]：4.11%–6.07%）。母猪的阳性率最高，尽管各生产阶段间的阳性率无显著差异。基于衣壳（Cap）基因的系统发育分析表明，大多数分析序列属于PCV4a亚型。探索性贝叶斯分析估计全球最近共同祖先时间（tMRCA）约为1996.6年，进化速率为4.76 × 10?4 substitutions/site/year。离散系统地理学重建识别出模型推断的迁移联系，而非经证实的传播路径。三个候选正选择位点（Ser27、Ser53和Val98）与预测的结构或表位相关区域重叠。总体而言，该检测方法为PCV4监测提供了实用的分子工具，流行病学和进化发现为未来的系统采样、重组感知基因组分析以及可能的垂直传播研究提出了假设。
论文解读文章
研究背景与目的
猪圆环病毒（PCV）是圆环病毒科圆环病毒属成员，目前已确认四种基因型：PCV1、PCV2、PCV3和PCV4。PCV4是一种新出现的圆环病毒，于2019年在中国首次发现，随后在韩国、泰国、马来西亚等亚洲国家以及美国和欧洲的野猪群中被检测到，显示其地理分布正在扩大。PCV4与其他PCV的核苷酸同源性较低（与PCV1为50.3%，PCV2为51.5%，PCV3为43.2%），而与貂源圆环病毒（MiCV）亲缘关系较近（66.9%）。回顾性研究表明，PCV4可能早在2011年就已在野生猪群中隐秘循环，在欧洲野生动物中的证据可追溯至20世纪90年代末和21世纪初。
目前PCV4的流行病学和进化动态仍不完全清楚。现有检测方法多针对衣壳（Cap）基因，但Cap基因编码主要结构蛋白，含有重要抗原决定簇，可能经历更强的免疫介导选择和更大的序列变异，导致引物或探针结合位点的突变影响检测效率，增加假阴性风险。相比之下，Rep基因通常更保守，受纯化选择约束。因此，本研究旨在开发一种靶向Rep基因保守区域的SYBR Green I实时荧光定量PCR（qPCR）方法，用于筛查中国临床猪样本，表征PCV4感染的时间、地理和生产阶段模式，并通过Cap基因系统发育、贝叶斯进化动力学和系统地理学分析揭示病毒的进化动态和空间传播格局。
研究方法概述
研究人员从2022年1月至2025年11月间，收集了覆盖中国29个省份的1920份猪组织样本（包括流产胎儿、仔猪、保育猪、育肥猪和母猪的组织样本），来自华中农业大学动物疾病诊断中心。主要关键技术方法包括：基于85条PCV4 Rep基因序列比对设计保守区引物，建立SYBR Green I qPCR方法并进行优化；制备标准质粒建立标准曲线；评估方法的特异性（交叉反应性测试）、灵敏度（probit回归估算LOD95%）和重复性（批内和批间变异系数）；与参考TaqMan qPCR进行一致性比较；应用该方法进行临床样本检测；对Cap基因进行最大似然法系统发育重建；采用贝叶斯系统动力学分析（严格分子钟、贝叶斯天际线合并先验）估计进化速率和tMRCA；进行离散系统地理学重建；利用多种方法（SLAC、FEL、MEME、FUBAR）分析选择压力；使用AlphaFold3预测Cap蛋白三维结构并映射候选正选择位点；预测线性B细胞表位和T细胞表位。
研究结果
3.1 PCV4 qPCR方法的分析性能：优化确定最佳退火温度为63℃，最佳引物浓度为400 nM。标准曲线方程为Y = ?3.268X + 38.41（R2 = 0.9991），扩增效率为102.3%。所有PCV4阳性反应产生单一的熔解峰，平均熔解温度（Tm）为85.5 ± 0.5℃。非PCV4病原对照（PCV2、PCV3、PRRSV、PEDV、CSFV、PPV、PRV）均未产生PCV4特异性熔解峰。
3.2 特异性、灵敏度和重复性：该方法对PCV4具有高度分析特异性，无交叉反应。Probit回归分析估计LOD95%为8.22 copies/μL（95% CI: 7.08–11.22 copies/μL）。在三个浓度水平（1.17 × 10?、1.17 × 10?、1.17 × 10? copies/μL）下，批内变异系数（CV）为0.46%–1.26%，批间CV为0.93%–0.97%，均低于1.5%，表明重复性和重现性良好。在接近LOD95%的浓度下，≥1× LOD95%时所有重复均检测到阳性，而0.5× LOD95%时检测一致性下降。
3.3 不同方法在临床样本中的一致性：开发的SYBR Green I qPCR方法与参考TaqMan qPCR方法在40份临床样本（30份阳性分高、中、低病毒载量组，10份阴性）中完全一致，总体符合率为100%（40/40），Cohen's kappa系数为1.00（p < 0.001）。两种方法的Ct值呈极强相关性（Spearman's ρ = 0.9978，p < 0.001），线性回归方程为Y = 1.0021X ? 0.0571（R2 = 0.9995）。Bland-Altman分析显示平均Ct差值为?0.011个循环，95%一致性界限为?0.197至0.174个循环，无比例偏倚证据。
3.4 PCV4临床阳性率及时空分布：1920份诊断样本的总PCV4阳性率为5.00%（96/1920；95% CI: 4.11%–6.07%）。2022至2025年的年度阳性率分别为3.07%、4.70%、5.09%和8.56%，呈现统计学显著上升趋势（χ2趋势 = 11.60，p = 0.00066）。七大地理区域中，东北地区阳性率最高（7.76%，9/116），其次是华南（7.60%，19/250）和华东（5.49%，18/328），华北最低（2.49%，5/201），华中为3.85%（22/572）。由于区域样本量不均，这些数据描述的是提交样本集的模式，不能确立区域群体流行率。
3.5 PCV4病毒载量及不同生产阶段的阳性率：病毒载量范围从3.39 × 102到1.73 × 10? copies/μL（均值4.64 × 10? copies/μL），其中12.5%（12/96）的阳性病例病毒载量超过1.0 × 10? copies/μL。各生产阶段阳性率：母猪最高（7.89%，6/76），其次为流产胎儿（6.27%，20/319）、保育猪（4.98%，14/281）、仔猪（4.68%，40/854）和育肥猪（4.10%，16/390），但差异无统计学意义（p > 0.05）。在有症状和无症状猪的比较中，有症状仔猪的病毒载量显著高于无症状仔猪（调整后p = 0.0006），有症状流产胎儿的病毒载量也显著高于无症状流产胎儿（调整后p = 0.0044）。
3.6 系统发育关系与宿主相关模式：基于15条本研究获得的Cap基因序列和82条参考菌株的最大似然树显示，大部分序列属于PCV4a基因型。本研究的15株新分离株中，13株聚集在PCV4a分支，2株（SXYL_220906N和GDYJ_220905N）位于PCV4b分支。系统发育分析还纳入了来自野猪、奶牛、狗、狐狸、猫和貉等不同宿主的序列，这些序列并未形成独立的宿主特异性单系群，而是分散在猪源菌株之间。国际菌株（西班牙、美国、韩国、马来西亚、泰国）与中国序列相互穿插，表明地理和宿主相关的混合。
3.7 重组筛选：RDP5检测到两个暂定的重组信号，但被归类为可能归因于非重组进化过程；GARD未检测到统计学支持的重组断点。因此，重组未被确认，所有序列保留在后续贝叶斯分析中。
3.8 PCV4 Cap基因进化动态：基于89条Cap基因序列的严格时钟贝叶斯分析估计进化速率为4.76 × 10?4 substitutions/site/year（95% HPD: 2.73 × 10?4 – 7.00 × 10?4），全球tMRCA约为1996.6年（95% HPD: 1980.5 – 2008.4）。贝叶斯天际线重建显示相对有效种群大小相对稳定，随后出现近期下降，但近期的不确定性增大。
3.9 国内系统地理学分析：省级离散系统地理学重建的MCC树显示地理上散布的谱系。BSSVS在选定阈值下保留了13个模型推断的地理联系，这些联系代表基于模型的系统发育关联，而非经证实的传播路径，其数量和方向受各省序列数量少且不均匀的影响。
3.10 国际系统地理学分析：全球Cap基因数据集中，中国序列在数量上占主导地位。马来西亚、西班牙、泰国、韩国和美国序列占据不同的系统发育位置。由于严重的空间采样不平衡，这些位置和BSSVS推导的联系需谨慎解读，不能确定地理起源、传播方向或主要传播枢纽。
3.11 选择压力与蛋白质结构映射：选择分析表明大多数Cap密码子受到纯化选择。MEME、FEL和FUBAR识别出三个候选正选择位点（Ser27、Ser53和Val98），而SLAC未检测到任何位点。AlphaFold3预测的Cap单体pTM评分为0.87。三个候选位点被映射到以PCV2 Cap结构（PDB ID: 6OLA）为模板生成的60聚体病毒样颗粒（VLP）模型上。
3.12 表位预测：三个候选正选择位点与预测的表位相关区域重叠：Ser27位于预测的线性B细胞表位（残基17–36），Ser53位于预测的SLA-I结合肽（残基51–60；排名<0.5%），Val98位于重叠的预测B细胞和SLA-I结合区域（残基93–100；排名<0.5%）。
讨论与结论总结
讨论部分：本研究开发的Rep靶向SYBR Green I qPCR方法具有良好的线性、效率和灵敏度。与TaqMan方法相比，SYBR Green I只需引物对，无需标记水解探针，减少了特异性寡核苷酸数量，便于实施。但SYBR Green I会与任何双链DNA产物结合产生荧光，非特异性扩增子和引物二聚体可能贡献信号；熔解曲线分析不能完全排除具有相似熔解温度的产物。共提取抑制剂可能导致低浓度时的假阴性。在评估的临床样本中观察到完全一致，但不能视为等效分析灵敏度或更广泛诊断准确性的证据，需要在更多低浓度标本和多样化临床基质中进一步验证。
将该方法应用于1920份诊断样本，总阳性率为5.00%。不同研究报道的阳性率差异很大（从2.2%到17%），反映了抽样设计、地理覆盖、研究群体、临床状态、生产阶段、样本类型和检测方法的影响。2022至2025年观察到的上升趋势同样受年度差异影响，不应解释为国家代表性患病率或群体水平增加的明确证据。各生产阶段的阳性率差异虽无统计学意义，但母猪和流产胎儿较高的阳性率值得进一步研究可能的垂直传播。
重组筛选方面，RDP5和GARD的不同假设和搜索策略可能在短的低多样性序列比对中产生不一致结果。本研究未确认重组，未来排除相关序列的分析可提供保守的稳健性评估。
Cap基因贝叶斯分析估计tMRCA约为1996.6年，进化速率为4.76 × 10?4 substitutions/site/year，低于PCV2的报道速率，与PCV3相当。tMRCA比先前PCV4估计值（1952.48）更晚，可能反映采样结构和数据集组成的差异。tMRCA描述的是采样序列共同祖先的时间尺度，而非病毒出现的精确日期。
采样结构影响系统地理学和人口统计推断。重建的地理联系应视为模型推断的关联而非经证实的传播路径。不均匀的地理代表性会影响祖先状态重建和根位置推断。近期相对有效种群大小的下降不应视为群体收缩的明确证据。
三个候选正选择位点与预测的B细胞和/或SLA-I结合区域重叠，提示PCV4 Cap蛋白可能存在免疫相关选择，与PCV2和PCV3的观察一致。但这些计算机模拟关联不能确立对抗原性、免疫逃逸、受体相互作用或病毒适应性的影响。
未来研究应解决这些候选位点的生物学意义以及系统地理学推断的稳健性，包括引入自然发生的单次替换及其组合单倍型到PCV4感染性克隆或复制子中，以及通过平衡采样、牲畜流动数据、重复降采样等方法改进系统地理学推断。
结论翻译：本研究开发的Rep靶向SYBR Green I实时荧光定量PCR方法为PCV4检测和监测提供了实用的分子工具。在测试样本中，PCV4阳性率为5.00%。母猪和流产胎儿中观察到的较高阳性率值得进一步研究可能的垂直传播。基于Cap基因的系统发育、分子钟、系统地理学、选择压力、表位和结构分析产生了关于PCV4进化历史、地理联系和潜在免疫相关变异的探索性假设。
