IF50.5｜绿色撒哈拉古DNA揭示北非祖先谱系

7000年前撒哈拉牧民基因组揭示一支此前未知的北非深部谱系，其与撒哈拉以南人群长期隔离。

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

研究背景

撒哈拉沙漠在距今14500至5000年前的非洲湿润期曾为绿色草原，湖泊河流广布，支持人类狩猎、采集与畜牧活动。此前西北非的古DNA研究显示，从至少15000年前到7500年前，当地存在一个稳定且遗传隔离的人群，即与Iberomaurusian石器文化相关的Taforalt洞穴采集者。约8300年前，家畜经西奈半岛和红海路线进入非洲，迅速扩散至东北非并抵达中撒哈拉。然而，由于DNA在干旱环境下难以保存，撒哈拉地区古人类基因组数据一直缺失，关于畜牧传播是人群迁徙还是文化扩散仍存争议。

科学问题

卡在哪里： 此前仅有来自Takarkori遗址的线粒体DNA数据，但线粒体DNA作为单一遗传位点，统计效力有限，无法解析群体动态。撒哈拉地区缺乏全基因组水平的古DNA数据，导致绿色撒哈拉时期人群的遗传来源、与西北非采集者的关系、以及畜牧传播机制等关键问题悬而未决。

本文要回答： 本文旨在通过获取Takarkori遗址两个约7000年前牧民个体的全基因组数据，揭示绿色撒哈拉人群的遗传构成，并探讨畜牧传播是文化扩散还是人群迁徙。

技术路线

作者主要用了：全基因组捕获测序、PCA、f3/f4统计、qpAdm、admixfrog、DATES、admixture graph。

作者从Takarkori岩棚遗址的两个女性个体（TKH001和TKH009）中提取DNA，由于内源DNA含量极低，采用Twist Ancient DNA和1240K捕获panel富集SNP位点，并对TKH001额外进行古人类混血SNP捕获。获得基因型数据后，通过PCA投影比较Takarkori与现代及古代非洲、近东人群的遗传位置；利用outgroup f3统计寻找与Takarkori共享遗传漂变最多的群体；通过f4统计检验Takarkori与Taforalt等群体的等位基因共享模式；使用qpAdm和admixture graph建模祖先成分；用admixfrog检测尼安德特人混血；用DATES估计混血时间。

图 1　Fig. 1。图1展示了晚更新世和全新世东北非和中撒哈拉的气候阶段与生计策略时间线，以及9000年前绿色撒哈拉时期和近代的生态区分布。Takarkori个体的碳十四年代用黑色菱形和圆形标出，位于中全新世畜牧时期。地图显示该遗址在绿色撒哈拉时期处于草原-稀树草原环境，而现代则为极端干旱区，凸显气候变迁对人群的影响。

核心亮点

亮点 1｜Takarkori代表未知北非谱系

PCA显示Takarkori个体位于西非和东非人群之间，与所有现代人群均不重叠。outgroup f3统计表明，Takarkori与摩洛哥旧石器时代晚期和早期新石器时代个体共享最多遗传漂变，特别是Ifri Ouberrid的Epipalaeolithic个体和15000年前的Taforalt采集者。admixture graph模型显示Takarkori约93%的祖先来自一个此前未知的北非人群，该人群与走出非洲人群关系密切，但在Mota之后分支，且长期保持隔离。

图 2　Fig. 2。图2的PCA投影显示Takarkori个体（合并数据）在PC1上介于西非和近东人群之间，但更接近西非；在PC2上位于萨赫勒/西非和东非人群之间，形成独立聚类。与Taforalt等西北非古代人群相比，Takarkori位置更偏东非方向，但仍与所有现代人群分离，表明其独特的遗传构成。

亮点 2｜与Taforalt的密切关系及模型修正

f4统计显示Takarkori与Taforalt的等位基因共享显著高于其他任何群体，且没有任何群体与Takarkori的亲和力超过Taforalt。qpAdm模型表明，Taforalt的祖先可由约60%的Natufian和40%的Takarkori相关成分解释，取代了此前约63.5% Natufian和36.5%撒哈拉以南非洲成分的模型。这一定位将Taforalt的非洲祖先成分从模糊的“撒哈拉以南”修正为具体的北非深部谱系。

图 3　Fig. 3。图3a的outgroup f3统计地图显示Takarkori与摩洛哥的Ifri Ouberrid、Taforalt和Ifri n’Amr o’Moussa个体共享遗传漂变最多（蓝色最深），而与撒哈拉以南非洲群体共享较少。图3b的f4统计显示所有测试群体与Takarkori的亲和力均不高于Taforalt，且撒哈拉以南群体与两者等距，表明Takarkori和Taforalt同属一个北非深部谱系，且与撒哈拉以南人群长期隔离。

亮点 3｜极低尼安德特人混血

admixfrog检测到Takarkori基因组中约0.15%的尼安德特人混血，片段长度超过0.05 cM。这一比例不到Taforalt和新石器时代摩洛哥个体（0.6-0.9%）的四分之一，比大多数走出非洲人群（1.4-2.36%）低约十倍，但显著高于撒哈拉以南非洲古人和现代人（几乎为零）。这表明Takarkori接受了少量来自走出非洲群体的基因流，但混血事件非常古老且不确定。

图 4　Fig. 4。图4a显示Takarkori的尼安德特人混血比例（约0.15%）显著低于Taforalt和黎凡特农民，但高于撒哈拉以南非洲人群。图4c的admixture graph模型将Takarkori建模为93%未知北非成分和7%黎凡特成分，Taforalt则为60% Natufian和40% Takarkori相关成分，与qpAdm结果一致。

亮点 4｜畜牧传播为文化扩散

admixture graph模型显示Takarkori仅有约7%的祖先来自黎凡特新石器时代人群，且DATES估计的混血时间非常古老。结合考古学证据，如物质文化的连续性和渐进变化，作者认为撒哈拉地区畜牧业的出现主要通过文化传播而非大规模人群迁徙。这与摩洛哥和东非地区新石器时代黎凡特血统显著增加的模式形成对比。

亮点 5｜撒哈拉长期遗传隔离

Takarkori与Taforalt个体与撒哈拉以南非洲谱系的遗传距离相等，f4统计未检测到任何撒哈拉以南群体与Takarkori有额外亲和力。这表明在非洲湿润期及之前的湿润时期，撒哈拉南北之间没有发生显著的基因交流。现代人群数据也显示撒哈拉两侧存在显著的遗传分化，支持撒哈拉作为长期基因流屏障的角色。

生信可带走

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

• 方法栈： 全基因组捕获测序、PCA、f3/f4统计、qpAdm、admixfrog、DATES、admixture graph
• 公开数据： 原文未给出公开组学登录号
• 代码： 原文未给出公开 GitHub/GitLab 仓库
• 谁该点开原文： 从事古DNA、群体遗传学、非洲史前史及畜牧传播研究的科研人员。

带走一句

古DNA捕获技术使极低内源DNA样本也能获得足够SNP进行群体遗传分析，为极端环境下的古基因组研究提供了可行路径。