《Journal of Fish Diseases》：High-Fidelity Whole-Genome Sequencing of White Spot Syndrome Virus From Infected Shrimp via Iodixanol Density Gradient Centrifugation and Enzymatic Host Depletion
编辑推荐：
白斑病（WSD）由白斑综合征病毒（WSSV；Whispovirus xiabaidian）引起，对全球虾类养殖造成灾难性损失。表征毒株水平的毒力变异需要比较基因组学；然而，由于缺乏连续甲壳类细胞系，受感染组织的直接测序因压倒性的宿主DNA背景而严重受限。本研究
白斑病（WSD）由白斑综合征病毒（WSSV；Whispovirus xiabaidian）引起，对全球虾类养殖造成灾难性损失。表征毒株水平的毒力变异需要比较基因组学；然而，由于缺乏连续甲壳类细胞系，受感染组织的直接测序因压倒性的宿主DNA背景而严重受限。本研究建立了一种高效的病毒颗粒纯化方案，将碘克沙醇密度梯度超速离心与酶切核酸酶消化相结合，以选择性去除宿主核DNA和线粒体DNA。所开发的方案在Illumina MiSeq和Oxford Nanopore MinION平台上均将目标WSSV读段富集至>70%，同时在非扩增文库中几乎消除了覆盖度偏差。逐步降采样分析表明，研究人员建立的参考引导组装流程可从仅5–10 Mbp的测序数据中可靠重建约300-kbp的完整基因组，实现每个分离株约10美元的低成本多重测序。应用该优化框架，研究人员成功解析了来自泰国（TH-W11）和菲律宾（PHR3-ZS-M101724，代表该国首个完整WSSV基因组）的两株新型田间分离株的完整基因组，揭示了独特的系统发育谱系和约5.6–5.8 kbp的大基因组缺失，并通过侧翼PCR进行了实验验证。该经济高效的平台为新兴WSSV分离株的高通量基因组监测和分子流行病学提供了标准化流程。
**研究背景与科学问题**
虾类是全球重要的水产养殖品种，东南亚和南美洲的养殖活动尤为活跃。白斑综合征病毒（WSSV；Whispovirus xiabaidian）是威胁虾类养殖最严重的病原之一，感染后3–10天内可导致高达100%的累积死亡率。WSSV是有囊膜的双链DNA病毒，基因组约280–315 kbp，属于已知最大的动物病毒之一。流行病学调查已证实全球虾类养殖区均存在WSSV，且不同分离株之间存在显著的毒力差异。全基因组流行病学分析表明，冷冻虾的跨境贸易是WSD全球传播的主要驱动力。然而，由于缺乏可连续传代的甲壳类细胞系，比较基因组学研究只能依赖从田间样本直接进行全基因组测序和一致性序列重建。传统方法中，从感染组织提取总DNA测序会产生大量宿主基因组读段，严重稀释病毒数据；而使用蔗糖或CsCl密度梯度超速离心的传统病毒颗粒纯化方案会产生高渗和高黏度条件，常损伤囊膜完整性。碘克沙醇密度梯度离心可在严格等渗条件下分离颗粒，有效保持病毒粒子结构完整性，但此前研究主要关注颗粒产量和形态纯度，其在高通量基因组测序中的实际应用和定量效果尚未被系统评估。因此，本研究旨在建立一套专为病毒全基因组测序优化的WSSV颗粒纯化方案，整合碘克沙醇密度梯度离心与酶切核酸酶消化以选择性去除宿主DNA背景，并系统评估全基因组扩增（MDA）相关的覆盖度偏差，最终应用于田间新分离株的完整基因组重建和特征分析。
**关键技术方法**
研究人员以WSSV-JP02株（GenBank：AP027280.1）人工感染日本囊对虾（Penaeus japonicus）制备实验样本，并从泰国和菲律宾的水产养殖场及当地零售市场采集WSSV阳性凡纳滨对虾（Penaeus vannamei）田间样本。主要技术方法包括：（1）碘克沙醇密度梯度超速离心联合Turbo Nuclease酶解纯化WSSV病毒颗粒；（2）qPCR绝对定量病毒基因组（vp28）、宿主核基因组（Hsp70）和线粒体基因组（MT-ND2）拷贝数；（3）Illumina MiSeq短读长和Oxford Nanopore MinION长读长双平台测序；（4）基于PVGA（Pore-guided Viral Genome Assembler）的参考引导基因组组装及逐步降采样分析；（5）侧翼PCR实验验证大基因组缺失；（6）基于核心基因组SNP的最大似然系统发育分析。
**研究结果**
**WSSV病毒粒子纯化与qPCR拷贝数定量**：通过碘克沙醇不连续密度梯度超速离心，在PBS/5%、5%/10%、10%/15%、15%/20%和20%/25%界面形成5条乳白色条带，收集13个连续组分。PCR筛选显示，经核酸酶处理组（Nuc+）在组分5和7–10中检测到vp28扩增信号，未处理组（Nuc?）在组分7–13中检测到信号。qPCR定量表明，组分7病毒富集度最高，核酸酶处理后WSSV/宿主核DNA比值显著提高，WSSV/线粒体DNA比值在组分7中无论是否酶处理均保持较高水平。酶处理后WSSV绝对拷贝数有所下降，推测系因匀浆和离心过程中囊膜或衣壳受损的病毒粒子释放的裸露病毒DNA被酶解所致，而结构完整的病毒粒子内部基因组受保护不被降解。这表明酶解核酸酶消化联合碘克沙醇密度梯度离心可成功消除宿主污染并保持高纯度病毒基因组模板。
**不同富集方案下的读段分类与病毒回收率**：将质量过滤后的读段比对至复合参考基因组进行分类。在ONT MinION平台上，优化条件（Nuc+/MDA+）产生14,115条总读段，其中78.9%比对至WSSV基因组，宿主核DNA污染降至1.7%；而未做核酸酶处理的对照（Nuc?/MDA+）病毒读段仅占8.2%，未比对读段高达73.7%。在Illumina MiSeq平台上，Nuc+/MDA+条件产生96,316条总读段，74.3%来源于WSSV，宿主核DNA污染仅1.2%；而未处理对照（Nuc?/MDA?）中宿主核DNA占81.1%，WSSV读段仅1.5%。Nuc+/MDA?条件（非扩增）病毒读段占8.5%，虽低于MDA+条件，但完全规避了酶促扩增偏差。这些比较明确表明，酶解宿主核酸是MDA前不可或缺的前提条件；优化预处理方案可在两个平台上可靠回收超过70%的总测序输出为目标病毒读段。
**10 Mbp标准化分析下的覆盖度偏差与基因组广度评估**：将非宿主数据标准化至10 Mbp进行降采样分析。Nuc+/MDA+条件在ONT和Illumina平台上分别产生36,658,743 bp和14,750,736 bp非宿主病毒碱基，标准化后两个平台平均比对深度几乎一致，均在≥1×深度下实现100%覆盖广度，在≥10×深度下均匀覆盖91%的WSSV参考基因组。Nuc+/MDA?条件虽非宿主产量仅6,929,129 bp，但覆盖均匀性最高（CVBias=0.43），表明省略MDA可规避phi29聚合酶的随机局部扩增偏差。相反，Nuc?/MDA+和Nuc?/MDA?条件下覆盖度偏差严重（CVBias分别为1.37、1.64和0.90），基因组覆盖广度受损。这表明核酸酶辅助富集方案不仅通过消除宿主基因组背景最大化测序通量，还显著减轻全基因组扩增过程中的覆盖度偏差。
**基因组组装与逐步数据缩减分析**：使用PVGA组装流程，优化文库（Nuc+/MDA+）在全数据集或10 Mbp标准化数据下均产生约305,000–311,000 bp的单一连续一致性基因组。即使输入数据降至7.5、5.0和2.5 Mbp，组装仍保持高结构完整性：Illumina Nuc+/MDA+条件从仅2.5 Mbp输入中恢复超过262,000 bp基因组，ONT Nuc+/MDA+条件在5.0 Mbp下保留超过238,000 bp核心基因组区域。相反，省略核酸酶预处理的条件在测序深度低于10 Mbp时组装严重碎片化，仅产生几千碱基至不足50 kb的碎片化contig。这表明核酸酶辅助MDA富集联合PVGA组装流程可从 minimal 序列产量（5.0–10 Mbp）中准确重建完整病毒基因组。
**侧翼PCR实验验证大基因组缺失**：比较基因组分析鉴定出新测序分离株中存在独特的约5.6–5.8 kbp大基因组缺失。侧翼PCR实验显示，非缺失日本分离株（JP02）产生约6.7-kbp完整扩增子，菲律宾分离株（PHR3-ZS-M101724）产生900–1000 bp之间的截短扩增子，泰国分离株（TH-W11）产生略高于1000 bp的条带，与in silico预测精确吻合。这证实这些大缺失是自然发生的真实结构变异（SVs），而非组装假象，验证了PVGA组装流程准确捕捉宏观基因组重排的保真性。
**系统发育分析**：基于39个WSSV全基因组序列的核心基因组SNP比对（309,286 bp，616个简约信息位点），最大似然系统发育树将全球分离株分为早期分化的历史亚洲谱系和衍生的当代东南亚谱系。泰国新分离株TH-W11以98%的超快自举支持率稳健聚类于当代东南亚谱系内，与泰国尖竹汶府和宋卡府近期分离株形成紧密单系群，表明TH-W11是活跃循环的国内谱系，已适应多种对虾宿主。菲律宾分离株PHR3-ZS-M101724则位于历史亚洲谱系的基部分支，具有异常长的终端分支长度，代表在地理隔离中积累了独特基因组替换的早期分化谱系，这是菲律宾首次报道的WSSV全基因组序列，填补了全球WSD基因组监测的重要地理空白。
**成本效益分析**：Illumina MiSeq Nano试剂盒（500 Mbp输出）零售价约76,500日元，单个MinION Flow Cell约159,900日元。基于约10 Mbp输入即可实现完整基因组覆盖，且优化富集方案确保近80%读段来源于WSSV，单个MiSeq Nano试剂盒理论上可多重测序多达50个分离株（500 Mbp/10 Mbp=50），每个病毒基因组测序成本约1530日元（约10美元）。该估算仅涵盖测序消耗试剂，不含上游组织解剖、DNA提取试剂盒、人工费用或设备折旧，但通过将所需数据量降至每基因组10 Mbp，显著降低了每样本运行成本。
**讨论与结论**
本研究建立了一套高效、高保真的WSSV全基因组测序框架，将酶切核酸酶消化与碘克沙醇密度梯度离心相结合。超速离心前用核酸酶预处理组织匀浆可大幅去除污染的宿主核DNA和线粒体DNA，使目标WSSV序列占总测序读段的70%以上。显著的宿主去除确保病毒基因组DNA在纯化模板中占绝对优势，与未处理对照相比有效减轻了深度波动和覆盖度偏差。值得注意的是，对于短读长测序，省略MDA（Nuc+/MDA?条件）通过规避phi29聚合酶相关的随机扩增偏差，产生了最高的覆盖均匀性，表明全基因组扩增应根据初始模板可获得性和所选测序平台的具体技术约束选择性使用。该方案的实践效果和保真性通过新全球田间分离株的完整全基因组重建得到了实验验证，新解析基因组中鉴定的大基因组缺失经侧翼PCR证实为真实的自然结构变异而非生物信息学组装假象。系统发育基因组学分析解析了精细传播动态，揭示了泰国对虾物种间的活跃宿主转换，并提供了菲律宾首个基因组序列以填补全球WSSV演化格局的关键空白。由于优化流程可从 minimal 测序产量（5.0–10 Mbp）中一致实现连续的全长一致性组装，在保持碱基级准确性的同时显著降低了测序成本。因此，该高通量、经济高效的方案为未来WSSV的大规模分子流行病学、接触追踪和演化监测提供了强有力的标准化平台。此外，由于该纯化方案产生的病毒颗粒宿主基因组污染可忽略不计，其在单病毒粒子基因组学中具有特殊前景，特别是通过基于微型琼脂糖凝胶的单扩增基因组（SAG）分析。在甲壳类连续传代细胞培养系统持续缺失的情况下，该方法有望成为直接从田间水生样本研究病毒种群动态、基因组可塑性和微演化途径的重要方法学改进。
