[论文解读] Interbreeding conditions for explaining Neandertal DNA in living humans: the nonneutral case

A. Neves|arXiv (Cornell University)|May 30, 2011

Pleistocene-Era Hominins and Archaeology参考文献 17被引用 3

一句话总结

本文通过引入适应度差异，扩展了先前关于智人与尼安德特人杂交的中性模型，表明即使智人仅具有1%的适应度优势，尼安德特人灭绝速度也过快，无法解释中东地区长期共存的考古证据。非中性模型预测，为达到非非洲人群中1–4%的尼安德特人DNA，需要更高的杂交率；但这种快速灭绝与以色列洞穴中尼安德特人与智人长期交替存在的记录相矛盾，因此中性模型在解释遗传与考古数据方面仍更可信。

ABSTRACT

We consider here an extension of a previous work by Neves and Serva, still unpublished, which estimates the amount of interbreeding between anatomically modern Africans and Neandertals necessary for explaining the experimental fact that 1 to 4% of the DNA in non-African living humans is of Neandertal origin. In that work we considered that Africans and Neandertals had the same fitness (neutral hypothesis) and Neandertal extinction was thus an event of fortune. In this work we consider that Africans had larger fitnesses. We show results for four values for the fitness difference: 1%, 5%, 10% and 20% and compare them with the corresponding neutral results. Some technical differences with respect to the neutral case appear. We conclude that even with 1% fitness difference Neandertals extinction comes up in too small a time, so the neutral model looks more suitable for explaining the known data on occupation of some caves in Israel for a very long time, alternately by Africans and Neandertals.

研究动机与目标

• 探究非中性模型（即智人相对于尼安德特人具有更高适应度）是否能解释非非洲人群中1–4%的尼安德特人DNA。
• 评估智人与尼安德特人之间的适应度差异是否与中东地区长期共存的考古记录相容。
• 比较非中性模型与先前研究的中性模型的预测，特别是关于所需杂交率和灭绝时间的差异。
• 评估非中性模型中使用的确定性近似在不同种群规模和适应度差异下的有效性。

提出的方法

• 本研究采用确定性、平均场模型，模拟两个子种群（智人与尼安德特人）在恒定大小为 $ N $ 的种群中的动态，世代不重叠。
• 引入适应度差异 $ h $，其中智人相对于尼安德特人具有相对适应度优势，并模拟智人种群中尼安德特人等位基因比例随时间的变化。
• 使用缩放后的杂交参数 $ \beta = \alpha / N $，其中 $ \alpha $ 为每代的杂交率，以确保结果与 $ N $ 无关。
• 在非中性假设下推导尼安德特人DNA贡献的概率分布，并与中性情况（$ h = 0 $）通过数值模拟和解析近似进行比较。
• 以20年为一代，确定性地计算灭绝时间，作为初始智人比例和适应度差异 $ h $ 的函数。
• 通过检查条件 $ \alpha \ll 1 / \log_2 N $ 来评估模型的有效性，以确保平均场近似成立。

实验结果

研究问题

• RQ1 引入智人适应度优势后，为在非非洲人群中实现1–4%的尼安德特人DNA，所需杂交率如何变化？
• RQ2 非中性模型能否解释在以色列洞穴中观察到的尼安德特人与智人长期交替存在的考古记录？
• RQ3 随着适应度差异 $ h $ 的变化，尼安德特人的灭绝时间如何变化？是否与斯库勒和卡夫济遗址观察到的13万年交替占据模式一致？
• RQ4 非中性模型中使用的确定性近似在真实种群规模和杂交率下，其有效性在多大程度上仍成立？
• RQ5 在调和遗传数据与共存考古证据方面，中性模型是否仍比非中性模型更合适？

主要发现

• 即使智人具有1%的适应度优势（$ h = 0.01 $），尼安德特人灭绝也仅需约600至1,000年，这一速度过快，无法匹配以色列洞穴中尼安德特人与智人长达13万年的交替占据。
• 当 $ h = 0.01 $ 时，所需杂交参数 $ \beta $ 约为 $ 9.3 \times 10^{-4} $，且概率密度的最大值出现在 $ \beta_{\mathrm{max}} \approx 5.0 \times 10^{-4} $，表明需要极高的杂交率。
• 随着适应度差异 $ h $ 增大，所需 $ \beta $ 增加，灭绝时间缩短，其中 $ h = 0.2 $ 时灭绝时间为600–1,000年，$ h = 0.01 $ 时为11,000–20,000年。
• 模型表明，当 $ N = 1000 $ 时，若 $ \beta \gg 10^{-4} $，确定性近似将失效，提示对较大 $ \beta $ 值的结果需谨慎解读；然而，较小的 $ N $ 可能提高准确性。
• 尽管非中性模型在数学上更简洁，但中性模型仍更符合长期共存的考古记录，因此更适合解释遗传与考古数据。