Archived DNA reveals fisheries and climate induced collapse of a major fishery

Unidirectional trans-Atlantic gene flow and a mixed spawning area shape the genetic connectivity of Atlantic bluefin tuna

The commercially important Atlantic bluefin tuna (Thunnus thynnus), a large migratory fish, has experienced notable recovery aided by accurate resource assessment and effective fisheries management efforts. Traditionally, this species has been perceived as consisting of eastern and western populations, spawning respectively in the Mediterranean Sea and the Gulf of Mexico, with mixing occurring throughout the Atlantic. However, recent studies have challenged this assumption by revealing weak genetic differentiation and identifying a previously unknown spawning ground in the Slope Sea used by Atlantic bluefin tuna of uncertain origin. To further understand the current and past population structure and connectivity of Atlantic bluefin tuna, we have assembled a unique dataset including thousands of genome-wide single-nucleotide polymorphisms (SNPs) from 500 larvae, young of the year and spawning adult samples covering the three spawning grounds and including individuals of other Thunnus species. Our analyses support two weakly differentiated but demographically connected ancestral populations that interbreed in the Slope Sea. Moreover, we also identified signatures of introgression from albacore (Thunnus alalunga) into the Atlantic bluefin tuna genome, exhibiting varied frequencies across spawning areas, indicating strong gene flow from the Mediterranean Sea towards the Slope Sea. We hypothesize that the observed genetic differentiation may be attributed to increased gene flow caused by a recent intensification of westward migration by the eastern population, which could have implications for the genetic diversity and conservation of western populations. Future conservation efforts should consider these findings to address potential genetic homogenization in the species.

Population Genetic Study on the European Flounder (Platichthys flesus) from the Southern Baltic Sea Using SNPs and Microsatellite Markers

The European flounder (Platichthys flesus), which is closely related to the recently discovered Baltic flounder (Platichthys solemdali), is currently the third most commercially fished species in the Baltic Sea. According to the available data from the Polish Fisheries Monitoring Center and fishermen's observations, the body condition indices of the species in the Baltic Sea have declined in recent years. The aim of the present study was to obtain information on the current patterns of genetic variability and the population structure of the European flounder and to verify whether the Baltic flounder is present in the southern Baltic Sea. Moreover, we aimed to verify whether the observed decline in the body condition indices of the species in the Baltic Sea might be associated with adaptive alterations in its gene pool due to increased fishing pressure. For this purpose, 190 fish were collected from four locations along the central coastline of Poland, i.e., Mechelinki, Władysławowo, the Vistula Lagoon in 2018, and the Słupsk Bank in 2020. The fish were morphologically analyzed and then genetically screened by the application of nineteen microsatellite DNA and two diagnostic SNP markers. The examined European flounder specimens displayed a high level of genetic diversity (PIC = 0.832-0.903, I = 2.579-2.768). A lack of significant genetic differentiation (Fst = 0.004, p > 0.05) was observed in all the examined fish, indicating that the European flounder in the sampled area constitutes a single genetic cluster. A significant deficiency in heterozygotes (Fis = 0.093, p < 0.05) and overall deviations from Hardy-Weinberg expectations (H-WE) were only detected in fish sampled from the Słupsk Bank. The estimated effective population size (Ne) among the sampled fish groups varied from 712 (Słupsk Bank) to 10,115 (Władysławowo and Mechelinki). However, the recorded values of the Garza-Williamson indicator (M = 0.574-0.600) and the lack of significant (p > 0.05) differences in Heq > He under the SMM model did not support the species' population size changes in the past. The applied SNP markers did not detect the presence of the Baltic flounder among the fish sampled from the studied area. The analysis of an association between biological traits and patterns of genetic diversity did not detect any signs of directional selection or density-dependent adaptive changes in the gene pool of the examined fish that might be caused by increased fishing pressure.

Hidden but revealed: After years of genetic studies behavioural monitoring combined with genomics uncover new insight into the population dynamics of Atlantic cod in Icelandic waters

Stock structure is of paramount importance for sustainable management of exploited resources. In that context, genetic markers have been used for more than two decades to resolve spatial structure of marine exploited resources and to fully fathom stock dynamics and interactions. While genetic markers such as allozymes and RFLP dominated the debate in the early era of genetics, technology advances have provided scientists with new tools every decade to better assess stock discrimination and interactions (i.e. gene flow). Here, we provide a review of genetic studies performed to understand stock structure of Atlantic cod in Icelandic waters, from the early allozyme approaches to the genomic work currently carried out. We further highlight the importance of the generation of a chromosome-anchored genome assembly together with whole-genome population data, which drastically changed our perception of the possible management units to consider. After nearly 60 years of genetic investigation of Atlantic cod structure in Icelandic waters, genetic (and later genomic) data combined with behavioural monitoring using Data Storage Tags shifted the attention from geographical population structures to behavioural ecotypes. This review also demonstrates the need for future research to further disentangle the impact of these ecotypes (and gene flow among them) on the population structure of Atlantic cod in Icelandic waters. It also highlights the importance of whole-genome data to unravel unexpected within-species diversity related to chromosomal inversions and associated supergenes, which are important to consider for future development of sustainable management programmes of the species within the North Atlantic.

Genotyping-in-Thousands by sequencing of archival fish scales reveals maintenance of genetic variation following a severe demographic contraction in kokanee salmon

Historical DNA analysis of archival samples has added new dimensions to population genetic studies, enabling spatiotemporal approaches for reconstructing population history and informing conservation management. Here we tested the efficacy of Genotyping-in-Thousands by sequencing (GT-seq) for collecting targeted single nucleotide polymorphism genotypic data from archival scale samples, and applied this approach to a study of kokanee salmon (Oncorhynchus nerka) in Kluane National Park and Reserve (KNPR; Yukon, Canada) that underwent a severe 12-year population decline followed by a rapid rebound. We genotyped archival scales sampled pre-crash and contemporary fin clips collected post-crash, revealing high coverage (> 90% average genotyping across all individuals) and low genotyping error (< 0.01% within-libraries, 0.60% among-libraries) despite the relatively poor quality of recovered DNA. We observed slight decreases in expected heterozygosity, allelic diversity, and effective population size post-crash, but none were significant, suggesting genetic diversity was retained despite the severe demographic contraction. Genotypic data also revealed the genetic distinctiveness of a now extirpated population just outside of KNPR, revealing biodiversity loss at the northern edge of the species distribution. More broadly, we demonstrated GT-seq as a valuable tool for collecting genome-wide data from archival samples to address basic questions in ecology and evolution, and inform applied research in wildlife conservation and fisheries management.

Ancient DNA SNP-panel data suggests stability in bluefin tuna genetic diversity despite centuries of fluctuating catches in the eastern Atlantic and Mediterranean

Atlantic bluefin tuna (Thunnus thynnus; BFT) abundance was depleted in the late 20th and early 21st century due to overfishing. Historical catch records further indicate that the abundance of BFT in the Mediterranean has been fluctuating since at least the 16th century. Here we build upon previous work on ancient DNA of BFT in the Mediterranean by comparing contemporary (2009–2012) specimens with archival (1911–1926) and archaeological (2nd century BCE–15th century CE) specimens that represent population states prior to these two major periods of exploitation, respectively. We successfully genotyped and analysed 259 contemporary and 123 historical (91 archival and 32 archaeological) specimens at 92 SNP loci that were selected for their ability to differentiate contemporary populations or their association with core biological functions. We found no evidence of genetic bottlenecks, inbreeding or population restructuring between temporal sample groups that might explain what has driven catch fluctuations since the 16th century. We also detected a putative adaptive response, involving the cytoskeletal protein synemin which may be related to muscle stress. However, these results require further investigation with more extensive genome-wide data to rule out demographic changes due to overfishing, and other natural and anthropogenic factors, in addition to elucidating the adaptive drivers related to these.

Intra-annual variation in feeding of Atlantic cod Gadus morhua: the importance of ephemeral prey bursts

Seasonal prey bursts are important for the life cycles and energy budgets of many predators. This study documents the diet and, especially, the importance of the ephemeral occurrence of capelin as prey for Atlantic cod (Gadus morhua) in Godthaabsfjord, west Greenland, over an annual cycle. The cod showed clear differences in diet composition on the 11 sampling dates resulting in a spring-summer, late summer-autumn and winter cluster. Moreover, a single sampling date, 12 May, was defined by cod gorge feeding on spawning capelin, which led to average stomach contents 4.3 times higher than the average for the remaining sampling dates. Changes in nitrogen stable isotope values from 22 April to 7 July in cod liver and muscle tissue were used to calculate the consumption of capelin. Based on this, the consumption of capelin varied between 538 and 658 g wet weight for a 1.3 kg cod. Using published consumption/biomass estimates and observed growth rates, the capelin intake corresponds to 10.1%-33.3% of the annual food consumption and accounts for 28.1%-34.5% of the annual growth of the cod. The present study documents the omnivorous feeding mode of Atlantic cod but highlights the utilization and importance of ephemeral prey bursts for the annual energy budget of the cod. It is hypothesized that access to capelin is critical for the postspawning recovery of Godthaabsfjord cod.

Genetic analyses reveal complex dynamics within a marine fish management area

Genetic data have great potential for improving fisheries management by identifying the fundamental management units-that is, the biological populations-and their mixing. However, so far, the number of practical cases of marine fisheries management using genetics has been limited. Here, we used Atlantic cod in the Baltic Sea to demonstrate the applicability of genetics to a complex management scenario involving mixing of two genetically divergent populations. Specifically, we addressed several assumptions used in the current assessment of the two populations. Through analysis of 483 single nucleotide polymorphisms (SNPs) distributed across the Atlantic cod genome, we confirmed that a model of mechanical mixing, rather than hybridization and introgression, best explained the pattern of genetic differentiation. Thus, the fishery is best monitored as a mixed-stock fishery. Next, we developed a targeted panel of 39 SNPs with high statistical power for identifying population of origin and analyzed more than 2,000 tissue samples collected between 2011 and 2015 as well as 260 otoliths collected in 2003/2004. These data provided high spatial resolution and allowed us to investigate geographical trends in mixing, to compare patterns for different life stages and to investigate temporal trends in mixing. We found similar geographical trends for the two time points represented by tissue and otolith samples and that a recently implemented geographical management separation of the two populations provided a relatively close match to their distributions. In contrast to the current assumption, we found that patterns of mixing differed between juveniles and adults, a signal likely linked to the different reproductive dynamics of the two populations. Collectively, our data confirm that genetics is an operational tool for complex fisheries management applications. We recommend focussing on developing population assessment models and fisheries management frameworks to capitalize fully on the additional information offered by genetically assisted fisheries monitoring.

Cryptic temporal changes in stock composition explain the decline of a flounder (Platichthys spp.) assemblage

Unobserved diversity, such as undetected genetic structure or the presence of cryptic species, is of concern for the conservation and management of global biodiversity in the face of threatening anthropogenic processes. For instance, unobserved diversity can lead to overestimation of maximum sustainable yields and therefore to overharvesting of the more vulnerable stock components within unrecognized mixed-stock fisheries. We used DNA from archival (otolith) samples to reconstruct the temporal (1976-2011) genetic makeup of two mixed-stock flounder fisheries in the Åland Sea (AS) and the Gulf of Finland (GoF). Both fisheries have hitherto been managed as a single stock of European flounders (Platichthys flesus), but were recently revealed to target two closely related species: the pelagic-spawning P. flesus and the newly described, demersal-spawning P. solemdali. While the AS and GoF fisheries were assumed to consist exclusively of P. solemdali, P. flesus dominated the GoF flounder assemblage (87% of total) in 1983, had disappeared (0%) by 1993, and remained in low proportions (10%-11%) thereafter. In the AS, P. solemdali dominated throughout the sampling period (>70%), and P. flesus remained in very low proportions after 1983. The disappearance of P. flesus from the GoF coincides in time with a dramatic (~60%) decline in commercial landings and worsening environmental conditions in P. flesus' northernmost spawning ground, the Eastern Gotland Basin, in the preceding 4-6 years. These results are compatible with the hypothesis that P. flesus in the GoF is a sink population relying on larval subsidies from southern spawning grounds and the cause of their disappearance is a cessation of larval supply. Our results highlight the importance of uncovering unobserved genetic diversity and studying spatiotemporal changes in the relative contribution of different stock components, as well as the underlying environmental causes, to manage marine resources in the age of rapid anthropogenic change.

Stable coexistence of genetically divergent Atlantic cod ecotypes at multiple spatial scales

Coexistence in the same habitat of closely related yet genetically different populations is a phenomenon that challenges our understanding of local population structure and adaptation. Identifying the underlying mechanisms for such coexistence can yield new insight into adaptive evolution, diversification and the potential for organisms to adapt and persist in response to a changing environment. Recent studies have documented cryptic, sympatric populations of Atlantic cod (Gadus morhua) in coastal areas. We analysed genetic origin of 6,483 individual cod sampled annually over 14 years from 125 locations along the Norwegian Skagerrak coast and document stable coexistence of two genetically divergent Atlantic cod ecotypes throughout the study area and study period. A "fjord" ecotype dominated in numbers deep inside fjords while a "North Sea" ecotype was the only type found in offshore North Sea. Both ecotypes coexisted in similar proportions throughout coastal habitats at all spatial scales. The size-at-age of the North Sea ecotype on average exceeded that of the fjord ecotype by 20% in length and 80% in weight across all habitats. Different growth and size among individuals of the two types might be one of several ecologically significant variables that allow for stable coexistence of closely related populations within the same habitat. Management plans, biodiversity initiatives and other mitigation strategies that do not account for the mixture of species ecotypes are unlikely to meet objectives related to the sustainability of fish and fisheries.

Don't bet against the natal homing abilities of marine fishes

Whether marine fishes are capable of homing to their natal areas has long been something of an enigma. For some estuarine species or sharks (which have extended nondispersal juvenile stages or are born as relatively large, fully formed juveniles), the answer is clearly 'yes' (Thorrold et al. ; Feldheim et al. ), but for most marine fishes, the issue is much more mysterious. Many species have free-floating eggs, and most have pelagic, passively dispersing larvae. It is challenging to imagine how adult fish might navigate to a region of the ocean they experienced only as eggs or larvae, and easier to assume that such dispersal leads inexorably to high gene flow, and even panmixia. One way to resolve the conundrum would be to track fish from hatching to reproduction, but for marine fishes with tiny eggs and drifting larvae, this is notoriously difficult to do (Bradbury & Laurel ). In this issue of Molecular Ecology, Bonanomi et al. () use a creative approach to solve this challenge for Atlantic cod (Gadus morhua) populations that mingle in the vicinity of Greenland. They show that cod that disperse more than a 1000 km away from Iceland as eggs and larvae, then spend years growing on the far side of Greenland, while mixing with two local populations, return as adults to spawning areas near Iceland - and further, that this behaviour has remained stable over more than six decades. They manage this feat with a clever use of historical cod tracking data, modern genomic data and genetic analysis of decades-old DNA obtained from archived materials. Their results have important implications for our view of the biocomplexity of marine fish populations, and how we should manage them.

Stocking activities for the Arctic charr in Lake Geneva: Genetic effects in space and time

Artificial stocking practices are widely used by resource managers worldwide, in order to sustain fish populations exploited by both recreational and commercial activities, but their benefits are controversial. Former practices involved exotic strains, although current programs rather consider artificial breeding of local fishes (supportive breeding). Understanding the complex genetic effects of these management strategies is an important challenge with economic and conservation implications, especially in the context of population declines. In this study, we focus on the declining Arctic charr (Salvelinus alpinus) population from Lake Geneva (Switzerland and France), which has initially been restocked with allochtonous fishes in the early eighties, followed by supportive breeding. In this context, we conducted a genetic survey to document the evolution of the genetic diversity and structure throughout the last 50 years, before and after the initiation of hatchery supplementation, using contemporary and historical samples. We show that the introduction of exotic fishes was associated with a genetic bottleneck in the 1980-1990s, a break of Hardy-Weinberg Equilibrium (HWE), a reduction in genetic diversity, an increase in genetic structure among spawning sites, and a change in their genetic composition. Together with better environmental conditions, three decades of subsequent supportive breeding using local fishes allowed to re-establish HWE and the initial levels of genetic variation. However, current spawning sites have not fully recovered their original genetic composition and were extensively homogenized across the lake. Our study demonstrates the drastic genetic consequences of different restocking tactics in a comprehensive spatiotemporal framework and suggests that genetic alteration by nonlocal stocking may be partly reversible through supportive breeding. We recommend that conservation-based programs consider local diversity and implement adequate protocols to limit the genetic homogenization of this Arctic charr population.

Genetic structure of West Greenland populations of lumpfish Cyclopterus lumpus

In this study, 11 microsatellite markers were used to determine the structure of West Greenlandic lumpfish Cyclopterus lumpus populations across six spawning locations spanning >1500 km and compared with neighbouring populations in Canada and Iceland. To evaluate whether data allow for identification of origin of C. lumpus in Greenlandic waters, genetic assignment analysis was performed for 86 C. lumpus sampled on a feeding migration. Significant structuring with isolation by distance was observed in the West Greenland samples and two major subpopulations, north and south, were suggested. Based on F_ST values, closer relationships were observed between Greenland and Canada, than Greenland and Iceland. Surprisingly, the North Greenland population showed more similarities with Canadian samples, than did the geographically closer south-west Greenland population. Origin could be assigned for a high proportion of non-spawning fish and demonstrated a marked east-west spatial separation of fish of Greenlandic and Icelandic genotypes.

Navigating the currents of seascape genomics: how spatial analyses can augment population genomic studies

Population genomic approaches are making rapid inroads in the study of non-model organisms, including marine taxa. To date, these marine studies have predominantly focused on rudimentary metrics describing the spatial and environmental context of their study region (e.g., geographical distance, average sea surface temperature, average salinity). We contend that a more nuanced and considered approach to quantifying seascape dynamics and patterns can strengthen population genomic investigations and help identify spatial, temporal, and environmental factors associated with differing selective regimes or demographic histories. Nevertheless, approaches for quantifying marine landscapes are complicated. Characteristic features of the marine environment, including pelagic living in flowing water (experienced by most marine taxa at some point in their life cycle), require a well-designed spatial-temporal sampling strategy and analysis. Many genetic summary statistics used to describe populations may be inappropriate for marine species with large population sizes, large species ranges, stochastic recruitment, and asymmetrical gene flow. Finally, statistical approaches for testing associations between seascapes and population genomic patterns are still maturing with no single approach able to capture all relevant considerations. None of these issues are completely unique to marine systems and therefore similar issues and solutions will be shared for many organisms regardless of habitat. Here, we outline goals and spatial approaches for landscape genomics with an emphasis on marine systems and review the growing empirical literature on seascape genomics. We review established tools and approaches and highlight promising new strategies to overcome select issues including a strategy to spatially optimize sampling. Despite the many challenges, we argue that marine systems may be especially well suited for identifying candidate genomic regions under environmentally mediated selection and that seascape genomic approaches are especially useful for identifying robust locus-by-environment associations.

Coupling Demographic and Genetic Variability from Archived Collections of European Anchovy (Engraulis encrasicolus)

It is well known that temporal fluctuations in small populations deeply influence evolutionary potential. Less well known is whether fluctuations can influence the evolutionary potentials of species with large census sizes. Here, we estimated genetic population parameters from as survey of polymorphic microsatellite DNA loci in archived otoliths from Adriatic European anchovy (Engraulis encrasicolus), a fish with large census sizes that supports numerous local fisheries. Stocks have fluctuated greatly over the past few decades, and the Adriatic fishery collapsed in 1987. Our results show a significant reduction of mean genetic parameters as a consequence of the population collapse. In addition, estimates of effective population size (Ne) are much smaller than those expected in a fishes with large population census sizes (Nc). Estimates of Ne indicate low effective population sizes, even before the population collapse. The ratio Ne/Ne ranged between 10-6 and 10-8, indicating a large discrepancy between the anchovy gene pool and population census size. Therefore, anchovy populations may be more vulnerable to fishery effort and environmental change than previously thought.

Historical DNA documents long-distance natal homing in marine fish

The occurrence of natal homing in marine fish remains a fundamental question in fish ecology as its unequivocal demonstration requires tracking of individuals from fertilization to reproduction. Here, we provide evidence of long-distance natal homing (>1000 km) over more than 60 years in Atlantic cod (Gadus morhua), through genetic analysis of archived samples from marked and recaptured individuals. Using a high differentiation single-nucleotide polymorphism assay, we demonstrate that the vast majority of cod tagged in West Greenland and recaptured on Icelandic spawning grounds belonged to the Iceland offshore population, strongly supporting a hypothesis of homing. The high degree of natal fidelity observed provides the evolutionary settings for development of locally adapted populations in marine fish and emphasize the need to consider portfolio effects in marine fisheries management strategies.