IF50.5｜古DNA揭示斯拉夫语系扩散伴随大规模东欧人群迁徙

公元6-8世纪，东欧人群大规模西迁，在德国东部、波兰、克罗地亚替换了超过80%的当地基因库。

这篇按生信解读号的读法拆：先看背景和问题，再看作者怎么做、做出了哪几条硬发现，最后把数据和代码收走。

研究背景

公元一千年下半叶，中欧与东欧经历了深刻的文化与政治变革，传统上认为与斯拉夫人的出现有关。考古学上，这一时期出现了以布拉格-科尔恰克文化为代表的、高度相似的简单物质文化，包括小型半地穴式房屋、火葬墓、手工素面陶器等。历史文献记载了斯拉夫群体在6世纪后从多瑙河下游向北、向西扩散。然而，关于这一文化变迁是源于大规模人口迁徙、当地人群的文化同化，还是两者结合，学界长期存在争议。此前基于现代和少量古代DNA的研究已发现东欧向巴尔干北部的基因流，但也有研究主张波兰地区存在人群连续性，因此斯拉夫扩张的规模、时序及其与物质文化的关联仍不清楚。

科学问题

卡在哪里： 此前古代DNA数据稀疏，尤其因为斯拉夫早期普遍实行火葬，缺乏可直接测序的遗骸。波兰地区的研究甚至提出本地起源说，认为斯拉夫基因库自青铜时代以来连续。因此，缺乏跨区域、高分辨率的古代基因组数据来直接检验大规模迁徙假说，也无法量化基因替换程度或评估性别偏向。

本文要回答： 本文旨在利用跨中欧、东欧的古代基因组时间断面，直接检验斯拉夫扩张是否伴随大规模人口迁徙，并量化基因替换比例、追溯来源地、评估社会结构变化。

技术路线

作者主要用了：全基因组捕获测序、PCA、F4统计、qpAdm、ADMIXTURE、ancIBD、DATES、MOBEST。

作者选取了三个研究断面：德国东部Elbe-Saale地区、西北巴尔干、波兰-西北乌克兰，并补充波罗的海和俄罗斯伏尔加-奥卡地区作为参考。对555个古代个体进行全基因组捕获测序，结合已发表数据共1840个古代个体。首先通过PCA和F统计量观察遗传结构变化，然后利用监督式ADMIXTURE分解祖先成分，用qpAdm量化SP祖先比例，用ancIBD检测IBD共享以追溯共同来源，用DATES估计混合时间，最后用MOBEST推断地理起源。

图 1　Fig. 1。图1展示了三个研究断面（西北巴尔干、德国东部、波兰-西北乌克兰）及参考区域（波罗的海、伏尔加-奥卡）的样本分布，共1840个古代个体。新测序550个，已发表723个。符号大小代表每个遗址的个体数。图b显示1243个新测序个体的时间序列，覆盖从青铜时代到中世纪，重点在迁移期和斯拉夫时期。

核心亮点

亮点 1｜SP人群遗传结构剧变

在德国东部、波兰和克罗地亚，公元600年前后的样本在PCA上从与北欧或南欧人群聚类，转变为与现今斯拉夫语人群（如波兰人、白俄罗斯人）紧密聚类。ADMIXTURE分析显示，东北欧祖先成分（BAL）在迁移期仅占5-7%，但在斯拉夫时期跃升至47-65%，成为最大成分。qpAdm估算显示，西北巴尔干、德国东部、波兰-西北乌克兰分别有82%、83%、93%的当地基因库被东欧移民替换。

图 2　Fig. 2。图2a是基于10528个现代欧洲人的PCA，古代样本投影其上。大符号代表各现代人群的平均PC1和PC2坐标。图b显示基于前10个主成分的马氏距离层次聚类，圆圈代表罗马时期和迁移期样本，方块代表斯拉夫时期样本。可见SP样本与MP样本明显分离，SP样本聚集在现代斯拉夫语人群附近。

亮点 2｜SP人群共享近期共同祖先

ancIBD分析发现，来自克罗地亚、德国东部和波兰-乌克兰的SP个体之间共享大量IBD片段，包括许多大于16cM的长片段，表明他们源自一个共同的、在几代前才向西、向南扩张的源人群。相比之下，他们与当地前斯拉夫人群几乎不共享IBD。基于IBD网络的层次聚类识别出一个包含大多数SP个体的大社区，内部又分为喀尔巴阡山以北和以南两个亚群，可能对应两波不同的扩张或不同的本地融合模式。

图 3　Fig. 3。图3a显示三个研究区域对东欧亲和力的时间轨迹，使用F4统计量（如F4(Han Chinese, test; Italy, Poland)）。误差线为2倍标准差。可见SP时期亲和力急剧上升。图b是监督式ADMIXTURE建模，将1344个古代和现代样本分解为12个祖先成分。BAL成分在SP时期显著增加，成为主导。

亮点 3｜SP祖先形成于约公元前1000年

DATES估计SP基因池的形成混合事件发生在约公元前1000年（Niederwünsch：972±250 BCE；Poland_EMA：906±362 BCE）。这一时间与波罗的-斯拉夫语系分化的贝叶斯语言年代学估计（约公元前2000-1000年）部分重叠。qpAdm建模显示，SP个体可被建模为约71%波罗的海青铜时代相关祖先和约29%早期欧洲农民（EEF）富集祖先的混合。MOBEST推断SP个体的地理起源最可能位于白俄罗斯南部和乌克兰北部。

图 4　Fig. 4。图4a显示MOBEST推断的SP个体祖先地理起源概率面，预测时间设为1950年前。平均概率面指向白俄罗斯南部和乌克兰北部。图b显示DATES估计的混合时间与语言分化时间的比较。图c显示qpAdm估算的SP祖先替换比例，在德国东部、波兰、克罗地亚均超过80%，且无性别偏向。

亮点 4｜无性别偏向的混合

通过比较常染色体和X染色体上的SP祖先比例，作者未发现任何SP人群（德国、克罗地亚、波兰、俄罗斯）存在性别偏向的混合（|z|<2）。这与迁移期德国东部发现的南欧相关基因流显著女性偏向形成对比。这表明斯拉夫扩张过程中，男性和女性移民对基因库的贡献大致相等，暗示家庭或群体性迁徙而非以男性为主的征服。

亮点 5｜SP社会结构转向父系聚居

在德国东部，SP墓地的遗传亲缘关系显著增强：平均每个个体有6.41个近亲，是迁移期的近6倍。大型父系谱系常见，发现一对夫妇有7个后代，其中4个存活至生育年龄。女性在墓地内的近亲数显著少于男性，但跨墓地共享更多IBD，表明女性外婚和父系居住制。SP墓地中遗传距离与空间距离显著相关（Mantel检验P=0.0001），说明墓地围绕父系亲属群规划。

生信可带走

这一块是给做分析的人用的，不是科普点缀。

• 方法栈： 全基因组捕获测序、PCA、F4统计、qpAdm、ADMIXTURE、ancIBD、DATES、MOBEST
• 公开数据： 原文未给出公开组学登录号
• 代码： 原文未给出公开 GitHub/GitLab 仓库
• 谁该点开原文： 关注群体遗传学、历史语言学、考古学交叉研究的学者，以及做古代DNA数据分析的生信人员。

带走一句

大规模古代基因组时间断面结合IBD共享和祖先成分建模，可以量化历史语言扩散背后的人口迁徙，并揭示社会结构变化。