乳糖耐受相关等位基因13910T与13915G的地理分布与欧亚大陆男性始祖事件紧密关联

《American Journal of Human Biology》：The Geographical Distribution of Lactose Tolerance-Associated Alleles 13910*T and 13915*G Is Strongly Linked to Male Founder Events in Eurasia

　　

乳糖耐受（lactase persistence, LP）是人类适应性进化的重要表现之一，其起源机制尚未完全阐明。本研究旨在分析欧亚大陆两种乳糖耐受相关等位基因——13910T与13915G的地理分布特征，并探讨其与主要Y染色体单倍群（Y haplogroups）及常染色体祖先成分的关联。研究人员收集了欧洲、近东及北非52个国家的13910T等位基因频率数据，以及30个国家中仅在近东与北非分布的13915G等位基因频率数据。结果显示，这两种等位基因的频率与所检测的遗传因素之间存在显著的线性相关性。13910T等位基因的出现与颜那亚（Yamnaya）祖先成分高度相关，支持其起源于欧亚草原颜那亚文化的假说。然而，该等位基因的后续扩散可归因于一系列区域性男性始祖事件及特定Y染色体单倍群的扩张，尤其是Y单倍群I1。相比之下，13915G等位基因的分布起源相对单一，与携带纳图夫（Natufian）祖先成分及Y单倍群J1的阿拉伯半岛牧民群体的地域性扩张密切相关。

该研究发表于《American Journal of Human Biology》，针对乳糖耐受这一人类快速进化性状的遗传基础与扩散机制展开深入解析。现有研究已明确乳糖耐受由位于MCM6基因调控区的单核苷酸多态性（SNP）驱动，其中欧亚大陆最常见的为13910T与13915G，但其起源时间与扩散路径仍存在争议，尤其缺乏结合父系遗传标记与全基因组祖先成分的系统性分析。为此，研究人员整合了多源基因频率数据，通过跨区域的遗传相关性分析，揭示了乳糖耐受等位基因传播背后的性别偏向性扩散模式及其与古代人群迁徙事件的耦合关系，为理解人类文化实践与遗传进化的相互作用提供了新证据。

关键技术方法包括：1）基于全球乳糖耐受关联数据库及文献补充的基因型数据，覆盖60个国家共5种LP相关SNP频率；2）收集25个Y染色体单倍群（含12个主类群及13个I、G、R1b亚支）的频率数据，样本量中位数为618人，范围60至169477人；3）采用三种常染色体祖先模型（k=7、k=11、k=12）解析人群遗传结构，其中k=7模型适用于欧洲全新世人群分析，k=11/k=12模型则纳入纳图夫、塔福拉尔特等关键祖先成分以覆盖近东与北非区域；4）结合青年男性平均身高数据进行表型关联分析，所有统计采用Pearson与Spearman相关性检验及因子分析。

研究结果如下：

3.1 乳糖耐受相关等位基因的分布

13910T在欧洲呈现北高南低的梯度分布，峰值见于爱尔兰（86.6%）与冰岛（85.3%），在东欧及西北非中等频率（15%-30%），其余地区低于10%；13915G仅分布于近东与北非，沙特阿拉伯（58.7%）与也门（54.9%）为高频区，埃及与叙利亚骤降至3.7%与2.9%。其他LP相关SNP（14009G、13907G、14010*C）因频率极低或分布局限未纳入后续分析。

3.2 乳糖耐受与Y染色体单倍群的关联

13910T在全样本中与斯堪的纳维亚Y单倍群I1相关性最强（r=0.80, p<0.001），联合日耳曼相关单倍群（I1、I2a-M223、R1b-U106）可使相关性提升至r=0.85；进一步整合Q、R1a、R1b-S116后达r=0.94（p<0.001）。东欧地区R1a与N单倍群对相关性有微弱正向贡献，而巴尔干特有单倍群（E1b-M78、J2、R1b-M269）呈显著负相关（r=-0.79）。13915*G在近东与北非仅与Y单倍群J1强相关（r=0.87, p<0.001），联合E1b-M123与T单倍群后相关性增至r=0.89，且其频率需在纳图夫相关单倍群累计占比>30%时才显著上升。

3.3 乳糖耐受与常染色体祖先成分的关联

13910T在欧洲唯一稳定正相关的祖先成分为颜那亚（Yamnaya^k=7，r=0.83），与维拉布鲁纳（Villabruna）、东欧狩猎采集者（EHG）的相关性存在东西欧异质性；近东地区其与高加索狩猎采集者（CHG）的正相关反映了颜那亚血统的渗入，北非地区则与主耕时代塔福拉尔特（Taforalt）成分相关（r=0.91）。13915G与纳图夫（Natufian）祖先成分呈极强正相关（近东r=0.96，近东北非r=0.93），且该关联强度与Y单倍群J1、E1b-M123、T的组合高度一致。

3.4 男性身高与乳糖耐受的关联

13910T与欧洲男性身高呈正相关（校正营养与儿童死亡率后仍r=0.56），该关联由颜那亚祖先成分介导；13915G则与身高呈负相关，对应阿拉伯相关Y单倍群（E1b-M123、J1、T）及纳图夫祖先成分的负向效应。

3.5 因子分析结果

欧洲群体中13910T与颜那亚祖先及日耳曼/颜那亚相关Y单倍群聚为一类，与巴尔干及近东Y单倍群呈拮抗；近东北非群体中13915G与J1及纳图夫祖先成分紧密聚集，与身高及安纳托利亚新石器成分拮抗。

讨论部分指出，13910T的扩散核心为北欧与西欧，Y单倍群I1（源自维拉布鲁纳集群）与R1b-S116（贝尔 beaker 文化标志）是关键驱动因素，其选择压力可能与贫瘠土壤环境下的饥荒适应相关，而非单纯乳制品摄入量。13915G的起源则局限于阿拉伯半岛，由J1单倍群下游分支（J1a-L858）与纳图夫祖先群体在特定干旱化时期经强烈遗传瓶颈扩张所致，与身高无正相关可能受伯格曼法则及资源限制影响。

结论重申，乳糖耐受等位基因的地理分布与父系始祖事件紧密耦合：13910T起源于颜那亚草原，经区域性男性奠基者效应在北欧加速扩散；13915G则源于阿拉伯半岛纳图夫相关群体的局部扩张，未来需更高分辨率的单倍群亚支数据以细化其演化细节。