澳大利亚与新西兰 OsHV-1「μVar」类变异毒株随时间的全基因组表征

《Virus Research》：Whole genome characterisation of Australian and New Zealand OsHV-1‘μVar’ like variants over time

　　

研究人员从澳大利亚和新西兰太平洋牡蛎（Magallana gigas）死亡事件中获集样本，组装出15个 OsHV-1 μVar-like 全基因组（澳大利亚12个、新西兰3个来自公开原始数据）。全基因组系统发育（whole-genome phylogeny）将澳新病毒归入一个区别于欧洲与亚太已知基因组的独立分支（clade）。研究人员检测到具有州别与时间结构的分化，包括一个塔斯马尼亚2024年谱系。澳新样本与 MF509813.1 共有约6.5 kb 缺失（deletion），重复丰富区域出现反复发生的插入缺失（indels）。基于标记基因（marker genes）如 ORF4/C-region 的基因树未能发现该分支分化，因此研究人员支持采用全基因组分型（genome-wide typing）与监测（surveillance）来识别 OsHV-1 变异与地理谱系。

研究背景方面，太平洋牡蛎（Magallana gigas）养殖在全球软体动物水产中占重要份额，而太平洋牡蛎大量死亡事件主要由牡蛎疱疹病毒1型（Ostreid herpesvirus 1，OsHV-1）引起。2008年欧洲暴发与 OsHV-1 微变异毒株（microvariant，μVar）相关，法国群体死亡率达80–100%。OsHV-1 属 Malacoherpesviridae 科 Ostreavirus 属，参考基因组约207 kbp，具有 U_L、U_S 及末端/内部反向重复区 TR_L/IR_L、TR_S/IR_S，并含独特 X 区。现有分类多依赖 ORF4 微卫星区突变，会掩盖删除易感区、替换富集区、宿主互作/复制/膜相关基因差异及全基因组多样性，因此难以反映地理谱系与自然选择压力。澳大利亚2010年乔治河、2013年霍克斯伯里河、2016年塔斯马尼亚、2018年南澳大利亚，以及新西兰2010–2011年夏季都报告过 μVar-like 死亡事件，但此前澳新研究多只测 ORF4，缺少全基因组层面证据。

本研究论文发表于《Virus Research》。研究人员对2010–2024年澳大拉西亚太平洋牡蛎死亡样本开展全基因组测序与系统发育分析，目的是判断澳新 OsHV-1 μVar-like 与欧洲、亚太 μVar 的进化邻近性，并识别可能改变关键基因表达的突变。

主要关键技术方法方面，研究人员使用2010–2024年澳大利亚新南威尔士、南澳大利亚、塔斯马尼亚 farmed/wild 太平洋牡蛎样本12个，以及新西兰2011年公开 Illumina HiSeq 双端原始数据3个；短读长经 fastp 过滤、比对太平洋牡蛎基因组去宿主、SPAdes 从头组装、DIAMOND blastx 鉴定 OsHV-1 来源；用 MAFFT 比对至 AY509253.2 与 MF509813.1 并重构重复区为单拷贝；对重复边界、低相似区和候选大缺失用长片段 PCR 与 Nanopore 长读长验证；全基因组用 MAFFT 比对、IQ-TREE 建最大似然树；用 GARD/HyPhy 检测系统发育信号变化；用 Liftoff 多参考转移并人工校正完成注释，筛选核心基因做分区系统发育；用 C-region 与五区域标记基因树比较单标记分类能力；用2,000 bp 分窗计算与参考基因组两两一致性。

研究结果方面，3.1 组装基因组概要：研究人员构建15个基因组，11个完整、4个部分；部分基因组因 ≥5 读深覆盖低于95%被定义为 partial，基因组长度最大相差133 bp，113个非片段且符合完整度标准的基因进入核心基因系统发育，约占完整基因组71%。

3.2 澳新与其他 OsHV-1 基因组系统发育关系：GARD 发现9个候选断点并将比对分为10个分区，但主要全基因组拓扑不被单一分区驱动；全基因组树与113个核心基因树都支持澳新序列成独立 clade，中国多数基因组最远，欧洲成另一支，MF509813.1 因共有约6.5 kb 缺失在全基因组树中更近；澳新内部新南威尔士2010–2013、南澳2018、塔斯马尼亚2016–2019 聚支，塔斯马尼亚2024三个野生样本另成一支并与新西兰2011样本共基部分支；核心基因树内澳新分支核苷酸分歧较全基因组树压缩约2–5倍。

3.3 OsHV-1 标记基因系统发育关系：C-region/ORF4 树未恢复澳新独占单系，澳样与单个中国序列 KM886595.1 的聚类无 >95% 支持，说明单标记信息不足；五区域树中澳新样本形成独立于中国组、姐妹于欧洲 clade 的分支，证明标记选择会扭曲变异分类。

3.4 两两一致性比较：澳新基因组逾90%区域与 MF509813.1、MW412420.1、AY509253.2、OM811597.1 近完全保守；低相似段集中在基因组末段 IR_L、IR_S、TR_S、U_S 与 X 区，并由100–5,000 bp 插缺造成；所有澳新样本与 MF509813.1 共有约6.5 kb 缺失，删除 IR_S 后 ORF120–123 所在区段，使二者整体相似性更高。

3.5 新 ORF 与已分类 ORF 偏差：每完整基因组见123个 ORF，ORF36、ORF48、ORF123 在参考或澳新组中缺失/删除；ORF38 缩短604 bp 见于除 MF509813.1 与 AY509253.2 外所有样本；ORF70/115/117/119 因移码或无义突变缩短；ORF5/32/33/62/63/65 稳定呈多片段预测，并非澳新特有新假基因化。功能已知/预测基因中，DNA聚合酶（ORF100）、RR1（ORF51）、RR2（ORF20）、核酸外切酶（ORF95）、IAP（ORF42、ORF87）、RING指蛋白（ORF9）、包膜膜蛋白（ORF41）等8个基因有27个编码序列变异；5个仅存于塔斯马尼亚2024，3个存于塔斯马尼亚/南澳2016–2019，其余遍在于澳新样本。

讨论部分总结方面，研究人员指出澳新 μVar-like 并非2008法国 μVar 简单扩散，中等大小 indels 勾勒区域结构，SNP 显示澳大利亚14年间时间分化；单靠 ORF4 会使欧洲 μVar 与澳新 μVar-like 显得过近，掩盖多分支结构。东亚被既往研究视为 μVar 可能起源地，但本全基因组数据不能独立证明起源或传播方向；欧洲 μVar 在 C-region 树中被嵌套的现象在信息更丰富树中不被支持。塔斯马尼亚2024野生牡蛎亚支可能与野生宿主选择压力、未采样谱系引入、批次/覆盖差异混淆，不能由现有数据区分，需同期 wild/farmed 抽样与形式化系统地理学分析。大片段 ORF 破坏均未发生在功能已知基因；澳新毒株应被视为有别于欧洲 μVar 的谱系。

研究结论部分翻译：研究人员对澳大利亚与新西兰死亡事件中的15个 OsHV-1 基因组进行表征；12个为新测序，3个来自既往研究原始数据。尽管所有澳新基因组被认为近似2008法国暴发的 OsHV-1 μVar，全基因组系统发育揭示其演化与地理分布比单一变异扩散更复杂。中等大小插入缺失界定区域遗传结构，个体 SNP 凸显澳大利亚过去14年基因组的时间分化；这些模式关联采样地与时期，但不能独立确立病毒移动方向与时间。虽观察到大片段 ORF 破坏，但 OsHV-1 极少基因有功能注释，且已知功能基因中未见大片段破坏。若仅用 ORF4 推断系统发育，欧洲 μVar 与澳新 μVar-like 会比全基因组分析显得更近，从而漏掉采样 OsHV-1 多样性的多分支结构。澳新基因组与欧洲 μVar 及中国基因组均不同，应视为不同于 μVar 的谱系；C-region 相似不等于全基因组等同。塔斯马尼亚2024亚支成因不能为当前数据所分辨，但其基因组独立性对育种计划与水产生物安全有意义。研究人员强调应以全基因组分型与监测替代单标记分类，以捕捉 OsHV-1 变异、地理谱系与潜在选择压力。