Abstract
Investigating gene flow between closely related species and its variation across the genome is important to understand how reproductive barriers shape genome divergence before speciation is complete. An efficient way to characterize differential gene flow is to study how the genetic interactions that take place in hybrid zones selectively filter gene exchange between species, leading to heterogeneous genome divergence. In the present study, genome-wide divergence and introgression patterns were investigated between two sole species, Solea senegalensis and Solea aegyptiaca, using restriction-associated DNA sequencing (RAD-Seq) to analyze samples taken from a transect spanning the hybrid zone. An integrative approach combining geographic and genomic clines methods with an analysis of individual locus introgression accounting for the demographic history of divergence was conducted. Our results showed that the two sole species have come into secondary contact postglacially, after experiencing a prolonged period (ca. 1.1 to 1.8 Myrs) of allopatric separation. Secondary contact resulted in the formation of a tension zone characterized by strong reproductive isolation, which only allowed introgression in a limited fraction of the genome. We found multiple evidence for a preferential direction of introgression in the S. aegyptiaca genetic background, indicating a possible recent or ongoing movement of the hybrid zone. Deviant introgression signals found in the opposite direction suggested that S. senegalensis could have possibly undergone adaptive introgression that has not yet spread throughout the entire species range. Our study thus illustrates the varied outcomes of genetic interactions between divergent gene pools that recently met after a long history of divergence.