IF50.5｜马格里布东部新石器时代人群的采集者祖先高度延续

古DNA揭示北非东部新石器人群以本地采集者血统为主，外来基因流影响远小于欧洲。

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

研究背景

地中海地区的古DNA研究已揭示与食物生产经济传播相关的长距离联系和人群转变。然而，与欧洲的丰富数据相比，北非这一关键转变期的数据此前仅来自马格里布西部（摩洛哥）。摩洛哥的古代个体显示，新石器时代人群主要源自本地“马格里布”采集者基因库，但也受到了来自伊比利亚的欧洲农民和来自黎凡特的牧人群体的基因流影响。马格里布东部（突尼斯和阿尔及利亚东北部）的考古记录表明其新石器化过程与西部不同，例如缺乏Cardial陶器，且农业出现较晚，暗示该地区可能有着不同的人口历史。

科学问题

卡在哪里： 此前缺乏马格里布东部新石器时代及更早时期的古DNA数据，无法检验该地区是否与西部有着相似的遗传转变，也无法评估欧洲农民和黎凡特人群对东部的影响程度。

本文要回答： 通过分析阿尔及利亚和突尼斯九个古代个体的全基因组数据，重建马格里布东部从晚石器时代到新石器时代的人群遗传历史，并比较其与西部及更广泛地中海地区的异同。

技术路线

作者主要用了：全基因组测序、PCA、qpAdm、f4统计、ALDER、DATES、hapROH、ancIBD。

作者对来自阿尔及利亚和突尼斯三个遗址的九个古代个体进行了全基因组测序，覆盖度约0.4–6.4x。首先通过主成分分析（PCA）将新个体投射到由现代地中海人群构建的遗传空间中，观察其与已发表古代人群的聚类关系。随后使用qpAdm对每个个体或群体进行祖先成分建模，以马格里布采集者、欧洲农民和黎凡特人群为代理来源，并通过f4统计量检验特定的等位基因共享模式。此外，利用ALDER和DATES推断混合事件的时间，用hapROH和ancIBD分析纯合性片段和IBD共享以评估群体大小和近期亲缘关系。

图 1　Figure 1:。图1展示了本研究中九个新报道个体的放射性碳测年结果和测序覆盖度，并与已发表个体对比。新个体来自阿尔及利亚的ABR（约15000–11000 BP）和突尼斯的Djebba（约8000 BP）、DEK（约7000–6350 BP）及Hergla（约5900 BP）。覆盖度从0.05x到6.4x不等，大部分个体覆盖度较高，适合群体遗传分析。地图显示这些遗址位于北非地中海沿岸，填补了马格里布东部的数据空白。

核心亮点

亮点 1｜马格里布祖先谱系时空分布广泛

来自阿尔及利亚Afalou Bou Rhummel（ABR，约15000–11000 BP）和突尼斯Djebba（约8000 BP）的早期个体，在PCA中与摩洛哥的Taforalt（TAF）和Ifri Ouberrid（OUB）等前新石器时代人群聚集在一起。qpAdm分析显示ABR个体（I13901）可完全由马格里布祖先成分拟合（p>0.18），Djebba个体之一（I20824）也可由100%马格里布祖先拟合（p=0.1）。这表明“马格里布”祖先谱系在晚更新世至早全新世期间广泛分布于北非西部至东部，具有较大的时空跨度。

图 2　Figure 2:。图2的主成分分析将新个体投射到由现代地中海人群构建的PC空间。早期个体（ABR和Djebba）与摩洛哥的TAF和OUB等前新石器时代人群聚集，而晚期新石器时代个体（DEK和Hergla）则向黎凡特和欧洲农民方向偏移，但仍处于马格里布祖先的三角区域内。这表明东部马格里布新石器时代人群保留了主要的本地祖先成分，同时吸收了少量外来基因流。

亮点 2｜Djebba个体携带西欧狩猎采集者血统

Djebba的另一位个体（I20825，约7900 BP）在qpAdm中无法由纯马格里布祖先拟合，但加入西欧狩猎采集者（WHG）相关成分后模型成功（5.7±1.1% WHG，p=0.36）。f4统计量显示I20825与WHG有显著的等位基因共享（Z值最高达5）。进一步比较不同欧洲狩猎采集者群体，发现其与西西里狩猎采集者的分化程度低于与西班牙和俄罗斯狩猎采集者，提示基因流可能来自西西里方向，与考古证据中跨西西里海峡的文化联系一致。

图 3　Figure 3:。图3的qpAdm结果对比了西部和东部马格里布人群的祖先成分。西部人群（如IAM、SKH、KEB）显示出较高的欧洲农民和黎凡特祖先比例，而东部人群（ABR、Djebba、DEK、Hergla）则以马格里布祖先为主，外来成分比例较低。例如，DEK1亚群约92%马格里布+8%欧洲农民，DEK2亚群约76%马格里布+24%黎凡特。这凸显了东西部在新石器化过程中遗传转变的差异。

亮点 3｜新石器时代东部马格里布人群以本地祖先为主

来自Doukanet el Khoutifa（DEK）和Hergla的新石器时代个体（约7000–5900 BP）在qpAdm中均需要马格里布祖先作为主要成分，比例约74–92%。DEK1亚群（约7000 BP）可由马格里布+欧洲农民祖先拟合（约92%和8%，p=0.16），而DEK2亚群（约6800 BP）则需要加入黎凡特祖先（约76%和24%，p=0.14）。Hergla个体（约5900 BP）的最佳模型包括马格里布、欧洲农民和黎凡特三种成分（约74%、13%和14%）。总体而言，东部马格里布新石器时代人群的外来祖先比例远低于西部（如KTG遗址欧洲农民祖先可达80%以上），显示出更强的本地连续性。

图 4　Figure 4:。图4展示了ROH和IBD分析结果。左侧显示东部马格里布个体的ROH分布，DEK1个体具有最长的ROH片段，提示其父母可能为二级表亲；而DEK2个体ROH很少，暗示较大的有效群体规模。右侧的IBD网络图显示DEK个体与西部马格里布人群（如SKH）共享中等长度片段，但与欧洲农民和黎凡特人群的共享较少，支持本地连续性为主、外来基因流有限的结论。

亮点 4｜混合事件时间与考古阶段对应

利用ALDER和DATES推断的混合时间显示，DEK个体的混合发生在其生前约5–25代之内，即混合事件发生在个体生活时间前数百年内。结合遗址的考古分期，DEK1亚群（约7400–7000 BP）的欧洲农民混合可能对应遗址早期阶段，而DEK2亚群（约7000–6300 BP）的黎凡特混合可能对应晚期阶段，该阶段驯化动物和典型新石器时代陶器及石器使用增加。

亮点 5｜东部马格里布群体规模较大且连续性高

ROH分析显示，东部马格里布个体的纯合性片段总长度低于西部，暗示其有效群体规模较大。例如，DEK2个体的有效群体大小估计至少约3500，而西部新石器时代个体则表现出广泛的ROH和低祖先群体大小。IBD分析发现DEK个体与西部马格里布人群（如SKH）共享中等长度片段（12–20 cM），表明东西部之间存在近期共同祖先，但未发现与黎凡特个体的IBD共享，提示黎凡特基因流可能通过间接方式或较大群体稀释了IBD信号。

生信可带走

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

• 方法栈： 全基因组测序、PCA、qpAdm、f4统计、ALDER、DATES、hapROH、ancIBD
• 公开数据： 原文未给出公开组学登录号
• 代码： 原文未给出公开 GitHub/GitLab 仓库
• 谁该点开原文： 从事古DNA、群体遗传学或北非考古学研究的学者，以及对新石器化过程感兴趣的研究生。

带走一句

古DNA分析需结合区域考古背景，避免将单一地区的模式推广到整个大陆；本地采集者群体的韧性和群体规模可能影响外来基因流的渗透程度。