Biotic and abiotic factors shaping the genome of cockle (Cerastoderma edule) in the Northeast Atlantic: a baseline for sustainable management of its wild resources

Knowledge on how environmental factors shape the genome of marine species is crucial for sustainable management of fisheries and wild populations. The edible cockle (Cerastoderma edule) is a marine bivalve distributed along the Northeast Atlantic coast of Europe and is an important resource from both commercial and ecological perspectives. We performed a population genomics screening using 2b-RAD genotyping on 9,309 SNPs localised in the cockles genome on a sample of 536 specimens pertaining to 14 beds in the Northeast Atlantic to ascertain its genetic structure regarding environmental variation. Larval dispersal modelling considering species behaviour and interannual variability in ocean conditions was carried out, as an essential background to compare genetic information with. Cockle populations in the Northeast Atlantic were shown to be panmictic and displayed low but significant geographical differentiation across populations (FST = 0.0240; P < 0.001), albeit not across generations. We identified 441 outlier SNPs related to divergent selection, sea surface temperature being the main environmental driver following a latitudinal axis. Two main genetic groups were identified, northwards and southwards of French Brittany, in accordance with our modelling, which demonstrated a barrier for larval dispersal linked to the Ushant front. Further genetic subdivision was observed using outlier loci and considering larval behaviour. The northern group was divided into the Irish/Celtic Seas and the English Channel/North Sea, while the southern group was divided into three subgroups. This information represents the baseline for management of cockles, designing conservation strategies, founding broodstock for depleted beds, and producing suitable seed for aquaculture production.