本方案介绍了一系列最佳实践操作流程,用于从五个不同的骨骼部位上八个推荐的解剖学取样位置(特定骨骼元素上的具体位置)采集中世纪个体的骨粉样本(放射性碳测年校准后2-sigma范围约为公元1040-1400年)。
本方案介绍了一系列最佳实践操作流程,用于从五个不同的骨骼部位上八个推荐的解剖学取样位置(特定骨骼元素上的具体位置)采集中世纪个体的骨粉样本(放射性碳测年校准后2-sigma范围约为公元1040-1400年)。
本文所述方法旨在最大限度提高从古代考古遗存中回收人类DNA的可能性,同时减少所需取样材料的用量。该方法通过靶向解剖学取样部位实现,这些部位在全骨骼范围内DNA回收的比较分析中已被证实可获得最丰富的古代DNA(aDNA)。先前研究表明,这些实验方案能够最大限度提高从考古遗存中成功回收古代人类及病原体DNA的概率。Parker等人于2020年对来自德国佩森附近废弃中世纪定居点克拉考尔贝格(Krakauer Berg)墓地的11名个体、多个骨骼部位的aDNA保存状况进行了广泛调查,评估了DNA产量。该墓地经放射性碳测年确定为中世纪(约公元1040–1400年,校准后2-sigma范围)。研究发现,涵盖五个骨骼部位(pars petrosa、恒磨牙、胸椎、远端指骨和距骨)的这八个取样位点成功获得了高质量的古代人类DNA,其产量显著高于所有骨骼部位和个体的总体平均水平。所得DNA产量足以满足大多数常见下游群体遗传学分析的需求。我们的研究结果支持在绝大多数涉及从考古遗存中分析古代人类DNA的研究中优先选用这些解剖学取样部位。采用这些方法将有助于最小化对珍贵考古标本的破坏。
为获取和分析DNA而对古代人类遗骸进行取样本质上具有破坏性1,2,3,4。样本本身是珍贵的标本,应尽可能保持其形态学完整性。因此,必须优化取样方法,以避免对不可替代的材料造成不必要的破坏,并最大限度地提高实验成功的概率。目前的最佳实践技术主要基于少量研究,这些研究仅限于法医学调查5,6、并非以优化取样为主要目标的古代标本研究7,或使用非人类遗骸8或仅针对极少数解剖学取样部位的研究(本文中“解剖学取样部位”指从骨骼单元的特定区域钻取骨粉,用于后续DNA分析)9,10。本文所介绍的取样方案是在首次针对多个个体的多种骨骼单元开展的大规模系统性DNA保存研究中优化而成的11。所有样本均来自德国萨克森-安哈尔特州佩森附近废弃中世纪定居点克拉考尔贝格(Krakauer Berg)教堂墓地发掘出的11名个体的骨骼单元(详细的样本人口学信息见表1),因此,若用于此地理和时间范围之外的样本,可能需要进行适当调整。
| 个体 | 性别 | 死亡时估计年龄 | 14C 测年(公元,日历年代 2-西格玛) |
| KRA001 | 男性 | 25–35 | 1058–1219 |
| KRA002 | 女性 | 20–22 | 1227–1283 |
| KRA003 | 男性 | 25 | 1059–1223 |
| KRA004 | 男性 | 15 | 1284–1392 |
| KRA005 | 男性 | 10–12 | 1170–1258 |
| KRA006 | 女性 | 30–40 | 1218–1266 |
| KRA007 | 女性 | 25–30 | 1167–1251 |
| KRA008 | 男性 | 20 | 1301–1402 |
| KRA009 | 男性 | 未知 | 1158–1254 |
| KRA010 | 男性 | 25 | 1276–1383 |
| KRA011 | 女性 | 30–45 | 1040–1159 |
表1:11例样本个体的遗传学性别、考古学推断的死亡年龄估计值以及放射性碳测年(14C 校准后2-sigma范围)。 本表格改编自 Parker, C. 等人,202011。
这些方案可从五个骨骼部位(包括pars petrosa)的八个解剖采样位置相对简便且高效地制备骨粉,同时将实验室引入的DNA污染控制在较低水平。在这五个骨骼部位中,已确定其中四个骨骼部位上的七个解剖采样位置可作为破坏性采集岩骨锥体的可行替代方案11,12。这些位置包括恒磨牙的牙骨质、牙本质和牙髓腔;从胸椎椎体及上椎切迹处采集的皮质骨;从远节指(趾)骨远端膨大部底面及骨干获取的皮质骨;以及距骨外侧部分致密的皮质骨。尽管目前已有多种广泛应用的方法用于pars petrosa4,12,13,14、牙本质及牙髓腔1,2,15的取样,但有关从牙骨质16、椎体、下椎切迹和距骨成功制备骨粉的已发表方法却较难获取。因此,本文展示了优化后的采样方案,涵盖岩骨锥体(步骤3.1);成人磨牙的牙骨质(步骤3.2.1)、牙本质(步骤3.2.2)和牙髓(步骤3.2.3);椎体皮质骨(步骤3.3.1)和上椎弓(步骤3.3.2);远节指(趾)骨(步骤3.4)以及距骨(步骤3.5),以促进这些骨骼部位在古DNA和法医学研究中的更广泛应用。
本文所展示的所有研究均遵循德国耶拿马克斯·普朗克人类历史科学研究所关于处理古代人类遗骸的指导方针。在执行本方案的任何步骤之前,请确保遵守您所在地区所有相关的本地/州/联邦伦理要求,包括获取开展科学研究以及对人类遗骸进行破坏性取样的许可。所有操作和化学品储存均应按照各机构的安全指南进行。
1. 样品处理前的注意事项
2. 预处理
3. 骨粉制备
注意:以下方案旨在用于遵循 Dabney 等人 2019 年方案26 进行的 DNA 提取。

图1:包含pars petrosa的颞骨。(A)切割前样本,显示岩骨锥体和sulcus petrosa的位置。(B)切割后岩骨部分,突出显示需钻取的致密区域。请点击此处查看该图的放大版本。

图 2:永久性磨牙采样前处理。(A)采样前的预处理磨牙,显示牙冠、牙骨质(牙根的淡黄色层)以及牙骨质-釉质交界处的切割位置。(B)采集牙骨质后的同一颗磨牙,显示在牙骨质-釉质交界处的切口位置。(C)切割并采样后的磨牙,显示牙冠内牙髓腔和牙本质的解剖学采样位置。请点击此处查看该图的放大版本。

图3:胸椎椎体及上位椎弓皮质骨的解剖采样位置。 请点击此处查看此图的放大版本。

图4:远节指骨,显示骨干及远端小梁下侧致密皮质骨的分布位置。请点击此处查看该图的放大版本。

图 5:距骨皮质骨回收的取样区域。请点击此处查看该图的放大版本。
注意:距骨的皮质骨极少(仅为一层薄的外层)。取材时不仅应采集表面,还应包括其下方的致密松质骨层。
在另一项研究中11从11名个体每个解剖采样部位获得的骨粉中,采用针对钙化组织中短片段DNA优化的标准DNA提取方案提取DNA2单链文库随后被构建28 在 HiSeq 4000 平台(75 bp 双端测序)上进行测序,每个样本测序深度约为 20,000,000 条读长。随后使用 EAGER 分析流程对获得的序列数据进行内源性人类 DNA 含量评估。29 (BWA 设置:种子长度为 32,错配罚分为 0.1,比对质量过滤阈值为 37)。所有代表性结果均采用与 Parker 等 2020 年研究相同的指标进行报告。11 为保持一致性。文库来自样品粉末部分的 岩部 平均而言,产生的内源性DNA高于所调查的其他23个解剖采样部位中的任何一个(图6A-B)。本方案中提供的其余七个解剖采样部位(恒磨牙的牙骨质、牙髓腔初次冲洗液及牙本质;胸椎椎体和上椎弓的皮质骨;远端指骨远端丛状突的皮质骨;以及距骨颈部的皮质骨)获得了次高的产量(这些解剖采样部位之间的产量差异无统计学意义; 图6A-B;补充文件1: 内源性DNA预建库)。这些替代取样部位均能稳定产生足以满足标准群体遗传学分析(如线粒体分析和单核苷酸多态性(SNP)分析)所需的DNA产量。来自所有解剖取样部位文库的重复率均较低(簇因子 < 1.2(平均值),计算方法为所有比对读段数与唯一比对读段数的比值, 表2;补充文件1:ClusterFactor),表明所筛选的所有文库均具有极高的复杂度。类似地,外源性人类DNA污染的平均估计值较低,平均为 < 2%(雄性样本中X染色体污染,n = 7,由ANGSD报告)30 除上椎弓外的所有解剖采样位置(平均估计污染率为2.11%,其中一份样本因属异常值被剔除;KRA005:19.52%,见 表2;补充文件1:污染)。平均片段长度(过滤去除所有读段后) < 30 bp)在从牙髓腔和牙本质采集的材料中最低,其他解剖采样部位之间无显著差异(分别为55.14 bp和60.22 bp,相较于平均中位数62.87,成对p值 < 0.019, 表2;补充文件1:平均片段长度)。此外,牙齿和胸椎各自包含多个解剖学取样位点,这些位点均观察到较高的内源性DNA得率,因此特别适合作为替代材料 岩部.

图6:所有筛选样本的人类DNA含量。 黑色线条表示总体均值,红色线条表示中位数(实线:人类DNA比例,虚线:每百万条测序读段中比对到人类基因组的唯一读段数)。在所有分析中,平均人类DNA比例高于总体均值(8.16%)的各个解剖采样位置均以颜色标注。(A)比对到hg19参考基因组的读段比例。蓝色虚线表示根据分析流程的比对参数所计算的理论最大值(使用Gargammel31模拟从hg19参考基因组随机分布的5,000,000条读段,并加入模拟损伤)。对于平均人类DNA比例高于总体均值的样本,图中报告了其个体均值(黑色X)和中位数(红色圆圈)。置信区间表示排除统计异常值后的上下限。(B)每百万测序读段(75 bp双端测序)中唯一比对到hg19参考基因组的读段数量。置信区间表示排除统计异常值后的上下限。本图改编自Parker, C. 等,202011。请点击此处查看此图的放大版本。
表2:所有解剖采样部位的平均重复水平(比对上的读段数/唯一读段数)、平均片段长度与中位片段长度,以及X染色体污染估计值。 误差以均值的标准误表示。本表改编自 Parker, C. 等,202011。
| 采样地点 | 平均重复因子(#比对上的读段数 / #唯一比对上的读段数) | 平均片段长度(bp) | X染色体污染的平均估计比例 |
| 岩骨锥体 | 1.188 ± 0.006 | 65.40 ± 1.36 | 0.000 ± 0.003 |
| 牙骨质 | 1.197 ± 0.028 | 67.28 ± 1.76 | 0.011 ± 0.003 |
| 牙本质 | 1.188 ± 0.061 | 60.22 ± 2.37 | 0.002 ± 0.007 |
| 牙髓 | 1.179 ± 0.024 | 55.14 ± 2.90 | 0.013 ± 0.006 |
| 远节指骨 | 1.191 ± 0.049 | 65.95 ± 1.08 | 0.013 ± 0.005 |
| 椎体 | 1.194 ± 0.037 | 66.14 ± 1.03 | 0.008 ± 0.003 |
| 上椎弓 | 1.19 ± 0.017 | 63.02 ± 1.23 | 0.021 ± 0.009* |
| 距骨 | 1.198 ± 0.010 | 68.20 ± 1.24 | 0.011 ± 0.003 |
| *样本 KRA005 因偏离值被移除,其数值为 0.1952 |
代码可用性
本研究中使用的所有分析程序和R模块均可从其相应作者处免费获取。所有自定义的R代码可根据要求提供。
数据可用性
用于计算代表性结果的所有原始数据均可在欧洲核苷酸档案库(European Nucleotide Archive ENA)中免费获取(登录号 PRJ-EB36983),或见Parker, C. 等人11的补充材料。
补充文件 1。 请点击此处下载该文件。
目前,在古人类群体遗传学研究中,优先选择从pars petrosa(步骤2.1)取样。然而,pars petrosa往往难以获取,因为它在多种骨骼评估中具有重要价值(例如,群体历史32、胎儿死亡年龄估算33以及性别鉴定34);此外,历史上为DNA分析而取样pars petrosa可能造成严重的破坏性3,4(包括本文所介绍的方案,尽管目前已有多种新型微创取样方案13,14被广泛采用,以缓解这一问题)。更为复杂的是,直到最近之前,学界尚未开展大规模、系统性的人体骨骼各部位DNA提取效果研究11,因此当无法获取岩骨锥体时,确定合适的取样策略极具挑战性。
本文提供的方案通过一系列优化的DNA取样操作流程,有助于缓解这一挑战,这些流程适用于考古学和法医学骨骼遗骸的DNA取样,包括pars petrosa以及另外四种骨骼元素上的七个替代解剖取样部位。所包含的关键步骤均旨在最大限度地减少因取样效率低下(步骤2.1.6和3.2.1.3)或在钻孔/切割过程中样品过热(步骤3.1.6)而导致的DNA丢失或损伤。此外,本方案中多次指出,可能需要修改或省略预处理步骤,以确保在高度降解的样本中获得最佳效果。还应注意的是,即使在本文所列的选定骨骼部位中,仍存在多种可能的替代取样技术(特别是针对pars petrosa13,14),并且本文提出的尚未被充分开发的解剖取样部位仍有很大的进一步优化空间(例如, 距骨:步骤2.5,以及椎骨:步骤2.3)。
还需要注意的是,这些实验方案的设计与测试均基于形态保存良好(形态学保存状况优良)的古代青少年至成年个体遗骸,旨在进行内源性人类DNA分析。所呈现的结果可能不适用于降解程度更高的材料、其他保存环境、婴儿遗骸、非人类遗骸,或病原体与微生物组的研究,因为这些方案在更多情境下的应用仍需进一步探索。此外,本文所介绍的替代骨骼样本(牙齿、椎骨、远端指骨和距骨)在混合遗骸中可能难以准确归属到单一个体,因此有必要从多个骨骼部位取样,以确保样本来源的单一性。尽管存在这些局限性,广泛公开这些实验方案有助于缓解样本选择与处理过程中存在的部分异质性问题,为未来各类人类遗骸的古DNA/法医学研究提供一个通用且经过定量优化的框架。
作者声明无任何利益冲突。
作者谨向马克斯·普朗克人类历史科学研究所的实验室工作人员致以诚挚感谢,感谢他们在本实验方案的开发与实施过程中提供的帮助。若无古多·布兰特博士、伊丽莎白·尼尔森博士、安特耶·维塞戈特和弗朗西斯卡·阿隆的贡献与辛勤工作,本研究将无法完成。本研究由马克斯·普朗克学会、欧洲研究理事会(ERC)在欧盟“地平线2020”研究与创新计划下资助,资助协议编号分别为771234 - PALEoRIDER(WH, ABR)和起始资助编号805268 CoDisEASe(授予KIB)。
| 姓名 | 公司 | 目录编号 | 评论 |
|---|---|---|---|
| #16 牙科钻头 | NTI | H1-016-HP | 示例钻头 |
| 0.6 mm 带锯条 | Fisher Scientific | 50-949-097 | 珠宝锯用锯条 |
| 22 mm 金刚石切割轮 | Kahla | SKU 806 104 358 514 220 | Dremel 切割附件 |
| 市售漂白剂 | Fisher Scientific | NC1818018 | |
| Control Company Ultra-Clean Supreme 铝箔 | Fisher Scientific | 15-078-29X | |
| DNA LoBind 管(2 mL) | Eppendorf | 22431048 | |
| Dremel 225-01 柔性轴附件 | Dremel | 225-01 | Dremel 柔性延长杆 |
| Dremel 4300 旋转工具 | Dremel | 4300 | 示例电钻 |
| Dremel 夹头和螺母套件 | Dremel | 4485 | 用于各种 Dremel 工具附件/钻头的转接器 |
| Eagle 33 加仑红色生物危害废物袋 | Fisher Scientific | 17-988-501 | |
| Eppendorf DNA LoBind 2 mL 微量离心管 | Fisher Scientific | 13-698-792 | |
| 乙醇(分子生物学级) | Millipore Sigma | 1.08543 | |
| FDA 批准的 2 级外科口罩 | Fisher Scientific | 50-206-0397 | 个人防护装备 |
| Fisherbrand 舒适型腈手套 | Fisher Scientific | 19-041-171X | 个人防护装备 |
| Fisherbrand 安全护目镜 | Fisher Scientific | 19-130-208X | 个人防护装备 |
| Granger 固定式台虎钳 | Fisher Scientific | NC1336173 | 台式虎钳 |
| Invitrogen UltraPure 无 DNase/RNase 蒸馏水 | Fisher Scientific | 10-977-023 | |
| 珠宝锯 | Fisher Scientific | 50-949-231 | |
| Kimwipes 纸巾 | Sigma-Aldritch | Z188956 | |
| Labconco Purifier Logic 生物安全柜 | Fisher Scientific | 30-368-1101 | |
| LookOut DNA Erase | Millipore Sigma | L9042-1L | |
| 中号称量舟 | Heathrow Scientific | HS120223 | |
| MSC 10 件套钳子/夹具组 | Fisher Scientific | 50-129-5352 | 杂项 夹具/虎钳,用于在钻孔/切割时牢固固定样品 |
| Sartorius Quintix 半微量天平 | Fisher Scientific | 14-560-019 | 封闭式天平 |
| 带连体头罩的 Tyvek 防护服 | Fisher Scientific | 01-361-7X | 个人防护装备 |
| 称量纸 | Heathrow Scientific | HS120116 |