《Frontiers in Microbiology》：Hybridization capture increases on-target nanopore sequencing of plant RNA tobamovirus- derived cDNA libraries
编辑推荐：
高通量测序（HTS）能够支持植物病毒监测，但宿主核酸常常降低靶向读段回收率。研究人员评估了一种靶向杂交捕获工作流程，其中由添加了冻干Tobamovirus阳性对照的植物RNA提取物生成的条形码双链cDNA（ds-cDNA）文库在Oxford Nanopore测
高通量测序（HTS）能够支持植物病毒监测，但宿主核酸常常降低靶向读段回收率。研究人员评估了一种靶向杂交捕获工作流程，其中由添加了冻干Tobamovirus阳性对照的植物RNA提取物生成的条形码双链cDNA（ds-cDNA）文库在Oxford Nanopore测序前进行了富集。生物素化探针靶向黄瓜绿斑驳花叶病毒（CGMMV，种名Tobamovirus viridimaculae）、辣椒轻斑驳病毒（PMMoV，种名Tobamovirus capsici）和烟草花叶病毒（TMV，种名Tobamovirus tabaci）的保守区域。对每种病毒的四个配对样本，捕获后相对靶标读段丰度从CGMMV的0.76 ± 0.33%增加到37.62 ± 15.72%，PMMoV从8.16 ± 3.86%增加到24.68 ± 12.34%，TMV从15.62 ± 10.40%增加到36.83 ± 30.33%。精确双侧Wilcoxon符号秩检验得出每种病毒P = 0.125；对于四个同向非零差异，这是可达到的最小双侧P值。基因组覆盖广度得以维持。回顾性双重qPCR支持CGMMV的病毒与18S比率增加，显示PMMoV反应可变，并显示TMV的病毒与18S比率下降，因为18S信号提前的程度与TMV信号相当或更大。这些发现为靶标依赖性文库富集提供了概念验证证据，但未确立分析灵敏度、诊断性能或田间有效性。需要使用自然感染、低滴度和混合感染样本进行验证，并与更简单的靶向工作流程进行比较。
论文解读：杂交捕获富集植物RNA Tobamovirus cDNA文库的靶向纳米孔测序研究
一、研究背景与目的
高通量测序（HTS）技术凭借单次运行即可回收多种病原信息的优势，在植物病毒监测、检疫和认证流程中得到日益广泛的应用。然而，植物核酸的鸟枪法测序常被宿主来源分子主导，导致病原读段比例低下，显著增加测序成本。同时，Oxford Nanopore测序因其便携式仪器、灵活的运行时间和长读段生成能力，在分散式或时间敏感型工作流程中具有独特价值，但靶向长读段工作流程仍需针对特定应用场景进行优化。Tobamoviruses是一类刚性杆状、正义单链RNA病毒，基因组约6.3–6.8 kb，可通过机械接触和种子相关污染传播，对生物安全计划具有重要意义。黄瓜绿斑驳花叶病毒（CGMMV）、辣椒轻斑驳病毒（PMMoV）和烟草花叶病毒（TMV）分别代表影响葫芦科、辣椒和茄科作物的典型烟草花叶病毒属成员。在人类和脊椎动物系统中，基于生物素化探针和链霉亲和素磁珠的杂交捕获已有效改善低丰度病毒序列的回收，但在植物病毒学中，特别是针对长读段测序准备的植物RNA病毒文库，探针捕获技术尚未得到充分开发。由于探针捕获本质上具有靶标偏向性，不适合作为新病毒发现的独立策略，其潜在价值在于已知病毒靶标的面板检测，以提高植物主导核酸背景中靶标序列的回收效率。为此，研究人员评估了杂交捕获作为测序前靶向富集步骤，以增加条形码文库中Tobamovirus来源ds-cDNA分子的回收率。
二、研究设计与技术方法
研究人员采用配对设计，比较三种Tobamoviruses捕获前后条形码ds-cDNA文库。健康植物RNA提取物中添加冻干阳性组织对照，经锚定逆转录、模板转换、PCR扩增、杂交捕获和Oxford Nanopore测序处理。每种病毒-宿主组合使用四个重复配对。捕获探针基于NCBI Virus数据库截至2024年4月的完整基因组多序列比对设计，分别比对CGMMV、PMMoV和TMV的161、48和94个完整基因组，从保守区域合成5′-生物素化探针。技术流程包括：使用携带病毒互补区的锚定逆转录引物启动逆转录并引入条形码；模板转换寡核苷酸引入反向PCR手柄；扩增的ds-cDNA文库经变性后与生物素化探针和封闭寡核苷酸杂交，通过链霉亲和素磁珠捕获、洗涤、再扩增和测序。同时，保留的ds-cDNA等分试样采用双重qPCR进行回顾性验证，以独立评估捕获后靶标与宿主信号比例的变化，使用18S作为宿主来源标记物，而非假定恒定的内参基因。生物信息学分析包括Dorado碱基调用、按条形码序列去复用、Filtlong长度和质量过滤、Cutadapt修剪、Minimap2比对和SAMtools处理，采用精确双侧Wilcoxon符号秩检验进行统计分析。
三、研究结果
1. 片段分析结果
片段分析显示，所有三种病毒-宿主组合在杂交捕获后均出现短片段的消耗。捕获前文库的平均大小估计在1.1至2.3 kb之间，捕获后文库在1.2至2.8 kb之间。
2. 靶向读段比例的配对增加
CGMMV的四个配对均显示捕获后相对丰度增加，平均值从0.76 ± 0.33%升至37.62 ± 15.72%，配对富集倍数范围45.16至53.78倍，单位点差值为38.79个百分点，精确Wilcoxon检验P = 0.125，配对t检验P = 0.017。PMMoV的四个配对同样增加，从8.16 ± 3.86%升至24.68 ± 12.34%，配对富集倍数范围1.71至5.14倍，中位配对差值为17.37个百分点，精确Wilcoxon检验P = 0.125，配对t检验P = 0.040。TMV的四个配对也增加，从15.62 ± 10.40%升至36.83 ± 30.33%，但幅度变化较大（1.06至3.90倍），中位配对差值16.82个百分点，精确Wilcoxon检验P = 0.125，配对t检验P = 0.139。由于四个同向非零配对差值的精确双侧Wilcoxon符号秩检验可达到的最小P值即为0.125，结果被描述为观察到的配对增加而非统计学显著增加。TMV的增加幅度较小且变异较大。
3. 双重qPCR验证结果
全部三个双重检测产生跨越五个数量级的近似线性六点标准曲线。CGMMV和PMMoV标准曲线显示R2值0.9995至0.9998，扩增效率94.8%至95.4%；TMV检测保持线性但效率较低（89.8%至90.5%）。CGMMV在全部五个配对文库中，捕获后病毒Ct值均提前，病毒/18S浓度比增加12.49至32.20倍。PMMoV在六个配对中五个显示相同方向，病毒/18S比率增加1.25至18.32倍。TMV在全部四个配对中绝对病毒Ct值均提前，但18S Ct值也提前且幅度相等或更大，导致TMV/18S比率在所有四个配对中下降（0.33至0.62倍）。18S Ct平均值从捕获前的18.68 ± 0.72变为捕获后的15.94 ± 1.13（CGMMV），从20.01 ± 1.63变为18.67 ± 2.47（PMMoV），从19.51 ± 1.13变为17.17 ± 1.44（TMV），表明检测的18S靶点捕获后显得更丰富而非消耗。该模式可能反映18S含分子的非特异性保留和捕获后扩增，或短18S扩增子的基因座特异性行为。
4. 基因组覆盖广度维持
CGMMV的平均覆盖率从捕获前的78.81 ± 17.17%变为捕获后的71.16 ± 6.46%，中位配对差值?4.66个百分点，精确Wilcoxon检验P = 0.625。PMMoV覆盖率在所有捕获前后文库均为100.00%，无法进行推断性检验。TMV的平均覆盖率从99.88 ± 0.19%变为98.93 ± 0.66%，中位配对差值?1.07个百分点，精确Wilcoxon检验P = 0.125。这些结果表明捕获未导致配对文库中基因组覆盖广度的显著观察损失。
四、讨论与结论
研究发现为杂交捕获可增加植物RNA背景下条形码ds-cDNA文库中靶标来源读段代表性提供了概念验证证据。CGMMV的基线靶标丰度最低但变化最显著，PMMoV增加较小且qPCR响应变异较大，TMV测序相对丰度增加但缺乏宿主归一化的qPCR支持。TMV的qPCR结果需特别考量：绝对Ct值提前表明TMV靶点回收增加，但18S Ct同样提前且幅度相等或更大，可能与18S含分子的非特异性保留和高效扩增有关。该方法定位于已知病毒靶向序列生成富集步骤，可补充RT-PCR、扩增子测序或非偏倚宏基因组测序。由于烟草花叶病毒属缺乏3′poly(A)尾，基因特异性锚定逆转录引物可能引入3′近端区域偏向。操作负担方面，杂交本身需18至20小时，总流程增加约一天，对少量已知病毒而言，多重RT-PCR结合扩增子测序通常更快更简单。研究局限包括使用冻干阳性对照而非自然感染样本、样本集仅限于三种病毒且各一种宿主背景、未进行稀释系列灵敏度测试、未评估混合感染、qPCR评价为回顾性且依赖单一18S扩增子、四个配对测序文库限制统计精度。未来工作应包括自然感染和低滴度样本、探针优化、与替代靶向工作流程的直接比较。研究结论为：杂交捕获在受控植物RNA背景下对三种Tobamoviruses产生了观察到的配对靶标读段相对丰度增加，CGMMV效应最大而TMV变异最大，基因组覆盖广度维持，但结果确立的是工作流程可行性而非分析灵敏度、诊断性能、成本效益或田间有效性。
