Search

Ostrea edulis, the European flat oyster, was once a widespread economically and ecologically important marine species, but has suffered dramatic declines over the past two centuries. Consequently, there has been a surge in European restoration efforts, many of which focus on restocking as a conservation measure. In this study, we used whole-genome sequencing (WGS) data to investigate the population structure, demographic history, and patterns of local adaptation of O. edulis across its natural distribution with increased sampling densities at Scandinavian localities. Results revealed seven distinct genetic clusters, including previously undescribed complex population structure in Norway, and evidence for introgression between genetic clusters in Scandinavia. We detected large structural variants (SVs) on three pseudo-chromosomes. These megabase long regions were characterised by strong linkage disequilibrium and clear geographical differentiation, suggestive of chromosomal inversions potentially associated with local adaptation. The results indicated that genomic traces of past translocations of non-native O. edulis were still present in some individuals, but overall, we found limited evidence of major impacts of translocations on the scale of contemporary population structure. Our findings highlight the importance of considering population structure and signatures of selection in the design of effective conservation strategies to preserve and restore wild native European flat oyster populations, and we provide direct knowledge safeguarding sustainable mitigation actions in this important species., (© 2024 The Author(s). Molecular Ecology published by John Wiley & Sons Ltd.)

Sustainable fisheries management is important for the continued harvest of the world's marine resources, especially as they are increasingly challenged by a range of climatic and anthropogenic factors. One of the pillars of sustainable fisheries management is the accurate identification of the biological units, i.e., populations. Here, we developed and implemented a genetic baseline for Atlantic herring harvested in the Norwegian offshore fisheries to investigate the validity of the current management boundaries. This was achieved by genotyping > 15,000 herring from the northern European seas, including samples of all the known populations in the region, with a panel of population-informative SNPs mined from existing genomic resources. The final genetic baseline consisted of ~1000 herring from 12 genetically distinct populations. We thereafter used the baseline to investigate mixed catches from the North and Norwegian Seas, revealing that each management area consisted of multiple populations, as previously suspected. However, substantial numbers (up to 50% or more within a sample) of herring were found outside of their expected management areas, e.g., North Sea autumn-spawning herring north of 62° N (average = 19.2%), Norwegian spring-spawning herring south of 62° N (average = 13.5%), and western Baltic spring-spawning herring outside their assumed distribution area in the North Sea (average = 20.0%). Based upon these extensive observations, we conclude that the assessment and management areas currently in place for herring in this region need adjustments to reflect the populations present. Furthermore, we suggest that for migratory species, such as herring, a paradigm shift from using static geographic stock boundaries towards spatial dynamic boundaries is needed to meet the requirements of future sustainable management regimes., Competing Interests: The authors declare no conflicts of interest., (© 2024 The Author(s). Evolutionary Applications published by John Wiley & Sons Ltd.)

Populations of anadromous brown trout, also known as sea trout, have suffered recent marked declines in abundance due to multiple factors, including climate change and human activities. While much is known about their freshwater phase, less is known about the species' marine feeding migrations. This situation is hindering the effective management and conservation of anadromous trout in the marine environment. Using a panel of 95 single nucleotide polymorphism markers we developed a genetic baseline, which demonstrated strong regional structuring of genetic diversity in trout populations around the English Channel and adjacent waters. Extensive baseline testing showed this structuring allowed high-confidence assignment of known-origin individuals to region of origin. This study presents new data on the movements of anadromous trout in the English Channel and southern North Sea. Assignment of anadromous trout sampled from 12 marine and estuarine localities highlighted contrasting results for these areas. The majority of these fisheries are composed predominately of stocks local to the sampling location. However, there were multiple cases of long-distance movements of anadromous trout, with several individuals originating from rivers in northeast England being caught in the English Channel and southern North Sea, in some cases more than 1000 km from their natal region. These results have implications for the management of sea trout in inshore waters around the English Channel and southern North Sea., Competing Interests: The authors declare no conflicts of interest., (© 2024 The Author(s). Evolutionary Applications published by John Wiley & Sons Ltd.)

The European sprat is a small plankton-feeding clupeid present in the northeastern Atlantic Ocean, in the Mediterranean Sea, and in the brackish Baltic Sea and Black Sea. This species is the target of a major fishery and, therefore, an accurate characterization of its genetic population structure is crucial to delineate proper stock assessments that aid ensuring the fishery's sustainability. Here, we present (i) a draft genome assembly, (ii) pooled whole genome sequencing of 19 population samples covering most of the species' distribution range, and (iii) the design and test of a single nucleotide polymorphism (SNP)-chip resource and use this to validate the population structure inferred from pooled sequencing. These approaches revealed, using the populations sampled here, three major groups of European sprat: Oceanic, Coastal, and Brackish with limited differentiation within groups even over wide geographical stretches. Genetic structure is largely driven by six large putative inversions that differentiate Oceanic and Brackish sprats, while Coastal populations display intermediate frequencies of haplotypes at each locus. Interestingly, populations from the Baltic and the Black Seas share similar frequencies of haplotypes at these putative inversions despite their distant geographic location. The closely related clupeids European sprat and Atlantic herring both show genetic adaptation to the brackish Baltic Sea, providing an opportunity to explore the extent of genetic parallelism. This analysis revealed limited parallelism because out of 125 independent loci detected in the Atlantic herring, three showed sharp signals of selection that overlapped between the two species and contained single genes such as PRLRA, which encodes the receptor for prolactin, a freshwater-adapting hormone in euryhaline species, and THRB, a receptor for thyroid hormones, important both for metabolic regulation and the development of red cone photoreceptors., (© The Author(s) 2024. Published by Oxford University Press on behalf of Society for Molecular Biology and Evolution.)

Introgression of non-native conspecifics changes the genetic composition of wild populations, potentially leading to loss of local adaptations and fitness declines. However, long-term data from wild populations are still relatively few. Here, we studied the effects of introgression in a Danish brown trout ( Salmo trutta , L.) population, subjected to intensive stocking with domesticated hatchery fish of non-native origin. We used wild-caught genetically wild and admixed trout as well as fish from the partly domesticated hatchery strain used for stocking the river up until ~15 years prior to this study, to produce 22 families varying in hatchery/wild admixture. Following a replicated common-garden experiment conducted in fish tanks from first feeding through 23 weeks at 7, 12, and 16°C, we observed a significant positive relationship between family admixture and fish size upon termination, an effect observed through all levels of admixture. Furthermore, the admixture effect was most distinct at the higher rearing temperatures. Although the hatchery strain used for stocking had been in culture for ~7 generations, it had not been deliberately selected for increased growth. These data thus demonstrate: (i) that growth had increased in the hatchery strain even in the absence of deliberate directional selection for this trait, (ii) that the increasing effect of admixture by temperature could represent inadvertent selection for performance in the hatchery strain at higher temperatures, and most significantly, (iii) that despite undergoing up to five generations of natural selection in the admixed wild population, the genetically increased growth potential was still detectable and thus persistent. Our findings suggest that altered growth patterns and potentially their cascading effects are of importance to the severity of hatchery/wild introgression, especially under changing-climate scenarios and are of general significance to conservation practitioners seeking to evaluate long-term effects of intra-specific hybridization including under recovery., Competing Interests: Authors declare no conflicts of interest., (© 2024 The Author(s). Evolutionary Applications published by John Wiley & Sons Ltd.)

Population genomics analysis holds great potential for informing conservation of endangered populations. We focused on a controversial case of European whitefish (Coregonus spp.) populations. The endangered North Sea houting is the only coregonid fish that tolerates oceanic salinities and was previously considered a species (C. oxyrhinchus) distinct from European lake whitefish (C. lavaretus). However, no firm evidence for genetic-based salinity adaptation has been available. Also, studies based on microsatellite and mitogenome data suggested surprisingly recent divergence (c. 2500 years bp) between houting and lake whitefish. These data types furthermore have provided no evidence for possible inbreeding. Finally, a controversial taxonomic revision recently classified all whitefish in the region as C. maraena, calling conservation priorities of houting into question. We used whole-genome and ddRAD sequencing to analyse six lake whitefish populations and the only extant indigenous houting population. Demographic inference indicated post-glacial expansion and divergence between lake whitefish and houting occurring not long after the Last Glaciation, implying deeper population histories than previous analyses. Runs of homozygosity analysis suggested not only high inbreeding (F ROH up to 30.6%) in some freshwater populations but also F ROH up to 10.6% in the houting prompting conservation concerns. Finally, outlier scans provided evidence for adaptation to high salinities in the houting. Applying a framework for defining conservation units based on current and historical reproductive isolation and adaptive divergence led us to recommend that the houting be treated as a separate conservation unit regardless of species status. In total, the results underscore the potential of genomics to inform conservation practices, in this case clarifying conservation units and highlighting populations of concern., (© 2024 The Authors. Molecular Ecology published by John Wiley & Sons Ltd.)

Background: The Atlantic mackerel, Scomber scombrus (Linnaeus, 1758) is a commercially valuable migratory pelagic fish inhabiting the northern Atlantic Ocean and the Mediterranean Sea. Given its highly migratory behaviour for feeding and spawning, several studies have been conducted to assess differentiation among spawning components to better define management units, as well as to investigate possible adaptations to comprehend and predict recent range expansion northwards., Methods: Here, a high-quality genome of S. scombrus was sequenced and annotated, as an increasing number of population genetic studies have proven the relevance of reference genomes to investigate genomic markers/regions potentially linked to differences at finer scale. Such reference genome was used to map Restriction-site-associated sequencing (RAD-seq) reads for SNP discovery and genotyping in more than 500 samples distributed along the species range. The resulting genotyping tables have been used to perform connectivity and adaptation analyses., Results: The assembly of the reference genome for S. scombrus resulted in a high-quality genome of 741 Mb. Our population genetic results show that the Atlantic mackerel consist of three previously known genetically isolated units (Northwest Atlantic, Northeast Atlantic, Mediterranean), and provide no evidence for genetically distinct spawning components within the Northwest or Northeast Atlantic., Conclusions: Therefore, our findings resolved previous uncertainties by confirming the absence of genetically isolated spawning components in each side of the northern Atlantic, thus rejecting homing behaviour and the need to redefine management boundaries in this species. In addition, no further genetic signs of ongoing adaptation were detected in this species., Competing Interests: No competing interests were disclosed., (Copyright: © 2024 Manuzzi A et al.)

Maintenance of genetic diversity in marine fishes targeted by commercial fishing is a grand challenge for the future. Most of these species are abundant and therefore important for marine ecosystems and food security. Here, we present a road map of how population genomics can promote sustainable fisheries. In these species, the development of reference genomes and whole genome sequencing is key, because genetic differentiation at neutral loci is usually low due to large population sizes and gene flow. First, baseline allele frequencies representing genetically differentiated populations within species must be established. These can then be used to accurately determine the composition of mixed samples, forming the basis for population demographic analysis to inform sustainably set fish quotas. SNP-chip analysis is a cost-effective method for determining baseline allele frequencies and for population identification in mixed samples. Finally, we describe how genetic marker analysis can transform stock identification and management.

Pelagic fish like herring, sardines, and mackerel constitute an essential and nutritious human food source globally. Their sustainable harvest is promoted by the application of precise, accurate, and cost-effective methods for estimating bycatch. Here, we experimentally test the new concept of using eDNA for quantitative bycatch assessment on the illustrative example of the Baltic Sea sprat fisheries with herring bycatch. We investigate the full pipeline from sampling of production water on vessels and in processing factories to the estimation of species weight fractions. Using a series of controlled mixture experiments, we demonstrate that the eDNA signal from production water shows a strong, seasonally consistent linear relationship with herring weight fractions, however, the relationship is influenced by the molecular method used (qPCR or metabarcoding). In four large sprat landings analyzed, despite examples of remarkable consistency between eDNA and visual reporting, estimates of herring bycatch biomass varied between the methods applied, with the eDNA-based estimates having the highest precision for all landings analyzed. The eDNA-based bycatch assessment method has the potential to improve the quality and cost effectiveness of bycatch assessment in large pelagic fisheries catches and in the long run lead to more sustainable management of pelagic fish as a precious marine resource., (© 2024. The Author(s).)

Lumpfish (Cyclopterus lumpus) is a transatlantic marine fish displaying large population sizes and a high potential for dispersal and gene-flow. These features are expected to result in weak population structure. Here, we investigated population genetic structure of lumpfish throughout its natural distribution in the North Atlantic using two approaches: I) 4,393 genome wide SNPs and 95 individuals from 10 locations, and II) 139 discriminatory SNPs and 1,669 individuals from 40 locations. Both approaches identified extensive population genetic structuring with a major split between the East and West Atlantic and a distinct Baltic Sea population, as well as further differentiation of lumpfish from the English Channel, Iceland, and Greenland. The discriminatory loci displayed ~2-5 times higher divergence than the genome wide approach, revealing further evidence of local population substructures. Lumpfish from Isfjorden in Svalbard were highly distinct but resembled most fish from Greenland. The Kattegat area in the Baltic transition zone, formed a previously undescribed distinct genetic group. Also, further subdivision was detected within North America, Iceland, West Greenland, Barents Sea, and Norway. Although lumpfish have considerable potential for dispersal and gene-flow, the observed high levels of population structuring throughout the Atlantic suggests that this species may have a natal homing behavior and local populations with adaptive differences. This fine-scale population structure calls for consideration when defining management units for exploitation of lumpfish stocks and in decisions related to sourcing and moving lumpfish for cleaner fish use in salmonid aquaculture., Competing Interests: The authors have declared that no competing interests exist., (Copyright: © 2023 Jansson et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.)

Atlantic herring in International Council for Exploration of the Sea (ICES) Divisions 6.a, 7.b-c comprises at least three populations, distinguished by temporal and spatial differences in spawning, which have until recently been managed as two stocks defined by geographical delineators. Outside of spawning the populations form mixed aggregations, which are the subject of acoustic surveys. The inability to distinguish the populations has prevented the development of separate survey indices and separate stock assessments. A panel of 45 single-nucleotide polymorphisms, derived from whole-genome sequencing, were used to genotype 3480 baseline spawning samples (2014-2021). A temporally stable baseline comprising 2316 herring from populations known to inhabit Division 6.a was used to develop a genetic assignment method, with a self-assignment accuracy greater than 90%. The long-term temporal stability of the assignment model was validated by assigning archive (2003-2004) baseline samples (270 individuals) with a high level of accuracy. Assignment of non-baseline samples (1514 individuals) from Divisions 6.a, 7.b-c indicated previously unrecognized levels of mixing of populations outside of the spawning season. The genetic markers and assignment models presented constitute a 'toolbox' that can be used for the assignment of herring caught in mixed survey and commercial catches in Division 6.a into their population of origin with a high level of accuracy., (© 2022 The Authors.)

Monitoring the distribution of marine nonindigenous species is a challenging task. To support this monitoring, we developed and validated the specificity of 12 primer-probe assays for detection of environmental DNA (eDNA) from marine species, all nonindigenous to Europe. The species include sturgeons, a Pacific red algae, oyster thief, a freshwater hydroid from the Black Sea, Chinese mitten crab, Pacific oyster, warty comb jelly, sand gaper, round goby, pink salmon, rainbow trout and North American mud crab. We tested all assays in the laboratory, on DNA extracted from both the target and non-target species to ensure that they only amplified DNA from the intended species. Subsequently, all assays were used to analyse water samples collected at 16 different harbours across two different seasons during 2017. We also included six previously published assays targeting eDNA from goldfish, European carp, two species of dinoflagellates of the genera Karenia and Prorocentrum, two species of the heterokont flagellate genus Pseudochattonella. Conventional monitoring was carried out alongside eDNA sampling but with only one sampling event over the one year. Because eDNA was relatively fast and easy to collect compared to conventional sampling, we sampled eDNA twice during 2017, which showed seasonal changes in the distribution of nonindigenous species. Comparing eDNA levels with salinity gradients did not show any correlation. A significant correlation was observed between number of species detected with conventional monitoring methods and number of species found using eDNA at each location. This supports the use of eDNA for surveillance of the distribution of marine nonindigenous species, where the speed and relative easy sampling in the field combined with fast molecular analysis may provide advantages compared to conventional monitoring methods. Prior validation of assays increases taxonomic precision, and laboratorial setup facilitates analysis of multiple samples simultaneously. The specific eDNA assays presented here can be implemented directly in monitoring programmes across Europe and potentially worldwide to infer a more precise picture of the dynamics in the distribution of marine nonindigenous species., Competing Interests: Declaration of competing interest The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (Copyright © 2022 The Authors. Published by Elsevier B.V. All rights reserved.)

Understanding how eco-evolutionary processes and environmental factors drive population differentiation and adaptation are key challenges in evolutionary biology of relevance for biodiversity protection. Differentiation requires at least partial reproductive separation, which may result from different modes of isolation such as geographic isolation (allopatry) or isolation by distance (IBD), resistance (IBR), and environment (IBE). Despite that multiple modes might jointly influence differentiation, studies that compare the relative contributions are scarce. Using RADseq, we analyse neutral and adaptive genetic diversity and structure in 11 pike (Esox lucius) populations from contrasting environments along a latitudinal gradient (54.9-63.6°N), to investigate the relative effects of IBD, IBE and IBR, and to assess whether the effects differ between neutral and adaptive variation, or across structural levels. Patterns of neutral and adaptive variation differed, probably reflecting that they have been differently affected by stochastic and deterministic processes. The importance of the different modes of isolation differed between neutral and adaptive diversity, yet were consistent across structural levels. Neutral variation was influenced by interactions among all three modes of isolation, with IBR (seascape features) playing a central role, wheares adaptive variation was mainly influenced by IBE (environmental conditions). Taken together, this and previous studies suggest that it is common that multiple modes of isolation interactively shape patterns of genetic variation, and that their relative contributions differ among systems. To enable identification of general patterns and understand how various factors influence the relative contributions, it is important that several modes are simultaneously investigated in additional populations, species and environmental settings., (© 2021 The Authors. Molecular Ecology published by John Wiley & Sons Ltd.)

Changing environmental conditions can lead to population diversification through differential selection on standing genetic variation. Structural variant (SV) polymorphisms provide examples of ancient alleles that in time become associated with novel environmental gradients. The European plaice (Pleuronectes platessa) is a marine flatfish showing large allele-frequency differences at two putative SVs associated with environmental variation. In this study, we explored the contribution of these SVs to population structure across the North East Atlantic. We compared genome-wide population structure using sets of RAD-sequencing SNPs with the spatial structure of the SVs. We found that in contrast to the rest of the genome, the SVs were only weakly associated with an isolation-by-distance pattern. Indeed, both SVs showed important variation in haplogroup frequencies, with the same haplogroup increasing both along the salinity gradient of the Baltic Sea, and found in high frequency in the northern-range margin of the Atlantic. Phylogenetic analyses suggested that the SV alleles are much older than the age of the Baltic Sea itself. These results suggest that the SVs are older than the age of the environmental gradients with which they currently co-vary. Altogether, our results suggest that the plaice SVs were shaped by evolutionary processes occurring at two time frames, firstly following their origin, ancient spread and maintenance in the ancestral populations, and secondly related to their current association with more recently formed environmental gradients such as those found in the North Sea-Baltic Sea transition zone.

Harvest of fish and wildlife, both commercial and recreational, is a selective force that can induce evolutionary changes to life history and behavior. Naturally selective forces may create countering selection pressures. Assessing natural fitness represents a considerable challenge in broadcast spawners. Thus, our understanding about the relative strength of natural and fisheries selection is slim. In the field, we compared the strength and shape of harvest selection to natural selection on body size over four years and behavior over one year in a natural population of a freshwater top predator, the northern pike ( Esox lucius ). Natural selection was approximated by relative reproductive success via parent-offspring genetic assignments over four years. Harvest selection was measured by comparing individuals susceptible to recreational angling with individuals never captured by this gear type. Individual behavior was measured by high-resolution acoustic telemetry. Harvest and natural size selection operated with equal strength but opposing directions, and harvest size selection was consistently negative in all study years. Harvest selection also had a substantial behavioral component independent of body length, while natural behavioral selection was not documented, suggesting the potential for directional harvest selection favoring inactive, timid fish. Simulations of the outcomes of different fishing regulations showed that traditional minimum size-based harvest limits are unlikely to counteract harvest selection without being completely restrictive. Our study suggests harvest selection may be inevitable and recreational fisheries may thus favor small, inactive, shy, and difficult-to-capture fish. Increasing fractions of shy fish in angling-exploited stocks would have consequences for stock assessment and all fisheries operating with hook and line., Competing Interests: The authors declare no competing interest.

Habitat changes represent one of the five most pervasive threats to biodiversity. However, anthropogenic activities also have the capacity to create novel niche spaces to which species respond differently. In 1880, one such habitat alterations occurred in Landvikvannet, a freshwater lake on the Norwegian coast of Skagerrak, which became brackish after being artificially connected to the sea. This lake is now home to the European sprat, a pelagic marine fish that managed to develop a self-recruiting population in barely few decades. Landvikvannet sprat proved to be genetically isolated from the three main populations described for this species; that is, Norwegian fjords, Baltic Sea, and the combination of North Sea, Kattegat, and Skagerrak. This distinctness was depicted by an accuracy self-assignment of 89% and a highly significant F ST between the lake sprat and each of the remaining samples (average of ≈0.105). The correlation between genetic and environmental variation indicated that salinity could be an important environmental driver of selection (3.3% of the 91 SNPs showed strong associations). Likewise, Isolation by Environment was detected for salinity, although not for temperature, in samples not adhering to an Isolation by Distance pattern. Neighbor-joining tree analysis suggested that the source of the lake sprat is in the Norwegian fjords, rather than in the Baltic Sea despite a similar salinity profile. Strongly drifted allele frequencies and lower genetic diversity in Landvikvannet compared with the Norwegian fjords concur with a founder effect potentially associated with local adaptation to low salinity. Genetic differentiation ( F ST ) between marine and brackish sprat is larger in the comparison Norway-Landvikvannet than in Norway-Baltic, which suggests that the observed divergence was achieved in Landvikvannet in some 65 generations, that is, 132 years, rather than gradually over thousands of years (the age of the Baltic Sea), thus highlighting the pace at which human-driven evolution can happen., Competing Interests: None declared., (© 2021 The Authors. Ecology and Evolution published by John Wiley & Sons Ltd.)

Atlantic herring is widespread in North Atlantic and adjacent waters and is one of the most abundant vertebrates on earth. This species is well suited to explore genetic adaptation due to minute genetic differentiation at selectively neutral loci. Here, we report hundreds of loci underlying ecological adaptation to different geographic areas and spawning conditions. Four of these represent megabase inversions confirmed by long read sequencing. The genetic architecture underlying ecological adaptation in herring deviates from expectation under a classical infinitesimal model for complex traits because of large shifts in allele frequencies at hundreds of loci under selection., Competing Interests: FH, MJ, MP, LS, AF, BD, DB, FB, MC, GD, EF, AF, LA No competing interests declared, (© 2020, Han et al.)

Environmental DNA (eDNA) is increasingly used for monitoring marine organisms; however, offshore sampling and time lag from sampling to results remain problematic. In order to overcome these challenges a robotic sampler, a 2nd generation Environmental Sample Processor (ESP), was tested for autonomous analysis of eDNA from four commercial fish species in a 4.5 million liter mesocosm. The ESP enabled in situ analysis, consisting of water collection, filtration, DNA extraction and qPCR analysis, which allowed for real-time remote reporting and archival sample collection, consisting of water collection, filtration and chemical preservation followed by post-deployment laboratory analysis. The results demonstrate that the 2G ESP was able to consistently detect and quantify target molecules from the most abundant species (Atlantic mackerel) both in real-time and from the archived samples. In contrast, detection of low abundant species was challenged by both biological and technical aspects coupled to the ecology of eDNA and the 2G ESP instrumentation. Comparison of the in situ analysis and archival samples demonstrated variance, which potentially was linked to diel migration patterns of the Atlantic mackerel. The study demonstrates strong potential for remote autonomous in situ monitoring which open new possibilities for the field of eDNA and marine monitoring.

Massive genotyping of single nucleotide polymorphisms (SNP) has opened opportunities for analyzing the way in which selection shapes genomes. Artificial or natural selection usually leaves genomic signatures associated with selective sweeps around the responsible locus. Strong selective sweeps are most often identified either by lower genetic diversity than the genomic average and/or islands of runs of homozygosity (ROHi). Here, we conducted an analysis of selective sweeps in turbot ( Scophthalmus maximus ) using two SNP datasets from a Northeastern Atlantic population (36 individuals) and a domestic broodstock (46 individuals). Twenty-six families (∼ 40 offspring per family) from this broodstock and three SNP datasets applying differing filtering criteria were used to adjust ROH calling parameters. The best-fitted genomic inbreeding estimate (F ROH ) was obtained by the sum of ROH longer than 1 Mb, called using a 21,615 SNP panel, a sliding window of 37 SNPs and one heterozygous SNP per window allowed. These parameters were used to obtain the ROHi distribution in the domestic and wild populations (49 and 0 ROHi, respectively). Regions with higher and lower genetic diversity within each population were obtained using sliding windows of 37 SNPs. Furthermore, those regions were mapped in the turbot genome against previously reported genetic markers associated with QTL (Quantitative Trait Loci) and outlier loci for domestic or natural selection to identify putative selective sweeps. Out of the 319 and 278 windows surpassing the suggestive pooled heterozygosity thresholds (ZHp) in the wild and domestic population, respectively, 78 and 54 were retained under more restrictive ZHp criteria. A total of 116 suggestive windows (representing 19 genomic regions) were linked to either QTL for production traits, or outliers for divergent or balancing selection. Twenty-four of them (representing 3 genomic regions) were retained under stricter ZHp thresholds. Eleven QTL/outlier markers were exclusively found in suggestive regions of the domestic broodstock, 7 in the wild population and one in both populations; one (broodstock) and two (wild) of those were found in significant regions retained under more restrictive ZHp criteria in the broodstock and the wild population, respectively. Genome mining and functional enrichment within regions associated with selective sweeps disclosed relevant genes and pathways related to aquaculture target traits, including growth and immune-related pathways, metabolism and response to hypoxia, which showcases how this genome atlas of genetic diversity can be a valuable resource to look for candidate genes related to natural or artificial selection in turbot populations., (Copyright © 2020 Aramburu, Ceballos, Casanova, Le Moan, Hemmer-Hansen, Bekkevold, Bouza and Martínez.)

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., Competing Interests: None declared, (© 2020 The Authors. Evolutionary Applications published by John Wiley & Sons Ltd.)

Background: Quantitative traits are typically considered to be under additive genetic control. Although there are indications that non-additive factors have the potential to contribute to trait variation, experimental demonstration remains scarce. Here, we investigated the genetic basis of growth in Atlantic salmon by exploiting the high level of genetic diversity and trait expression among domesticated, hybrid and wild populations., Results: After rearing fish in common-garden experiments under aquaculture conditions, we performed a variance component analysis in four mapping populations totaling ~ 7000 individuals from six wild, two domesticated and three F1 wild/domesticated hybrid strains. Across the four independent datasets, genome-wide significant quantitative trait loci (QTLs) associated with weight and length were detected on a total of 18 chromosomes, reflecting the polygenic nature of growth. Significant QTLs correlated with both length and weight were detected on chromosomes 2, 6 and 9 in multiple datasets. Significantly, epistatic QTLs were detected in all datasets., Discussion: The observed interactions demonstrated that the phenotypic effect of inheriting an allele deviated between half-sib families. Gene-by-gene interactions were also suggested, where the combined effect of two loci resulted in a genetic effect upon phenotypic variance, while no genetic effect was detected when the two loci were considered separately. To our knowledge, this is the first documentation of epistasis in a quantitative trait in Atlantic salmon. These novel results are of relevance for breeding programs, and for predicting the evolutionary consequences of domestication-introgression in wild populations.

Sandeels are an ecologically important group of fishes; they are a key part of the food chain serving as food for marine mammals, seabirds and fish. Sandeels are further targeted by a large industrial fishery, which has led to concern about ecosystem effects. In the North Sea, the lesser sandeel Ammodytes marinus is by far the most prevalent species of sandeel in the fishery. Management of sandeel in the North Sea plus the Kattegat is currently divided into seven geographical areas, based on subtle differences in demography, population dynamics and results from simulations of larval dispersal. However, little is known about the underlying genetic population structure. In this study, we used 2,522 SNPs derived from restriction site-associated DNA sequencing (RADseq) typed in 429 fish representing four main sandeel management areas. Our main results showed (a) a lack of a clear spatially defined genetic structure across the majority of genetic markers and (b) the existence of a group of at least 13 SNPs under strong linkage disequilibrium which together separate North Sea sandeel into three haplotype clusters, suggestive of one or more structural variants in the genome. Analyses of the spatial distribution of these putative structural variants suggest at least partial reproductive isolation of sandeel in the western management area along the Scottish coast, supporting a separate management. Our results highlight the importance of the application of a large number of markers to be able to detect weak patterns of differentiation. This study contributes to increasing the genetic knowledge of this important exploited species, and results can be used to improve our understanding of population dynamics and stock structure., Competing Interests: None declared., (© 2019 The Authors. Evolutionary Applications published by John Wiley & Sons Ltd.)

The salmonid fish Brown trout is iconic as a model for the application of conservation genetics to understand and manage local interspecific variation. However, there is still scant information about relationships between local and large-scale population structure, and to what extent geographical and environmental variables are associated with barriers to gene flow. We used information from 3,782 mapped SNPs developed for the present study and conducted outlier tests and gene-environment association (GEA) analyses in order to examine drivers of population structure. Analyses comprised >2,600 fish from 72 riverine populations spanning a central part of the species' distribution in northern Europe. We report hitherto unidentified genetic breaks in population structure, indicating strong barriers to gene flow. GEA loci were widely spread across genomic regions and showed correlations with climatic, abiotic and geographical parameters. In some cases, individual loci showed consistent GEA across the geographical regions Britain, Europe and Scandinavia. In other cases, correlations were observed only within a sub-set of regions, suggesting that locus-specific variation was associated with local processes. A paired-population sampling design allowed us to evaluate sampling effects on detection of outlier loci and GEA. Two widely applied methods for outlier detection ( pcadapt and bayescan ) showed low overlap in loci identified as statistical outliers across sub-sets of data. Two GEA analytical approaches (LFMM and RDA) showed good correspondence concerning loci associated with specific variables, but LFMM identified five times more statistically significant associations than RDA. Our results emphasize the importance of carefully considering the statistical methods applied for the hypotheses being tested in outlier analysis. Sampling design may have lower impact on results if the objective is to identify GEA loci and their population distribution. Our study provides new insights into trout populations, and results have direct management implications in serving as a tool for identification of conservation units., Competing Interests: None declared., (© 2019 The Authors. Evolutionary Applications published by John Wiley & Sons Ltd.)

The evolutionary process that occurs when a species colonizes a new environment provides an opportunity to explore the mechanisms underlying genetic adaptation, which is essential knowledge for understanding evolution and the maintenance of biodiversity. Atlantic herring has an estimated total breeding stock of about 1 trillion (10 12 ) and has colonized the brackish Baltic Sea within the last 10,000 y. Minute genetic differentiation between Atlantic and Baltic herring populations at selectively neutral loci combined with this rapid adaptation to a new environment facilitated the identification of hundreds of loci underlying ecological adaptation. A major question in the field of evolutionary biology is to what extent such an adaptive process involves selection of novel mutations with large effects or genetic changes at many loci, each with a small effect on phenotype (i.e., selection on standing genetic variation). Here we show that a missense mutation in rhodopsin (Phe261Tyr) is an adaptation to the red-shifted Baltic Sea light environment. The transition from phenylalanine to tyrosine differs only by the presence of a hydroxyl moiety in the latter, but this results in an up to 10-nm red-shifted light absorbance of the receptor. Remarkably, an examination of the rhodopsin sequences from 2,056 species of fish revealed that the same missense mutation has occurred independently and been selected for during at least 20 transitions between light environments across all fish. Our results provide a spectacular example of convergent evolution and how a single amino acid change can have a major effect on ecological adaptation., Competing Interests: The authors declare no conflict of interest., (Copyright © 2019 the Author(s). Published by PNAS.)

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., Competing Interests: None declared.

This work investigated whether multiple freshwater populations of three-spined stickleback Gasterosteus aculeatus in different freshwater catchments in the Jutland Peninsula, Denmark, derived from the same marine populations show repeated adaptive responses. A total of 327 G. aculeatus collected at 13 sampling locations were screened for genetic variation using a combination of 70 genes putatively under selection and 26 neutral genes along with a marker linked to the ectodysplasin gene (eda), which is strongly correlated with plate armour morphs in the species. A highly significant genetic differentiation was found that was higher among different freshwater samples than between marine-freshwater samples. Tests for selection between marine and freshwater populations showed a very low degree of parallelism and no single nucleotide polymorphism was detected as outlier in all freshwater-marine pairwise comparisons, including the eda. This suggests that G. aculeatus is not necessarily the prime example of parallel local adaptation suggested in much of the literature and that important exceptions exist (i.e. the Jutland Peninsula). While marine populations in the results described here showed a high phenotype-genotype correlation at eda, a low association was found for most of the freshwater populations. The most extreme case was found in the freshwater Lake Hald where all low-plated phenotypes were either homozygotes for the allele supposed to be associated with completely plated morphs or heterozygotes, but none were homozygotes for the putative low-plated allele. Re-examination of data from seven G. aculeatus studies agrees in showing a high but partial association between phenotype-genotype at eda in G. aculeatus freshwater populations and that mismatches occur everywhere in the European regions studied (higher in some areas, i.e. Denmark). This is independent of the eda marker used., (© 2017 The Fisheries Society of the British Isles.)

Background: Distinct hybrid zones and phenotypic and genomic divergence is often observed between marine and freshwater threespine sticklebacks (Gasterosteus aculeatus). Nevertheless, cases also exist where marine-freshwater divergence is diffuse despite seemingly similar environmental settings. In order to assess what characterizes these highly different outcomes, we focused on the latter kind of system in the Odder River, Denmark. Here, a previous study based on RAD (Restriction site Associated DNA) sequencing found non-significant genome-wide differentiation between marine and freshwater sticklebacks. In the present study, we analyzed samples on a finer geographical scale. We assessed if the system should be regarded as panmictic, or if fine-scale genetic structure and local selection was present but dominated by strong migration. We also asked if specific population components, that is the two sexes and different lateral plate morphs, contributed disproportionally more to dispersal., Results: We assessed variation at 96 SNPs and the Eda gene that affects lateral plate number, conducted molecular sex identification, and analyzed morphological traits. Genetic differentiation estimated by F ST was non-significant throughout the system. Nevertheless, spatial autocorrelation analysis suggested fine scale genetic structure with a genetic patch size of 770 m. There was no evidence for sex-biased dispersal, but full-plated individuals showed higher dispersal than low- and partial-plated individuals. The system was dominated by full-plated morphs characteristic of marine sticklebacks, but in the upstream part of the river body shape and frequency of low-plated morphs changed in the direction expected for freshwater sticklebacks. Five markers including Eda were under possible diversifying selection. However, only subtle clinal patterns were observed for traits and markers., Conclusions: We suggest that gene flow from marine sticklebacks overwhelms adaptation to freshwater conditions, but the short genetic patch size means that the effect of gene flow on the most upstream region must be indirect and occurs over generations. The occurrence of both weak unimodal and strong bimodal hybrid zones within the same species is striking. We suggest environmental and demographic factors that could determine these outcomes, but also highlight the possibility that long-term population history and the presence or absence of genomic incompatibilities could be a contributing factor.

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., (© 2016 The Fisheries Society of the British Isles.)

Herring, Clupea harengus, is one of the ecologically and commercially most important species in European northern seas, where two distinct ecotypes have been described based on spawning time; spring and autumn. To date, it is unknown if these spring and autumn spawning herring constitute genetically distinct units. We assessed levels of genetic divergence between spring and autumn spawning herring in the Baltic Sea using two types of DNA markers, microsatellites and Single Nucleotide Polymorphisms, and compared the results with data for autumn spawning North Sea herring. Temporally replicated analyses reveal clear genetic differences between ecotypes and hence support reproductive isolation. Loci showing non-neutral behaviour, so-called outlier loci, show convergence between autumn spawning herring from demographically disjoint populations, potentially reflecting selective processes associated with autumn spawning ecotypes. The abundance and exploitation of the two ecotypes have varied strongly over space and time in the Baltic Sea, where autumn spawners have faced strong depression for decades. The results therefore have practical implications by highlighting the need for specific management of these co-occurring ecotypes to meet requirements for sustainable exploitation and ensure optimal livelihood for coastal communities.

Background: Species showing complex life cycles provide excellent opportunities to study the genetic associations between life cycle stages, as selective pressures may differ before and after metamorphosis. The European eel presents a complex life cycle with two metamorphoses, a first metamorphosis from larvae into glass eels (juvenile stage) and a second metamorphosis into silver eels (adult stage). We tested the hypothesis that different genes and gene pathways will be under selection at different life stages when comparing the genetic associations between glass eels and silver eels., Results: We used two sets of markers to test for selection: first, we genotyped individuals using a panel of 80 coding-gene single nucleotide polymorphisms (SNPs) developed in American eel; second, we investigated selection at the genome level using a total of 153,423 RAD-sequencing generated SNPs widely distributed across the genome. Using the RAD approach, outlier tests identified a total of 2413 (1.57%) potentially selected SNPs. Functional annotation analysis identified signal transduction pathways as the most over-represented group of genes, including MAPK/Erk signalling, calcium signalling and GnRH (gonadotropin-releasing hormone) signalling. Many of the over-represented pathways were related to growth, while others could result from the different conditions that eels inhabit during their life cycle., Conclusions: The observation of different genes and gene pathways under selection when comparing glass eels vs. silver eels supports the adaptive decoupling hypothesis for the benefits of metamorphosis. Partitioning the life cycle into discrete morphological phases may be overall beneficial since it allows the different life stages to respond independently to their unique selection pressures. This might translate into a more effective use of food and niche resources and/or performance of phase-specific tasks (e.g. feeding in the case of glass eels, migrating and reproducing in the case of silver eels).

Hatching early in the season is often assumed to elevate fitness, particularly in cannibalistic fish in which size-dependent predation mortality is a major selective force. While the importance of the thermal environment for the growth of fish is undisputed, the relevance of maternal effects for offspring growth in the wild is largely unknown. Otoliths of 366 age-0 pike (Esox lucius L.) were sampled in a natural lake over three seasons. All offspring were assigned to more than 330 potential mothers using 16 informative microsatellites. We found temperature and past maternal environment (as represented by juvenile growth rate), but not female total length, to jointly contribute to explain within- and among-season size variation in juvenile pike. While there was no statistical evidence for maternal effects on offspring growth rate, fast female juvenile growth positively correlated with the offspring length in early summer. One mechanism could be related to fast-growing females spawning somewhat earlier in the season. However, the more likely mechanism emerging in our study was that fast-growing females could have been in better condition prior to spawning, in turn possibly producing higher numbers of high-quality eggs. Our study is among the few to reveal carry-over effects related to past maternal environments on offspring performance in a naturally reproducing fish stock. At the same time, our study underscores recent arguments that size-dependent maternal effects may not be expressed in the wild and that early hatching does not generally produce size advantages in light of stochastically varying temperature conditions.

Background: The threespine stickleback (Gasterosteus aculeatus) has become an important model species for studying both contemporary and parallel evolution. In particular, differential adaptation to freshwater and marine environments has led to high differentiation between freshwater and marine stickleback populations at the phenotypic trait of lateral plate morphology and the underlying candidate gene Ectodysplacin (EDA). Many studies have focused on this trait and candidate gene, although other genes involved in marine-freshwater adaptation may be equally important. In order to develop a resource for rapid and cost efficient analysis of genetic divergence between freshwater and marine sticklebacks, we generated a low-density SNP (Single Nucleotide Polymorphism) array encompassing markers of chromosome regions under putative directional selection, along with neutral markers for background., Results: RAD (Restriction site Associated DNA) sequencing of sixty individuals representing two freshwater and one marine population led to the identification of 33,993 SNP markers. Ninety-six of these were chosen for the low-density SNP array, among which 70 represented SNPs under putatively directional selection in freshwater vs. marine environments, whereas 26 SNPs were assumed to be neutral. Annotation of these regions revealed several genes that are candidates for affecting stickleback phenotypic variation, some of which have been observed in previous studies whereas others are new., Conclusions: We have developed a cost-efficient low-density SNP array that allows for rapid screening of polymorphisms in threespine stickleback. The array provides a valuable tool for analyzing adaptive divergence between freshwater and marine stickleback populations beyond the well-established candidate gene Ectodysplacin (EDA).

The two North Atlantic eel species, the European eel (Anguilla anguilla) and the American eel (Anguilla rostrata), spawn in partial sympatry in the Sargasso Sea, providing ample opportunity to interbreed. In this study, we used a RAD (Restriction site Associated DNA) sequencing approach to identify species-specific diagnostic single-nucleotide polymorphisms (SNPs) and design a low-density array that combined with screening of a diagnostic mitochondrial DNA marker. Eels from Iceland (N=159) and from the neighboring Faroe Islands (N=29) were genotyped, along with 94 larvae (49 European and 45 American eel) collected in the Sargasso Sea. Our SNP survey showed that the majority of Icelandic eels are pure European eels but there is also an important contribution of individuals of admixed ancestry (10.7%). Although most of the hybrids were identified as F1 hybrids from European eel female × American eel male crosses, backcrosses were also detected, including a first-generation backcross (F1 hybrid × pure European eel) and three individuals identified as second-generation backcrosses originating from American eel × F1 hybrid backcrosses interbreeding with pure European eels. In comparison, no hybrids were observed in the Faroe Islands, the closest bodies of land to Iceland. It is possible that hybrids show an intermediate migratory behaviour between the two parental species that ultimately brings hybrid larvae to the shores of Iceland, situated roughly halfway between the Sargasso Sea and Europe. Only two hybrids were observed among Sargasso Sea larvae, both backcrosses, but no F1 hybrids, that points to temporal variation in the occurrence of hybridization.

Local adaptation and its underlying molecular basis has long been a key focus in evolutionary biology. There has recently been increased interest in the evolutionary role of plasticity and the molecular mechanisms underlying local adaptation. Using transcriptome analysis, we assessed differences in gene expression profiles for three brown trout (Salmo trutta) populations, one resident and two anadromous, experiencing different temperature regimes in the wild. The study was based on an F2 generation raised in a common garden setting. A previous study of the F1 generation revealed different reaction norms and significantly higher QST than FST among populations for two early life-history traits. In the present study we investigated if genomic reaction norm patterns were also present at the transcriptome level. Eggs from the three populations were incubated at two temperatures (5 and 8 degrees C) representing conditions encountered in the local environments. Global gene expression for fry at the stage of first feeding was analysed using a 32k cDNA microarray. The results revealed differences in gene expression between populations and temperatures and population × temperature interactions, the latter indicating locally adapted reaction norms. Moreover, the reaction norms paralleled those observed previously at early life-history traits. We identified 90 cDNA clones among the genes with an interaction effect that were differently expressed between the ecologically divergent populations. These included genes involved in immune- and stress response. We observed less plasticity in the resident as compared to the anadromous populations, possibly reflecting that the degree of environmental heterogeneity encountered by individuals throughout their life cycle will select for variable level of phenotypic plasticity at the transcriptome level. Our study demonstrates the usefulness of transcriptome approaches to identify genes with different temperature reaction norms. The responses observed suggest that populations may vary in their susceptibility to climate change.

Shallow population structure is generally reported for most marine fish and explained as a consequence of high dispersal, connectivity and large population size. Targeted gene analyses and more recently genome-wide studies have challenged such view, suggesting that adaptive divergence might occur even when neutral markers provide genetic homogeneity across populations. Here, 381 SNPs located in transcribed regions were used to assess large- and fine-scale population structure in the European hake (Merluccius merluccius), a widely distributed demersal species of high priority for the European fishery. Analysis of 850 individuals from 19 locations across the entire distribution range showed evidence for several outlier loci, with significantly higher resolving power. While 299 putatively neutral SNPs confirmed the genetic break between basins (F(CT) = 0.016) and weak differentiation within basins, outlier loci revealed a dramatic divergence between Atlantic and Mediterranean populations (F(CT) range 0.275-0.705) and fine-scale significant population structure. Outlier loci separated North Sea and Northern Portugal populations from all other Atlantic samples and revealed a strong differentiation among Western, Central and Eastern Mediterranean geographical samples. Significant correlation of allele frequencies at outlier loci with seawater surface temperature and salinity supported the hypothesis that populations might be adapted to local conditions. Such evidence highlights the importance of integrating information from neutral and adaptive evolutionary patterns towards a better assessment of genetic diversity. Accordingly, the generated outlier SNP data could be used for tackling illegal practices in hake fishing and commercialization as well as to develop explicit spatial models for defining management units and stock boundaries., (© 2013 John Wiley & Sons Ltd.)

The genomic architecture underlying ecological divergence and ecological speciation with gene flow is still largely unknown for most organisms. One central question is whether divergence is genome-wide or localized in 'genomic mosaics' during early stages when gene flow is still pronounced. Empirical work has so far been limited, and the relative impacts of gene flow and natural selection on genomic patterns have not been fully explored. Here, we use ecotypes of Atlantic cod to investigate genomic patterns of diversity and population differentiation in a natural system characterized by high gene flow and large effective population sizes, properties which theoretically could restrict divergence in local genomic regions. We identify a genomic region of strong population differentiation, extending over approximately 20 cM, between pairs of migratory and stationary ecotypes examined at two different localities. Furthermore, the region is characterized by markedly reduced levels of genetic diversity in migratory ecotype samples. The results highlight the genomic region, or 'genomic island', as potentially associated with ecological divergence and suggest the involvement of a selective sweep. Finally, we also confirm earlier findings of localized genomic differentiation in three other linkage groups associated with divergence among eastern Atlantic populations. Thus, although the underlying mechanisms are still unknown, the results suggest that 'genomic mosaics' of differentiation may even be found under high levels of gene flow and that marine fishes may provide insightful model systems for studying and identifying initial targets of selection during ecological divergence., (© 2013 Blackwell Publishing Ltd.)