古DNA样本采集与测序
获取高分辨率古基因组以用于群体遗传学分析
从考古遗址采集骨骼样本,提取DNA,构建双索引文库,进行全基因组测序(Illumina HiSeq3000)
The genomic origins of the world's first farmers.
包含15个新测序古基因组(>10x)和10个已发表古基因组,以及现代参考基因组;进行古DNA提取、测序、生物信息学分析和群体遗传学建模(MDS、ADMIXTURE、f-statistics、MSMC2、fastsimcoal2)等多种验证。
15个新测序的古人类个体(来自欧洲和安纳托利亚的早期农民和狩猎采集者) 10个已发表的古人类基因组 77个现代人基因组(SGDP)
古DNA的污染评估:通过blank controls和ContamMix/ANGSD估计污染率 样本的放射性碳定年(95.4%概率区间) 基因组覆盖度均>10X 群体遗传学分析(MDS、ADMIXTURE、f-statistics)验证群体结构 人口统计学建模(fastsimcoal2)的似然比较和参数估计
提炼研究问题、关键发现与证据,快速把握文章的核心贡献。
按研究推进顺序梳理实验设计、验证步骤与关键观察。
获取高分辨率古基因组以用于群体遗传学分析
从考古遗址采集骨骼样本,提取DNA,构建双索引文库,进行全基因组测序(Illumina HiSeq3000)
获得高质量的基因组变异数据用于下游分析
使用ATLAS管道进行 reads 修剪、比对(bwa)、局部重比对(GATK)、基因型调用(ATLAS)、过滤和数据集组装
描述群体结构和遗传关系
进行MDS、ADMIXTURE、f-statistics、杂合度分析、ROH分析等
推断群体分化、混合和有效群体大小等历史人口事件
使用fastsimcoal2进行最大似然估计,拟合多群体SFS,并通过模拟验证模型
验证推断模型的可靠性
将模拟数据与观察数据进行比较(MDS、ADMIXTURE、f-statistics)
按研究目的归类文中使用的方法,便于定位所需技术。
| 实验环节 | 名称 | 品牌 | 货号 |
|---|---|---|---|
| DNA提取 | EDTA | Ambion/Applied Biosystems, Life Technologies | AM9262 |
| 蛋白酶K | Roche Diagnostics, Mannheim, Germany | 3115828001 | |
| 月桂酰肌氨酸钠盐 | Merck Millipore, Merck KGaA, Darmstadt, Germany | 428010 | |
| MinElute PCR纯化试剂盒 | QIAGEN, Hilden, Germany | 28006 | |
| Amicon Ultra-4离心过滤装置,30kDa | Merck Millipore, Darmstadt, Germany | UFC803096 | |
| 文库构建 | USER酶 | New England Biolabs GmbH | M5505L |
| NEBNext末端修复酶混合物 | New England Biolabs GmbH | E6050L | |
| T4 DNA连接酶 | Thermo Scientific | EL0011 | |
| AccuPrime Pfx SuperMix | Invitrogen | 12344040 | |
| Agencourt AMPure XP磁珠 | Beckmann Coulter | A63880 | |
| MSB Spin PCRapace | Invitek Molecular GmbH, Berlin, Germany | 1020220400 | |
| 文库定量 | Qubit荧光定量和dsDNA HS检测试剂盒 | Invitrogen | Q32854 |
| Agilent 2100 Expert生物分析仪系统和高灵敏度DNA分析试剂盒 | Agilent Technologies | 5067-4626 | |
| 样本准备 | 羟基磷灰石 | Sigma-Aldrich | 21223 |
| 特殊刚玉(EW60/250μ) | Harnisch+Rieth, Winterbach, Germany | 75250 | |
| 特殊刚玉(30B/50μ) | Harnisch+Rieth, Winterbach, Germany | 75308 |
汇总复现实验时建议重点确认的条件及原文阅读提示。
| 环节 | 核对要点 |
|---|---|
| 古DNA样本测序 | 测序深度:平均>10X;测序平台:Illumina HiSeq3000;测序方式:单端100循环或双端150循环 阅读提示:见STAR Methods中“Whole-genome sequencing”小节 |
| 生物信息学分析 | 过滤标准:DP≥8,GQ≥30,等位基因平衡p≤0.1;参考基因组:hs37d5;变异检测软件:ATLAS 阅读提示:见STAR Methods中“Data filtering”小节 |
| 群体遗传学分析 | 中性位点集定义:重组率>1 cM/Mb,非CpG,A↔T和G↔C突变;SNP数量:473,834 阅读提示:见STAR Methods中“Dataset assembling”小节 |
| 人口统计学建模 | fastsimcoal2运行参数:50次独立运行,每次50个ECM循环,500,000次溯祖模拟;置信区间通过100次参数自举估计 阅读提示:见STAR Methods中“fastsimcoal2”小节 |
| 放射性碳定年 | 校准曲线:IntCal20;校正软件:OxCal 4.4.2;淡水水库效应校正方法:Cook et al. (2001) 和 Bonsall et al. (2015) 阅读提示:见STAR Methods中“Radiocarbon dating and stable isotope analysis”小节 |