世界首批农民的基因组起源

The genomic origins of the world's first farmers.

作者信息Nina Marchi, Laura Winkelbach, Ilektra Schulz, Maxime Brami, Zuzana Hofmanová, Jens Blöcher, Carlos S Reyna-Blanco, Yoan Diekmann, Alexandre Thiéry, Adamandia Kapopoulou, Vivian Link, Valérie Piuz, Susanne Kreutzer, Sylwia M Figarska, Elissavet Ganiatsou, Albert Pukaj, Travis J Struck, Ryan N Gutenkunst, Necmi Karul, Fokke Gerritsen, Joachim Pechtl, Joris Peters, Andrea Zeeb-Lanz, Eva Lenneis, Maria Teschler-Nicola, Sevasti Triantaphyllou, Sofija Stefanović, Christina Papageorgopoulou, Daniel Wegmann, Joachim Burger, Laurent Excoffier
PMID35561686
期刊Cell
发布时间2022-05-26
DOI10.1016/j.cell.2022.04.008

实验完整度

高

包含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分析和进化建模已实现对西南亚和欧洲新石器时代人群的祖先溯源,以解析世界上首批定居农民的遗传起源。

实验结论

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

研究问题
世界首批农民的遗传起源是什么?欧洲和西南亚新石器时代人群的祖先来源及分化过程如何?
核心机制
欧洲早期农民起源于西南亚人群与西部狩猎采集者在大约14.2 kya和12.9 kya的两次混合,随后经历强烈遗传漂变(有效群体大小约620人)而分化;欧洲狩猎采集者在末次盛冰期(LGM)经历严重瓶颈导致遗传多样性降低。
主要证据
对25个古基因组(15个新测序)进行人口基因组学建模,结合MDS、ADMIXTURE、f-statistics和fastsimcoal2分析,模拟数据重现了观察到的群体关系。
研究意义
阐明了新石器转型背后的人口统计学过程,并提供了晚更新世和早全新世西南亚和欧洲人群历史的空间解释;表明首批农民并非单一文化起源,而是经历了与欧洲狩猎采集者的重复混合。

研究路径

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

1

古DNA样本采集与测序

获取高分辨率古基因组以用于群体遗传学分析

从考古遗址采集骨骼样本,提取DNA,构建双索引文库,进行全基因组测序(Illumina HiSeq3000)

2

生物信息学数据处理与变异检测

获得高质量的基因组变异数据用于下游分析

使用ATLAS管道进行 reads 修剪、比对(bwa)、局部重比对(GATK)、基因型调用(ATLAS)、过滤和数据集组装

3

群体遗传学分析

描述群体结构和遗传关系

进行MDS、ADMIXTURE、f-statistics、杂合度分析、ROH分析等

4

人口统计学建模

推断群体分化、混合和有效群体大小等历史人口事件

使用fastsimcoal2进行最大似然估计,拟合多群体SFS,并通过模拟验证模型

5

模型验证与模拟数据比较

验证推断模型的可靠性

将模拟数据与观察数据进行比较(MDS、ADMIXTURE、f-statistics)

研究方法

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

产品清单

实验环节名称品牌货号
EDTAAmbion/Applied Biosystems, Life TechnologiesAM9262
蛋白酶KRoche Diagnostics, Mannheim, Germany3115828001
月桂酰肌氨酸钠盐Merck Millipore, Merck KGaA, Darmstadt, Germany428010
MinElute PCR纯化试剂盒QIAGEN, Hilden, Germany28006
Amicon Ultra-4离心过滤装置,30kDaMerck Millipore, Darmstadt, GermanyUFC803096
USER酶New England Biolabs GmbHM5505L
NEBNext末端修复酶混合物New England Biolabs GmbHE6050L
T4 DNA连接酶Thermo ScientificEL0011
AccuPrime Pfx SuperMixInvitrogen12344040
Agencourt AMPure XP磁珠Beckmann CoulterA63880
MSB Spin PCRapaceInvitek Molecular GmbH, Berlin, Germany1020220400
Qubit荧光定量和dsDNA HS检测试剂盒InvitrogenQ32854
Agilent 2100 Expert生物分析仪系统和高灵敏度DNA分析试剂盒Agilent Technologies5067-4626
羟基磷灰石Sigma-Aldrich21223
特殊刚玉(EW60/250μ)Harnisch+Rieth, Winterbach, Germany75250
特殊刚玉(30B/50μ)Harnisch+Rieth, Winterbach, Germany75308

关键环节

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

环节核对要点
古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”小节