黄河中下游地区数千年遗传连续性与跨区域混合的古基因组研究

《Cell Reports》：Genomic insights into millennia of genetic continuity and trans-regional admixture in the middle-lower Yellow River basin

　　

黄河中下游地区对东亚早期农业发展及复杂社会的出现至关重要，但其从新石器时代至历史时期的精细尺度人群历史仍未完全解析。研究人员呈现来自河南和山东的61例古代个体（约4,500–200年前（before present, BP））的全基因组数据。分析显示，河南大河村（Dahecun）中晚新石器时代人群遗传同质，并与同期黄河粟作农业人群（Yellow River millet-farming populations）密切相关。至约3,000年前，山东人群向新石器时代晚期中黄河相关祖源谱（ancestry profile）转变，这与中黄河农业人群东向扩张一致。该祖源谱持续至历史时期，表明山东长期遗传连续性。来自山东西部的两例唐宋时期个体携带约15%西欧亚相关祖源（West Eurasian-related ancestry），表明有限的跨区域遗传互动。这些发现细化了黄河中下游地区的人群历史。

论文解读：黄河中下游地区人群遗传连续性与跨区域混合的古基因组研究

研究背景：黄河中下游是东亚早期农业出现、社会复杂化和早期文明发展的关键区域。约10,000年前，该地区出现作物驯化，形成以黍和粟为中心的农业体系。中原地区从新石器时代至历史时期具有连续文化序列，包括仰韶文化（Yangshao culture，约7,000–5,000 BP，中全新世粟作农业文化）、龙山文化（Longshan culture，约5,000–4,000 BP，晚全新世文化）、二里头文化（约3,900–3,500 BP）、商文化（约3,600–3,000 BP）及后续王朝。关于中黄河流域在史前文化发展中的作用及其与周边互动长期存在争论：Su Bingqi的“Man Tian Xing Dou”（“full constellation”）模型强调多区域文化传统并行出现；Yan Wenming的“diversity in unity”（petal-like structure）模型认为约5,000年前逐渐形成以中原为核心、周边为外围的层级文化结构。考古证据不能完全说明文化互动是否伴随人群扩散。已有古DNA研究显示，以仰韶相关群体为代表的新石器时代中期黄河人群（middle Neolithic Yellow River populations, YR_MN）与以龙山相关群体为代表的新石器时代晚期黄河人群（late Neolithic Yellow River populations, YR_LN）之间存在显著遗传连续性，后者携带更多南方东亚相关祖源。但河南为中心的中原史前人群是否存在遗传亚结构、仰韶至龙山转变期的区域文化互动是否伴随可检测基因流，以及约3,000 BP以来山东建立的遗传结构是否长期连续并贡献现代山东汉族，仍不清楚。为此，研究人员分析来自河南大河村（Dahecun）、山东西高（Xigao）和朝城三里营南（Chaocheng Sanliyingnan）的61例古代个体（约4,500–200 BP）全基因组数据，以刻画黄河中下游精细人群遗传史。

研究人员开展的研究、结论及意义：研究人员整合新测与已发表古基因组数据，发现河南中晚新石器时代大河村人群遗传同质，并与同期黄河粟作农业人群关系密切；约3,000 BP山东人群转向晚新石器时代中黄河相关祖源谱，符合中黄河农业人群东扩；该祖源谱持续至历史时期，表明山东长期遗传连续；两例唐宋山东西部个体携带约15%西欧亚相关祖源，提示有限跨区域遗传互动。这些发现细化了黄河中下游地区人群历史，为理解农业扩散、文化互动与人群迁徙关系提供古基因组证据。

关键技术方法（≤250字）：样本队列来自河南大河村、山东西高和朝城三里营南，共采集112例古代人骨，时间约4,500–200 BP；提取古DNA并建库测序，保留61例全基因组数据。采用加速器质谱（AMS）放射性碳定年、死后损伤（postmortem damage, PMD）与污染质控、1240K单核苷酸多态性（SNP）面板伪单倍体分型、线粒体与Y染色体单倍群分析、READv.2亲缘分析。群体遗传方法包括主成分分析（PCA）、无监督ADMIXTURE、qpWave、f₄统计量、Outgroup-f₃统计量、qpAdm、TreeMix和admixfrog。

研究结果：

Ancient genome-wide data from the middle and lower Yellow River basin：研究人员从河南和山东三个遗址的112例人骨中提取DNA并建库，获得61例古代个体全基因组数据，平均覆盖度为0.002×–1.6×。放射性碳定年显示，大河村个体约4,800–4,500 cal. BP，属中晚新石器时代；西高个体约3,400–1,800 cal. BP，并延伸至清代。经PMD和污染质控后，线粒体单倍群包括D、B、M、G、C、N、F，父系以O为主，另有C、N和1例R。READv.2鉴定15对二级及以上亲缘个体，并去除低覆盖个体。qpWave分组显示，大河村仰韶和龙山个体分别同质，西高六个历史时期组同质，朝城三里营南发现CC32和CC35两个离群个体，命名为Shandong_CC_Tang_Song_o1，其余合并为Shandong_CC_Tang_Song。

Close genetic affinities among populations in the middle-lower Yellow River basin：PCA显示东亚人群呈东北亚、东南亚和青藏高原相关三向遗传梯度。大河村仰韶和龙山人群紧密聚集，并与河南仰韶/YR_MN人群大量重叠，提示中黄河流域中晚新石器时代遗传连续。山东史前人群形成HG、大汶口期和龙山期三个主要遗传簇，历史时期山东人群向YR_LN和YR_LBIA亲和并紧密聚集；两个CC离群个体向中亚和西欧亚相关人群偏移。ADMIXTURE在K=3时支持上述模式。

Genetic continuity in the middle-to-late Neolithic Central Plains：大河村类型为晚仰韶文化区域变体。qpWave显示其与中原仰韶人群整体遗传同质；Outgroup-f₃显示其与中黄河流域人群和现代汉族遗传漂变最高；f₄未发现大河村仰韶与YR_MN存在显著等位基因共享偏差，qpAdm支持单源YR_MN。龙山期大河村两个已发表个体保留YR_MN相关祖源，与新样本无显著差异，合并为Dahecun_Longshan。Outgroup-f₃显示其与Dahecun_Yangshao和现代汉族遗传漂变最高；f₄未检测到ANA（Northeast Asian-related ancestry，东北亚相关祖源）、SEA（southern East Asian coastal ancestry，中国南方沿海相关祖源）或山东狩猎采集者（Shandong hunter-gatherers, Shandong_HG）额外等位基因共享，qpAdm可用Dahecun_Yangshao单源建模，TreeMix无外部基因流证据。结论为从仰韶到龙山期大河村人群显著遗传连续。

Persistence of Yellow River-related ancestry in Shandong across historical periods：约3,000 BP山东人群向YR_LN/YR_LBIA相关祖源谱转变，可能反映中原农业人群东扩。qpWave显示商代西高人群与YR_LN、YR_LBIA遗传同质，且与后续山东人群同质，提示黄河农民相关祖源长期连续。Outgroup-f₃显示各历史时期西高人群与黄河相关人群和现代汉族遗传漂变最高；f₄未发现显著偏离；以商代西高为局部参考时，后续历史时期无显著额外等位基因共享；qpAdm可用商代西高单源建模所有历史时期人群，现代山东汉族也可用历史西高人群单源建模。

Subtle Western Eurasian-related genetic signals in Tang-Song Shandong populations：聊城唐福大门遗址多数人群与黄河农业人群连续，少数个体约5%中亚或西欧亚相关祖源，与丝绸之路互动及“Hu”人移民一致。新测唐宋聊城CC人群多数与黄河农业人群聚集，两个离群个体向中亚和西欧亚偏移。f₄显示CC主组相对福大门无显著偏差；qpWave支持唐时期两邻遗址同质；qpAdm主CC可用单源YR_LBIA建模。o1组被qpWave拒绝与山东历史人群同质；f₄显示其与中亚和西欧亚参考共享更多等位基因；qpAdm两向模型显示约85%黄河相关祖源和约15%中亚或西欧亚相关祖源；仅转换SNP重复仍支持该信号。一例o1个体Y染色体单倍群暂定为R1b，但覆盖度低，需谨慎。

Archaic introgression signals in ancient Yellow River populations：admixfrog分析41例信息位点充足的古代个体。两例Shandong_CC_Tang_Song_o1的丹尼索瓦相关和尼安德特相关渗入比例与同期唐宋组无明确偏离，标准化Z分数分别为?0.57和?0.61。时间组比较中，丹尼索瓦相关渗入比例在仰韶组高于唐宋和清组；尼安德特相关渗入比例在仰韶和商组高于唐宋和清组；其余比较无显著差异。77个基因与至少20例个体的古渗入片段重叠，14个基因在30例以上个体中检出；TBC1D1最常见，在41例中36例重叠，该基因与葡萄糖转运和能量代谢相关，曾报道可能参与“oceanian”和高海拔人群适应性渗入。

讨论部分总结：中原新石器时代总体遗传同质但存在局部混合。考古上仰韶文化中心东移并形成多区域“花瓣”结构；大河村处于中黄河与下黄河之间，与 Dawenkou 和 Qujialing 文化互动频繁，但遗传结果未检测到周边基因流，说明文化接触未必伴随大规模人群移动；而南部朱岗遗址因更靠近长江流域而检测到Qujialing相关基因流，提示地理邻近和接触强度影响混合。山东在商、齐、秦汉、魏晋南北朝、宋元等历史变迁中经历政治整合和人口流动，但西高各历史阶段均可由YR_LN/YR_LBIA相关祖源解释，并延续至现代山东汉族，与临淄、淄博古基因组一致。研究人员提出机制：定居汉族农业人口规模大产生稀释效应；游牧政权政治主导不一定导致大规模人群迁移，农牧生活方式差异可能构成生态屏障；移民可数代内文化和遗传同化。由于样本量和时空覆盖有限，尤其集中于山东中部，未来需更大队列和更广地理横断面验证。唐宋时期东西祖源贡献复杂。丝绸之路和大运河使聊城成为交通节点；CC少数个体约15%外源祖源，qpAdm最优代理指向中亚青铜至铁器时代草原相关人群，如Kazakhstan_MLBA、Russia_MLBA_Krasnoyarsk/Sintashta和Uzbekistan_Dzharkutan_BA，与粟特商人和景教传教士历史背景一致。仅转换SNP排除PMD偏差，同地区福大门人群验证该信号；当地景教教堂遗迹及平民墓地中检出该祖源，提示中亚相关个体已融入基层社会。Y染色体R1b结果需高覆盖数据验证。古渗入分析中，o1外源祖源比例较小，未明显改变古渗入比例；历史组古渗入低于史前组可能反映祖源稀释，也受重组、覆盖度、样本量和古参考基因组限制。TBC1D1等基因重叠需进一步功能验证。局限包括样本数量和地理分布有限、性别相关全基因组祖源差异未系统评估、Y染色体标记稀疏、古渗入估计受覆盖度和参考基因组影响。

研究结论翻译：这些发现细化了黄河中下游地区的人群历史。