IF50.5｜阿根廷中部八千年未间断的未知遗传谱系

阿根廷中部发现一个延续八千年的未知遗传谱系，重塑南美洲人群历史。

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

研究背景

南美洲是人类最晚定居的大陆之一，但其中部南锥体区域（Central Southern Cone, CSC）的古DNA研究长期不足。此前基于少量样本的研究表明，约9000BP后南美洲至少存在三个主要遗传簇：中安第斯、热带/亚热带森林（包括亚马逊）、以及智利中部/巴塔哥尼亚/潘帕斯。然而，由于采样稀疏，这些谱系的分化时间、相互关系以及CSC内部的人群动态仍不清楚。考古学证据显示该地区自晚更新世以来就有人类活动，并在中晚全新世出现了多样化的文化群体，但遗传学证据一直缺乏。

科学问题

卡在哪里： CSC地区古DNA样本极度匮乏，尤其是早中全新世个体，导致无法确定该地区早期人群的遗传构成、谱系分化时间以及与其他南美人群的基因交流历史。此前研究多基于现代人或零散古代个体，无法揭示区域内部长期的人口连续性和迁徙事件。

本文要回答： 通过大规模古DNA测序，重建CSC地区过去一万年的人群遗传结构，识别主要谱系，并探究其分化、迁徙和混合历史。

技术路线

作者主要用了：全基因组捕获测序、f4统计、qpAdm、qpWave、ROH分析、MDS。

作者对来自阿根廷中部及邻近地区的344份骨骼/牙齿样本进行DNA提取和文库构建，利用1240k或Twist1.4M探针进行全基因组SNP捕获测序，最终获得238个古代个体的基因组数据（中位覆盖659,011个SNP）。结合已发表的588个古代美洲个体和现代美洲原住民SNP芯片数据，通过f4统计、qpAdm/qpWave混合建模、outgroup-f3距离矩阵和MDS等方法，系统解析了CSC地区从早全新世到殖民接触前的人群遗传结构、谱系分化及基因流动事件。

图 1　Figure 1.。图1展示了本研究新报道的古代个体地理分布（a）和时间分布（b）。地图上标出了阿根廷不同生态区（西北、中部、大查科、巴拉那河、潘帕斯）及巴拉圭Pantanal的采样点，时间分布图显示样本跨越10000BP至150BP，其中深色部分表示有直接放射性碳测年的个体比例。

核心亮点

亮点 1｜10000BP个体揭示南锥体早期分化

最古老的个体LagunadelosPampas_10000BP（潘帕斯地区，10000BP）与中全新世南锥体个体（包括阿根廷潘帕斯7700-6800BP、南巴塔哥尼亚5100-7300BP）共享更多等位基因，而与中安第斯（9000-8600BP）和中东巴西（10400-6800BP）个体显著不同（f4统计Z值显著）。这表明南锥体人群与中安第斯和中东巴西人群的分化在10000BP前已经开始。此外，该个体与Anzick相关谱系有微弱但可检测的亲和力，暗示其可能混合了基底南锥体谱系和基底Anzick相关谱系，支持早期美洲人群结构呈梯度而非完全替代的模型。

图 2　Figure 2.。图2通过f4统计量量化了LagunadelosPampas_10000BP与早中全新世南美人群的亲和力。顶部六个显著非零的统计量表明该个体与中全新世南锥体个体（南巴塔哥尼亚、阿根廷潘帕斯）共享更多等位基因，而与中安第斯和中东巴西个体较少，支持南锥体谱系在10000BP前已分化。

亮点 2｜中部阿根廷谱系延续八千年

在科尔多瓦山区的JesusMaria_8500BP个体代表了一个此前未取样的深部谱系，该谱系在距今8500BP已具有独特的遗传漂变。从4200BP到150BP的中部阿根廷个体均与该个体形成clade，f4统计显示显著的正向偏斜（Z<5.54），表明该谱系在区域内持续存在至少八千年，未受大规模跨区域迁徙干扰。现代中部阿根廷混血个体仍显示出与晚期中部阿根廷古代个体的亲和力，说明该谱系是当地主要的原住民遗传成分。

图 3　Figure 3.。图3是基于outgroup-f3距离矩阵的MDS图（第1和第3主成分）。大多数新样本形成一个包含JesusMaria_8500BP的簇，水平轴区分中安第斯（右）与中部阿根廷（左），垂直轴区分南巴塔哥尼亚（上）与中东巴西（下）。西北阿根廷个体向中安第斯偏移，大查科和Pantanal个体向巴西偏移，反映了不同的混合历史。

亮点 3｜三次独立的基因流动事件

中部阿根廷谱系参与了三次不同的基因流动：1）约3300BP开始向南渗入潘帕斯，与当地中全新世谱系混合，并在800BP后成为潘帕斯的主要成分（qpAdm估计中部阿根廷成分在3300BP潘帕斯个体中占58±10%）；2）在西北阿根廷与中安第斯谱系混合，形成沿安第斯山脉的遗传梯度；3）在大查科地区与热带/亚热带森林谱系混合，现代大查科人群（如Wichí、Toba）可建模为中部阿根廷与亚马逊相关来源的混合。

图 4　Figure 4.。图4展示了qpAdm祖先成分估计和中部南锥体内部的精细结构。a图显示不同区域群体的祖先成分比例，中部阿根廷个体主要由中部阿根廷谱系构成，西北阿根廷和潘帕斯显示混合。b图基于outgroup-f3距离矩阵的MDS揭示中部阿根廷内部存在两条cline，分别反映与中安第斯和中全新世潘帕斯的混合。c图显示潘帕斯个体中中部阿根廷成分随时间增加。

亮点 4｜西北阿根廷近亲婚配率上升

通过ROH分析发现，西北阿根廷个体（≤3000BP）的>20cM ROH累积长度显著高于中部阿根廷（p<0.01）和潘帕斯（p<0.03），表明该地区在1000BP左右近亲婚配率增加。这一模式与中安第斯地区在Wari和Tiwanaku社会衰落后出现的现象平行，可能反映了类似ayllu社会组织（基于亲属关系的社区单位）的采用或强化。

亮点 5｜巴拉那河地区Guaraní相关个体

一个400BP、具有Guaraní考古学关联的个体在遗传上聚类于巴西群体，与Guaraní扩张的时间相符。然而，由于数据稀疏，无法进行祖先成分建模。此外，未发现巴拉那河地区群体与已知Arawak相关古代或现代群体之间的特定亲和力，不支持大规模Arawak迁徙，但文化传播或小规模迁移仍有可能。

生信可带走

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

• 方法栈： 全基因组捕获测序、f4统计、qpAdm、qpWave、ROH分析、MDS
• 公开数据： 原文未给出公开组学登录号
• 代码： 原文未给出公开 GitHub/GitLab 仓库
• 谁该点开原文： 从事古DNA、群体遗传学、南美考古学或人类演化研究的学者。

带走一句

大规模古DNA采样结合严格的统计建模，能揭示长期被忽视的区域人群连续性和局部基因流动，提示在解读现代人群遗传结构时需考虑深部谱系的贡献。