IF50.5｜大麦驯化的单倍型演化史

用单倍型追溯大麦驯化，揭示其马赛克祖先来源与复杂演化路径。

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

研究背景

大麦是最古老的栽培作物之一，其驯化历史复杂。早期研究已否定单一驯化中心，支持马赛克基因组祖先模型，即驯化大麦由多个野生大麦群体贡献遗传物质。随着全基因组测序成本下降，研究者得以在更高分辨率下解析作物演化。PSMC等方法可推断历史有效群体大小，IntroBlocker等工具则能定义祖先单倍型组（AHG），在单倍型水平而非全基因组水平推断局部祖先。这些进展使研究者不仅能追问大麦在何处被驯化，还能探究驯化大麦基因组的不同区域如何追溯至其野生祖先。

科学问题

卡在哪里： 此前研究多依赖简化代表性测序或在驯化系中确定的标记，分辨率有限，难以精确回答驯化大麦的单倍型究竟源自哪些野生群体，以及不同野生群体是否对特定驯化谱系贡献不均。

本文要回答： 本文旨在利用单倍型方法，重建驯化大麦的单倍型水平祖先，明确哪些野生群体贡献了哪些基因组区域，以及单倍型如何通过驯化和驯化后基因流被重塑。

技术路线

作者主要用了：全基因组测序、古代DNA、IntroBlocker、PSMC、GEVA。

作者对682份种质资源（380份野生大麦、302份驯化大麦）进行全基因组测序，并结合23份考古大麦样本的古代DNA。首先利用ADMIXTURE和PCA将野生大麦划分为五个地理群体，然后使用IntroBlocker定义祖先单倍型组，比较野生与驯化大麦的单倍型共享情况。通过设置不同阈值对应不同分化时间区间，推断单倍型进入驯化基因库的时间和空间来源。同时利用PSMC重建群体历史，用GEVA估计关键驯化基因单倍型年龄，并用古代DNA验证现代样本的发现。

图 1　Fig. 1。图1展示野生大麦的采集地点和群体结构。ADMIXTURE分析（K=5）将143个野生大麦划分为五个地理群体：南黎凡特、叙利亚沙漠、北黎凡特、北美索不达米亚和中亚。图中仅显示未混合样本（主要祖先成分≥0.85）。不同颜色代表不同祖先群体，地理分布大致沿南黎凡特至中亚的路径。

核心亮点

亮点 1｜野生大麦的深层单倍型分化

在野生大麦中，远端区域的分化时间分布呈单峰，峰值约60万年前，对应有效群体大小的低谷。而近端区域（着丝粒附近）的分化时间分布呈多峰，且不同染色体间存在差异。特别地，在5H染色体上发现一个约90万年前分化的‘孑遗’单倍型，远深于其他野生群体间的分化。该单倍型在驯化大麦中占主导，并检测到选择性清除信号，提示其可能因适应性优势而被保留。

图 2　Fig. 2。图2a显示南黎凡特和北黎凡特群体间序列分化（SNPs/Mb）在远端、间质和近端区域的分布。右侧灰色阴影标记5H染色体上由于深度分化单倍型导致的最高分化。图2b展示所有样本对在远端区域的分化分布（上）以及PSMC推断的野生大麦有效群体大小历史轨迹（红）和全球平均表面温度（灰）。

亮点 2｜驯化大麦的马赛克单倍型结构

驯化大麦的单倍型分化呈双峰分布，近期峰对应约8000年前的瓶颈，早期峰对应驯化前分化。通过IntroBlocker，55.9%的驯化单倍型在至少一个野生样本中存在，而仅7.0%的野生单倍型与驯化大麦共享。约91%的驯化单倍型与野生对应单倍型在32,000至8,000年前分化，表明马赛克结构在驯化早期即已形成。所有五个野生群体均对驯化大麦有贡献，但比例不同：南黎凡特、北黎凡特和中亚各贡献20-27%，叙利亚沙漠和北美索不达米亚分别贡献16.4%和12.9%。

图 3　Fig. 3。图3a显示野生与驯化大麦间的序列分化分布。蓝线标记400 SNPs/Mb阈值，用于区分驯化前和驯化后单倍型。绿色阴影表示从8 ka bp开始分裂的创始群体的持续存在。近期峰约300 SNPs对应25 ka bp，标记祖先驯化大麦与现存野生大麦的分化。图3b展示1H染色体上使用5-Mb窗口推断的AHG。

亮点 3｜驯化基因的单倍型年龄与起源

利用GEVA估计关键驯化基因单倍型年龄：btr1约27,000年，btr2约15,000年，vrs1.a1约25,000年，vrs1.a2约8,000年，vrs1.a4约7,000年。这些单倍型的最接近野生对应物分布在不同野生群体中：btr1、nud和vrs1.a3来自南黎凡特，btr2和vrs1.a2来自北黎凡特，vrs1.a1和vrs1.a4来自北美索不达米亚和中亚。这支持多中心驯化模型，即不同驯化性状的等位基因源自不同野生群体。

图 4　Fig. 4。图4a展示EU-THS群体在2H染色体上20-Mb窗口的单倍型时空起源。彩色条高度与单倍型在特定时期（左）或来自特定野生群体（右）的概率成正比。黄色表示来源未知。图4b展示所有驯化群体中单倍型的时间起源分布，约91%的单倍型在32-8 ka bp间与野生对应物分化。

亮点 4｜古代DNA揭示持久群体结构

23份古代大麦样本（6000-2000 cal bp）均与栽培类型聚类，携带驯化型btr1Btr2单倍型。大多数古代样本携带与现代驯化大麦相同的长着丝粒单倍型，但7H染色体上存在私有单倍型。遗传多样性分析显示，从Yoram Cave（6 ka bp）到Abi'or Cave（2 ka bp）多样性逐渐增加，D统计检测到Timna 34和Abi'or Cave样本中存在来自西方群体的基因流，表明长距离贸易和人口流动对以色列大麦遗传多样性有重要贡献。

图 5　Fig. 5。图5a展示不同驯化大麦群体间序列分化分布的小提琴图，虚线标记分布峰值（分裂时间）。多峰分布可能源于基因流事件。图5b示意15个驯化大麦群体和西藏青稞的谱系分化时间。图5c地图显示驯化大麦从新月沃地起源中心向外扩散的时间和路线。

亮点 5｜驯化大麦的谱系分化与扩散

利用驯化后单倍型的分化时间推断群体分裂时间：西方、东方和埃塞俄比亚大麦约8,500年前分化，反映农业从新月沃地扩散。西方大麦随后在约7,500年前分化为近东、欧洲二棱和欧洲六棱三个谱系，与考古记录一致。欧洲二棱和六棱大麦的分化可能受形态和物候选择驱动，而非地理隔离。

生信可带走

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

• 方法栈： 全基因组测序、古代DNA、IntroBlocker、PSMC、GEVA
• 公开数据： 原文未给出公开组学登录号
• 代码： 原文未给出公开 GitHub/GitLab 仓库
• 谁该点开原文： 从事作物演化、群体遗传学或单倍型分析的生信研究者。

带走一句

单倍型水平的分析揭示了大麦驯化的马赛克本质，提示在作物演化研究中，局部祖先推断比全基因组平均更能捕捉真实历史。