Population genetic analysis of Euro-Arctic polar cod Boreogadus saida suggests fjord and oceanic structuring

Genomic architecture and population structure of Boreogadus saida from Canadian waters

The polar cod, Boreogadus saida, is an abundant and ubiquitous forage fish and a crucial link in Arctic marine trophic dynamics. Our objective was to unravel layers of genomic structure in B. saida from Canadian waters, specifically screening for potential hybridization with the Arctic cod, Arctogadus glacialis, large chromosomal inversions, and sex-linked regions, prior to interpreting population structure. Our analysis of 53,384 SNPs in 522 individuals revealed hybridization and introgression between A. glacialis and B. saida. Subsequent population level analyses of B. saida using 12,305 SNPs in 511 individuals revealed three large (ca. 7.4-16.1 Mbp) chromosomal inversions, and a 2 Mbp region featuring sex-linked loci. We showcase population structuring across the Western and Eastern North American Arctic, and subarctic regions ranging from the Hudson Bay to the Canadian Atlantic maritime provinces. Genomic signal for the inferred population structure was highly aggregated into a handful of SNPs (13.8%), pointing to potentially important adaptive evolution across the Canadian range. Our study provides a high-resolution perspective on the genomic structure of B. saida, providing a foundation for work that could be expanded to the entire circumpolar range for the species.

Thyasirid species composition (Bivalvia: Thyasiridae) and genetic connectivity of Parathyasira equalis (A. E. Verrill & K. J. Bush, 1898) in deep basins of sub-Arctic fjords

BACKGROUND

Thyasirid bivalves are often recorded as a dominant component of macrobenthic infaunal communities in depositional environments such as fjord basins. Fjord basins comprise patchy soft-bottom habitats bounded by steep walls and sills; however, little is known how this semi-isolated nature of fjords affects benthic populations. Accordingly, data on the composition and population connectivity of thyasirids can provide valuable information on the ecology of these ecosystems.

RESULTS

The species composition of thyasirid bivalves has been studied in the basins of three sub-Arctic fjords (Nordland, Northern Norway). Overall, six thyasirid species were recorded: Parathyasira equalis, Parathyasira dunbari, Mendicula ferruginosa, Genaxinus eumyarius, Thyasira sarsii, and Thyasira obsoleta. The species composition remained stable within the basins during the sampling period (2013-2020) and suggested the importance of local reproduction over advection of individuals for population dynamics. Only one species, Parathyasira equalis, was common in all fjords. We have further investigated the population genetics of this species by combining two types of genetic markers: a 579 bp fragment of the cytochrome c oxidase subunit I (COI) gene and 4043 single-nucleotide polymorphisms (SNPs) generated by genotyping-by-sequencing. The latter provided a more in-depth resolution on the population genetics of this species and revealed a weak but significant differentiation of populations within fjords, further indicating limited connectivity between basins.

CONCLUSION

Based on our findings, we conclude that limited dispersal between the basin communities results in weakly connected populations and might be an important structuring factor for macrobenthic communities.

An evolutionarily distinct ringed seal in the Ilulissat Icefjord

The Earth's polar regions are low rates of inter- and intraspecific diversification. An extreme mammalian example is the Arctic ringed seal (Pusa hispida hispida), which is assumed to be panmictic across its circumpolar Arctic range. Yet, local Inuit communities in Greenland and Canada recognize several regional variants; a finding supported by scientific studies of body size variation. It is however unclear whether this phenotypic variation reflects plasticity, morphs or distinct ecotypes. Here, we combine genomic, biologging and survey data, to document the existence of a unique ringed seal ecotype in the Ilulissat Icefjord (locally 'Kangia'), Greenland; a UNESCO World Heritage site, which is home to the most productive marine-terminating glacier in the Arctic. Genomic analyses reveal a divergence of Kangia ringed seals from other Arctic ringed seals about 240 kya, followed by secondary contact since the Last Glacial Maximum. Despite ongoing gene flow, multiple genomic regions appear under strong selection in Kangia ringed seals, including candidate genes associated with pelage coloration, growth and osmoregulation, potentially explaining the Kangia seal's phenotypic and behavioural uniqueness. The description of 'hidden' diversity and adaptations in yet another Arctic species merits a reassessment of the evolutionary processes that have shaped Arctic diversity and the traditional view of this region as an evolutionary freezer. Our study highlights the value of indigenous knowledge in guiding science and calls for efforts to identify distinct populations or ecotypes to understand how these might respond differently to environmental change.

Multiple exposure of the Boreogadus saida from bessel fjord (NE Greenland) to legacy and emerging pollutants

This work investigates the occurrence of OCPs, such as hexachlorocyclohexane (α-, β-, γ- and δ-HCH) isomers, dichlorodiphenyltrichloroethane (p,p'-DDT) and its metabolite dichlorodiphenyldichloroethylene (p,p'-DDE), endosulfan (α- and β-EDS) isomers, chlorpyrifos (CPF), dacthal (DAC) and phenolic compounds, such as 4-nonylphenol (4-NP) and its precursors nonylphenol polyethoxylates (NP1EO and NP2EO) and bisphenol A (BPA), in polar cod sampled in and outside Bessel Fjord (NE Greenland). Linear regressions between target contaminants and morphological parameters (age, length, weight, gonad- and hepato-somatic indices and Fulton K) have been also evaluated. Polar cod collected at shelf had higher average concentrations of BPA, NP1EO, NP2EO and 4-NP (muscle: 6.2, 13.2, 8.9 and 1.9 ng/g w.w., respectively; liver: 5.8, 7.5, 5.2 and 0.9 ng/g w.w. respectively), than fjord's specimens (muscle: 3.5, 9.1, 3.9 and 1.0 ng/g w.w., respectively; liver: 2.4, 5.3, 2.9 and 1.1 ng/g w.w. respectively). ΣHCHs, ΣEDSs, ΣDDTs, CPF and DAC, were more accumulated in the polar cod from the fjord (average amount in muscle: 9.1, 4.8, 7.9, 3.8 and 2.8 ng/g w.w., respectively; average amount in the liver: 11.2, 9.0, 3.8, 5.9 and 4.9 ng/g w.w., respectively) than shelf's ones (average amount in muscle 3.9, 4.5, 4.2, 0.9 and 1.2 ng/g w.w., respectively; average amount in liver 7.8, 6.3, 2.1, 3.4 and 2.5 ng/g w.w., respectively). The comparison between the concentration of target contaminants and morphologic parameters suggested a different exposure of polar cod occupying the fjord and shelf habitats, due to a combination of genetic and dietary differences, climate change effects and increased human activities.

Distinct genetic clustering in the weakly differentiated polar cod, Boreogadus saida Lepechin, 1774 from East Siberian Sea to Svalbard

The cold-adapted polar cod Boreogadus saida, a key species in Arctic ecosystems, is vulnerable to global warming and ice retreat. In this study, 1257 individuals sampled in 17 locations within the latitudinal range of 75–81°N from Svalbard to East Siberian Sea were genotyped with a dedicated suite of 116 single-nucleotide polymorphic loci (SNP). The overall pattern of isolation by distance (IBD) found was driven by the two easternmost samples (East Siberian Sea and Laptev Sea), whereas no differentiation was registered in the area between the Kara Sea and Svalbard. Eleven SNP under strong linkage disequilibrium, nine of which could be annotated to chromosome 2 in Atlantic cod, defined two genetic groups of distinct size, with the major cluster containing seven-fold larger number of individuals than the minor. No underlying geographic basis was evident, as both clusters were detected throughout all sampling sites in relatively similar proportions (i.e. individuals in the minor cluster ranging between 4 and 19% on the location basis). Similarly, females and males were also evenly distributed between clusters and age groups. A differentiation was, however, found regarding size at age: individuals belonging to the major cluster were significantly longer in the second year. This study contributes to increasing the population genetic knowledge of this species and suggests that an appropriate management should be ensured to safeguard its diversity.

High gene flow in polar cod (Boreogadus saida) from West-Svalbard and the Eurasian Basin

The current and projected environmental change of the Arctic Ocean contrasts sharply with the limited knowledge of its genetic biodiversity. Polar cod Boreogadus saida (Lepechin, 1774) is an abundant circumpolar marine fish and ecological key species. The central role of polar cod in the Arctic marine food web warrants a better understanding of its population structure and connectivity. In this study, the genetic population structure of 171 juveniles, collected from several fjords off West-Svalbard (Billefjorden, Hornsund and Kongsfjorden), the northern Sophia Basin and the Eurasian Basin of the Arctic Ocean, was analysed using nine DNA microsatellite loci. Genetic analyses indicated moderate to high genetic diversity, but absence of spatial population structure and isolation-by-distance, suggesting ongoing gene flow between the studied sampling regions. High levels of connectivity may be key for polar cod to maintain populations across wide spatial scales. The adaptive capacity of the species will be increasingly important to face challenges such as habitat fragmentation, ocean warming and changes in prey composition. In view of a limited understanding of the population dynamics and evolution of polar cod, a valuable next step to predict future developments should be an integrated biological evaluation, including population genomics, a life-history approach, and habitat and biophysical dispersal modelling.

Distribution and ecology of polar cod (Boreogadus saida) in the eastern Barents Sea: A review of historical literature

The polar cod (Boreogadus saida) has a circumpolar distribution and is the most abundant planktivorous fish in the Arctic. Declining sea-ice coverage impacts polar cod directly and also facilitates expansion of human activities in the region leading to increasing anthropogenic pressures on biota. Here we summarize current data and knowledge on polar cod from the Russian sector of the Barents Sea and discuss knowledge needs for the management of polar cod under changing environmental conditions and anthropogenic impacts. We review 36 Russian historical (1935 - 2020) sources of data and knowledge largely unknown to western researchers, in addition to sources already published in the English language. This effort allowed for digitalization and visualization of 69 separate datasets on polar cod ecology, including maturation, fertility, feeding intensity, diet, lipid content, length-weight relationships and seasonal variation in larval size. Our review suggests that polar cod abundances are particularly large in the eastern Barents Sea and adjacent waters. Here, we identify and discuss key knowledge gaps. The review of polar cod in the eastern Barents Sea revealed 1) major variation in the timing and area of polar cod spawning, 2) uncertainty as to what degree the polar cod is dependent on sea ice, 3) deficient knowledge of juvenile (e.g., 0-group) distributions, particularly in the north-eastern Barents Sea, 4) deficient knowledge of the species' genetic structure and spatio-temporal distributions, and 5) insufficient understanding as to whether ongoing environmental change may induce phenological changes affecting the availability of potential food items for polar cod larvae and their match in space and time. Filling these knowledge gaps would provide an important step towards the reliable knowledge base needed in order to perform well-founded management and impact assessment under environmental changes and increasing anthropogenic impacts.

Seasonal shifts in feeding patterns: Individual and population realized specialization in a high Arctic fish

Species with a broad and flexible diet may be at an advantage in a rapidly changing environment such as in today's Arctic ecosystems. Polar cod (Boreogadus saida), an abundant and ecologically important circumpolar Arctic fish, is often described as a zooplankton generalist feeder, which suggests that it may cope successfully with changes in prey composition. This description is justified based on the relatively broad diet of polar cod across sites and seasons. In this case study, we used polar cod dietary data from fall and winter and from two distinct environments, dominated either by Arctic or Atlantic water masses in Svalbard. Our results point to the importance of time and space when drawing conclusions on dietary plasticity and degree of specialization. Polar cod diet differed significantly between fall and the winter and between Arctic and Atlantic domains. Polar cod from Arctic domains displayed a strong realized population specialization on Themisto libellula in fall, and the larger dietary niche width observed in the winter was the product of realized individual specialization, with increased feeding on fish prey. Overall, we did not observe a generalized feeding behavior. If dietary niche width is to inform conservation management, we argue it must be recognized that populations from a single species may adopt seasonally contrasting degrees of dietary specialization and that these populations may differ in their vulnerability to climate-induced changes in prey community composition.