《Canadian Journal of Microbiology》：Coincidental observation of a parvovirus during a targeted bacterial survey of cattle suffering from septic arthritis
编辑推荐：
基于杂交探针富集的DNA捕获测序（CapSeq）实现了高通量、靶向的基因至微生物检测。尽管CapSeq平台设计了特异性，非靶向DNA仍代表了额外的基因与目标微生物。为此，从头宏基因组组装（de novo metagenome assembly）CapSeq数据
基于杂交探针富集的DNA捕获测序（CapSeq）实现了高通量、靶向的基因至微生物检测。尽管CapSeq平台设计了特异性，非靶向DNA仍代表了额外的基因与目标微生物。为此，从头宏基因组组装（de novo metagenome assembly）CapSeq数据集提供了发现机会。本研究报道了一个实例。对来自患有化脓性关节炎（septic arthritis, SA）牛的样本进行CapSeq，结果检测并组装了两个有蹄四细小病毒1型（ungulate tetraparvovirus 1, UTPV1）基因组。这些偶然的检测与组装来源于一个旨在表征细菌而非病毒病原体的平台。获取匹配的鸟枪法宏基因组数据集进一步证实了该结果。重要的是，SA是肉牛跛行的重要诱因，目前缺乏有效且获批的治疗方案，其病原体清单尚不完整。此次宏基因组观察提示病毒因素可能参与牛SA的发病机制，而该病目前被视为细菌性疾病。该结果例证了无偏倚的CapSeq数据挖掘在支持临床发现方面的潜在应用价值。
研究背景
牛化脓性关节炎（septic arthritis, SA）是导致肉牛跛行的主要原因之一，但目前缺乏有效且经批准的治疗方案，其已知病原体清单亦不完整。传统上，SA被视为一种细菌性疾病。然而，人类关节炎的临床研究表明，包括人细小病毒B19在内的病毒性病因同样值得关注。与此同时，基于杂交探针富集的DNA捕获测序（CapSeq）技术因其高通量和高保真度而被广泛应用于微生物学领域，例如检测抗菌药物耐药性（antimicrobial resistance, AMR）基因或鉴定蝙蝠冠状病毒。尽管CapSeq设计用于靶向富集特定DNA，但非靶向DNA也会被同时测序，这为从看似“背景”的数据中发现意外微生物提供了机会。在此背景下，研究人员在利用CapSeq平台对加拿大西部育肥场牛源样本进行细菌病原体和AMR基因调查的过程中，偶然发现了病毒序列的存在。
研究目的
鉴于上述背景，本研究旨在通过无偏倚地挖掘CapSeq宏基因组数据，验证能否从专门针对细菌设计的靶向测序数据中检测并组装出完整的病毒基因组，从而探索病毒因素在牛SA中可能扮演的角色，并强调充分利用宏基因组数据进行科学发现的潜力。
主要关键技术方法
- 1.
宏基因组组装（Metagenomic assembly）：使用SqueezeMeta流程（v1.7.0）对51份CapSeq数据集进行宏基因组组装，其中重叠群由MEGAHIT软件拼接，并对长度大于1000 bp的重叠群进行人工核查。
- 2.
靶向病毒检测（Targeted detection of UTPV1）：使用fastp对CapSeq数据集和一份鸟枪法宏基因组数据进行预处理和质量控制，随后通过bowtie2比对去除宿主（Bos taurus）序列。处理后的读数以95%的同一性阈值映射到相关病毒序列集合（NCBI Virus taxid: 1096141）。
- 3.
病毒序列组装与分析（Assembly and analysis of UTPV1 sequences）：将获得的推定UTPV1读数使用Geneious组装器进行从头组装，并基于覆盖度和读段保守性修剪得到最终草图基因组序列。利用PHyML软件和TIM3ef+G替代模型，基于19条独特UTPV1序列的非结构蛋白（NS）基因起始端至衣壳蛋白（VP）基因末端的4868个保守位点构建最大似然树。
研究结果
UTPV1基因组的检测与组装
通过对51份CapSeq数据集的宏基因组组装重叠群进行检查，研究人员最初从一份死后关节液样本（PMSJ3）中检测到一个长度为2796 bp的重叠群。该重叠群与源自中国牛蜱的牛hokovirus 2（即UTPV1）同源。进一步的靶向评估显示，第二份关节液样本（PMSJ49）也含有UTPV1 DNA。两份样本分别获得1015和17425条推定属于UTPV1的读数。对于PMSJ3，可用的鸟枪法宏基因组数据也识别出16681条与UTPV1高度同源的读数，证实了CapSeq的发现。组装后得到两个UTPV1草图基因组，长度分别为5368 bp（PMSJ3）和5230 bp（PMSJ49），且读段覆盖度分布均匀，表明其检测依赖于样本中的高丰度而非探针富集。这两个基因组共享99.2%的核苷酸序列一致性，仅在非结构蛋白上有三个氨基酸位点差异，而结构蛋白完全保守。
UTPV1存在于患病牛中且观察依赖于DNA
目前尚无研究建立UTPV1与疾病的因果关系，也未能在实验室成功培养该病毒。本研究中检测到的UTPV1序列增加了该病毒在不同动物来源材料中的少数公开报告。PMSJ3和PMSJ49的UTPV1基因组与美国呼吸道疾病相关的病毒序列亲缘关系最近，且均来自加拿大阿尔伯塔省不同育肥场的死亡牛只。两份样本的细菌群落组成及其推定的SA主要致病菌不同：PMSJ3以支原体属（Mycoplasmopsis）为主，而PMSJ49的败血症推定由嗜组织菌属（Histophilus）引起。UTPV1是否促成SA或其他牛病仍是一个悬而未决的问题。在牛呼吸道疾病（常与SA共病）中，病毒感染被认为会启动复杂的疾病进程，进而涉及多种细菌。此外，人类关节炎的病毒病因学（如人类细小病毒B19）也为这一方向提供了参考。
最大化宏基因组数据的价值
靶向宏基因组数据可作为支持健康、疾病和微生物生态学新兴假说的信息资源库，包括那些在设计实验阶段被认为是“超出范围”的假说。未来，开发专门算法（例如从低丰度数据恢复和生成宏基因组组装基因组的方法）将进一步推动此类努力。本研究首次在患SA的牛中报告UTPV1，且来自于靶向细菌而非病毒的CapSeq实验，突显了从通常被忽视的读段数据中进行发现的潜力。更广泛地说，该结果呼吁重视对所有宏基因组数据集（无论靶向与否）进行无偏倚的宏基因组组装基因组分析所带来的益处。
研究结论讨论
讨论部分指出，虽然尚未建立UTPV1与疾病之间的因果联系，但在牛SA样本中检测到该病毒，结合此前在呼吸系统疾病和跛行牛中的相关报道，提示有必要开展更多研究以明确这些病毒究竟是旁观者还是疾病（尤其是犊牛疾病）的贡献者。本研究的结果强调了从靶向宏基因组数据中进行无偏倚数据挖掘以支持临床发现的潜在应用价值。
