内蒙古北狄族群在龙山晚期的南下是否影响到豫东的有虞氏族群？

在距今4300年前（龙山时代晚期），内蒙古岱海地区老虎山文化突然消失，其时间与使全新世大暖期结束的寒冷事件相吻合。此后，老虎山文化因素渗入到中原地区。韩建业教授认为老虎山文化族群与古籍中的北狄系统族群相关。

不久前发表的研究，证实河南仰韶村龙山时代人群，相对于本地仰韶文化先民，呈现出显著的类似石峁族群（老虎山文化向陕北延伸）。但，这种现象并没有出现在此前报道的豫东同时代的造律台类型等中原龙山族群。其中，豫东龙山文化的平粮台古城族群，“可能就是传说中的有虞氏文化”（李伯谦《论造律台类型》），舜是有虞部落的首领。

NO.样本编号遗址年代推测C14测年（BP）采样部位

14PLTM310平粮台新石器晚期（龙山）4,225-3,995颞骨

15PLTM311平粮台新石器晚期（龙山）4,151-3,974颞骨

16PLTM312平粮台新石器晚期（龙山）4,085-3,889颞骨

17PLTM313平粮台新石器晚期（龙山）4,068-3,844颞骨

Science主刊，最近发表了最新混合估算工具Twigstats（https://www.23mofang.com/community/6791cf1e0ac10427faf53f05），Twigstats 可以准确且无偏差地推断混合比例，排除远古基因流的干扰信号。趁着假期有空，部署了Twigstats以及上下游工具包，以此工具，针对性考察，消失于内蒙古的老虎山文化北狄系统族群的南下，是否影响到豫东的有虞氏华夏系统族群。这里以老虎山文化本地更古老的庙子沟-海生不浪文化族群（pop2，MZGM）作为代理人群，汪沟遗址（pop1，WGM）作为河南本地土著的代理人群。外群（popO）为非洲约鲁巴人，参考人群（popI，若干西北喇家遗址样本），以下为4个平粮台古城样本（popX，PLTM）的计算结果如下，其中，标准差se<0.1作为阀值，t_value代表的世系为50-1000代。可以观察到，Twigstats的确可以获得较为理想的se值，结果可靠性大幅提升。平梁台古城样本的本地仰韶文化汪沟组血统接近100%，并没有显著的庙子沟组血统信号。这样看来，在新石器时代，太行山脉很可能很大程度上影响了东亚地区种群的基因流动。

popO  popI  pop1  pop2  popX  val   se    t_value

sample_30   sample_3    WGM20 MZGM10      PLTM310     1.002204708 0.007513335 50

sample_30   sample_9    WGM20 MZGM10      PLTM310     0.998453594 0.048169718 50

sample_30   sample_6    WGM20 MZGM10      PLTM310     0.997593137 0.051258654 50

sample_30   sample_3    WGM20 MZGM10      PLTM310     0.98362124  0.070972986 100

sample_30   sample_9    WGM20 MZGM10      PLTM310     0.970671555 0.07577214  100

sample_30   sample_5    WGM20 MZGM10      PLTM310     1.421006768 0.080929389 1000

popO  popI  pop1  pop2  popX  val   se    t_value

sample_30   sample_3    WGM20 MZGM10      PLTM311     1.001969358 0.007039124 50

sample_30   sample_8    WGM20 MZGM10      PLTM311     0.998716114 0.034184659 100

sample_30   sample_9    WGM20 MZGM10      PLTM311     0.998484913 0.047198234 50

sample_30   sample_6    WGM20 MZGM10      PLTM311     0.997712502 0.048733666 50

sample_30   sample_9    WGM20 MZGM10      PLTM311     0.963782838 0.094912486 100

sample_30   sample_3    WGM20 MZGM10      PLTM311     0.968182885 0.098662286 100

popO  popI  pop1  pop2  popX  val   se    t_value

sample_30   sample_7    WGM20 MZGM10      PLTM312     -4.15E-17   1.35E-16    50

sample_30   sample_3    WGM20 MZGM10      PLTM312     0.999087511 0.006576457 50

sample_30   sample_5    WGM20 MZGM10      PLTM312     0.534820238 0.038471016 200

sample_30   sample_4    WGM20 MZGM10      PLTM312     1.02789811  0.05842565  500

sample_30   sample_9    WGM20 MZGM10      PLTM312     0.564050393 0.058961366 1000

sample_30   sample_5    WGM20 MZGM10      PLTM312     0.280735456 0.06230908  1000

sample_30   sample_6    WGM20 MZGM10      PLTM312     0.58059661  0.063343314 200

sample_30   sample_5    WGM20 MZGM10      PLTM312     0.595604374 0.068129192 300

popO  popI  pop1  pop2  popX  val   se    t_value

sample_30   sample_7    WGM20 MZGM10      PLTM313     -5.55E-15   1.62E-14    50

sample_30   sample_3    WGM20 MZGM10      PLTM313     1.001367688 0.004909636 50

sample_30   sample_3    WGM20 MZGM10      PLTM313     0.991538594 0.041124502 100

sample_30   sample_5    WGM20 MZGM10      PLTM313     0.994694834 0.081011252 50

sample_30   sample_6    WGM20 MZGM10      PLTM313     0.965501899 0.091373329 200