The planktonic pteropod Limacina helicina is increasingly studied as a bioindicator of climate-driven changes in polar marine ecosystems. Although broadly distributed across the Arctic Basin and the North Pacific, its population structure and dispersal pathways remain poorly understood, especially in the Siberian Arctic. We analyzed mitochondrial COI sequences from populations sampled in the Barents, Kara, Laptev, East Siberian, and White Seas, as well as adjacent Pacific regions. Three major haplogroups (H1, H2, H3) were identified with distinct spatial patterns. H1 is widespread, occurring across the Pacific and most Arctic seas except the White Sea. H2 is confined to the western Arctic shelves (Barents–Kara–Laptev), and H3 is unique to the White Sea. We found a pronounced genetic discontinuity corresponding to hydrographic barriers, particularly the strong freshwater inflow from the Lena River, which restricts eastward dispersal of H2 from the Laptev to the East Siberian Sea. These patterns suggest postglacial expansions from geographically separated populations that survived the Last Glacial Maximum in isolated marine regions. The White Sea population is highly isolated and genetically distinct. Our results highlight how both glacial history and modern oceanography shape Arctic plankton diversity and define biogeographic boundaries in a rapidly changing climate.

Plankton communities often consist of cosmopolitan species with an extensive gene flow between populations. Nevertheless, populations of some plankton species are genetically structured, owing to various barriers such as ocean currents, hydrological fronts, and continents. Drivers that could explain the genetic structures of most mesopelagic species remain unknown on an ocean-basin scale, and our study aims to analyze the genetic and morphological differences between populations of a cosmopolitan mesopelagic shrimp, Systellaspis debilis, from the Southern and Northern Atlantic Ocean, and the Southwest Indian Ocean. We analyzed the ITS-1 and COI markers of 75 specimens and assessed the genetic integrity and within-species variability of these genes. We also coded 32 morphological characteristics in 73 specimens, analyzed their variability, and assessed the correlation between morphological and genetic characteristics using a Redundancy analysis and Mantel test. Systellaspis debilis was genetically cohesive across the whole Atlantic and Southwest Indian Oceans, which is possibly a result of an intensive gene flow through ecological barriers, the resistance of species to hydrological gradients, a purifying selection of mitochondrial genes, etc. In contrast, we found significant morphological differences between populations from different regions, which mirrors morphological diversification and calls for further genomic approaches in order to understand the basis of these variations and uncover potential local adaptations.