Genomic analysis reveals neutral and adaptive patterns that challenge the current management regime for East Atlantic cod Gadus morhua L

Putting Structural Variants Into Practice: The Role of Chromosomal Inversions in the Management of Marine Environments

Major threats to marine species and ecosystems include overfishing, invasive species, pollution and climate change. The changing climate not only imposes direct threats through the impacts of severe marine heatwaves, cyclones and ocean acidification but also complicates fisheries and invasive species management by driving species range shifts. The dynamic nature of these threats means that the future of our oceans will depend on the ability of species to adapt. This has led to calls for genetic interventions focussed on enhancing species' adaptive capacity, including translocations, restocking and selective breeding. Assessing the benefits and risks of such approaches requires an improved understanding of the genetic architecture of adaptive variation, not only in relation to climate-resilient phenotypes but also locally adapted populations and the fitness of hybrids. Large structural genetic variants such as chromosomal inversions play an important role in local adaptation by linking multiple adaptive loci. Consequently, inversions are likely to be particularly important when managing for adaptive capacity. However, under some circumstances, they also accumulate deleterious mutations, potentially increasing the risk of inbreeding depression. Genetic management that takes account of these dual roles on fitness is likely to be more effective at ensuring population persistence. We summarise evolutionary factors influencing adaptive and deleterious variation of inversions, review inversions found in marine taxa, and provide a framework to predict the consequences of ignoring inversions in key management scenarios. We conclude by describing practical methods to bridge the gap between evolutionary theory and practical application of inversions in conservation.

First Report of 'Soft Flesh' Induced by the Parasite Kudoa thyrsites (Myxosporea) in Commercial Codfish From Norway

Kudoa thyrsites is a myxosporean parasite that infects the skeletal muscle of various teleost fish species globally. Severe infections lead to 'soft flesh' in fish fillets, resulting in food spoilage and subsequent discard. While K. thyrsites has previously been identified in migratory Atlantic mackerel in the northern Northeast Atlantic Ocean, it has not been observed in resident or farmed fish species in Norwegian waters until now. This study presents the first evidence of K. thyrsites infection and the associated 'soft flesh' condition in resident commercially important gadoid species from Norwegian waters, including Norwegian coastal cod (NCC), Northeast Arctic cod (NEA) and tusk. Molecular analyses confirmed the parasitic infection in 'soft flesh'-affected fish sampled from multiple coastal locations in Norway. The life cycle of Kudoa remains unknown but likely involves an alternating annelid host as in other myxosporeans. These findings in resident hosts suggest that the parasite completes its life cycle also at higher latitudes, in northern Norway. Consequently, there is a risk for the Norwegian fishing industry, as the effect of the parasite on fish fillet texture can occasionally occur and impact both consumer acceptance and industry revenues.

Reduced-Representation Sequencing Detects Trans-Arctic Connectivity and Local Adaptation in Polar Cod (Boreogadus saida)

Information on connectivity and genetic structure of marine organisms remains sparse in frontier ecosystems such as the Arctic Ocean. Filling these knowledge gaps becomes increasingly urgent, as the Arctic is undergoing rapid physical, ecological and socio-economic changes. The abundant and widely distributed polar cod (Boreogadus saida) is highly adapted to Arctic waters, and its larvae and juveniles live in close association with sea ice. Through a reduced-representation sequencing approach, this study explored the spatial genetic structure of polar cod at a circum-Arctic scale. Genomic variation was partitioned into neutral and adaptive components to respectively investigate genetic connectivity and local adaptation. Based on 922 high-quality single nucleotide polymorphism (SNP) markers genotyped in 611 polar cod, broad-scale differentiation was detected among three groups: (i) Beaufort -Chukchi seas, (ii) all regions connected by the Transpolar Drift, ranging from the Laptev Sea to Iceland, including the European Arctic and (iii) West Greenland. Patterns of neutral genetic structure suggested broadscale oceanographic and sea ice drift features (i.e., Beaufort Gyre and Transpolar Drift) as important drivers of connectivity. Genomic variation at 35 outlier loci indicated adaptive divergence of the West Greenland and the Beaufort-Chukchi Seas populations, possibly driven by environmental conditions. Sea ice decline and changing ocean currents can alter or disrupt connectivity between polar cod from the three genetic groups, potentially undermining their resilience to climate change, even in putative refugia, such as the Central Arctic Ocean and the Arctic Archipelago.

Haemoglobin revisited: delineating population structure with the world's first molecular genetic marker used in fisheries research

When haemoglobin genotyping was implemented in the early 1960s to investigate population genetic structure in Atlantic cod (Gadus morhua), it became one of the first molecular genetic markers deployed in fisheries research worldwide. However, its suitability was questioned due to its potential for selection. While the issue of neutrality concerned the first population geneticists, markers under selection are now routinely used to study population genetic structure. Here, we revisited haemoglobin genotyping half a decade later to analyse >6000 mature Atlantic cod from 73 spawning locations throughout Norway's approximately 2500 km coastline. A latitudinal gradient in allele frequencies, with a decrease in the HbI-2 allele towards the south, was observed. Our observed HbI-2 frequencies were consistently slightly lower than data from the 1960s, potentially reflecting adaptive changes to increasing sea temperatures. However, despite this difference, the observed north-south pattern in allele frequencies observed here and in the historical studies overlapped, aligning with current knowledge of population genetic structure in this species. We therefore conclude that this once questioned marker, which provided the first molecular genetic insights into genetic structure in Atlantic cod, provides knowledge consistent with the isolation by distance pattern revealed through decades of research in this species in this region.

Genetics in the Ocean's Twilight Zone: Population Structure of the Glacier Lanternfish Across Its Distribution Range

The mesopelagic zone represents one of the few habitats that remains relatively untouched from anthropogenic activities. Among the many species inhabiting the north Atlantic mesopelagic zone, glacier lanternfish (Benthosema glaciale) is the most abundant and widely distributed. This species has been regarded as a potential target for a dedicated fishery despite the scarce knowledge of its population genetic structure. Here, we investigated its genetic structure across the North Atlantic and into the Mediterranean Sea using 121 SNPs, which revealed strong differentiation among three main groups: the Mediterranean Sea, oceanic samples, and Norwegian fjords. The Mediterranean samples displayed less than half the genetic variation of the remaining ones. Very weak or nearly absent genetic structure was detected among geographically distinct oceanic samples across the North Atlantic, which contrasts with the low motility of the species. In contrast, a longitudinal gradient of differentiation was observed in the Mediterranean Sea, where genetic connectivity is known to be strongly shaped by oceanographic processes such as current patterns and oceanographic discontinuities. In addition, 12 of the SNPs, in linkage disequilibrium, drove a three clusters' pattern detectable through Principal Component Analysis biplot matching the genetic signatures generally associated with large chromosomal rearrangements, such as inversions. The arrangement of this putative inversion showed frequency differences between open-ocean and more confined water bodies such as the fjords and the Mediterranean, as it was fixed in the latter for the second most common arrangement of the fjord's samples. However, whether genetic differentiation was driven by local adaptation, secondary contact, or a combination of both factors remains undetermined. The major finding of this study is that B. glaciale in the North Atlantic-Mediterranean is divided into three major genetic units, information that should be combined with demographic properties to outline the management of this species prior to any eventual fishery attempt.

Tracing 600 years of long-distance Atlantic cod trade in medieval and post-medieval Oslo using stable isotopes and ancient DNA

Marine resources have been important for the survival and economic development of coastal human communities across northern Europe for millennia. Knowledge of the origin of such historic resources can provide key insights into fishing practices and the spatial extent of trade networks. Here, we combine ancient DNA and stable isotopes (δ13C, δ15N, non-exchangeable δ2H and δ34S) to investigate the geographical origin of archaeological cod remains in Oslo from the eleventh to seventeenth centuries CE. Our findings provide genetic evidence that Atlantic cod was obtained from different geographical populations, including a variety of distant-water populations like northern Norway and possibly Iceland. Evidence for such long-distance cod trade is already observed from the eleventh century, contrasting with archaeological and historical evidence from Britain and other areas of Continental Europe around the North and Baltic Seas, where such trade increased during the thirteenth to fourteenth centuries. The genomic assignments of specimens to different populations coincide with significantly different δ13C values between those same specimens, indicating that multiple Atlantic cod populations living in different environments were exploited. This research provides novel information about the exploitation timeline of specific Atlantic cod stocks and highlights the utility of combining ancient DNA (aDNA) methods and stable isotope analysis to describe the development of medieval and post-medieval marine fisheries.

Integration of genomic and ecological methods inform management of an undescribed, yet highly exploited, sardine species

Assessing genetic diversity within species is key for conservation strategies in the context of human-induced biotic changes. This is important in marine systems, where many species remain undescribed while being overfished, and conflicts between resource-users and conservation agencies are common. Combining niche modelling with population genomics can contribute to resolving those conflicts by identifying management units and understanding how past climatic cycles resulted in current patterns of genetic diversity. We addressed these issues on an undescribed but already overexploited species of sardine of the genus Harengula. We find that the species distribution is determined by salinity and depth, with a continuous distribution along the Brazilian mainland and two disconnected oceanic archipelagos. Genomic data indicate that such biogeographic barriers are associated with two divergent intraspecific lineages. Changes in habitat availability during the last glacial cycle led to different demographic histories among stocks. One coastal population experienced a 3.6-fold expansion, whereas an island-associated population contracted 3-fold, relative to the size of the ancestral population. Our results indicate that the island population should be managed separately from the coastal population, and that a Marine Protected Area covering part of the island population distribution can support the viability of this lineage.

Adaptation to seasonal reproduction and environment-associated factors drive temporal and spatial differentiation in northwest Atlantic herring despite gene flow

Understanding how marine organisms adapt to local environments is crucial for predicting how populations will respond to global climate change. The genomic basis, environmental factors and evolutionary processes involved in local adaptation are however not well understood. Here we use Atlantic herring, an abundant, migratory and widely distributed marine fish with substantial genomic resources, as a model organism to evaluate local adaptation. We examined genomic variation and its correlation with environmental variables across a broad environmental gradient, for 15 spawning aggregations in Atlantic Canada and the United States. We then compared our results with available genomic data of northeast Atlantic populations. We confirmed that population structure lies in a fraction of the genome including likely adaptive genetic variants of functional importance. We discovered 10 highly differentiated genomic regions distributed across four chromosomes. Nine regions show strong association with seasonal reproduction. One region, corresponding to a known inversion on chromosome 12, underlies a latitudinal pattern discriminating populations north and south of a biogeographic transition zone on the Scotian Shelf. Genome-environment associations indicate that winter seawater temperature best correlates with the latitudinal pattern of this inversion. The variation at two so-called 'islands of divergence' related to seasonal reproduction appear to be private to the northwest Atlantic. Populations in the northwest and northeast Atlantic share variation at four of these divergent regions, simultaneously displaying significant diversity in haplotype composition at another four regions, which includes an undescribed structural variant approximately 7.7 Mb long on chromosome 8. Our results suggest that the timing and geographic location of spawning and early development may be under diverse selective pressures related to allelic fitness across environments. Our study highlights the role of genomic architecture, ancestral haplotypes and selection in maintaining adaptive divergence in species with large population sizes and presumably high gene flow.

Ancient DNA evidence for the ecological globalization of cod fishing in medieval and post-medieval Europe

Understanding the historical emergence and growth of long-range fisheries can provide fundamental insights into the timing of ecological impacts and the development of coastal communities during the last millennium. Whole-genome sequencing approaches can improve such understanding by determining the origin of archaeological fish specimens that may have been obtained from historic trade or distant water. Here, we used genome-wide data to individually infer the biological source of 37 ancient Atlantic cod specimens (ca 1050-1950 CE) from England and Spain. Our findings provide novel genetic evidence that eleventh- to twelfth-century specimens from London were predominantly obtained from nearby populations, while thirteenth- to fourteenth-century specimens were derived from distant sources. Our results further suggest that Icelandic cod was indeed exported to London earlier than previously reported. Our observations confirm the chronology and geography of the trans-Atlantic cod trade from Newfoundland to Spain starting by the early sixteenth century. Our findings demonstrate the utility of whole-genome sequencing and ancient DNA approaches to describe the globalization of marine fisheries and increase our understanding regarding the extent of the North Atlantic fish trade and long-range fisheries in medieval and early modern times.

Geographic variation in gene flow from a genetically distinct migratory ecotype drives population genetic structure of coastal Atlantic cod (Gadus morhua L.)

Identifying how physical and biotic factors shape genetic connectivity among populations in time and space is essential to our understanding of the evolutionary trajectory as well as the management of marine species. Atlantic cod is a widespread and commercially important marine species displaying several ecotypes with different life history strategies. Using three sets of SNPs: neutral, informative, and genome-inversion linked, we studied population genetic structure of ~2500 coastal Atlantic cod (CC) from 40 locations along Norway's 2500 km coastline, including nine fjords. We observed: (1) a genetic cline, suggesting a mechanism of isolation by distance, characterized by a declining F ST between CC and North East Arctic Cod (NEAC-genetically distinct migratory ecotype) with increasing latitude, (2) that in the north, samples of CC from outer-fjord areas were genetically more similar to NEAC than were samples of CC from their corresponding inner-fjord areas, (3) greater population genetic differentiation among CC sampled from outer-fjord areas along the coast, than among CC sampled from their corresponding inner-fjord areas, (4) genetic differentiation among samples of CC from both within and among fjords. Collectively, these results permit us to draw two main conclusions. First, that differences in the relative presence of the genetically highly distinct, migratory ecotype NEAC, declining from north to south and from outer to inner fjord, plays the major role in driving population genetic structure of the Norwegian CC. Second, that there is limited connectivity between CC from different fjords. These results suggest that the current management units implemented for this species in Norway should be divided into smaller entities. Furthermore, the situation where introgression from one ecotype drives population genetic structure of another, as is the case here, may exist in other species and geographical regions, thus creating additional challenges for sustainable fisheries management.

Phylogenomics and species delimitation of the economically important Black Basses (Micropterus)

Informed management and conservation efforts are vital to sustainable recreational fishing and biodiversity conservation. Because the taxonomic rank of species is important in conservation and management strategies, success of these efforts depends on accurate species delimitation. The Black Basses (Micropterus) are an iconic lineage of freshwater fishes that include some of the world’s most popular species for recreational fishing and world's most invasive species. Despite their popularity, previous studies to delimit species and lineages in Micropterus suffer from insufficient geographic coverage and uninformative molecular markers. Our phylogenomic analyses of ddRAD data result in the delimitation of 19 species of Micropterus, which includes 14 described species, the undescribed but well-known Altamaha, Bartram’s, and Choctaw basses, and two additional undescribed species currently classified as Smallmouth Bass (M. dolomieu). We provide a revised delimitation of species in the Largemouth Bass complex that necessitates a change in scientific nomenclature: Micropterus salmoides is retained for the Florida Bass and Micropterus nigricans is elevated from synonymy for the Largemouth Bass. The new understanding of diversity, distribution, and systematics of Black Basses will serve as important basis for the management and conservation of this charismatic and economically important clade of fishes.

Ancient DNA reveals a southern presence of the Northeast Arctic cod during the Holocene

Climate change has been implicated in an increased number of distributional shifts of marine species during the last century. Nonetheless, it is unclear whether earlier climatic fluctuations had similar impacts. We use ancient DNA to investigate the long-term spawning distribution of the Northeast Arctic cod (skrei) which performs yearly migrations from the Barents Sea towards spawning grounds along the Norwegian coast. The distribution of these spawning grounds has shifted northwards during the last century, which is thought to be associated with food availability and warming temperatures. We genetically identify skrei specimens from Ruskeneset in west Norway, an archaeological site located south of their current spawning range. Remarkably, ¹⁴C analyses date these specimens to the late Holocene, when temperatures were warmer than present-day conditions. Our results either suggest that temperature is not the only driver influencing the spawning distribution of Atlantic cod, or could be indicative of uncertainty in palaeoclimate reconstructions in this region. Regardless, our findings highlight the utility of aDNA to reconstruct the historical distribution of economically important fish populations and reveal the complexity of long-term ecological interactions in the marine environment.

Stabilizing selection on Atlantic cod supergenes through a millennium of extensive exploitation

Ecological disruption due to human impacts is evident worldwide, and a key to mitigation lies in characterizing the underlying mechanisms of species and ecosystem stability. Here we show that three extensive “supergenes” are maintained in Atlantic cod by stabilizing selection, tying these genes to the persistence of a keystone species distributed across the northern Atlantic Ocean. Removal of this species has caused severe ecosystem reshuffling in several areas of its range. Genomic inference of historic stock sizes further shows that cod has been under pressure in the North Sea system since the Viking period, in line with zooarchaeological records. Expansion of fisheries in Northern Europe through the past millennium is well documented and supports the inferred long-term declines.

Life on Earth has been characterized by recurring cycles of ecological stasis and disruption, relating biological eras to geological and climatic transitions through the history of our planet. Due to the increasing degree of ecological abruption caused by human influences many advocate that we now have entered the geological era of the Anthropocene, or “the age of man.” Considering the ongoing mass extinction and ecosystem reshuffling observed worldwide, a better understanding of the drivers of ecological stasis will be a requisite for identifying routes of intervention and mitigation. Ecosystem stability may rely on one or a few keystone species, and the loss of such species could potentially have detrimental effects. The Atlantic cod (Gadus morhua) has historically been highly abundant and is considered a keystone species in ecosystems of the northern Atlantic Ocean. Collapses of cod stocks have been observed on both sides of the Atlantic and reported to have detrimental effects that include vast ecosystem reshuffling. By whole-genome resequencing we demonstrate that stabilizing selection maintains three extensive “supergenes” in Atlantic cod, linking these genes to species persistence and ecological stasis. Genomic inference of historic effective population sizes shows continued declines for cod in the North Sea–Skagerrak–Kattegat system through the past millennia, consistent with an early onset of the marine Anthropocene through industrialization and commercialization of fisheries throughout the medieval period.

Supergene origin and maintenance in Atlantic cod

Supergenes are sets of genes that are inherited as a single marker and encode complex phenotypes through their joint action. They are identified in an increasing number of organisms, yet their origins and evolution remain enigmatic. In Atlantic cod, four megabase-scale supergenes have been identified and linked to migratory lifestyle and environmental adaptations. Here we investigate the origin and maintenance of these four supergenes through analysis of whole-genome-sequencing data, including a new long-read-based genome assembly for a non-migratory Atlantic cod individual. We corroborate the finding that chromosomal inversions underlie all four supergenes, and we show that they originated at different times between 0.40 and 1.66 million years ago. We reveal gene flux between supergene haplotypes where migratory and stationary Atlantic cod co-occur and conclude that this gene flux is driven by gene conversion, on the basis of an increase in GC content in exchanged sites. Additionally, we find evidence for double crossover between supergene haplotypes, leading to the exchange of an ~275 kilobase fragment with genes potentially involved in adaptation to low salinity in the Baltic Sea. Our results suggest that supergenes can be maintained over long timescales in the same way as hybridizing species, through the selective purging of introduced genetic variation.

Atlantic cod carries four supergenes linked to migratory lifestyle and environmental adaptations. Using whole-genome sequencing, the authors show that the genome inversions that underlie the supergenes originated at different times and show gene flux between supergene haplotypes.

Population and seascape genomics of a critically endangered benthic elasmobranch, the blue skate Dipturus batis

The blue skate (Dipturus batis) has a patchy distribution across the North-East Atlantic Ocean, largely restricted to occidental seas around the British Isles following fisheries-induced population declines and extirpations. The viability of remnant populations remains uncertain and could be impacted by continued fishing and by-catch pressure, and the projected impacts of climate change. We genotyped 503 samples of D. batis, obtained opportunistically from the widest available geographic range, across 6 350 single nucleotide polymorphisms (SNPs) using a reduced-representation sequencing approach. Genotypes were used to assess the species' contemporary population structure, estimate effective population sizes and identify putative signals of selection in relation to environmental variables using a seascape genomics approach. We identified genetic discontinuities between inshore (British Isles) and offshore (Rockall and Faroe Island) populations, with differentiation most pronounced across the deep waters of the Rockall Trough. Effective population sizes were largest in the Celtic Sea and Rockall, but low enough to be of potential conservation concern among Scottish and Faroese sites. Among the 21 candidate SNPs under positive selection was one significantly correlated with environmental variables predicted to be affected by climate change, including bottom temperature, salinity and pH. The paucity of well-annotated elasmobranch genomes precluded us from identifying a putative function for this SNP. Nevertheless, our findings suggest that climate change could inflict a strong selective force upon remnant populations of D. batis, further constraining its already-restricted habitat. Furthermore, the results provide fundamental insights on the distribution, behaviour and evolutionary biology of D. batis in the North-East Atlantic that will be useful for the establishment of conservation actions for this and other critically endangered elasmobranchs.

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.

Charting a course for genetic diversity in the UN Decade of Ocean Science

The health of the world's oceans is intrinsically linked to the biodiversity of the ecosystems they sustain. The importance of protecting and maintaining ocean biodiversity has been affirmed through the setting of the UN Sustainable Development Goal 14 to conserve and sustainably use the ocean for society's continuing needs. The decade beginning 2021-2030 has additionally been declared as the UN Decade of Ocean Science for Sustainable Development. This program aims to maximize the benefits of ocean science to the management, conservation, and sustainable development of the marine environment by facilitating communication and cooperation at the science-policy interface. A central principle of the program is the conservation of species and ecosystem components of biodiversity. However, a significant omission from the draft version of the Decade of Ocean Science Implementation Plan is the acknowledgment of the importance of monitoring and maintaining genetic biodiversity within species. In this paper, we emphasize the importance of genetic diversity to adaptive capacity, evolutionary potential, community function, and resilience within populations, as well as highlighting some of the major threats to genetic diversity in the marine environment from direct human impacts and the effects of global climate change. We then highlight the significance of ocean genetic diversity to a diverse range of socioeconomic factors in the marine environment, including marine industries, welfare and leisure pursuits, coastal communities, and wider society. Genetic biodiversity in the ocean, and its monitoring and maintenance, is then discussed with respect to its integral role in the successful realization of the 2030 vision for the Decade of Ocean Science. Finally, we suggest how ocean genetic diversity might be better integrated into biodiversity management practices through the continued interaction between environmental managers and scientists, as well as through key leverage points in industry requirements for Blue Capital financing and social responsibility.

The pantophysin gene and its relationship with survival in early life stages of Atlantic cod

Genetic markers are widely used in fisheries management around the world. While the genetic structure and markers selected are usually based on samples from the wild, very few controlled experiments have been carried out to investigate possible differences in influence on traits between markers. Here we examine the bi-allelic gene pantophysin (Pan I), widely used in the management of Atlantic cod, in a series of in vitro crosses under a range of temperatures. It has been proposed that this gene, or another tightly linked gene, may be under strong divergent selection. Resolving this issue is essential in order to interpret results when using this gene marker for stock management. We found no evidence of departure from the expected 1 : 2 : 1 Mendelian ratio for any of the three genotypes during the egg stage, while both the 6 and 12°C temperature regimes in tank experiments favoured the survival of the Pan IAA genotype. No difference in genotype survival was, however, found in a more natural mesocosm environment. Collectively, these results suggest that for the early life stages of Atlantic cod, and under the current experimental conditions, there is no strong consistent influence of Pan I genotype on survival. The results also emphasize the importance of varied experimental studies to verify the importance of environmental factors influencing genotype selection.