The date of interbreeding between Neandertals and modern humans

Abstract: Comparisons of DNA sequences between Neandertals and present-day humans have shown that Neandertals share more genetic variants with non-Africans than with Africans. This could be due to interbreeding between Neandertals and modern humans when the two groups met subsequent to the emergence of modern humans outside Africa. However, it could also be due to population structure that antedates the origin of Neandertal ancestors in Africa. We measure the extent of linkage disequilibrium (LD) in the genomes of present-day Europeans and find that the last gene flow from Neandertals (or their relatives) into Europeans likely occurred 37,000-86,000 years before the present (BP), and most likely 47,000-65,000 years ago. This supports the recent interbreeding hypothesis, and suggests that interbreeding may have occurred when modern humans carrying Upper Paleolithic technologies encountered Neandertals as they expanded out of Africa.

• The paper introduces a novel LD decay-based statistic to accurately estimate the timing of gene flow from Neandertals into modern human populations.
• Simulation results validate the method under varied demographic models while highlighting biases from ancient structure and population bottlenecks.
• The study enhances paleogenetic reconstructions by providing actionable insights into human evolutionary dynamics and historical migration patterns.

Summary of Insights from the Research on Inferring Gene Flow Timing Using LD Decay

This paper addresses the challenging task of estimating the timing of ancient gene flow events between populations, particularly focusing on gene flow from Neandertals to modern non-African human populations. The research introduces a robust statistical method grounded in the decay pattern of linkage disequilibrium (LD) due to admixture. The primary objective is to estimate extremely old gene flow dates, like those older than 10,000 years BP, which necessitates managing uncertainties such as recombination rates.

Methodological Approach and Statistical Framework

The authors propose a novel statistic to estimate the gene flow timing, $t_{GF}$, utilizing patterns of LD decay in European (CEU) populations with reference to Neandertal-derived alleles. The method is structured around three main assumptions:

1. The ascertainment of SNPs carrying derived alleles from Neandertals, being polymorphic within Europeans, and exhibiting specific allele frequency thresholds.
2. The expectation that the negative exponential decay of LD with genetic distance—where the decay rate corresponds to $t_{GF}$—remains consistent across certain demographic shifts, though ascertainment imperfections can introduce variability.
3. The use of coalescent simulations to test the robustness of this method against various demographic models, such as those simulating recent gene flows, ancient structures, and no gene flow scenarios.

Simulation Results and Numerical Findings

The paper reports on extensive simulations that validate the proposed framework under various demographic histories:

• Robust Estimation Under Recent Gene Flow: The statistic accurately traces the true timing of gene flow when applied to recent exchanges, even amid demographic fluctuations post-admixture.
• Ancient Structure and No Gene Flow Models: In contrast, models without gene flow or with ancient population structure suggested by extended LD distortions demonstrate potential biases in estimation, indicating the necessity of careful demographic scenario considerations when applying this model.
• Impact of Population Size Changes and Sampling Schemes: Key results showed downward biases in cases exhibiting significant population bottlenecks or ancient structures, providing insights into ascertainment strategy impacts on the accuracy of gene flow timing estimates.

Implications and Future Directions

The study bridging genetic data analysis through LD decay estimation holds substantial implications for genetic anthropology and evolutionary biology:

• Theoretical Contributions: This research enriches understanding of the temporal nuances in human evolution, particularly in tracing the extent of interbreeding events facilitated by admixture.
• Practical Application: The framework offers a tool for more accurate paleogenetic reconstructions that could inform on historical population dynamics and migration patterns.
• Further Development: Future work may focus on refining SNP ascertainment strategies, integrating additional genetic maps, and correcting potential map-induced biases to enhance the predictive power of LD decay-based gene flow estimation.

Overall, this paper contributes to the methodological arsenal available to researchers aiming to resolve chronologies of historical gene flow events by combining robust statistical approaches with high-resolution genetic data.