优化溶液中富集超过一百万个古代人类SNP

Optimized in-solution enrichment of over a million ancient human SNPs

作者信息Roberta Davidson, Xavier Roca-Rada, Shyamsundar Ravishankar, Leonard Taufik, Christian Haarkötter, Evelyn Collen, Matthew P Williams, Peter Webb, M Irfan Mahmud, Erlin Novita Idje Djami, Gludhug A Purnomo, Cristina Santos, Assumpció Malgosa, Linda R Manzanilla, Ana Maria Silva, Sofia Tereso, Vítor Matos, Pedro C Carvalho, Teresa Fernandes, Anne-France Maurer, João C Teixeira, Raymond Tobler, Lars Fehren-Schmitz, Bastien Llamas
PMID40611166
发布时间2025-07-03
DOI10.1186/s13059-025-03622-6

实验完整度

高

研究包含24个古代人类样本的深度鸟枪测序、一轮和两轮富集、多种内源DNA含量、合并与非合并富集,以及等位基因偏差和索引跳跃等验证,证据层级包括体外实验、统计分析和样本验证,功能验证明确。

主要模型

24个古代人类骨骼样本(来自伊比利亚、中美洲和东南亚) 测序文库(分别进行鸟枪测序和Twist富集)

重点核对

样本内源DNA含量范围(mappable endo% 0.1-44%) 富集轮数(1轮和2轮) 文库合并(单样本、2个合并、4个合并) 测序平台和试剂(NovaSeq 6000,Twist Ancient DNA试剂) 等位基因偏差检测(f4统计)

摘要

背景:溶液中杂交富集遗传标记是古基因组学研究中的一种首选方法,因为目标DNA通常高度碎片化并受到环境DNA污染,且获取个体间可比较的遗传数据具有挑战性。在此,我们使用来自不同地理来源、内源DNA含量从低到高(0.1%–44%)的古代人类样本的测序文库,对Twist Biosciences的商业化“Twist Ancient DNA”试剂进行基准测试。对每个文库,我们测试了一轮和两轮富集,并评估了与深度鸟枪法测序相比的性能。结果:我们发现“Twist Ancient DNA”检测对约120万个目标SNP提供了稳健的富集,且未引入可能干扰下游群体遗传学分析的等位基因偏差。此外,我们表明,对于内源DNA含量低于27%的文库,合并多达4个测序文库并进行两轮富集既可靠又具有成本效益。当内源含量高于38%时,建议最多进行一轮富集,以保证成本效益并保持文库复杂度。结论:总之,我们为人类古基因组学领域的研究人员提供了对不同测序和富集策略的优势与局限性的全面理解,我们的结果为优化实验方案提供了实用指导。

实验结论

提炼研究问题、关键发现与证据,快速把握文章的核心贡献。

研究问题
商业化的“Twist Ancient DNA”试剂在多种内源DNA含量和合并策略下进行一轮或两轮富集后,与深度鸟枪测序相比,是否能够提供高效、无偏且成本有效的古代人类SNP数据?
核心机制
Twist富集检测未引入显著的等位基因偏差,其效率受内源DNA含量影响,高内源含量下两轮富集可能导致偏好性再捕获扩增分子,降低文库复杂度。
主要证据
对24个古代人类样本进行深度鸟枪测序、一轮和两轮Twist富集,比较SNP产量、内源含量和成本;使用f4统计检测等位基因偏差,使用多种方法评估索引跳跃。
研究意义
为人类古基因组学研究人员提供了不同测序和富集策略的优势与局限性,并为优化实验方案提供了实用指导,强调合并文库和两轮富集在低内源含量样本中的成本效益。

研究路径

按研究推进顺序梳理实验设计、验证步骤与关键观察。

1

样本收集和DNA提取

获取用于后续文库制备和测序的古代人类DNA。

从伊比利亚、中美洲和东南亚的24个古代人类骨骼样本中提取DNA,使用针对降解DNA优化的方法和部分UDG处理。

2

文库制备和筛选

制备双索引双链DNA文库,并通过鸟枪测序筛选评估文库质量和内源DNA含量。

生成部分UDG处理的双索引双链DNA文库,使用Qubit和TapeStation进行质量控制,对文库进行浅层鸟枪测序以确定内源DNA含量和复杂度。

3

文库合并和过扩增

合并文库以达到富集所需的输入DNA量,并降低每文库成本。

根据总DNA质量将文库合并为六组:6个单文库反应、3个含2个低内源文库的合并、3个含4个高内源文库的合并,然后进行过扩增以达到每个反应1000 ng DNA。

4

目标富集(一轮和两轮)

使用Twist Ancient DNA试剂对靶标SNP进行富集,并比较一轮与两轮富集的效率和成本。

按照制造商方案,对每个合并或单文库进行一轮和两轮富集,包括杂交和捕获,并进行后富集PCR。

5

测序和数据处理

对富集和鸟枪文库进行测序,并处理原始数据以生成遗传数据。

在NovaSeq 6000平台进行测序,数据使用nf-core/eager流程处理,包括合并reads、比对到参考基因组、去除重复和进行假单倍体变异调用。

6

评估富集效率和成本效益

比较不同方法(鸟枪、一轮富集、两轮富集)的SNP产量、内源含量和成本。

计算不同测序深度下的SNP数量、内源DNA百分比(多种定义)和每SNP成本,并使用统计模型(t检验、广义线性模型)比较方法间的差异。

7

评估等位基因偏差

检测Twist富集是否引入等位基因偏差。

计算f4统计量(形式:f4(Mbuti, Pop1.MethodA; Pop2.MethodA, Pop2.MethodB)),比较不同方法间的差异。

8

评估索引跳跃和污染

评估合并文库中是否存在导致交叉污染的索引跳跃。

计算单索引跳跃率,使用多种方法(HaploCheck、ANGSD、HapConX、DICE)估计线粒体、X染色体和常染色体污染。

研究方法

按研究目的归类文中使用的方法,便于定位所需技术。

产品清单

实验环节名称品牌货号
QubitThermo Fisher--
TapeStationAgilent--
KAPA HiFi HotStart 预混液Roche--
AmpureXP 磁珠----
Twist Ancient DNA 试剂Twist Biosciences--
NovaSeq 6000 系统Illumina--

关键环节

汇总复现实验时建议重点确认的条件及原文阅读提示。

环节核对要点
DNA提取和文库制备
骨粉用量,样本来源,提取方法,UDG处理
阅读提示:Methods: DNA extraction and library preparation
文库扩增
PCR反应成分(聚合酶、引物、循环数)
阅读提示:Methods: DNA extraction and library preparation
富集
输入DNA量,富集轮数,合并策略
阅读提示:Methods: Library enrichment
测序
测序平台,读取长度,流动槽类型
阅读提示:Methods: Sequencing
数据分析
比对参数,重复去除,变异调用,内源含量定义
阅读提示:Methods: Data processing, Tables 1 and 2