Limited Parallelism in Genetic Adaptation to Brackish Water Bodies in European Sprat and Atlantic Herring.

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.

Genetic response to human-induced habitat changes in the marine environment: A century of evolution of European sprat in Landvikvannet, Norway.

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.

Genetic analysis redraws the management boundaries for the European sprat.

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

Comparative biology and population mixing among local, coastal and offshore Atlantic herring (Clupea harengus) in the North Sea, Skagerrak, Kattegat and western Baltic.

The population structure of Atlantic herring (Clupea harengus) from 13 local, coastal and offshore areas of the North Sea, Skagerrak, Kattegat and western Baltic (northeast Atlantic) was studied using biological and environmental data from 1970-2015. The objective was to identify distinct populations by comparing variability in the temporal and spatial phenotypic characteristics and evaluate the potential for mixing of populations in time and space. The populations varied in biological characteristics such as mean vertebral counts (VS), growth and maturity ogives. Generalized additive models indicated temporally stable VS in the North Sea and western Baltic, whereas intra-annual temporal variation of VS occurred in other areas. High variability of VS within a population was not affected by environmental factors such as temperature and salinity. Consequently, seasonal VS variability can be explained by the presence or absence of herring populations as they migrate between areas. The three main populations identified in this paper correspond to the three managed stocks in this area: Norwegian spring spawners (NSS), western Baltic spring spawners (WBSS) and North Sea autumn spawners (NSAS). In addition, several local populations were identified in fjords or lakes along the coast, but our analyses could not detect direct mixing of local populations with the three main populations. Our results highlight the importance of recognizing herring dynamics and understanding the mixing of populations as a challenge for management of herring.

Seasonal dynamics of Atlantic herring (Clupea harengus L.) populations spawning in the vicinity of marginal habitats.

Gillnet sampling and analyses of otolith shape, vertebral count and growth indicated the presence of three putative Atlantic herring (Clupea harengus L.) populations mixing together over the spawning season February-June inside and outside an inland brackish water lake (Landvikvannet) in southern Norway. Peak spawning of oceanic Norwegian spring spawners and coastal Skagerrak spring spawners occurred in March-April with small proportions of spawners entering the lake. In comparison, spawning of Landvik herring peaked in May-June with high proportions found inside the lake, which could be explained by local adaptations to the environmental conditions and seasonal changes of this marginal habitat. The 1.85 km(2) lake was characterized by oxygen depletion occurring between 2.5 and 5 m depth between March and June. This was followed by changes in salinity from 1-7‰ in the 0-1 m surface layer to levels of 20-25‰ deeper than 10 m. In comparison, outside the 3 km long narrow channel connecting the lake with the neighboring fjord, no anoxic conditions were found. Here salinity in the surface layer increased over the season from 10 to 25‰, whereas deeper than 5 m it was stable at around 35‰. Temperature at 0-5 m depth increased significantly over the season in both habitats, from 7 to 14 °C outside and 5 to 17 °C inside the lake. Despite differences in peak spawning and utilization of the lake habitat between the three putative populations, there was an apparent temporal and spatial overlap in spawning stages suggesting potential interbreeding in accordance with the metapopulation concept.