Gene flow erodes spatial genomic clines while selection preserves adaptive loci over hundreds of generations

Abstract

Spatial genetic clines suggest local adaptation, but whether they persist or erode over evolutionary time remains largely unknown. Here, we use a unique genomic time series collected over several decades to quantify how gene flow and selection jointly shape spatial genomic structure. We analyzed pooled whole-genome sequencing from 16 natural populations of Drosophila melanogaster sampled along the North American east coast between 1997 and 2023, spanning hundreds of generations. We find pronounced homogenization of populations over time, accompanied by a genome-wide loss of clinal variants. Polymorphisms that remained clinal were less steep on chromosome 3R, reflecting declines in the frequency and clinality of major chromosomal inversions. Genomic context modulates this process: clinal SNPs are preferentially lost in regions of high recombination, whereas remaining clines are increasingly concentrated in coding and splice-associated sites. Despite widespread erosion of spatial structure, we detect signatures of selection shared across space at detoxification loci, including cytochrome P450 genes linked to insecticide resistance. Together, our results show widespread homogenization of neutral spatial variation consistent with gene flow, while spatially varying selection has preserved a subset of functionally important loci.