Genomic signatures associated with population divergence, speciation and the evolutionary mechanisms responsible for these are key research topics in evolutionary biology.Evolutionary radiations and parallel evolution have offered opportunities to study the role of the environment by providing replicates of ecologically driven speciation. Here, we apply an extension of the parallel evolution framework to study replicates of ecological speciation where multiple species went through a process of population divergence during the colonization of a common environmental gradient. We used the conditions offered by the North Sea -Baltic Sea environmental transition zone and found clear evidence of population structure linked to the Baltic Sea salinity gradient in four flatfish species. We found highly heterogeneous signatures of population divergence within and between species, and no evidence of parallel genomic architecture across species associated with the divergence.Analyses of demographic history suggest that Baltic Sea lineages are older than the age of the Baltic Sea itself. In most cases, divergence appears to involve reticulated demography through secondary contact, and our analyses revealed that genomic patterns of divergence were likely the result of a combination of effects from past isolation and subsequent adaptation to a new environment. In one case, we identified two large structural variants associated with the environmental gradient, where populations were inferred to have diverged in the presence of gene flow. Our results highlight the heterogeneous genomic effects associated with complex interplays of evolutionary forces, and stress the importance of genomic background for studies of parallel evolution. 2006), but so far the comparative framework has not been extended with the use of higher genome coverage.Here, we used a comparative population genomic approach to study the evolutionary processes involved during the colonization of the Baltic Sea. We studied population structure of four closely related (but genetically isolated) flatfish species: turbot (Scophthalmus maximus), common dab (Limanda limanda), European plaice (Pleuronectes platessa) andEuropean flounder (Platichthys flesus). These species were considered replicates of population divergence and were analysed independently. Our analyses showed heterogeneous genetic patterns of North Sea (NS) -Baltic Sea (BS) differentiation across species, both for magnitudes and genomic patterns of differentiation. Most of the genetic differences found were strongly associated with the environmental gradient of the NBTZ.However, the majority of the genomic regions underlying the NS-BS divergence were not shared between species. The timing of divergence between NS-BS populations was inferred to be five to ten times older than the age of the Baltic Sea. In three species, population divergence was inferred to involve past isolation followed by a secondary contact phase, and our analyses revealed that highly diverged genomic regions were likely the result of ...