Being at the western fringe of Europe, Iberia had a peculiar prehistory and a complex pattern of Neolithization. A few studies, all based on modern populations, reported the presence of DNA of likely African origin in this region, generally concluding it was the result of recent gene flow, probably during the Islamic period. Here, we provide evidence of much older gene flow from Africa to Iberia by sequencing whole genomes from four human remains from northern Portugal and southern Spain dated around 4000 years BP (from the Middle Neolithic to the Bronze Age). We found one of them to carry an unequivocal sub-Saharan mitogenome of most probably West or West-Central African origin, to our knowledge never reported before in prehistoric remains outside Africa. Our analyses of ancient nuclear genomes show small but significant levels of sub-Saharan African affinity in several ancient Iberian samples, which indicates that what we detected was not an occasional individual phenomenon, but an admixture event recognizable at the population level. We interpret this result as evidence of an early migration process from Africa into the Iberian Peninsula through a western route, possibly across the Strait of Gibraltar.

Domesticated crops have been disseminated by humans over vast geographic areas. Common bean (Phaseolus vulgaris L.) was introduced in Europe after 1492. Here, by combining whole-genome profiling, metabolic fingerprinting and phenotypic characterisation, we show that the first common bean cultigens successfully introduced into Europe were of Andean origin, after Francisco Pizarro’s expedition to northern Peru in 1529. We reveal that hybridisation, selection and recombination have shaped the genomic diversity of the European common bean in parallel with political constraints. There is clear evidence of adaptive introgression into the Mesoamerican-derived European genotypes, with 44 Andean introgressed genomic segments shared by more than 90% of European accessions and distributed across all chromosomes except PvChr11. Genomic scans for signatures of selection highlight the role of genes relevant to flowering and environmental adaptation, suggesting that introgression has been crucial for the dissemination of this tropical crop to the temperate regions of Europe.

The Arctic climate was warmer than today at the last interglacial and the Holocene thermal optimum. To reveal the impact of past climate-warming events on the demographic history of an Arctic specialist, we examined both mitochondrial and nuclear genomic variation in the collared lemming (Dicrostonyx torquatus, Pallas), a keystone species in tundra communities, across its entire distribution in northern Eurasia. The ancestral phylogenetic position of the West Beringian group and divergence time estimates support the hypothesis of continental range contraction to a single refugial area located in West Beringia during high-magnitude warming of the last interglacial, followed by westward recolonization of northern Eurasia in the last glacial period. The West Beringian group harbors the highest mitogenome diversity and its inferred demography indicates a constantly large effective population size over the Late Pleistocene to Holocene. This suggests that northward forest expansion during recent warming of the Holocene thermal optimum did not affect the gene pool of the collared lemming in West Beringia but reduced genomic diversity and effective population size in all other regions of the Eurasian Arctic. Demographic inference from genomic diversity was corroborated by species distribution modeling showing reduction in species distribution during past climate warming. These conclusions are supported by recent paleoecological evidence suggesting smaller temperature increases and moderate northward forest advances in the extreme northeast of Eurasia during the Late Pleistocene-to-Holocene warming events. This study emphasizes the importance of West Beringia as a potential refugium for cold-adapted Arctic species under ongoing climate warming.

SARS-CoV-2 epidemics quickly propagated worldwide, sorting virus genomic variants in newly established propagules of infections. Stochasticity in transmission within and between countries or an actual selective advantage could explain the global high frequency reached by some genomic variants. Using statistical analyses, demographic reconstructions, and molecular dynamics simulations, we show that the globally invasive G614 spike variant 1) underwent a significant demographic expansion in most countries explained neither by stochastic effects nor by overrepresentation in clinical samples, 2) increases the spike S1/S2 furin-like site conformational plasticity (short-range effect), and 3) modifies the internal motion of the receptor-binding domain affecting its cross-connection with other functional domains (long-range effect). Our results support the hypothesis of a selective advantage at the basis of the spread of the G614 variant, which we suggest may be due to structural modification of the spike protein at the S1/S2 proteolytic site, and provide structural information to guide the design of variant-specific drugs.