Population divergence manifested by genomic rearrangements in a keystone Arctic species with high gene flow

Genomic rearrangements have in recent years gained attention due to their evolutionary role in processes related to adaptation to local environmental conditions as well as diversification and speciation. In this study, we report on genomic rearrangements in the cold-water adapted polar cod (Boreogadus saida), a keystone Arctic fish species. By taking advantage of a new chromosome-level genome assembly in combination with whole-genome population sequencing data from specimens across the northern Barents Sea and adjacent regions, we identified a substantial number of larger chromosomal inversions (n=20) and characterized the previously identified chromosomal fusions (n=5). These genomic features -- encompassing over 20% of the genome -- exhibited genetic divergence, strong internal linkage disequilibrium, and signals of selection. Two of the identified inversions were associated with the two previously described hemoglobin clusters, while a third chromosomal region was found to differentiate between males and females. Moreover, clustering analyses on genotype frequencies of inversions revealed sub- structuring according to five geographic sub-groups suggesting sub-populations and/or the existence of cryptic ecotypes. These results provide novel insights into the impact of genomic rearrangements in population divergence and thus, potentially local adaptation, especially in species with high gene flow.