迷你型文库制备：微型化是被子植物DNA靶向测序中降低单样本成本的有效方法

《Applications in Plant Sciences》：Fun-sized library prep: Miniaturization is a valid method for per-sample cost reduction in targeted sequencing of angiosperm DNA

　　

前提：群体结构的基因组分析对于理解植物遗传多样性在群体间和群体内非常重要。获取大数据集的成本限制了遗传学在保护及相关领域的更广泛应用。靶向测序和低容量自动化液体处理装置可以减少文库制备反应体积和成本。
方法：研究人员使用Angiosperms353进行靶向测序，以0.5×和0.1×反应体积获得768个样本的数据，其中18个样本是相同的。研究人员计算了质量和数量控制统计数据，以比较组织年龄和文库反应体积对靶向核基因和脱靶质体（plastid）基因测序结果的影响。
结果：文库微型化至0.1×降低了成本，并且总体上与0.5×文库表现相当。在完整数据集中，仅十分之一体积（Tenth-Volume Only）文库显示出较小的插入片段大小和较少的已映射序列基因，但映射读数没有减少。在重叠集（Overlap Set）中，0.1×文库具有相同或更好的性能，测序效率没有显著下降。组织类型之间的差异可能反映了采样变异。
讨论：微型化至0.1×显著降低了单样本成本，同时在新鲜和标本被子植物DNA中保持了相当的测序质量。总体而言，微型化为靶向测序提供了一种可靠、经济有效的方法，增加了使用标本馆收藏的可行性，并使得群体水平基因组研究更广泛可及。

**论文解读：文库微型化在被子植物靶向测序中的成本降低与数据质量验证**

**研究背景**
群体遗传学分析是理解植物遗传多样性在群体内和群体间分布的关键手段，在保护遗传学中，遗传多样性直接影响短期适应性和长期进化潜力。然而，获取大规模数据集的高昂成本限制了遗传学在保护管理及其他领域的广泛应用。靶向测序（targeted sequencing）通过选择性地富集目标基因区域，相比全基因组测序可降低约三分之二的成本，但仍需进一步降低单样本支出，以支持大规模群体研究。标本馆标本（herbarium specimens）是重要的历史记录，但其DNA常存在降解问题，传统方法难以高效利用。文库微型化（library miniaturization）通过减少反应体积来降低试剂和耗材成本，此前在植物组织（尤其是标本DNA）和靶向测序中的适用性尚未得到系统验证。本研究旨在回答三个问题：（1）文库微型化是否影响被子植物DNA靶向测序后的数据质量和数量；（2）该方法是否同时适用于新鲜组织和标本组织；（3）微型化是否为有效的成本降低手段。论文发表在《Applications in Plant Sciences》。

**主要技术方法**
研究人员利用Angiosperms353（一组被子植物通用的353个核基因探针）进行靶向捕获测序，样本来源包括新鲜组织（590份）和标本组织（928份），涵盖243个物种、176个属、63个科。通过SPT Labtech Dragonfly和Mosquito低容量自动化液体处理机器人，在0.5×和0.1×反应体积下制备文库。数据处理使用fastp（去除poly-G尾并估计插入片段大小）和HybPiper（组装目标基因并提取外显子与旁侧非编码区序列）。统计分析采用Welch双样本t检验、方差分析（ANOVA）及Tukey诚实显著性差异（HSD）检验，并进行了多变量方差分析（MANOVA）和效能分析。

**研究结果**

**完整数据集分析**
在768个0.5×文库与768个0.1×文库的完整数据集中，0.1×文库的平均插入片段大小显著减小（t检验：t_1196.6 = -12.402, P < 0.001），具有映射序列的基因数也显著减少（t检验：t_1464.4 = -6.82, P < 0.001），但映射读数总数无显著差异（t检验：t_1247.4 = 1.556, P = 0.12）。

**重叠集直接比较**
为直接比较反应体积的影响，研究人员对18个相同的标本DNA样本分别以0.5×和0.1×制备文库并测序（Overlap Set）。结果显示，0.1×文库的插入片段大小无显著下降（t检验：t_32.7 = 0.952, P = 0.348），基因恢复数无显著减少（t检验：t_28.7 = -0.934, P = 0.358），总“超级重叠群”（supercontig）碱基恢复数无显著差异（t检验：t_28.4 = 0.247, P = 0.807），而映射读数总数反而显著增加（t检验：t_21.4 = 2.55, P = 0.0184）。热图显示，0.1×与0.5×文库在单个基因恢复百分比上无系统性差异。

**家族与体积效应**
MANOVA分析表明，文库反应体积对综合测序指标有显著影响（Pillai’s trace = 0.460, F_4,25 = 5.33, P = 0.003），分类学家族也有显著影响（Pillai’s trace = 0.685, F_12,81 = 1.996, P = 0.035），但二者交互作用不显著（Pillai’s trace = 0.393, F_12,81 = 1.018, P = 0.44）。单变量ANOVA表明，仅文库体积对映射读数有显著影响（F_1,28 = 6.43, P = 0.017）。

**组织类型与体积交互效应**
在完整数据集中，插入片段大小、映射读数和基因恢复数均存在显著差异（ANOVA：P < 0.05），但Tukey HSD检验显示，并非所有组织类型与体积组合间均有显著差异。例如，0.1×标本组织与0.1×新鲜组织文库的插入片段大小无显著差异；0.5×新鲜组织文库与0.1×标本组织或0.1×新鲜组织文库的基因恢复数也无显著差异。双因素ANOVA表明，映射读数的差异主要由组织类型（而非文库体积或交互作用）驱动。

**质体基因回收**
对子集样本的质体（plastid）基因回收分析显示，映射至质体目标的读数比例最高仅为1.3%，44.4%的0.5×样本和16.7%的0.1×样本未恢复任何质体基因。当映射读数超过2000时，基因恢复更稳定，但体积间存在个体差异。

**成本比较**
0.1×文库的单样本总成本（耗材、试剂、质控）为5.51美元，0.5×为23.26美元，1.0×为36.75美元。0.1×成本仅为1.0×的15%，且手工操作时间从14小时降至3小时。

**讨论与结论**
研究人员总结道，微型化至0.1×反应体积能显著降低单样本成本，同时在新鲜和标本被子植物DNA中维持可比的测序质量。虽然完整数据集中0.1×文库的插入片段显著减小、基因恢复数略降，但重叠集直接比较未发现显著负面效应，且映射读数有所增加。组织类型差异主要源于采样变异而非方法本身。标本组织文库在0.1×体积下表现与0.5×新鲜组织文库相当，表明DNA降解程度的影响大于反应体积。研究还指出，质体基因回收率较低可能与靶向捕获效率提高导致脱靶读数减少有关，建议需要细胞器数据的研究采用“加标”（spiking）方法。初始设备成本较高（约19万美元），但高样本量可在20次测序运行（每次960样本）后收回成本。本研究确认了文库微型化对于新鲜和标本组织的靶向测序是一种可靠、经济有效的方法，可扩大标本馆标本在遗传研究中的应用，降低群体水平基因组研究的门槛，有利于保护遗传学、系统发育等领域的发展。