Epistatic regulation of growth in Atlantic salmon revealed: a QTL study performed on the domesticated-wild interface.

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.

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

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

Outlier SNP markers reveal fine-scale genetic structuring across European hake populations (Merluccius merluccius).

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.

Introgression affects Salmo trutta juvenile life-history traits generations after stocking with non-native strains.

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.

A genomic island linked to ecotype divergence in Atlantic cod.

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.

Microevolution in time and space: SNP analysis of historical DNA reveals dynamic signatures of selection in Atlantic cod.

Little is known about how quickly natural populations adapt to changes in their environment and how temporal and spatial variation in selection pressures interact to shape patterns of genetic diversity. We here address these issues with a series of genome scans in four overfished populations of Atlantic cod (Gadus morhua) studied over an 80-year period. Screening of >1000 gene-associated single-nucleotide polymorphisms (SNPs) identified 77 loci that showed highly elevated levels of differentiation, likely as an effect of directional selection, in either time, space or both. Exploratory analysis suggested that temporal allele frequency shifts at certain loci may correlate with local temperature variation and with life history changes suggested to be fisheries induced. Interestingly, however, largely nonoverlapping sets of loci were temporal outliers in the different populations and outliers from the 1928 to 1960 period showed almost complete stability during later decades. The contrasting microevolutionary trajectories among populations resulted in sequential shifts in spatial outliers, with no locus maintaining elevated spatial differentiation throughout the study period. Simulations of migration coupled with observations of temporally stable spatial structure at neutral loci suggest that population replacement or gene flow alone could not explain all the observed allele frequency variation. Thus, the genetic changes are likely to at least partly be driven by highly dynamic temporally and spatially varying selection. These findings have important implications for our understanding of local adaptation and evolutionary potential in high gene flow organisms and underscore the need to carefully consider all dimensions of biocomplexity for evolutionarily sustainable management.

The memory remains: application of historical DNA for scaling biodiversity loss.

Few species worldwide have attracted as much attention in relation to conservation and sustainable management as Pacific salmon. Most populations have suffered significant reductions, many have disappeared, and even entire evolutionary significant units (ESUs) are believed to have been lost. Until now, no 'smoking gun' in terms of direct genetic evidence of the loss of a salmon ESU has been produced. In this issue of Molecular Ecology, Iwamoto et al. (2012) use microsatellite analysis of historical scale samples of Columbia River sockeye salmon (Oncorhynchus nerka) from 1924 (Fig. 1) to ask the pertinent question: Do the historical samples contain salmon from extirpated populations or ESUs? They identified four genetic groups in the historical samples of which two were almost genetically identical to contemporary ESUs in the river, one showed genetic relationship with a third ESU, but one group was not related to any of the contemporary populations. In association with ecological data, the genetic results suggest that an early migrating Columbia River headwater sockeye salmon ESU has been extirpated. The study has significant importance for conservation and reestablishment of sockeye populations in the Columbia River, but also underpins the general significance of shifting baselines in conservation biology, and how to assess loss of genetic biodiversity. The results clearly illustrate the huge and versatile potential of using historical DNA in population and conservation genetics. Because of the extraordinarily plentiful historical samples and rapid advances in fish genomics, fishes are likely to spearhead future studies of temporal ecological and population genomics in non-model organisms. [Figure: see text].

Monitoring adaptive genetic responses to environmental change.

Widespread environmental changes including climate change, selective harvesting and landscape alterations now greatly affect selection regimes for most organisms. How animals and plants can adapt to these altered environments via contemporary evolution is thus of strong interest. We discuss how to use genetic monitoring to study adaptive responses via repeated analysis of the same populations over time, distinguishing between phenotypic and molecular genetics approaches. After describing monitoring designs, we develop explicit criteria for demonstrating adaptive responses, which include testing for selection and establishing clear links between genetic and environmental change. We then review a few exemplary studies that explore adaptive responses to climate change in Drosophila, selective responses to hunting and fishing, and contemporary evolution in Daphnia using resurrected resting eggs. We further review a broader set of 44 studies to assess how well they meet the proposed criteria, and conclude that only 23% fulfill all criteria. Approximately half (43%) of these studies failed to rule out the alternative hypothesis of replacement by a different, better-adapted population. Likewise, 34% of the studies based on phenotypic variation did not test for selection as opposed to drift. These shortcomings can be addressed via improved experimental designs and statistical testing. We foresee monitoring of adaptive responses as a future valuable tool in conservation biology, for identifying populations unable to evolve at sufficiently high rates and for identifying possible donor populations for genetic rescue. Technological advances will further augment the realization of this potential, especially next-generation sequencing technologies that allow for monitoring at the level of whole genomes.

Environmental selection on transcriptome-derived SNPs in a high gene flow marine fish, the Atlantic herring (Clupea harengus).

High gene flow is considered the norm for most marine organisms and is expected to limit their ability to adapt to local environments. Few studies have directly compared the patterns of differentiation at neutral and selected gene loci in marine organisms. We analysed a transcriptome-derived panel of 281 SNPs in Atlantic herring (Clupea harengus), a highly migratory small pelagic fish, for elucidating neutral and selected genetic variation among populations and to identify candidate genes for environmental adaptation. We analysed 607 individuals from 18 spawning locations in the northeast Atlantic, including two temperature clines (5-12 °C) and two salinity clines (5-35‰). By combining genome scan and landscape genetic analyses, four genetically distinct groups of herring were identified: Baltic Sea, Baltic-North Sea transition area, North Sea/British Isles and North Atlantic; notably, samples exhibited divergent clustering patterns for neutral and selected loci. We found statistically strong evidence for divergent selection at 16 outlier loci on a global scale, and significant correlations with temperature and salinity at nine loci. On regional scales, we identified two outlier loci with parallel patterns across temperature clines and five loci associated with temperature in the North Sea/North Atlantic. Likewise, we found seven replicated outliers, of which five were significantly associated with low salinity across both salinity clines. Our results reveal a complex pattern of varying spatial genetic variation among outlier loci, likely reflecting adaptations to local environments. In addition to disclosing the fine scale of local adaptation in a highly vagile species, our data emphasize the need to preserve functionally important biodiversity.

Retrospective genomics highlights changes in genetic composition of tiger sharks (Galeocerdo cuvier) and potential loss of a south-eastern Australia population.

Over the last century, many shark populations have declined, primarily due to overexploitation in commercial, artisanal and recreational fisheries. In addition, in some locations the use of shark control programs also has had an impact on shark numbers. Still, there is a general perception that populations of large ocean predators cover wide areas and therefore their diversity is less susceptible to local anthropogenic disturbance. Here we report on temporal genomic analyses of tiger shark (Galeocerdo cuvier) DNA samples that were collected from eastern Australia over the past century. Using Single Nucleotide Polymorphism (SNP) loci, we documented a significant change in genetic composition of tiger sharks born between ~1939 and 2015. The change was most likely due to a shift over time in the relative contribution of two well-differentiated, but hitherto cryptic populations. Our data strongly indicate a dramatic shift in the relative contribution of these two populations to the overall tiger shark abundance on the east coast of Australia, possibly associated with differences in direct or indirect exploitation rates.

Hierarchical genetic structure in an evolving species complex: Insights from genome wide ddRAD data in Sebastes mentella.

The diverse biology and ecology of marine organisms may lead to complex patterns of intraspecific diversity for both neutral and adaptive genetic variation. Sebastes mentella displays a particular life-history as livebearers, for which existence of multiple ecotypes has been suspected to complicate the genetic population structure of the species. Double digest restriction-site associated DNA was used to investigate genetic population structure in S. mentella and to scan for evidence of selection. In total, 42,288 SNPs were detected in 277 fish, and 1,943 neutral and 97 tentatively adaptive loci were selected following stringent filtration. Unprecedented levels of genetic differentiation were found among the previously defined 'shallow pelagic', 'deep pelagic' and 'demersal slope' ecotypes, with overall mean FST = 0.05 and 0.24 in neutral and outlier SNPs, respectively. Bayesian computation estimated a concurrent and historical divergence among these three ecotypes and evidence of local adaptation was found in the S. mentella genome. Overall, these findings imply that the depth-defined habitat divergence of S. mentella has led to reproductive isolation and possibly adaptive radiation among these ecotypes. Additional sub-structuring was detected within the 'shallow' and 'deep' pelagic ecotypes. Population assignment of individual fish showed more than 94% agreement between results based on SNP and previously generated microsatellite data, but the SNP data provided a lower estimate of hybridization among the ecotypes than that by microsatellite data. We identified a SNP panel with only 21 loci to discriminate populations in mixed samples based on a machine-learning algorithm. This first SNP based investigation clarifies the population structure of S. mentella, and provides novel and high-resolution genomic tools for future investigations. The insights and tools provided here can readily be incorporated into the management of S. mentella and serve as a template for other exploited marine species exhibiting similar complex life history traits.

Remote, autonomous real-time monitoring of environmental DNA from commercial fish.

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.

Genomic signatures of local directional selection in a high gene flow marine organism; the Atlantic cod (Gadus morhua).

BACKGROUND: Marine fishes have been shown to display low levels of genetic structuring and associated high levels of gene flow, suggesting shallow evolutionary trajectories and, possibly, limited or lacking adaptive divergence among local populations. We investigated variation in 98 gene-associated single nucleotide polymorphisms (SNPs) for evidence of selection in local populations of Atlantic cod (Gadus morhua L.) across the species distribution. RESULTS: Our global genome scan analysis identified eight outlier gene loci with very high statistical support, likely to be subject to directional selection in local demes, or closely linked to loci under selection. Likewise, on a regional south/north transect of central and eastern Atlantic populations, seven loci displayed strongly elevated levels of genetic differentiation. Selection patterns among populations appeared to be relatively widespread and complex, i.e. outlier loci were generally not only associated with one of a few divergent local populations. Even on a limited geographical scale between the proximate North Sea and Baltic Sea populations four loci displayed evidence of adaptive evolution. Temporal genome scan analysis applied to DNA from archived otoliths from a Faeroese population demonstrated stability of the intra-population variation over 24 years. An exploratory landscape genetic analysis was used to elucidate potential effects of the most likely environmental factors responsible for the signatures of local adaptation. We found that genetic variation at several of the outlier loci was better correlated with temperature and/or salinity conditions at spawning grounds at spawning time than with geographic distance per se. CONCLUSION: These findings illustrate that adaptive population divergence may indeed be prevalent despite seemingly high levels of gene flow, as found in most marine fishes. Thus, results have important implications for our understanding of the interplay of evolutionary forces in general, and for the conservation of marine biodiversity under rapidly increasing evolutionary pressure from climate and fisheries induced changes in local environments.

Population genomics of marine fishes: identifying adaptive variation in space and time.

Studies of adaptive evolution have experienced a recent revival in population genetics of natural populations and there is currently much focus on identifying genomic signatures of selection in space and time. Insights into local adaptation, adaptive response to global change and evolutionary consequences of selective harvesting can be generated through population genomics studies, allowing the separation of the effects invoked by neutral processes (drift-migration) from those due to selection. Such knowledge is important not only for improving our basic understanding of natural as well as human-induced evolutionary processes, but also for predicting future trajectories of biodiversity and for setting conservation priorities. Marine fishes possess a number of features rendering them well suited for providing general insights into adaptive genomic evolution in natural populations. These include well-described population structures, substantial and rapidly developing genomic resources and abundant archived samples enabling temporal studies. Furthermore, superior possibilities for conducting large-scale experiments under controlled conditions, due to the economic resources provided by the large and growing aquaculture industry, hold great promise for utilizing recent technological developments. Here, we review achievements in marine fish genomics to date and highlight potential avenues for future research, which will provide both general insights into evolution in high gene flow species, as well as specific knowledge which can lead to improved management of marine organisms.

Author Correction: Gene-associated markers provide tools for tackling illegal fishing and false eco-certification.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

Interpopulation differences in expression of candidate genes for salinity tolerance in winter migrating anadromous brown trout (Salmo trutta L.).

BACKGROUND: Winter migration of immature brown trout (Salmo trutta) into freshwater rivers has been hypothesized to result from physiologically stressful combinations of high salinity and low temperature in the sea. RESULTS: We sampled brown trout from two Danish populations entering different saline conditions and quantified expression of the hsp70 and Na/K-ATPases alpha 1b genes following acclimation to freshwater and full-strength seawater at 2 degrees C and 10 degrees C. An interaction effect of low temperature and high salinity on expression of both hsp70 and Na/K-ATPase alpha 1b was found in trout from the river entering high saline conditions, while a temperature independent up-regulation of both genes in full-strength seawater was found for trout entering marine conditions with lower salinities. CONCLUSION: Overall our results support the hypothesis that physiologically stressful conditions in the sea drive sea-run brown trout into freshwater rivers in winter. However, our results also demonstrate intra-specific differences in expression of important stress and osmoregulative genes most likely reflecting adaptive differences between trout populations on a regional scale, thus strongly suggesting local adaptations driven by the local marine environment.

Geographic extent of introgression in Sebastes mentella and its effect on genetic population structure.

Genetic population structure is often used to identify management units in exploited species, but the extent of genetic differentiation may be inflated by geographic variation in the level of hybridization between species. We identify the genetic population structure of Sebastes mentella and investigate possible introgression within the genus by analyzing 13 microsatellites in 2,562 redfish specimens sampled throughout the North Atlantic. The data support an historical divergence between the "shallow" and "deep" groups, beyond the Irminger Sea where they were described previously. A third group, "slope," has an extended distribution on the East Greenland Shelf, in addition to earlier findings on the Icelandic slope. Furthermore, S. mentella from the Northeast Arctic and Northwest Atlantic waters are genetically different populations. In both areas, interspecific introgression may influence allele frequency differences among populations. Evidence of introgression was found for almost all the identified Sebastes gene pools, but to a much lower extent than suggested earlier. Greenland waters appear to be a sympatric zone for many of the genetically independent Sebastes groups. This study illustrates that the identified groups maintain their genetic integrity in this region despite introgression.

Population structure and connectivity of tiger sharks (Galeocerdo cuvier) across the Indo-Pacific Ocean basin.

Population genetic structure using nine polymorphic nuclear microsatellite loci was assessed for the tiger shark (Galeocerdo cuvier) at seven locations across the Indo-Pacific, and one location in the southern Atlantic. Genetic analyses revealed considerable genetic structuring (FST > 0.14, p < 0.001) between all Indo-Pacific locations and Brazil. By contrast, no significant genetic differences were observed between locations from within the Pacific or Indian Oceans, identifying an apparent large, single Indo-Pacific population. A lack of differentiation between tiger sharks sampled in Hawaii and other Indo-Pacific locations identified herein is in contrast to an earlier global tiger shark nDNA study. The results of our power analysis provide evidence to suggest that the larger sample sizes used here negated any weak population subdivision observed previously. These results further highlight the need for cross-jurisdictional efforts to manage the sustainable exploitation of large migratory sharks like G. cuvier.

The use of archived tags in retrospective genetic analysis of fish.

Collections of historical tissue samples from fish (e.g. scales and otoliths) stored in museums and fisheries institutions are precious sources of DNA for conducting retrospective genetic analysis. However, in some cases, only external tags used for documentation of spatial dynamics of fish populations have been preserved. Here, we test the usefulness of fish tags as a source of DNA for genetic analysis. We extract DNA from historical tags from cod collected in Greenlandic waters between 1950 and 1968. We show that the quantity and quality of DNA recovered from tags is comparable to DNA from archived otoliths from the same individuals. Surprisingly, levels of cross-contamination do not seem to be significantly higher in DNA from external (tag) than internal (otolith) sources. Our study therefore demonstrates that historical tags can be a highly valuable source of DNA for retrospective genetic analysis of fish.

Gene-associated markers provide tools for tackling illegal fishing and false eco-certification.

Illegal, Unreported and Unregulated fishing has had a major role in the overexploitation of global fish populations. In response, international regulations have been imposed and many fisheries have been 'eco-certified' by consumer organizations, but methods for independent control of catch certificates and eco-labels are urgently needed. Here we show that, by using gene-associated single nucleotide polymorphisms, individual marine fish can be assigned back to population of origin with unprecedented high levels of precision. By applying high differentiation single nucleotide polymorphism assays, in four commercial marine fish, on a pan-European scale, we find 93-100% of individuals could be correctly assigned to origin in policy-driven case studies. We show how case-targeted single nucleotide polymorphism assays can be created and forensically validated, using a centrally maintained and publicly available database. Our results demonstrate how application of gene-associated markers will likely revolutionize origin assignment and become highly valuable tools for fighting illegal fishing and mislabelling worldwide.