Temporally stable genetic structure of heavily exploited Atlantic herring (Clupea harengus) in Swedish waters

Genetic Stock Identification Reveals Mismatches Between Management Areas and Population Genetic Structure in a Migratory Pelagic Fish

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.

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.

Microsatellites as Molecular Markers with Applications in Exploitation and Conservation of Aquatic Animal Populations

A large number of species and taxa has been studied for genetic polymorphism. Microsatellites have been known as hypervariable neutral molecular markers with the highest resolution power in comparison with any other markers. However, the discovery of a new type of molecular marker—single nucleotide polymorphism (SNP) has put the existing applications of microsatellites to the test. To ensure good resolution power in studies of populations and individuals, a number of microsatellite loci from 14 to 20 was often used, which corresponds to about 200 independent alleles. Recently, these numbers have tended to be increased by the application of genomic sequencing of expressed sequence tags (ESTs), and the choice of the most informative loci for genotyping depends on the aims of research. Examples of successful applications of microsatellite molecular markers in aquaculture, fisheries, and conservation genetics in comparison with SNPs have been summarized in this review. Microsatellites can be considered superior markers in such topics as kinship and parentage analysis in cultured and natural populations, the assessment of gynogenesis, androgenesis and ploidization. Microsatellites can be coupled with SNPs for mapping QTL. Microsatellites will continue to be used in research on genetic diversity in cultured stocks, and also in natural populations as an economically advantageous genotyping technique.

Ecological adaptation in Atlantic herring is associated with large shifts in allele frequencies at hundreds of loci

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.

Time matters: genetic composition and evaluation of effective population size in temperate coastal fish species

Background

Extensive knowledge on the genetic characterization of marine organisms has been assembled, mainly concerning the spatial distribution and structuring of populations. Temporal monitoring assesses not only the stability in genetic composition but also its trajectory over time, providing critical information for the accurate forecast of changes in genetic diversity of marine populations, particularly important for both fisheries and endangered species management. We assessed fluctuations in genetic composition among different sampling periods in the western Portuguese shore in three fish species.

Methods

White seabream Diplodus sargus, sand smelt Atherina presbyter and shanny Lipophrys pholis were chosen, because of their genetic patterns in distinct ecological environments, insight into historical and contemporary factors influencing population effective size (N_e), and degree of commercial exploitation. Samples were obtained near Lisbon between 2003 and 2014 and screened for genetic variation with mitochondrial and nuclear markers. Analyses included genealogies, genetic diversities, temporal structures and contemporary N_e.

Results

For mtDNA no temporal structure was detected, while for nDNA significant differences were recorded between some sampling periods for the shanny and the sand smelt. Haplotype networks revealed deep genealogies, with various levels of diversification. The shanny revealed a smaller N_e/generation when compared to the other species, which, in turn, revealed no evidence of genetic drift for most study periods. These results highlight the fact that temporal variations in genetic pool composition should be considered when evaluating the population structure of fish species with long distance dispersal, which are more vulnerable to recruitment fluctuations.

A chromosome-level assembly of the Atlantic herring genome-detection of a supergene and other signals of selection

The Atlantic herring is a model species for exploring the genetic basis for ecological adaptation, due to its huge population size and extremely low genetic differentiation at selectively neutral loci. However, such studies have so far been hampered because of a highly fragmented genome assembly. Here, we deliver a chromosome-level genome assembly based on a hybrid approach combining a de novo Pacific Biosciences (PacBio) assembly with Hi-C-supported scaffolding. The assembly comprises 26 autosomes with sizes ranging from 12.4 to 33.1 Mb and a total size, in chromosomes, of 726 Mb, which has been corroborated by a high-resolution linkage map. A comparison between the herring genome assembly with other high-quality assemblies from bony fishes revealed few inter-chromosomal but frequent intra-chromosomal rearrangements. The improved assembly facilitates analysis of previously intractable large-scale structural variation, allowing, for example, the detection of a 7.8-Mb inversion on Chromosome 12 underlying ecological adaptation. This supergene shows strong genetic differentiation between populations. The chromosome-based assembly also markedly improves the interpretation of previously detected signals of selection, allowing us to reveal hundreds of independent loci associated with ecological adaptation.

Temporal stability and assignment power of adaptively divergent genomic regions between herring (Clupea harengus) seasonal spawning aggregations

Atlantic herring (Clupea harengus), a vital ecosystem component and target of the largest Northwest Atlantic pelagic fishery, undergo seasonal spawning migrations that result in elusive sympatric population structure. Herring spawn mostly in fall or spring, and genomic differentiation was recently detected between these groups. Here we used a subset of this differentiation, 66 single nucleotide polymorphisms (SNPs) to analyze the temporal dynamics of this local adaptation and the applicability of SNP subsets in stock assessment. We showed remarkable temporal stability of genomic differentiation corresponding to spawning season, between samples taken a decade apart (2005 N = 90 vs. 2014 N = 71) in the Gulf of St. Lawrence, and new evidence of limited interbreeding between spawning components. We also examined an understudied and overexploited herring population in Bras d'Or lake (N = 97); using highly reduced SNP panels (N SNPs > 6), we verified little-known sympatric spawning populations within this unique inland sea. These results describe consistent local adaptation, arising from asynchronous reproduction in a migratory and dynamic marine species. Our research demonstrates the efficiency and precision of SNP-based assessments of sympatric subpopulations; and indeed, this temporally stable local adaptation underlines the importance of such fine-scale management practices.

Decadal stability in genetic variation and structure in the intertidal seaweed Fucus serratus (Heterokontophyta: Fucaceae)

BACKGROUND

The spatial distribution of genetic diversity and structure has important implications for conservation as it reveals a species' strong and weak points with regard to stability and evolutionary capacity. Temporal genetic stability is rarely tested in marine species other than commercially important fishes, but is crucial for the utility of temporal snapshots in conservation management. High and stable diversity can help to mitigate the predicted northward range shift of seaweeds under the impact of climate change. Given the key ecological role of fucoid seaweeds along rocky shores, the positive effect of genetic diversity may reach beyond the species level to stabilize the entire intertidal ecosystem along the temperate North Atlantic. In this study, we estimated the effective population size, as well as temporal changes in genetic structure and diversity of the seaweed F. serratus using 22 microsatellite markers. Samples were taken across latitudes and a range of temperature regimes at seven locations with decadal sampling (2000 and 2010).

RESULTS

Across latitudes, genetic structure and diversity remained stable over 5-10 generations. Stable small-scale structure enhanced regional diversity throughout the species' range. In accordance with its biogeographic history, effective population size and diversity peaked in the species' mid-range in Brittany (France), and declined towards its leading and trailing edge to the north and south. At the species' southern edge, multi-locus-heterozygosity displayed a strong decline from 1999 to 2010.

CONCLUSION

Temporally stable genetic structure over small spatial scales is a potential driver for local adaptation and species radiation in the genus Fucus. Survival and adaptation of the low-diversity leading edge of F. serratus may be enhanced by regional gene flow and 'surfing' of favorable mutations or impaired by the accumulation of deleterious mutations. Our results have clear implications for the conservation of F. serratus at its genetically unique southern edge in Northwest Iberia, where increasing temperatures are likely the major cause for the decline not only of F. serratus, but also other intertidal and subtidal macroalgae. We expect that F. serratus will disappear from Northwest Iberia by 2100 if genetic rescue is not induced by the influx of genetic variation from Brittany.

Whole-genome sequencing approaches for conservation biology: Advantages, limitations and practical recommendations

Whole-genome resequencing (WGR) is a powerful method for addressing fundamental evolutionary biology questions that have not been fully resolved using traditional methods. WGR includes four approaches: the sequencing of individuals to a high depth of coverage with either unresolved or resolved haplotypes, the sequencing of population genomes to a high depth by mixing equimolar amounts of unlabelled-individual DNA (Pool-seq) and the sequencing of multiple individuals from a population to a low depth (lcWGR). These techniques require the availability of a reference genome. This, along with the still high cost of shotgun sequencing and the large demand for computing resources and storage, has limited their implementation in nonmodel species with scarce genomic resources and in fields such as conservation biology. Our goal here is to describe the various WGR methods, their pros and cons and potential applications in conservation biology. WGR offers an unprecedented marker density and surveys a wide diversity of genetic variations not limited to single nucleotide polymorphisms (e.g., structural variants and mutations in regulatory elements), increasing their power for the detection of signatures of selection and local adaptation as well as for the identification of the genetic basis of phenotypic traits and diseases. Currently, though, no single WGR approach fulfils all requirements of conservation genetics, and each method has its own limitations and sources of potential bias. We discuss proposed ways to minimize such biases. We envision a not distant future where the analysis of whole genomes becomes a routine task in many nonmodel species and fields including conservation biology.

Temporal genetic stability in natural populations of the waterflea Daphnia magna in response to strong selection pressure

Studies monitoring changes in genetic diversity and composition through time allow a unique understanding of evolutionary dynamics and persistence of natural populations. However, such studies are often limited to species with short generation times that can be propagated in the laboratory or few exceptional cases in the wild. Species that produce dormant stages provide powerful models for the reconstruction of evolutionary dynamics in the natural environment. A remaining open question is to what extent dormant egg banks are an unbiased representation of populations and hence of the species' evolutionary potential, especially in the presence of strong environmental selection. We address this key question using the water flea Daphnia magna, which produces dormant stages that accumulate in biological archives over time. We assess temporal genetic stability in three biological archives, previously used in resurrection ecology studies showing adaptive evolutionary responses to rapid environmental change. We show that neutral genetic diversity does not decline with the age of the population and it is maintained in the presence of strong selection. In addition, by comparing temporal genetic stability in hatched and unhatched populations from the same biological archive, we show that dormant egg banks can be consulted to obtain a reliable measure of genetic diversity over time, at least in the multidecadal time frame studied here. The stability of neutral genetic diversity through time is likely mediated by the buffering effect of the resting egg bank.

The genetic basis for ecological adaptation of the Atlantic herring revealed by genome sequencing

Ecological adaptation is of major relevance to speciation and sustainable population management, but the underlying genetic factors are typically hard to study in natural populations due to genetic differentiation caused by natural selection being confounded with genetic drift in subdivided populations. Here, we use whole genome population sequencing of Atlantic and Baltic herring to reveal the underlying genetic architecture at an unprecedented detailed resolution for both adaptation to a new niche environment and timing of reproduction. We identify almost 500 independent loci associated with a recent niche expansion from marine (Atlantic Ocean) to brackish waters (Baltic Sea), and more than 100 independent loci showing genetic differentiation between spring- and autumn-spawning populations irrespective of geographic origin. Our results show that both coding and non-coding changes contribute to adaptation. Haplotype blocks, often spanning multiple genes and maintained by selection, are associated with genetic differentiation.

Outlier Loci Detect Intraspecific Biodiversity amongst Spring and Autumn Spawning Herring across Local Scales

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.

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.

Temporal stability of genetic variability and differentiation in the three-spined stickleback (Gasterosteus aculeatus)

Temporal variation in allele frequencies, whether caused by deterministic or stochastic forces, can inform us about interesting demographic and evolutionary phenomena occurring in wild populations. In spite of the continued surge of interest in the genetics of three-spined stickleback (Gasterosteus aculeatus) populations, little attention has been paid towards the temporal stability of allele frequency distributions, and whether there are consistent differences in effective size (Ne) of local populations. We investigated temporal stability of genetic variability and differentiation in 15 microsatellite loci within and among eight collection sites of varying habitat type, surveyed twice over a six-year time period. In addition, Nes were estimated with the expectation that they would be lowest in isolated ponds, intermediate in larger lakes and largest in open marine sites. In spite of the marked differences in genetic variability and differentiation among the study sites, the temporal differences in allele frequencies, as well as measures of genetic diversity and differentiation, were negligible. Accordingly, the Ne estimates were temporally stable, but tended to be lower in ponds than in lake or marine habitats. Hence, we conclude that allele frequencies in putatively neutral markers in three-spined sticklebacks seem to be temporally stable - at least over periods of few generations - across a wide range of habitat types differing markedly in levels of genetic variability, effective population size and gene flow.

Conservation genetics of threatened Hippocampus guttulatus in vulnerable habitats in NW Spain: temporal and spatial stability of wild populations with flexible polygamous mating system in captivity

This study was focused on conservation genetics of threatened Hippocampus guttulatus on the Atlantic coast of NW Iberian Peninsula. Information about spatial structure and temporal stability of wild populations was obtained based on microsatellite markers, and used for monitoring a captive breeding program firstly initiated in this zone at the facilities of the Institute of Marine Research (Vigo, Spain). No significant major genetic structure was observed regarding the biogeographical barrier of Cape Finisterre. However, two management units under continuous gene flow are proposed based on the allelic differentiation between South-Atlantic and Cantabrian subpopulations, with small to moderate contemporary effective size based on single-sample methods. Temporal stability was observed in South-Atlantic population samples of H. guttulatus for the six-year period studied, suggesting large enough effective population size to buffer the effects of genetic drift within the time frame of three generations. Genetic analysis of wild breeders and offspring in captivity since 2009 allowed us to monitor the breeding program founded in 2006 in NW Spain for this species. Similar genetic diversity in the renewed and founder broodstock, regarding the wild population of origin, supports suitable renewal and rearing processes to maintain genetic variation in captivity. Genetic parentage proved single-brood monogamy in the wild and in captivity, but flexible short- and long-term mating system under captive conditions, from strict monogamy to polygamy within and/or among breeding seasons. Family analysis showed high reproductive success in captivity under genetic management assisted by molecular relatedness estimates to avoid inbreeding. This study provides genetic information about H. guttulatus in the wild and captivity within an uncovered geographical range for this data deficient species, to be taken into account for management and conservation purposes.

Genetic similarity of island populations of tent caterpillars during successive outbreaks

Cyclic or fluctuating populations experience regular periods of low population density. Genetic bottlenecks during these periods could give rise to temporal or spatial genetic differentiation of populations. High levels of movement among increasing populations, however, could ameliorate any differences and could also synchronize the dynamics of geographically separated populations. We use microsatellite markers to investigate the genetic differentiation of four island and one mainland population of western tent caterpillars, Malacosoma californicum pluviale, in two periods of peak or pre-peak density separated by 8 years. Populations showed high levels of genetic variation and little genetic differentiation either temporally between peaks or spatially among sites. Mitochondrial haplotypes were also shared between one island population and one mainland population in the two years studied. An isolation-by-distance analysis showed the FST values of the two geographically closest populations to have the highest level of differentiation in both years. We conclude that high levels of dispersal among populations maintain both synchrony of population dynamics and override potential genetic differentiation that might occur during population troughs. As far we are aware, this is the first time that genetic similarity between temporally separated population outbreaks in insects has been investigated. A review of genetic data for both vertebrate and invertebrate species of cyclic animals shows that a lack of spatial genetic differentiation is typical, and may result from high levels of dispersal associated with fluctuating dynamics.

Oceanographic connectivity and environmental correlates of genetic structuring in Atlantic herring in the Baltic Sea

Marine fish often show little genetic structuring in neutral marker genes, and Atlantic herring (Clupea harengus) in the Baltic Sea are no exception; historically, very low levels of population differentiation (F ST ≍ 0.002) have been found, despite a high degree of interpopulation environmental heterogeneity in salinity and temperature. Recent exome sequencing and SNP studies have however shown that many loci are under selection in this system. Here, we combined population genetic analyses of a large number of transcriptome-derived microsatellite markers with oceanographic modelling to investigate genetic differentiation and connectivity in Atlantic herring at a relatively fine scale within the Baltic Sea. We found evidence for weak but robust and significant genetic structuring (F ST = 0.008) explainable by oceanographic connectivity. Genetic differentiation was also associated with site differences in temperature and salinity, with the result driven by the locus Her14 which appears to be under directional selection (F ST = 0.08). The results show that Baltic herring are genetically structured within the Baltic Sea, and highlight the role of oceanography and environmental factors in explaining this structuring. The results also have implications for the management of herring fisheries, the most economically important fishery in the Baltic Sea, suggesting that the current fisheries management units may be in need of revision.

Genetic and life-history changes associated with fisheries-induced population collapse

Over the recent years, growing number of studies suggests that intensive size-selective fishing can cause evolutionary changes in life-history traits in the harvested population, which can have drastic negative effects on populations, ecosystems and fisheries. However, most studies to date have overlooked the potential role of immigration of fish with different phenotypes as an alternative plausible mechanism behind observed phenotypic trends. Here, we investigated the evolutionary consequences of intensive fishing simultaneously at phenotypic and molecular level in Eurasian perch (Perca fluviatilis L.) population in the Baltic Sea over a 24-year period. We detected marked changes in size- and age-distributions and increase in juvenile growth rate. We also observed reduction of age at sexual maturity in males that has frequently been considered to support the hypothesis of fisheries-induced evolution. However, combined individual-based life-history and genetic analyses indicated increased immigration of foreign individuals with different life-history patterns as an alternative mechanism behind the observed phenotypic change. This study demonstrates the value of combining genetic and phenotypic analyses and suggests that replacement or breakdown of locally adapted gene complexes may play important role in impeding the recovery of fish populations.

Whole mitochondrial genome scan for population structure and selection in the Atlantic herring

BACKGROUND

Marine fish, such as the Atlantic herring (Clupea harengus), often show a low degree of differentiation over large geographical regions. Despite strong environmental gradients (salinity and temperature) in the Baltic Sea, population genetic studies have shown little genetic differentiation among herring in this area, but some evidence for environmentally-induced selection has been uncovered. The mitochondrial genome is a likely target for selection in this system due to its functional role in metabolism.

RESULTS

We sequenced whole mitochondrial genomes for herring from throughout the Baltic region (n=98) in order to investigate evidence for geographical structuring, selection, and associations between genetic and environmental variation. Three well-supported clades that predate the formation of the Baltic Sea were identified, but geographic structuring of this variation was weak (ΦST = 0.036). There was evidence for significant positive selection, particularly in the ND2, ND4 and ND5 genes, and amino acids under significant selection in these genes explained some of the clade formation. Despite uncovering evidence for selection, correlations between genetic diversity or differentiation with environmental factors (temperature, salinity, latitude) were weak.

CONCLUSIONS

The results indicate that most of the current mtDNA diversity in herring predates the formation of the Baltic Sea, and that little structuring has evolved since. Thus, fisheries management units in this region cannot be determined on the basis of mtDNA variability. Preliminary evidence for selection underlying clade formation indicates that the NADH complex may be useful for examining adaptation and population structuring at a broader geographical scale.

Population-scale sequencing reveals genetic differentiation due to local adaptation in Atlantic herring

The Atlantic herring (Clupea harengus), one of the most abundant marine fishes in the world, has historically been a critical food source in Northern Europe. It is one of the few marine species that can reproduce throughout the brackish salinity gradient of the Baltic Sea. Previous studies based on few genetic markers have revealed a conspicuous lack of genetic differentiation between geographic regions, consistent with huge population sizes and minute genetic drift. Here, we present a cost-effective genome-wide study in a species that lacks a genome sequence. We first assembled a muscle transcriptome and then aligned genomic reads to the transcripts, creating an "exome assembly," capturing both exons and flanking sequences. We then resequenced pools of fish from a wide geographic range, including the Northeast Atlantic, as well as different regions in the Baltic Sea, aligned the reads to the exome assembly, and identified 440,817 SNPs. The great majority of SNPs showed no appreciable differences in allele frequency among populations; however, several thousand SNPs showed striking differences, some approaching fixation for different alleles. The contrast between low genetic differentiation at most loci and striking differences at others implies that the latter category primarily reflects natural selection. A simulation study confirmed that the distribution of the fixation index F(ST) deviated significantly from expectation for selectively neutral loci. This study provides insights concerning the population structure of an important marine fish and establishes the Atlantic herring as a model for population genetic studies of adaptation and natural selection.

Temporal genetic stability and high effective population size despite fisheries-induced life-history trait evolution in the North Sea sole

Heavy fishing and other anthropogenic influences can have profound impact on a species' resilience to harvesting. Besides the decrease in the census and effective population size, strong declines in mature adults and recruiting individuals may lead to almost irreversible genetic changes in life-history traits. Here, we investigated the evolution of genetic diversity and effective population size in the heavily exploited sole (Solea solea), through the analysis of historical DNA from a collection of 1379 sole otoliths dating back from 1957. Despite documented shifts in life-history traits, neutral genetic diversity inferred from 11 microsatellite markers showed a remarkable stability over a period of 50 years of heavy fishing. Using simulations and corrections for fisheries induced demographic variation, both single-sample estimates and temporal estimates of effective population size (N(e) ) were always higher than 1000, suggesting that despite the severe census size decrease over a 50-year period of harvesting, genetic drift is probably not strong enough to significantly decrease the neutral diversity of this species in the North Sea. However, the inferred ratio of effective population size to the census size (N(e) /N(c) ) appears very small (10(-5) ), suggesting that overall only a low proportion of adults contribute to the next generation. The high N(e) level together with the low N(e) /N(c) ratio is probably caused by a combination of an equalized reproductive output of younger cohorts, a decrease in generation time and a large variance in reproductive success typical for marine species. Because strong evolutionary changes in age and size at first maturation have been observed for sole, changes in adaptive genetic variation should be further monitored to detect the evolutionary consequences of human-induced selection.

The future of Baltic Sea populations: local extinction or evolutionary rescue?

Environmental change challenges local and global survival of populations and species. In a species-poor environment like the Baltic Sea this is particularly critical as major ecosystem functions may be upheld by single species. A complex interplay between demographic and genetic characteristics of species and populations determines risks of local extinction, chances of re-establishment of lost populations, and tolerance to environmental changes by evolution of new adaptations. Recent studies show that Baltic populations of dominant marine species are locally adapted, have lost genetic variation and are relatively isolated. In addition, some have evolved unusually high degrees of clonality and others are representatives of endemic (unique) evolutionary lineages. We here suggest that a consequence of local adaptation, isolation and genetic endemism is an increased risk of failure in restoring extinct Baltic populations. Additionally, restricted availability of genetic variation owing to lost variation and isolation may negatively impact the potential for evolutionary rescue following environmental change.

Detecting population structure in a high gene-flow species, Atlantic herring (Clupea harengus): direct, simultaneous evaluation of neutral vs putatively selected loci

In many marine fish species, genetic population structure is typically weak because populations are large, evolutionarily young and have a high potential for gene flow. We tested whether genetic markers influenced by natural selection are more efficient than the presumed neutral genetic markers to detect population structure in Atlantic herring (Clupea harengus), a migratory pelagic species with large effective population sizes. We compared the spatial and temporal patterns of divergence and statistical power of three traditional genetic marker types, microsatellites, allozymes and mitochondrial DNA, with one microsatellite locus, Cpa112, previously shown to be influenced by divergent selection associated with salinity, and one locus located in the major histocompatibility complex class IIA (MHC-IIA) gene, using the same individuals across analyses. Samples were collected in 2002 and 2003 at two locations in the North Sea, one location in the Skagerrak and one location in the low-saline Baltic Sea. Levels of divergence for putatively neutral markers were generally low, with the exception of single outlier locus/sample combinations; microsatellites were the most statistically powerful markers under neutral expectations. We found no evidence of selection acting on the MHC locus. Cpa112, however, was highly divergent in the Baltic samples. Simulations addressing the statistical power for detecting population divergence showed that when using Cpa112 alone, compared with using eight presumed neutral microsatellite loci, sample sizes could be reduced by up to a tenth while still retaining high statistical power. Our results show that the loci influenced by selection can serve as powerful markers for detecting population structure in high gene-flow marine fish species.