“Modern humans were essentially like waves crashing on a beach, slowly but steadily eroding the beach away.” With this vivid simile, Princeton geneticist Joshua Akey describes a process that, until now, was shrouded in the mist of deep time: the incorporation, not extermination, of Neanderthals into the gene pool of modern humans. The new science, driven by innovations in ancient DNA recovery and machine learning, is not only rewriting the history of our closest extinct cousins it’s questioning the very concept of extinction in human evolution.

For decades, the Neanderthal narrative was one of vanishing, an extinct line displaced by the inexorable expansion of Homo sapiens. But since geneticists have sequenced the DNA of recent and ancient humans, a much more complex portrait has emerged. Employing genomes from 2,000 contemporary humans, three Neanderthals, and a single Denisovan, Akey’s team sought to chart the genetic transactions that formed our species in the last 250,000 years. Their ace in the hole: IBDmix, a machine learning technology that can identify buried archaic segments of DNA missed by traditional methods.
The IBDmix methodology represents a quantum leap over earlier approaches, which typically depended on “reference populations” assumed to be Neanderthal- and Denisovan-free. It so happens that even these reference populations frequently drawn from groups thousands of miles distant from recognized Neanderthal locations are contaminated with infinitesimal amounts of Neanderthal DNA, perhaps introduced by ancient travelers or their offspring. With the aid of machine learning, IBDmix finds segments of identity-by-descent (IBD) that mark common descent between modern and archaic humans, even if those segments are scattered or low in frequency.
This technological leap has resulted in a dazzling discovery: there was not one, but at least three separate waves of recent human–Neanderthal admixture in the last 250,000 years. The earliest was between 200,000 and 250,000 years ago, the second one between 100,000 and 120,000 years ago, and the biggest and most recent at around 50,000 to 60,000 years ago. These results overturn the long-standing perception that modern humans were confined in Africa for 200,000 years until a single dispersal event settled the remainder of the world. Rather, the genetic data suggest a dynamic history of migration, return, and repeated contact a tale of “dispersal, that modern humans have been moving around and encountering Neanderthals and Denisovans much more than we previously recognized,” as Akey explained to Princeton University News.
The consequences reach further than the timing and frequency of the admixture. By emphasizing the existence of modern human DNA within Neanderthal genomes instead of the better-known quest for Neanderthal DNA in modern humans Akey’s team revealed that modern human DNA could have comprised 2.5% to 3.7% of the Neanderthal genome, their analysis in Science as reported by Live Science. This is about equivalent to one in thirty contemporary human parents in the ancient Neanderthal population. Such evidence indicates that the genetic inheritance of contemporary humans within Neanderthals was not an isolated occurrence, but rather an intermittent aspect of their common heritage.
Drawing these conclusions from ancient bones is no small achievement. The ancient DNA extraction protocols, as outlined in Nature, include doing all the work in purpose-built clean-room facilities to avoid contamination, drilling bone powder with sterile instruments, and applying buffer solutions that have been optimized for the broken, delicate DNA of ancient specimens. The DNA is then transformed into single-stranded libraries and sequenced on high-throughput platforms like HiSeq4000 and NextSeq500. Quality control procedures like mapping onto the human reference genome with parameters optimized for ancient DNA, filtering for sequence length and mapping quality, and estimating contamination using software like AuthentiCT are necessary to guarantee authenticity of the data.
Sequenced, the problem then shifts to determining regions of archaic ancestry against the enormous backdrop of human genetic variation. For this, machine learning software such as IBDmix and ArchaicSeeker 2.0 are utilized. Simulation experiments have demonstrated that ArchaicSeeker 2.0 is highly accurate (approximately 93%) and has a very low rate of false positives (0.14%) in identifying archaic segments even under realistic admixture scenarios. These methods use linkage disequilibrium patterns, segment lengths, and allele frequencies to separate true introgressed segments from noise, enabling scientists to infer not only the timing and amount of admixture events but also their geographic and population-specific footprints.
The strength of these approaches lies in the new estimates of Neanderthal population size. Classic genetic modeling estimated population size using gene diversity, but IBDmix uncovered that much of the observed diversity in Neanderthals was actually derived from introgressed modern human DNA. The effective population size of Neanderthals was therefore corrected downward from around 3,400 to around 2,400 breeding individuals. “This doesn’t sound like a large difference, but given that Neanderthals were already estimated to have a fairly small population size, the fact that it was even smaller is an important insight,” Akey said in Live Science.
The genetic data now tend firmly towards the “assimilation model” of Neanderthal extinction, originally formulated by anthropologist Fred Smith in 1989. Instead of a catastrophic extinction, Neanderthals were progressively assimilated into growing modern human groups, their uniqueness eroded by repeated gene flow and population pressure. As co-author Liming Li remarked: “This is the first time that geneticists have identified multiple waves of modern human-Neanderthal admixture.”
The tale does not stop with Neanderthals. The very same analytical paradigms are being used to study Denisovans and other Eurasian and Oceanean archaic hominins, and a pattern of integrated ancestries across Eurasia and Oceania is emerging. For example, ArchaicSeeker 2.0 analyses have identified repeated episodes of both Neanderthal and Denisovan introgression in contemporary human populations with regional patterns of archaic ancestry and functional impacts on traits like immunity, pigmentation, and metabolism, as reported in Nature Communications.
These developments are not technological achievements only; they are redefining our sense of what it means to be human. The distinctions between “us” and “them” are breaking down, replaced by a view of human evolution as a history of contact, exchange, and shared heritage one in which the echoes of old encounters reverberate in our genomes today.

