The origin and remolding of genomic islands of differentiation in the European sea bass

Speciation is a complex process that leads to the progressive establishment of reproductive isolation barriers between diverging populations. Genome-wide comparisons between closely related species have revealed the existence of heterogeneous divergence patterns, dominated by genomic islands of increased divergence supposed to contain reproductive isolation loci. However, this divergence landscape only provides a static picture of the dynamic process of speciation, during which confounding mechanisms unlinked to speciation can interfere. Here, we used haplotype-resolved whole-genome sequences to identify the mechanisms responsible for the formation of genomic islands between Atlantic and Mediterranean sea bass lineages. We show that genomic islands first emerged in allopatry through the effect of linked selection acting on a heterogeneous recombination landscape. Upon secondary contact, preexisting islands were strongly remolded by differential introgression, revealing variable fitness effects among regions involved in reproductive isolation. Interestingly, we found that divergent regions containing ancient polymorphisms conferred the strongest resistance to introgression.