样本收集与DNA提取
为后续遗传分析获取古代DNA。
从Sakhtysh遗址的43个个体中采集84块骨骼样本（包括岩骨和牙齿），进行DNA提取和部分尿嘧啶-DNA糖基化酶文库构建。
关于丁香通
公司信息
个人用户
企业机构
无忧采购 轻松科研
- 移动端
古代基因组学研究揭示东方狩猎采集者HLA多样性较低
Ancient genomics study reveals low HLA diversity in eastern hunter-gatherers
实验完整度
高研究包含多个独立证据层级：全基因组群体遗传分析（PCA、f-statistics、ADMIXTURE）、亲缘关系分析、病原体筛查及HLA基因分型，并包含功能验证（中性检验、选择信号分析）。
主要模型
Sakhtysh古代人群（Lyalovo–Volosovo文化） 外部比较人群（新石器早期农民、晚期农民、现代欧洲人群）
重点核对
样本纳入标准：1240K面板基因分型SNP数>20000 分组：按文化时期分组（Lyalovo、Early Volosovo、Volosovo、Late Volosovo、Volosovo-transitional） HLA分型方法：TARGT管线，OptiType确认 病原体筛查方法：MALTv0.4.1，最小一致性90%，数据库含27,730细菌和10,543病毒基因组
摘要
背景
东方狩猎采集者（EHG）为许多现代欧洲人贡献了不同比例的遗传祖先。然而，这种贡献如何塑造人类白细胞抗原（HLA）基因多样性仍未被探索。俄罗斯西部的Sakhtysh是最大的EHG相关埋葬群之一。该遗址的持续使用为研究狩猎采集者群体的遗传史、免疫基因多样性和病原体暴露提供了难得的机会。
结果
来自36个个体的全基因组数据（跨越Lyalovo–Volosovo文化传统，6800–4800 cal BP）将Sakhtysh群体置于EHG聚类中，与西北俄罗斯和东波罗的海的同时代群体具有高度的遗传亲和性。这种遗传组成在两千多年中保持稳定，没有与邻近的西方狩猎采集者或农民群体混合的证据。我们报告了27个Sakhtysh个体的HLA基因型，代表了古代狩猎采集者中首批真正的HLA分型。Sakhtysh群体的HLA多样性较低，可能是由于其较小的种群规模。特定等位基因如B*27:05以高频率（约50%）被观察到，尽管在现代人群中相对罕见，表明选择压力发生了变化。病原体筛查确定了有限的病原体谱，包括可能的猪红斑丹毒丝菌人畜共患病以及乙型肝炎病毒和细小病毒B19的慢性感染。
结论
Sakhtysh社区在文化转变中表现出遗传连续性。该群体表现出有限的HLA多样性和某些等位基因的意外高频率。这些模式可能受到低病原体多样性和对环境因素特定适应的可能影响。
实验结论
提炼研究问题、关键发现与证据，快速把握文章的核心贡献。
- 研究问题
- 东方狩猎采集者（EHG）对现代欧洲人的遗传贡献如何影响HLA基因多样性？
- 核心机制
- 文中未明确说明
- 主要证据
- 基于Sakhtysh遗址36个个体的全基因组数据和27个个体的HLA基因型，通过PCA、f-statistics、ADMIXTURE等群体遗传分析，发现该群体遗传连续且属于EHG聚类；通过Shannon多样性指数、FST和选择中性检验，发现HLA多样性低且存在高频率等位基因。
- 研究意义
- 该研究提供了古代狩猎采集者首批真正的HLA基因型，表明EHG群体HLA多样性低且特定等位基因频率高，可能受低病原体多样性和环境适应影响；强调需要更多史前狩猎采集者HLA数据以理解过去免疫遗传变异。
研究路径
按研究推进顺序梳理实验设计、验证步骤与关键观察。
全基因组测序与群体遗传分析
评估Sakhtysh人群的遗传亲和性、连续性和混合情况。
对43个个体进行全基因组鸟枪法测序，结合已发表的4个个体，共36个个体用于群体遗传分析。通过PCA、outgroup f3、f4统计、ADMIXTURE和qpAdm建模分析遗传结构和连续性。
亲缘关系分析
检测个体间的生物学亲缘关系。
使用READv2对合并后的BAM文件进行亲缘关系分析，限制到二级亲缘。
病原体筛查
检测古代样本中是否存在细菌或病毒DNA。
使用MALT v0.4.1对43个样本的宏基因组测序读段进行筛查，阳性样本通过参考基因组比对验证。
HLA靶向富集与基因分型
获取古代样本的HLA基因型，评估HLA多样性。
对31个有岩骨的样本进行HLA区域DNA靶向富集，27个样本成功进行HLA基因分型。使用TARGT管线确定HLA-A,B,C,DPB1,DQB1,DRB1等位基因，并用OptiType确认I类基因。
研究方法
按研究目的归类文中使用的方法，便于定位所需技术。
产品清单
| 实验环节 | 名称 | 品牌 | 货号 |
|---|---|---|---|
| 全基因组测序 | Illumina HiSeq 6000 | -- | -- |
| DNA提取和文库构建 | 部分尿嘧啶-DNA糖基化酶文库转化 | -- | -- |
| 病原体筛查 | MALT v0.4.1 | -- | -- |
| MEGAN6 v6.25.9 | -- | -- | |
| 比对 | BWA 0.7.15版 | -- | -- |
| 处理测序读数 | ClipAndMerge 1.7.7版 | -- | -- |
| 重复去除 | DeDup 0.12.1版 | -- | -- |
| 损伤分析 | DamageProfiler 1.1版 | -- | -- |
| 污染评估 | Schmutzi 1.5.5.5版 | -- | -- |
| ANGSD 0.935版 | -- | -- | |
| 线粒体单倍群 | HaploGrep2 | -- | -- |
| Y染色体单倍群 | yHaplo | -- | -- |
| 亲缘关系分析 | READv2 | -- | -- |
| 基因分型 | SequenceTools 1.2.2版 | -- | -- |
| 群体遗传分析 | smartpca | -- | -- |
| ADMIXTURE v1.3.0 | -- | -- | |
| qp3Pop、qpDstat、qpAdm | -- | -- | |
| HLA分析 | TARGT管线 | -- | -- |
| OptiType | -- | -- | |
| 统计检验 | Pypop | -- | -- |
| FST计算 | hierfstat包 | -- | -- |
关键环节
汇总复现实验时建议重点确认的条件及原文阅读提示。
| 环节 | 核对要点 |
|---|---|
| 样本采集与DNA提取 | 样本数量和类型（岩骨和牙齿），个体数（43），DNA提取和文库转化方法（部分UDG处理） 阅读提示：Methods: Sampling; DNA extraction and library preparation |
| 全基因组测序与数据处理 | 测序平台（Illumina HiSeq 6000）、读长（2x100bp）、比对参数（BWA -n 0.01 -l 300）、重复去除（DeDup） 阅读提示：Methods: Processing of reads and mapping to the human reference; Contamination assessment |
| 群体遗传分析 | 分析面板（1240K和HO）、样本纳入标准（>20000 SNPs）、分组（文化时期）、统计方法（PCA、f3/f4、ADMIXTURE、qpAdm） 阅读提示：Methods: Genotyping and merging; PCA; Unsupervised admixture; F-statistics |
| 病原体筛查 | 筛查工具（MALT v0.4.1）、数据库（细菌和病毒基因组数）、比对参数（最少一致性90%），阳性标准（>1000 reads） 阅读提示：Methods: Metagenomic pathogen screening; Results: Pathogen screening |
| HLA基因分型 | 富集样本数（31）、成功分型数（27）、分型工具（TARGT和OptiType）、分析分辨率（第一和第二字段） 阅读提示：Methods: HLA analysis; Results: HLA diversity |