Hidden genomic structure and widespread structural polymorphism across environmental gradients in the spiny sea star Marthasterias glacialis

Genomic regions of reduced recombination can preserve linkage among co-adapted alleles, facilitating local adaptation despite high connectivity. Such regions-often generated by chromosomal inversions-may be especially important in highly dispersive marine taxa yet remain poorly documented in echinoderms. Here, we combined a chromosome-level reference genome with genome-wide ddRAD-seq from 296 Marthasterias glacialis individuals across 19 Atlantic-Mediterranean locations to quantify population structure and scan for recombination-suppressed haploblocks. Genome-wide neutral markers showed significant population differentiation together with evidence of high connectivity, revealed by the presence of inter-ecoregion migrants. Additionally, we identified 16 polymorphic haploblocks with patterns consistent with putative chromosomal inversions spanning 18.6% of the genome. Haploblock haplotypes were strongly environmentally and geographically structured and contained genes with key functions in stress response, osmoregulation and thermal tolerance. Haplotype distributions also paralleled previously described mitochondrial lineages despite nuclear gene flow, consistent with a model of ancient divergence followed by secondary contact. Overall, our results suggest a role for widespread structural polymorphism in adaptive differentiation in Echinodermata, providing a framework for linking echinoderm genome rearrangements to ecological divergence. Marthasterias glacialis thus emerges as a promising system to explore how structural variation contributes to adaptation and genome evolution in highly dispersive organisms.