Assessing SNP-markers to study population mixing and ecological adaptation in Baltic cod

Diverse Transcriptome Responses to Salinity Change in Atlantic Cod Subpopulations

Adaptation to environmental variation caused by global climate change is a significant aspect of fisheries management and ecology. A reduction in ocean salinity is visible in near-shore areas, especially in the Baltic Sea, where it is affecting the Atlantic cod population. Cod is one of the most significant teleost species, with high ecological and economical value worldwide. The population of cod in the Baltic Sea has been traditionally divided into two subpopulations (western and eastern) existing in higher- and lower-salinity waters, respectively. In recent decades, both Baltic cod subpopulations have declined massively. One of the reasons for the poor condition of cod in the Baltic Sea is environmental factors, including salinity. Thus, in this study, an oligonucleotide microarray was applied to explore differences between Baltic cod subpopulations in response to salinity fluctuations. For this purpose, an exposure experiment was conducted consisting of salinity elevation and reduction, and gene expression was measured in gill tissue. We found 400 differentially expressed genes (DEGs) involved in the immune response, metabolism, programmed cell death, cytoskeleton, and extracellular matrix that showed a subpopulation-dependent pattern. These findings indicate that osmoregulation in Baltic cod is a complex process, and that western and eastern Baltic cod subpopulations respond differently to salinity changes.

Evidence of hybridization between genetically distinct Baltic cod stocks during peak population abundance(s)

Range expansions can lead to increased contact of divergent populations, thus increasing the potential of hybridization events. Whether viable hybrids are produced will most likely depend on the level of genomic divergence and associated genomic incompatibilities between the different entities as well as environmental conditions. By taking advantage of historical Baltic cod (Gadus morhua) otolith samples combined with genotyping and whole genome sequencing, we here investigate the genetic impact of the increased spawning stock biomass of the eastern Baltic cod stock in the mid 1980s. The eastern Baltic cod is genetically highly differentiated from the adjacent western Baltic cod and locally adapted to the brackish environmental conditions in the deeper Eastern basins of the Baltic Sea unsuitable for its marine counterparts. Our genotyping results show an increased proportion of eastern Baltic cod in western Baltic areas (Mecklenburg Bay and Arkona Basin)-indicative of a range expansion westwards-during the peak population abundance in the 1980s. Additionally, we detect high frequencies of potential hybrids (including F1, F2 and backcrosses), verified by whole genome sequencing data for a subset of individuals. Analysis of mitochondrial genomes further indicates directional gene flow from eastern Baltic cod males to western Baltic cod females. Our findings unravel that increased overlap in distribution can promote hybridization between highly divergent populations and that the hybrids can be viable and survive under specific and favourable environmental conditions. However, the observed hybridization had seemingly no long-lasting impact on the continuous separation and genetic differentiation between the unique Baltic cod stocks.

Living apart together: Long-term coexistence of Baltic cod stocks associated with depth-specific habitat use

Coexistence of fish populations (= stocks) of the same species is a common phenomenon. In the Baltic Sea, two genetically divergent stocks of Atlantic cod (Gadus morhua), Western Baltic cod (WBC) and Eastern Baltic cod (EBC), coexist in the Arkona Sea. Although the relative proportions of WBC and EBC in this area are considered in the current stock assessments, the mixing dynamics and ecological mechanisms underlying coexistence are not well understood. In this study, a genetically validated otolith shape analysis was used to develop the most comprehensive time series of annual stock mixing data (1977–2019) for WBC and EBC. Spatio-temporal mixing analysis confirmed that the two stocks coexist in the Arkona Sea, albeit with fluctuating mixing proportions over the 43-year observation period. Depth-stratified analysis revealed a strong correlation between capture depth and stock mixing patterns, with high proportions of WBC in shallower waters (48–61% in <20m) and increasing proportions of EBC in deeper waters (50–86% in 40-70m). Consistent depth-specific mixing patterns indicate stable differences in depth distribution and habitat use of WBC and EBC that may thus underlie the long-term coexistence of the two stocks in the Arkona Sea. These differences were also reflected in significantly different proportions of WBC and EBC in fisheries applying passive gears in shallower waters (more WBC) and active gears in deeper waters (more EBC). This highlights the potential for fishing gear-specific exploitation of different stocks, and calls for stronger consideration of capture depth and gear type in stock assessments. This novel evidence provides the basis for improved approaches to research, monitoring and management of Baltic cod stocks.

The effect of growth rate on otolith-based discrimination of cod (Gadus morhua) ecotypes

Otolith shape has previously been used to identify ecotypes within the Icelandic cod (Gadus morhua) stock, using DST profiles to validate the results. Fish otolith shape variation has repeatedly been found to be largely determined by growth rate. To examine the effect of growth rate on the relationship between otolith shape and cod ecotypes (using the Pan I genotype as a proxy for ecotype), 826 archived sagittal otoliths collected over a 58 year sampling period were retrieved, the individual growth rate calculated, and otolith shape described using both Normalized Elliptic Fourier transform and Discrete Wavelet transform. Discriminant functions of otolith shape successfully classified ecotype, whether using Fourier or Wavelet descriptors, but only when excluding a heterozygous genotype from the analysis. The otolith shape variability of this genotype lowered the classification success, while otolith shape, in turn, was significantly affected by growth rate and cohort. Growth rate differences previously reported for the ecotypes were present, but were less marked than expected and indeed, growth rate variance attributable to ecotype identity was dwarfed by cohort- and location-related variance in growth. Such a strong effect of growth rate suggests that cod ecotype discrimination based on otolith shape is sensitive to both temporal and spatial variations in growth, which can mask the effect of ecotype-related growth rate differences on otolith shape.

A star is torn-molecular analysis divides the Mediterranean population of Poli's stellate barnacle, Chthamalus stellatus (Cirripedia, Chtamalidae)

Poli's stellate barnacle, Chthamalus stellatus Poli, populates the Mediterranean Sea, the North-Eastern Atlantic coasts, and the offshore Eastern Atlantic islands. Previous studies have found apparent genetic differences between the Atlantic and the Mediterranean populations of C. stellatus, suggesting possible geological and oceanographic explanations for these differences. We have studied the genetic diversity of 14 populations spanning from the Eastern Atlantic to the Eastern Mediterranean, using two nuclear genes sequences revealing a total of 63 polymorphic sites. Both genotype-based, haplotype-based and the novel SNP distribution population-based methods have found that these populations represent a geographic cline along the west to east localities. The differences in SNP distribution among populations further separates a major western cluster into two smaller clusters, the Eastern Atlantic and the Western Mediterranean. It also separates the major eastern cluster into two smaller clusters, the Mid-Mediterranean and Eastern Mediterranean. We suggested here environmental conditions like surface currents, water salinity and temperature as probable factors that have formed the population structure. We demonstrate that C. stellatus is a suitable model organism for studying how geological events and hydrographic conditions shape the fauna in the Mediterranean Sea.

Tipping point realized in cod fishery

Understanding tipping point dynamics in harvested ecosystems is of crucial importance for sustainable resource management because ignoring their existence imperils social-ecological systems that depend on them. Fisheries collapses provide the best known examples for realizing tipping points with catastrophic ecological, economic and social consequences. However, present-day fisheries management systems still largely ignore the potential of their resources to exhibit such abrupt changes towards irreversible low productive states. Using a combination of statistical changepoint analysis and stochastic cusp modelling, here we show that Western Baltic cod is beyond such a tipping point caused by unsustainable exploitation levels that failed to account for changing environmental conditions. Furthermore, climate change stabilizes a novel and likely irreversible low productivity state of this fish stock that is not adapted to a fast warming environment. We hence argue that ignorance of non-linear resource dynamics has caused the demise of an economically and culturally important social-ecological system which calls for better adaptation of fisheries systems to climate change.

The forgotten feeding ground: patterns in seasonal and depth-specific food intake of adult cod Gadus morhua in the western Baltic Sea

This study presents the diet composition of western Baltic cod Gadus morhua based on 3150 stomachs sampled year-round between 2016 and 2017 using angling, gillnetting and bottom trawling, which enhanced the spatio-temporal coverage of cod habitats. Cod diet composition in shallow areas (<20 m depth) was dominated by benthic invertebrate species, mainly the common shore crab Carcinus maneas. Compared to historic diet data from the 1960s and 1980s (limited to depth >20 m), the contribution of herring Clupea harengus decreased and round goby Neogobius melanostomus occurred as a new prey species. Statistical modelling revealed significant relationships between diet composition, catch depth, fish length and season. Generalized additive modelling identified a negative relationship between catch depth and stomach content weight, suggesting reduced food intake in winter when cod use deeper areas for spawning and during peak summer when cod tend to avoid high water temperatures. The results of this study highlight the importance of shallow coastal areas as major feeding habitats of adult cod in the western Baltic Sea, which were previously unknown because samples were restricted to deeper trawlable areas. The results strongly suggest that historic stomach analyses overestimated the role of forage fish and underestimated the role of invertebrate prey. Eventually, this study shows the importance of a comprehensive habitat coverage for unbiased stomach sampling programmes to provide a more reliable estimation of top predator diet, a key information for food web analyses and multispecies models.

Identification of Arctic Food Fish Species for Anthropogenic Contaminant Testing Using Geography and Genetics

The identification of food fish bearing anthropogenic contaminants is one of many priorities for Indigenous peoples living in the Arctic. Mercury (Hg), arsenic (As), and persistent organic pollutants including polychlorinated biphenyls (PCBs) are of concern, and these are reported, in some cases for the first time, for fish sampled in and around King William Island, located in Nunavut, Canada. More than 500 salmonids, comprising Arctic char, lake trout, lake whitefish, and ciscoes, were assayed for contaminants. The studied species are anadromous, migrating to the ocean to feed in the summers and returning to freshwater before sea ice formation in the autumn. Assessments of muscle Hg levels in salmonids from fishing sites on King William Island showed generally higher levels than from mainland sites, with mean concentrations generally below guidelines, except for lake trout. In contrast, mainland fish showed higher means for As, including non-toxic arsenobetaine, than island fish. Lake trout were highest in As and PCB levels, with salmonid PCB congener analysis showing signatures consistent with the legacy of cold-war distant early warning stations. After DNA-profiling, only 4–32 Arctic char single nucleotide polymorphisms were needed for successful population assignment. These results support our objective to demonstrate that genomic tools could facilitate efficient and cost-effective cluster assignment for contaminant analysis during ocean residency. We further suggest that routine pollutant testing during the current period of dramatic climate change would be helpful to safeguard the wellbeing of Inuit who depend on these fish as a staple input to their diet. Moreover, this strategy should be applicable elsewhere.

Genetic analysis redraws the management boundaries for the European sprat

Sustainable fisheries management requires detailed knowledge of population genetic structure. The European sprat is an important commercial fish distributed from Morocco to the Arctic circle, Baltic, Mediterranean, and Black seas. Prior to 2018, annual catch advice on sprat from the International Council for the Exploration of the Sea (ICES) was based on five putative stocks: (a) North Sea, (b) Kattegat-Skagerrak and Norwegian fjords, (c) Baltic Sea, (d) West of Scotland-southern Celtic Seas, and (e) English Channel. However, there were concerns that the sprat advice on stock size estimates management plan inadequately reflected the underlying biological units. Here, we used ddRAD sequencing to develop 91 SNPs that were thereafter used to genotype approximately 2,500 fish from 40 locations. Three highly distinct and relatively homogenous genetic groups were identified: (a) Norwegian fjords; (b) Northeast Atlantic including the North Sea, Kattegat-Skagerrak, Celtic Sea, and Bay of Biscay; and (c) Baltic Sea. Evidence of genetic admixture and possibly physical mixing was detected in samples collected from the transition zone between the North and Baltic seas, but not between any of the other groups. These results have already been implemented by ICES with the decision to merge the North Sea and the Kattegat-Skagerrak sprat to be assessed as a single unit, thus demonstrating that genetic data can be rapidly absorbed to align harvest regimes and biological units.

SNP genotyping reveals substructuring in weakly differentiated populations of Atlantic cod (Gadus morhua) from diverse environments in the Baltic Sea

Atlantic cod (Gadus morhua) is one of the most important fish species in northern Europe for several reasons including its predator status in marine ecosystems, its historical role in fisheries, its potential in aquaculture and its strong public profile. However, due to over-exploitation in the North Atlantic and changes in the ecosystem, many cod populations have been reduced in size and genetic diversity. Cod populations in the Baltic Proper, Kattegat and North Sea have been analyzed using a species specific single nucleotide polymorphism (SNP) array. Using a subset of 8,706 SNPs, moderate genetic differences were found between subdivisions in three traditionally delineated cod management stocks: Kattegat, western and eastern Baltic. However, an FST measure of population differentiation based on allele frequencies from 588 outlier loci for 2 population groups, one including 5 western and the other 4 eastern Baltic populations, indicated high genetic differentiation. In this paper, differentiation has been demonstrated not only between, but also within western and eastern Baltic cod stocks for the first time, with salinity appearing to be the most important environmental factor influencing the maintenance of cod population divergence between the western and eastern Baltic Sea.

Effects of freezing on length and mass measurements of Atlantic cod Gadus morhua in the Baltic Sea

An aggregated sample of 925 Atlantic cod Gadus morhua collected by four countries in different regions of the Baltic Sea during different seasons were measured (total length, L_T = 161-890 mm and weighed (mass, M = 45-6900 g) both before freezing and after defrosting. The cod were found to decrease significantly in both L_T and M following death and frozen storage. There was an average (±SD) change in L_T of -2.91% (±0.05%) following freezing, independent of starting L_T . Total M changed by -2.65% (±0.14%), independent of starting mass. Shrinkage of L_T and M did not differ significantly between 1 and 4 months frozen storage, though L_T shrinkage was significantly greater after 1 or 4 months in the freezer compared with after 5 days. There was significant variation in L_T and M shrinkage between regions of capture. A significant negative relationship between condition of cod and L_T or M change was also observed. Equations to back-calculate fresh L_T and M from thawed L_T , M and standard length (L_S ), gutted L_T , gutted L_S and gutted M are provided.