《Viruses》：Genomic Epidemiology and National Seroprevalence Reveal the Widespread Distribution of Palyam Viruses in China
编辑推荐：
分类物种帕利亚姆病毒（Palyam virus，PALV）包含一组广泛分布、库蠓（Culicoides）传播的虫媒病毒（arbovirus），与牛繁殖障碍相关，但其遗传多样性和中国的血清流行率尚未被充分表征。为填补这一空白，研究人员从中国南方哨兵牛中获得29株
分类物种帕利亚姆病毒（Palyam virus，PALV）包含一组广泛分布、库蠓（Culicoides）传播的虫媒病毒（arbovirus），与牛繁殖障碍相关，但其遗传多样性和中国的血清流行率尚未被充分表征。为填补这一空白，研究人员从中国南方哨兵牛中获得29株PALV分离株，对15株代表性毒株进行全基因组测序，并利用新开发的竞争酶联免疫吸附测定（c-ELISA）对中国各地4660头牛开展全国血清调查。这些基因组数据与全球数据集整合后进行全面的系统发育（phylogenetic）和系统地理学（phylogeographic）分析。通过建立全球VP2（病毒蛋白2）框架，中国PALV分离株被归为Chuzan、Bunyip Creek和D’Aguilar基因群（genogroup）。在全球背景下，中国PALV毒株与日本毒株遗传亲缘关系最近，而系统地理学重构提示至少在20世纪80年代和90年代发生过两次来自日本的独立引入事件。调查揭示牛总体血清流行率达46.5%（95% CI：44.7–47.5%），表现出明显的纬度梯度，在32.5° N存在尖锐的生态阈值。该病毒在中国湿润南方为高地方性（hyperendemic），血清流行率为33.0%–88.8%，而在北方地区减弱，血清流行率低于8.0%。这些发现将PALV重新定义为东亚虫媒病毒生态系统中广泛的“沉默威胁”，强调需要协调跨境监测。
该研究发表于《Viruses》。研究背景方面，环状病毒属（Orbivirus，呼肠孤病毒科Sedoreoviridae）包含多种对畜牧业造成重大经济损失的虫媒病毒，如蓝舌病毒（Bluetongue virus，BTV）、非洲马瘟病毒（African horse sickness virus，AHSV）、流行性出血病病毒（Epizootic hemorrhagic disease virus，EHDV）等，但除主要病原体外，同属其他研究较少的成员仍对家畜健康存在潜在威胁，帕利亚姆病毒（Palyam virus，PALV）即属此类被忽视的病原。PALV主要由库蠓（Culicoides）传播，广泛分布于亚洲、非洲、澳洲的热带亚热带地区；其Chuzan病毒（Chuzan virus，CHUV）血清型曾在日本引发大规模犊牛积水性无脑－小脑发育不全（hydranencephaly–cerebellar hypoplasia，HCH）综合征暴发，证实其致病潜力，但PALV整体研究严重不足，缺乏标准化诊断工具，限制了对其分布、循环动力学和致病性的理解。PALV基因组含10个双链RNA（double-stranded RNA，dsRNA）节段（Seg-1至Seg-10），编码7个结构蛋白（VP1–VP7）和3个非结构蛋白（NS1–NS3）；外层衣壳蛋白VP2（Seg-2）和VP5（Seg-6）变异度最高，VP2是决定中和特异性和血清型的主要因子，历史上有多达13个PALV血清型被提出，但分类界限不清，部分血清型VP2序列同一性超过98%，存在分类冗余。既往全球PALV比较基因组分析将毒株分为非洲支系和亚洲－澳洲支系，近期日本15株的全面基因组表征提出日本、中国、印度毒株构成共享亚洲基因池，偶有与澳洲谱系的重配（reassortment），但这些基础研究所用数据集中中国序列缺失或代表性不足，东亚微观演化轨迹、跨境扩散路径、节段特异性重配动态完全未被刻画。中国作为东亚虫媒病毒生态系统关键组成部分，拥有庞大反刍动物种群和关键地理流行病学位置，但PALV的分子与血清学流行病学高度碎片化且大多未绘制；虽自2012年起在云南无症状哨兵牛中首次分离到CHUV，后续在广西牛、青海高原牦牛中也有分离，提示广泛地理足迹，长期哨兵监测还获得CHUV及BCV（Bunyip Creek virus）、DAV（D’Aguilar virus）分离株，但仅限病毒分离和Seg-2、Seg-7部分测序；中国分离株缺乏完整高可信度全基因组序列，加上无大规模血清学调查，导致中国大陆PALV的真实遗传架构、演化来源、综合血清流行率图景严重不明。为此，研究人员在2012–2020年从云南、广东哨兵牛获得29株PALV分离株，对15株代表性CHUV、BCV、DAV进行全基因组测序；整合全球数据集以明确中国分离株分子分型，重构东亚PALV谱系具体演化轨迹、局部重配动态和跨境扩散路径；同时开发群特异性竞争酶联免疫吸附测定（competitive ELISA，c-ELISA）并在2016–2018年对中国15省4660头牛开展横断面血清调查，综合勾画全国暴露格局及其地理空间决定因素，这是中国大陆首次针对流行PALV的大规模基因组与血清流行病学评估。
为开展研究，研究人员用到几个主要关键技术方法：样本队列方面包括2012–2020年云南（4个站点：石钟、江城、芒市、景洪）和广东（汕头1个站点）哨兵牛（云南方为6–12月龄本地黄牛和水牛，广东方为荷斯坦奶牛）采血获得的病毒分离源样本队列，以及2016–2018年全国15省119县129个采样组（按县－年度定义）共4660份临床健康牛血清构成的血清调查队列；关键技术方法有建立并验证基于纯化CHUV全病毒的c-ELISA（诊断截断值PI≥50.0%，敏感性93.3%、特异性97.3%）、以血清转化触发回顾选取抗凝全血做靶向Seg-7（群特异性）和Seg-2（血清型特异性）qRT-PCR筛查并盲传C6/36细胞和BHK-21细胞进行病毒分离与型特异性RT-PCR鉴定、采用FLAC法富集dsRNA并转为cDNA构建文库进行Illumina高通量测序获得全基因组、用MAFFT比对后通过IQ-TREE做最大似然（Maximum Likelihood，ML）系统发育推断、用RDP5筛查重组后经TempEst评估时间信号再以BEAST做贝叶斯（Bayesian）时间尺度系统地理学重构（离散特征连续时间马尔可夫链模型辅以贝叶斯随机搜索变量选择BSSVS划分中国、日本、印度、澳洲、非洲5个地理状态）、全国血清数据用混合效应logistic回归（省份和年份为固定效应，采样组为随机截距）计算调整比值比（adjusted Odds Ratio，aOR），并按1°纬度分箱用广义相加模型（Generalized Additive Model，GAM）惩罚回归样条分析纬度非线性阈值。
结果部分如下。
3.1 中国PALV血清流行率的地理空间分布与多变量风险分析：通过对全国4660份血清检测，总体血清流行率为46.5%（95% CI：44.7–47.5%），除吉林外所有调查省均检出阳性；北方PALV传播明显减弱，河北最低1.4%（95% CI：0.5–3.5%），内蒙古最高8.0%（95% CI：3.5–15.0%）；中南部九个省流行率超过30%，广西达88.8%（95% CI：82.9–92.9%），为高地方性区，集中于湿润亚热带和热带季风区、海拔一般低于1500 m；西藏虽属南方纬度但因高海拔（3000–5000 m）极端高原高山气候，流行率仅0.8%（95% CI：0.2–3.0%）。混合效应logistic回归以内蒙古为参照，北方其他省风险显著更低（吉林aOR=0.1，河北aOR=0.2），中南部风险显著升高（湖北aOR=17.0，广西aOR=89.5，所有p<0.001，aOR均>10.0）。地理分析界定高发地方性廊道（>20%流行率）约束在17.5° N–32.5° N、97.5° E–115° E窗口内，对应湿润中低海拔亚热带／热带季风气候区；Spearman分析显示组水平流行率与纬度显著负相关（ρ=?0.67，p<0.001），GAM进一步验证在约32.5° N存在锐利生态阈值：以南流行率维持高地方性水平，以北传播风险骤降；该纬度模型适用于低海拔区，但西藏因垂直高海拔超越纬度影响。
3.2 中国南方三种PALV血清型共循环：2012–2020年5个哨点共采集2496份血样，获得29株PALV分离株；血清型特异性RT-PCR鉴定为CHUV（n=17）、BCV（n=7）、DAV（n=5）；云南24株（CHUV、BCV、DAV均有），广东5株（BCV四、DAV一）；CHUV为最常见且见于所有4个云南哨点，石钟县占云南分离株75.0%（18/24，含CHUV 14、BCV 2、DAV 2）；分离有明显季节性，72.4%（21/29）在6–9月获得。
3.3 中国PALV分离株全基因组测序：同一血清型同一年同一地点的分离株前期筛显示Seg-2和Seg-7遗传同质性高（相似性>99.0%），故选15株代表性分离株（CHUV、BCV、DAV来源）做全基因组测序；de novo组装获得各株10个节段完整编码序列，平均比对深度超8000×（范围8500–38000×），严格多数规则一致性序列定型保证准确性，次要变异（<5%读段占比）不纳入一致性；序列已上传GenBank。
3.4 建立VP2和VP5系统发育框架以精确分子分型中国PALV分离株：用全球VP2序列建最大似然（ML）系统发育树，全局VP2分为7个高支持（bootstrap≥90%）单系基因群（genogroup）；结合一致性热图分析设定VP2氨基酸同一性73.7%为分子分型边界；据此将15株中国分离株明确归入Chuzan、Bunyip Creek、D’Aguilar三个基因群；中国分离株在各基因群内与日本毒株紧密共聚类，基因群内成对核苷酸同一性≥94.2%、氨基酸同一性≥95.4%；VP5（Seg-6）系统发育聚类模式与VP2框架高度一致，支持分型结果。
3.5 PALV毒株八个内部基因组节段的系统发育分析：对八个内部节段（Seg-1、Seg-3至Seg-5、Seg-7至Seg-10）做ML系统发育，全球序列分为亚洲、澳洲、非洲三大洲支系（bootstrap≥90%），但存在洲际间重配打破此拓扑；中国与日本序列紧密聚类（成对核苷酸同一性87.4%–96.8%）；多数内部节段（Seg-1、Seg-3、Seg-5、Seg-8至Seg-10）中中国CHUV、BCV、DAV与日本对应血清型不形成血清型单系群，而呈交织亚谱系，不论外壳血清型归类；Seg-4和Seg-7则将几乎所有中国分离株归为一个主导簇且嵌套少数日本分离株；节段间系统发育不一致强暗示两国循环谱系间活跃遗传交换。
3.6 PALV基因组节段内外基因重组筛查与检测：用RDP5在全基因组比对中筛查，检测到多个假定内基因重组信号；东亚病毒种群中，日本BCV株（ON-14/E/17）的Seg-3有统计强支持的重组事件，主要亲本为日本DAV株（ON-3/E/17）、次要亲本为新建序中国BCV株（V256）；证实东亚区域循环毒株间存在真实内基因重组；所有检出重组的节段严格排除出后续贝叶斯系统地理学分析。
3.7 PALV谱系的贝叶斯进化时间尺度与离散系统地理学重构：对非重组节段评估时间信号（R2>0.2适用时间尺度贝叶斯推断），平均进化速率为2.70–3.84×10?4替换／位点／年，可比BTV；因复杂历史重配全球最晚共同祖先时间（time to most recent common ancestor，tMRCA）差异大，Seg-1早至约1032 CE（95% HPD：547–1462），Seg-9晚至约1663 CE（95% HPD：1481–1832）；选用时间信号最强的Seg-3（R2=0.62）做高分辨系统地理学重构；时间尺度最大 clade 可信（Maximum Clade Credibility，MCC）树重现亚洲、澳洲、非洲三大地理支系，并重构亚洲支系内传输路径：日本为亚洲PALV谱系最可能源 reservoir（后验概率PP=0.69）；从日本源池向中国至少两次独立传入：较早引入（中国Clade 2）约1982 CE（95% HPD：1974–1990），较近引入（中国Clade 1）约1990 CE（95% HPD：1981–1999）；两起日本至中国地理转移事件均具极高后验概率（PP=0.99）。
讨论部分总结如下。研究人员指出，虽然血清中和试验（Serum Neutralization Test，SNT）是国际病毒分类委员会（International Committee on Taxonomy of Viruses，ICTV）认可的Orbivirus血清分型金标准，但受限于参考试剂稀缺；本研究建立基于序列的框架（VP2氨基酸同一性73.7%阈值）对齐BTV、EHDV等主流orbivirus标准，将中国新分离株明确归为Chuzan、Bunyip Creek、D’Aguilar基因群，理顺历史分类不一致，为中国分子诊断和监测提供可扩展基础。分段dsRNA病毒中内基因重组 historically 被认为罕见，但本研究检出多个重组信号；需注意检出信号 disproportionally 集中在单一历史数据集，短读次世代测序（Next-Generation Sequencing，NGS）de novo组装中嵌合序列易造成计算伪影，但Seg-3（VP3）中V256株重组经Sanger测序严格验证，排除组装错误，证实东亚PALV种群真实发生内基因重组。遗传重配是分段RNA病毒标志进化机制，可快速产生新基因型改变抗原性／传播性；除宏观洲际重配模式外，中国与日本毒株多节段系统发育拓扑严重不一致提供持续跨境扩散和局部重配证据，说明区域间遗传交换是东亚PALV进化的主要驱动，支撑“东亚虫媒病毒生态系统”概念，中日为互联流行病学网络节点；活跃区域重配不断重塑中国南方循环毒株遗传多样性，并可能催生毒力或宿主范围改变的新变种，需持续全基因组监测。系统地理学重构明确日本为中国家系主要源 reservoir，两国间频发的跨境传播反映在毒株高度遗传亲缘上；此定向扩散模式不独见PALV，亦见于BTV、EHDV、AKAV（Akabane virus）及新区域orbivirus如TIBOV（Tibet orbivirus）、YOUV（Yunnan orbivirus）、GXOUV（Guangxi orbivirus）等，重申“东亚虫媒病毒生态系统”统一网络假说。解释日本至中国定向扩散机制：人为层面，80–90年代中国经济改革开放带来密集区域牲畜贸易，可能传入带毒活牛；气候层面，日本南部（九州、琉球）与中国南方同属湿润亚热带气候，库蠓媒介全年活跃、病毒高地方性构成稳定演化 reservoir；向量生物学层面，夏季－秋季虫媒高峰季东亚季风和强低空急流常伴台风路径，形成高效大气通道可将感染库蠓跨东海远距离输运；耦合环境与人因驱动下，中日循环虫媒病毒遗传亲缘强，凸显东亚国家协调监测网络必要性。大规模血清调查揭示中国牛PALV高地方性（46.5%），南北二分法与BTV、EHDV模式一致；空间分层明确高强度传播廊道（17.5° N–32.5° N，97.5° E–115° E）内感染风险呈强纬度依赖，反映库蠓媒介基本生态约束（北方干冷冬季抑制，南方暖湿季风利于）；西藏极低流行（0.8%）虽纬度偏南但因高海拔寒冷限制媒介存活，构成生态反证。南方高地方性形成流行病学悖论：广泛循环却几乎无记载临床暴发；研究人员认为多因素交织：当前循环毒株内在致病性可能偏低（类比中国EHDV常有 endemic 循环而无临床病），高血清流行率建立“地方病稳定”（enzootic stability）机制——配种前广泛自然暴露使母牛获保护免疫，后续妊娠时胎儿受护；诊断忽视也掩盖疾病负担，现有牛繁殖障碍鉴别诊断面板优先BTV、EHDV、AKAV、布鲁氏菌、BVD（Bovine Viral Diarrhea）、IBR（Infectious Bovine Rhinotracheitis）等，PALV系统排除在外；无论如何广泛未监测传播代表“沉默风险”，病毒演化或稳定性破坏可触发暴发，需主动监测。研究局限：遗传学分析仅基于滇粤分离株（哨点后勤限制，其他高流行区无分离株），可能低估中国PALV遗传多样性与演化复杂度；自建c-ELISA基于全CHUV病毒（因其与日本临床关联及初期分离优势），未来用保守重组蛋白和单抗可提升标准化；系统地理学重构受公共数据库序列有限且有偏限制（印度、澳洲PALV基因组稀缺），引入祖先节点定年和扩散路径的不确定性。结论：本研究首次全面表征中国PALV，将其从被忽视共生病原重新定义为东亚虫媒病毒生态系统内广泛分布的“沉默威胁”；整合基因组与血清学证据勾画中国南方媒介驱动传播廊道，追溯病毒起源为多次自日本引入；提出将历史血清型合并为七个基于VP2序列的基因群框架，为全球监测提供合理基础；未来必须从被动观察转向主动风险评估，优先扩宽基因组监测解决演化盲区、实验量化中国毒株致病潜力、明确库蠓特异性媒介种，从而建立协调跨国预警网络以减缓PALV及相关orbivirus在东亚虫媒病毒生态系统中的新发。
