Genomic signatures of local directional selection in a high gene flow marine organism; the Atlantic cod (Gadus morhua)
- Einar E Nielsen†1Email author,
- Jakob Hemmer-Hansen†1,
- Nina A Poulsen1, 2,
- Volker Loeschcke2,
- Thomas Moen3,
- Torild Johansen4,
- Christian Mittelholzer5, 8,
- Geir-Lasse Taranger4,
- Rob Ogden6 and
- Gary R Carvalho7
© Nielsen et al; licensee BioMed Central Ltd. 2009
Received: 8 July 2009
Accepted: 1 December 2009
Published: 1 December 2009
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.
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.
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.
When genotype by environment interactions take place, populations can evolve traits that confer a Darwinian fitness advantage in their local habitat. The process and the resulting patterns are termed "local adaptation" . For adaptive divergence of populations to take place, the evolutionary force of directional selection should be stronger than random genetic drift and the homogenising effect of migration among populations . Accordingly, local adaptations are expected to be rare or absent in species with extensive gene flow, such as many marine fishes . Even though some groups of marine fishes, such as coral reef fishes, have been found to be highly genetically structured , more than 40 years of population genetic research has typically demonstrated low levels of genetic differentiation among local demes compared to other fish species and terrestrial organisms .
It is notoriously difficult to separate genetic from environmental effects when traits of presumed adaptive value are measured directly in the wild . Furthermore, disentangling genetic from environmental effects on phenotypes is a challenge in relatively large aquatic organisms such as many marine fish species, where experimental designs such as common garden and reciprocal transplantation require large facilities. Alternatively, genome scan approaches using a high number of genetic markers in natural populations can provide a powerful shortcut for demonstrating local selection pressures by allowing the identification of outlier loci with divergent levels of genetic differentiation [7, 8]. In particular, targeted approaches investigating gene-associated rather than random DNA marker variation holds great promise . Targeted genome scans have recently been used in natural populations of various organisms such as trees (white spruce)  and fish (sticklebacks) .
The Atlantic cod (Gadus morhua L.) has been the subject of numerous population genetic studies and thus represents one of the best studied marine fish species within the field. Population genetic studies have demonstrated low, albeit significant genetic differentiation on various hierarchical levels, ranging from clear transatlantic differences [12–14] to microgeographical population structure of demes separated by a few tens of kilometres [15–17]. At the same time there is growing evidence that commonly observed variation in life history traits among putative natural populations of Atlantic cod is at least partly determined by genetic rather than environmental differences, see  and references therein. Earlier studies have indicated that functional genes may be under selection in natural populations of Atlantic cod. Already in the 1960's the classical studies of haemoglobin polymorphisms in cod by Sick  identified a locus (Hb I) most likely under environmental selection. Recently, Andersen et al.  identified the molecular basis as well as the different oxygen binding properties of the haemoglobin alleles, thereby providing the important link between genotypes and phenotypes of adaptive importance. Another example is the identification of the Pantophysin (Pan I) locus in cod  which has been shown to display markedly elevated levels of genetic differentiation compared to presumed neutral markers, patterns which are potentially driven by temperature selection [22, 23]. A more extensive survey of genes was recently conducted by Moen et al. , who reported on the first genome scan approach in the species. They assessed variation in 318 Single nucleotide polymorphisms (SNPs) in a selection of samples from two cod populations (northeast Arctic and Norwegian coastal cod) and identified a number of outlier loci as candidates for being under adaptive evolution. However, in summation, these studies have focused on single genes and/or a restricted geographical area not using random population samples. Consequently, even though Atlantic cod represents one of the most frequently targeted marine species in population genetic research, the genomic architecture as well as the geographical scale and distribution of adaptive divergence in this and other high gene flow marine fishes is still relatively poorly known (see also review by Nielsen et al. ).
In this study, we tested the hypothesis that populations of Atlantic cod are adapted to local environmental conditions and therefore display genomic signatures of divergent selection. We applied genome scan methods to identify outlier loci potentially associated with adaptive population divergence on global, regional and local scales as well as an explorative landscape genetic analysis focusing on environmental factors suspected to be involved in local adaptation.
Genome scan approach
Atlantic cod samples analysed in the present study.
54.51° N 15.28° E
54.53° N 13.33° E
Western Baltic Sea3
54.56° N 12.28° E
55.11° N 10.28° E
57.15° N 11.35° E
Central North Sea1,2,3
55.17° N 03.39° E
Southern North Sea2
54.29° N 0.02° E
50.47° N 0.29° E
Northern North Sea2
58.00° N 03.00° W
Northeast Arctic cod, Lofoten1,2
68.35° N 12.13° E
Norwegian coastal cod, Lofoten1,2
68.12° N 14.44° E
62.53° N 06.18° W
60.56° N 08.52° W
63.49° N 21.05° W
66.17° N 15.45° W
66.35° N 53.32° W
Western Atlantic (Canada)1
47.39° N 55.24° W
Significant associations between candidate loci for local adaptation and environmental variables.
This study suggests that loci subject to directional selection can be found on various geographical scales in Atlantic cod. In addition, the results suggest that a large number of loci could be subject to directional selection on a genomic scale among local populations. Local directional selection appears to be general and relatively widespread and can be found on a number of geographical scales since the global or regional outliers identified are not exclusively dependent on one or a few particularly divergent cod populations. Furthermore, although it is difficult to completely disentangle the effects of geographic and environmental distance, a higher proportion of the identified outlier loci were associated with environmental parameters on spawning grounds than with geographic variables. This suggests, but does not establish, these environmental factors as potentially responsible for adaptive divergence among cod populations.
The remaining "neutral" SNP markers generally exhibited very low levels of genetic differentiation and, accordingly, limited resolution of patterns of population relationships, except for the highly divergent western Atlantic and Baltic Sea cod populations. This is somewhat in contrast with evidence from microsatellite markers (e.g. see [12, 14, 27]), showing higher levels of differentiation but also high variance among loci. Although great care should be taken when comparing levels of genetic structure between markers with different mutational properties [28, 29], this difference may suggest that estimates of genetic differentiation from a limited number of microsatellite loci may be even more inflated by hitch-hiking selection than previously assumed. Alternatively, levels of SNP differentiation may be depressed due to the presence of SNPs, which are monomorphic in a number of population samples. This in turn may also explain why we find no evidence for balancing selection in cod. Simulations have shown that scenarios with low levels of population structure provide very low power for detecting loci under balancing selection .
The apparent high number of loci with strong statistical support for selection may be somewhat unexpected given the very low level of neutral genetic differentiation among cod populations. However, neutral divergence has in general been shown to be a poor predictor of adaptive trait variation . The proportion of loci subject to selection resembles the findings using genome scan methods for highly structured Atlantic salmon (Salmo salar) populations . The result that the North Sea/Baltic Sea cod population comparison only yielded four outliers likely to be subject to directional selection may be surprising given the unique low saline, low oxygen environment for cod in the Baltic Sea. However, although adaptations in Baltic Sea cod may be essential, they may be distributed over few genes, also considering the relatively short history of the Baltic Sea . Alternatively, we cannot exclude that ascertainment bias might have affected results of analyses including Baltic Sea samples (see also below). More detailed studies of the genetic architecture of adaptive traits in Baltic Sea cod are warranted to elucidate these issues. It should be stressed, however, that results from this study should not be taken as a direct estimate of the proportion of genes under divergent (and balancing) selection in different cod populations. Ascertainment bias  could potentially have affected the geographical distribution of outliers, since the majority of loci were ascertained from EST libraries consisting of only Norwegian coastal cod  and the candidate loci originate from a different and wider ascertainment sample. This could potentially result in higher variability, different levels of genetic differentiation and a different frequency of outliers for comparisons including Norwegian coastal cod and/or our candidate loci. However, marine fishes have been shown to display relatively shallow genetic structure compared to other organisms including other fish species , suggesting that ascertainment bias may be less of a problem here . This theoretical consideration is supported by the fact that estimates of heterozygosity and proportion of monomorphic loci were not markedly different in the sample of Norwegian coastal cod (see additional file 2) compared to samples from other populations, except for the most differentiated populations of Canadian cod and Baltic Sea cod. This phenomenon of shared polymorphisms could be caused by the fact that Norwegian coastal cod is closely related to the ancestral population of cod  and therefore serve as a superior ascertainment sample (see  for an example for Drosophila). Furthermore, overall genetic differentiation of "neutral" loci for the EST derived loci and our own candidate loci was almost identical (illustrated in Figure 2), suggesting limited systematic difference between loci from the two ascertainment panels. In conclusion, we only think that ascertainment bias is likely to have significantly reduced the identified number of loci under selection for the very divergent Canadian cod. However, even if the proportion of outliers in the most divergent populations within the remaining cluster (i.e. northeast Arctic cod and Baltic Sea cod) might be mildly underestimated, we expect ascertainment bias to play a relatively minor role on regional and local scales based on theoretical expectations and our own evaluation of data.
Another potential concern in relation to establishing outliers as products of directional selection is the high number of loci tested (98). With an alpha level of 0.05 we would expect five significant outliers, which in fact is not so far from the present results. However, the statistical support for individual outliers identified here is generally much stronger. We should expect only one outlier by chance with a Bayes factors above 2 (log10), but the computed values are almost exclusively above 2 (log10), which is generally interpreted as "decisive". Four loci have log10 values as high as 5, which was the value ascribed to posterior probabilities of 1 (Bayes factor is infinity). Accordingly, we are still convinced that the majority of these outliers are affected by directional selection and not just outliers by chance. We think that our evidence for selection is as strong as possibly achievable given the statistical method and the number of loci at hand.
As evidenced above, the identified outlier genes originate both from randomly chosen gene associated SNPs and from specifically targeted SNPs in candidate genes of known functions. Interestingly, a relatively high proportion of candidate genes showed evidence of being subject to selection compared to the randomly chosen SNPs (app. 30% versus app. 6%). The identified outlier genes code for proteins with highly divergent functions (additional file 4), suggesting that directional selection for local adaptations in cod occur along multiple environmental dimensions. For example, heat-shock protein genes, represented by Hsp90, have been shown to be involved in various stress responses in fish including temperature, salinity and pollution . Aromatase plays a regulatory role in sex determination, gametogenesis, central nervous system development and sex behaviour  and rhodopsin is a pigment gene involved in the formation of photoreceptor cells and the perception of light . Alternatively, since temperature is important for proper functioning of all proteins and physiological processes [40, 41] temperature differences among geographical locations alone could be driving the evolutionary response for all genes in concert. The landscape genetic analysis partly supports this hypothesis, since almost all outlier loci appear to have alleles significantly associated with temperature. However, salinity also emerges as a likely driver of directional selection for local adaptation in Atlantic cod. Both temperature and salinity have been associated with phenotypic variation in cod. For instance, body shape variation was observed among juvenile cod originating from different Canadian spawning populations when reared under different temperature conditions . Likewise, adaptive trait variation in response to salinity has been strongly suggested for cod in the Baltic Sea, showing divergent egg buoyancy and sperm mobility . Still, the landscape genetic analysis should only be considered explorative. We have investigated a small subset of environmental factors potentially affecting allele frequencies at outliers. Likewise, most parameters are correlated to some extent, leading to significant associations with both geography and environment. This is illustrated by the positive correlations with environment for a number of genes not identified as outliers. It is not a simple task to unambiguously establish the link between genomic signatures of local adaptation and specific environmental conditions, particularly if only correlation analyses are applied. More research is needed on cod reared under controlled conditions in order to establish the genetic architecture of traits subject to local directional selection leading to fitness difference between native and non-native fish .
Previously, genetic evidence of local adaptation at the DNA level in marine fishes has been inferred from single genes [20, 44, 45], or from restricted geographical areas . To our knowledge, the present targeted genome scan approach represents the first attempt to elucidate genomic signatures of directional selection in natural populations of a marine fish on various geographical scales across its range. Even though our understanding of the genetic architecture of adaptive evolution in marine fishes is not as advanced as in other well studied fish species (e.g. sticklebacks  or whitefish ), our findings strongly suggest that marine fishes are not only isolated into local populations, even on relatively small geographical scales . Despite being connected by variable levels of gene flow, these populations can indeed follow semi-independent adaptive evolutionary trajectories shaped by selection by their local environments. Thus, adaptive population divergence seems to be possible and may even be prevalent despite seemingly high levels of gene flow often found in marine fishes. Our findings have implications for sustainable management of marine fishes by underpinning the local population as the focal management unit. Thus, we cannot rely on immigration to rehabilitate declining populations, since non-native individuals will have lower fitness and therefore are less likely to be successful in their new environment. Moreover, the extirpation of local populations may represent irreversible changes to the gene pool and associated adaptive evolutionary potential . Recent assertions of global warming causing recruitment failure in cold-adapted North Sea cod emphasize the importance of maintaining the adaptive potential of exploited species . Likewise, predictions on future distribution and abundance of marine fishes should not rely on simplified spatial models viewing whole species as single units with global bioclimatic niches. In order to improve the predictive power of future responses to environmental change, information on the dynamics of locally adapted populations, as well as on population-specific adaptive trait values, will be required.
In total 708 adult individuals were analysed from 17 spatial samples and one temporal sample, all consisting of approximately 40 individuals (Table 1 and Figure 1). Spatial analyses were conducted using samples from three different hierarchical levels; global (CAN, GRE, ICS, NCC, NEAC, ECH, CNS and BAS), a regional south-north transect of central and north-eastern Atlantic population samples (ECH, SNC, CNS, NNS, FBA, FPL, NCC NEAC, ICS, ICN and GRE) and a North Sea - Baltic Sea transect (CNS, KAT, BES, WBA, ARK and BAS). The temporal comparison was conducted on samples from the Faeroe Bank in the eastern Atlantic (FBA78 and FBA).
Novel SNPs analysed in this study
genbank accession no.
Follicle stimulating hormone
Haemoglobin alpha subunit
Haemoglobin beta subunit
Heat shock protein 90
Lactate Dehydrogenase B
Ribosomal protein S2
Ribosomal protein S2
Observed and expected levels of heterozygosity and tests for conformance to Hardy-Weinberg Equilibrium following  were calculated for each locus and population in ARLEQUIN ver. 3.1 . Summary statistics for each locus are presented in additional file 2. Overall and locus-specific FST values were estimated by Weir and Cockerham's θ  in GENEPOP 4.0.10 . Pairwise single locus FST between samples used in the global analyses were estimated for each of the loci identified as outliers in BAYESCAN. In addition, pairwise multi-locus FST were estimated for all loci, excluding the outlier loci (see additional file 1). Likewise, multi-locus "neutral" pairwise FST were estimated between all samples and used to generate a MDS plot (Figure 3) using the program Vista 5.6.3. .
To detect signatures of natural selection we used the Bayesian likelihood method implemented via reversible jump Markov Chain Monte Carlo in BAYESCAN . The Bayesian regression approach implemented in BAYESCAN has several advantages over the widely used approach based on summary statistics implemented in FDIST2 . Firstly, in BAYESCAN FST is modelled using a logistic regression model implementing a locus effect and a population effect, while simulations in FDIST2 are based on an Island Model. Thus, BAYESCAN allow population specific FST in contrast to the symmetrical Island Model in FDIST2, and therefore BAYESCAN should be appropriate for Atlantic cod, which is expected to display non-symmetrical patterns of gene flow among natural populations. Secondly, and supporting this theoretical reasoning, analyses of simulated data have shown that the Bayesian regression approach performs slightly better than FDIST2 when scenarios deviate from the standard Island Model [56, 57]. Thirdly, FDIST2 simulates a distribution of FST from the empirical estimate, which is itself influenced by loci potentially under selection. This problem can be alleviated by iteratively removing outlier loci before running the programme [30, 56]. However, this approach is subject to a certain degree of subjectivity, which is avoided in the Bayesian approach in BAYESCAN. BAYESCAN directly estimates the posterior probability that a locus is under selection (see below) in contrast to an earlier method build on the same basic regression model, which provided an approximated p-value . In short, BAYESCAN estimates the probability that a locus is under selection by calculating a Bayes factor, which is simply the ratio of the posterior probabilities of two models (selection/neutral) given the data. A Bayes factor between 32 and 100 (log10 = 1.5 - 2) is "very strong evidence" of different statistical support for the two models and corresponds to a posterior probability between 0.97-0.99. For Bayes factors above 100 (log10 > 2) the evidence is interpreted as "decisive" and correspond to posterior probabilities between0.99 and 1.
In BAYESCAN, monomorphic loci were removed from each run. Following 10 pilot runs of 5000 iterations and an additional burn-in of 50000 iterations, we used 100000 iterations (sample size of 5000 and thinning interval of 20) to identify loci under selection from locus specific Bayes factors. A Bayes factor of infinity, corresponding to a posterior probability of 1, was assigned a log10 value of 5. In order to alleviate potential problems arising from allele frequency correlations between samples [58, 59], we ran a global analysis with eight representative samples covering the distributional area of the species rather than including all available samples. In addition, we ran regional analyses with more samples from restricted geographical areas.
A spatial analysis, using the program SAM , was performed in order to identify associations between alleles and environmental variables. Two environmental parameters, temperature and salinity, as well as latitude and longitude were assessed (see below). We chose to use temperature and salinity at spawning time (but see also additional file 3 for analyses applying a wider set of closely related environmental variables) since this is as close as possible to the most critical time of the cod life-cycle, i.e. when the life stages (i.e. egg and larvae) have little option for actively avoiding unfavourable conditions (see also discussion in ). Since the early life stages of Atlantic cod are pelagic, we would expect a relatively strong correlation between surface/near surface data and the actual environmental conditions experienced by these life stages at most sampling locations. Temperature data were obtained from NOAAs Optimal Interpolation version 2 monthly SST analyses (OI.v2 monthly Sea Surface Temperature, available from the IRI/LDEO Climate Data Library at http://ingrid.ldgo.columbia.edu). Salinity data were collected from NOAAs Global Ocean Data Assimilation System (GODAS monthly below sea level salinity at a depth of 5 m, http://ingrid.ldgo.columbia.edu) and from the oceanographic database managed by the International Council for the Exploration of the Sea (ICES surface salinity data, available from http://www.ices.dk). Environmental data were integrated in a Geographic Information System (ArcGIS 9 from ESRI) in order to extract data covering a latitude/longitude grid of 0.5° *0.5° to 1° *1° around each sample position. Environmental variables used in SAM were means of a sampling period from 1990 to 2007. The null hypothesis of no association between an allele and an environmental variable was rejected if the examined variable explained the observed allelic distribution better than a model with a constant only. We changed the critical P value from 0.05 to 0.01 and 0.001 for the implemented Wald test to demonstrate robustness of significant associations. Since only two alleles were present at each locus, leaving only one independent variable, we present the results as associations between locus and environmental variable.
The authors thank the participants in the European Union networks and projects, "Marine Genomics Europe" (Grant number: GOCE-CT-2004-505403) "Seafood Plus" (Grant number: 506359) and "FinE" (Grant number: 044276) for input and fruitful discussions. We are grateful to Eydfinn Magnussen and Christophe Pampoulie for providing samples and to Thomas Damm Als for help with graphics.
- Williams GC: Adaptation and Natural Selection. 1966, Princeton University Press, PrincetonGoogle Scholar
- Kawecki TJ, Ebert D: Conceptual issues in local adaptation. Ecol Lett. 2004, 7: 1225-1241. 10.1111/j.1461-0248.2004.00684.x.View ArticleGoogle Scholar
- Conover DO, Clarke LM, Munch SB, Wagner GN: Spatial and temporal scales of adaptive divergence in marine fishes and the implications for conservation. J Fish Biol. 2006, 6: 21-47. 10.1111/j.1095-8649.2006.01274.x.View ArticleGoogle Scholar
- Palumbi SR, Warner RR: Why gobies are like hobbits. Science. 2003, 299: 51-52. 10.1126/science.1080775.View ArticlePubMedGoogle Scholar
- Waples RS: Separating the wheat from the chaff: patterns of genetic differentiation in high gene flow species. J Hered. 1998, 89: 439-450. 10.1093/jhered/89.5.438.View ArticleGoogle Scholar
- Endler JA: Natural selection in the wild. 1986, Princeton University Press, PrincetonGoogle Scholar
- Beaumont MA: Adaptation and speciation: what can F-st tell us?. Trends Ecol Evol. 2005, 20: 435-440. 10.1016/j.tree.2005.05.017.View ArticlePubMedGoogle Scholar
- Storz JF: Using genome scans of DNA polymorphism to infer adaptive population divergence. Mol Ecol. 2005, 14: 671-688. 10.1111/j.1365-294X.2005.02437.x.View ArticlePubMedGoogle Scholar
- Bonin A: Population genomics: a new generation of genome scans to bridge the gap with functional genomics. Mol Ecol. 2008, 17: 3583-3584. 10.1111/j.1365-294X.2008.03854.x.View ArticlePubMedGoogle Scholar
- Namroud MC, Beaulieu J, Juge N, Laroche J, Bousquet J: Scanning the genome for gene single nucleotide polymorphisms involved in adaptive population differentiation in white spruce. Mol Ecol. 2008, 17: 3599-3613. 10.1111/j.1365-294X.2008.03840.x.PubMed CentralView ArticlePubMedGoogle Scholar
- Mäkinen HS, Cano JM, Merilä J: Identifying footprints of directional and balancing selection in marine and freshwater three-spined stickleback (Gasterosteus aculeatus) populations. Mol Ecol. 2008, 17: 3565-3582. 10.1111/j.1365-294X.2008.03714.x.View ArticlePubMedGoogle Scholar
- Bentzen P, Taggart CT, Ruzzante DE, Cook D: Microsatellite polymorphism and the population structure of Atlantic cod (Gadus morhua) in the northwest Atlantic. Can J Fish Aquat Sci. 1996, 53: 2706-2721. 10.1139/cjfas-53-12-2706.View ArticleGoogle Scholar
- Arnason E: Mitochondrial cytochrome b DNA variation in the high-fecundity Atlantic cod: Trans-atlantic clines and shallow gene genealogy. Genetics. 2004, 166: 1871-1885. 10.1534/genetics.166.4.1871.PubMed CentralView ArticlePubMedGoogle Scholar
- Nielsen EE, Hansen MM, Meldrup D: Evidence of microsatellite hitch-hiking selection in Atlantic cod (Gadus morhua L.): Implications for inferring population structure in non-model organisms. Mol Ecol. 2006, 15: 3219-3229. 10.1111/j.1365-294X.2006.03025.x.View ArticlePubMedGoogle Scholar
- Ruzzante DE, Taggart CT, Cook DA: Nuclear DNA basis for shelf- and bank scale population structure in northwest Atlantic cod (Gadus morhua): Labrador to George Bank. Mol Ecol. 1998, 7: 1663-1680. 10.1046/j.1365-294x.1998.00497.x.View ArticleGoogle Scholar
- Knutsen H, Jorde PE, André C, Stenseth NC: Fine-scaled geographical population structuring in highly mobile marine species: the Atlantic cod. Mol Ecol. 2003, 12: 385-394. 10.1046/j.1365-294X.2003.01750.x.View ArticlePubMedGoogle Scholar
- Nielsen EE, Wright PJ, Hemmer-Hansen J, Poulsen NA, Gibb IM, Meldrup D: Population structure of cod (Gadus morhua L.) in the North Sea and west of Scotland; the role of sampling loci and individuals for inferring micro-geographical differentiation. Mar Ecol Prog Ser. 2008, 276: 213-225.View ArticleGoogle Scholar
- Hutchings JA, Fraser DJ: The nature of fisheries- and farming-induced evolution. Mol Ecol. 2008, 17: 294-313. 10.1111/j.1365-294X.2007.03485.x.View ArticlePubMedGoogle Scholar
- Sick K: Haemoglobin polymorphism of cod in the Baltic and the Danish Belt Sea. Hereditas. 1965, 54: 49-73. 10.1111/j.1601-5223.1965.tb02005.x.View ArticlePubMedGoogle Scholar
- Andersen Ø, Wetten OF, De Rosa MC, André C, Alinovi CC, Colafranceschi M, Brix O, Colosimo A: Haemoglobin polymorphisms affect the oxygen-binding properties in Atlantic cod populations. Proc R Soc Lond B Bio. 2009, 276: 833-841. 10.1098/rspb.2008.1529.View ArticleGoogle Scholar
- Pogson GH: Nucleotide polymorphism and natural selection at the pantophysin (Pan I) locus in the Atlantic cod, Gadus morhua (L.). Genetics. 2001, 157: 317-330.PubMed CentralPubMedGoogle Scholar
- Case RAJ, Hutchinson WF, Hauser L, Van Oosterhout C, Carvalho GR: Macro- and micro-geographic variation in pantophysin (Pan I) allele frequencies in NE Atlantic cod Gadus morhua . Mar Ecol Prog Ser. 2005, 301: 267-278. 10.3354/meps301267.View ArticleGoogle Scholar
- Pampoulie C, Ruzzante DE, Chosson V, Jorundsdottir TD, Taylor L, Thorsteinsson V, Danielsdottir AK, Marteinsdottir G: The genetic structure of Atlantic cod (Gadus morhua) around Iceland: insight from microsatellites, the Pan I locus, and tagging experiments. Can J Fish Aquat Sci. 2006, 63: 2660-2674. 10.1139/F06-150.View ArticleGoogle Scholar
- Moen T, Hayes B, Nilsen F, Delghandi M, Fjalestad KT, Fevolden SE, Berg PR, Lien S: Identification and characterisation of novel SNP markers in Atlantic cod: Evidence for directional selection. BMC Genet. 2008, 9: 18-10.1186/1471-2156-9-18.PubMed CentralView ArticlePubMedGoogle Scholar
- Nielsen EE, Hemmer-Hansen J, Larsen PF, Bekkevold D: Population genomics of marine fishes: identifying adaptive variation in space and time. Mol Ecol. 2009, 18: 3128-3150. 10.1111/j.1365-294X.2009.04272.x.View ArticlePubMedGoogle Scholar
- Johannesson K, André C: Life on the margin: genetic isolation and diversity loss in a peripheral marine ecosystem, the Baltic Sea. Mol Ecol. 2006, 15: 2013-2029. 10.1111/j.1365-294X.2006.02919.x.View ArticlePubMedGoogle Scholar
- O'Leary DB, Coughlan J, Dillane E, McCarthy TV, Cross TF: Microsatellite variation in cod Gadus morhua throughout its geographic range. J Fish Biol. 2007, 70: 310-335. 10.1111/j.1095-8649.2007.01451.x.View ArticleGoogle Scholar
- Hedrick PW: A standardized genetic differentiation measure. Evolution. 2005, 59: 1633-1638.View ArticlePubMedGoogle Scholar
- Jost L: G(ST) and its relatives do not measure differentiation. Mol Ecol. 2008, 17: 4015-4026. 10.1111/j.1365-294X.2008.03887.x.View ArticlePubMedGoogle Scholar
- Foll M, Gaggiotti O: A Genome-scan method to identify selected loci appropriate for both dominant and codominant markers: A Bayesian perspective. Genetics. 2008, 180: 977-993. 10.1534/genetics.108.092221.PubMed CentralView ArticlePubMedGoogle Scholar
- Leinonen T, O'hara RB, Cano JM, Merila J: Comparative studies of quantitative trait and neutral marker divergence: a meta-analysis. J Evolution Biol. 2008, 21: 1-17.Google Scholar
- Vasemägi A, Nilsson J, Primmer CR: Expressed sequence tag-linked microsatellites as a source of gene-associated polymorphisms for detecting signatures of divergent selection in Atlantic salmon (Salmo salar L.). Mol Biol Evol. 2005, 22: 1067-1076. 10.1093/molbev/msi093.View ArticlePubMedGoogle Scholar
- Clark AG, Hubisz MJ, Bustamante CD, Williamson SH, Nielsen R: Ascertainment bias in studies of human genome-wide polymorphism. Genome Res. 2005, 15: 1496-1502. 10.1101/gr.4107905.PubMed CentralView ArticlePubMedGoogle Scholar
- Rosenblum EB, Novembre J: Ascertainment bias in spatially structured populations: A case study in the eastern fence lizard. J Hered. 2007, 98: 331-336. 10.1093/jhered/esm031.View ArticlePubMedGoogle Scholar
- Bigg GR, Cunningham CW, Ottersen G, Pogson GH, Wadley MR, Williamson P: Ice-age survival of Atlantic cod: agreement between palaeoecology models and genetics. Proc R Soc B-Biol Sci. 2008, 275: 163-173. 10.1098/rspb.2007.1153.View ArticleGoogle Scholar
- Schlötterer C, Harr B: Single nucleotide polymorphisms derived from ancestral populations show no evidence for biased diversity estimates in Drosophila melanogaster. Mol Ecol. 2002, 11: 947-950. 10.1046/j.1365-294X.2002.01491.x.View ArticlePubMedGoogle Scholar
- Basu N, Todgham AE, Ackerman PA, Bibeau MR, Nakano K, Schulte PM, Iwama GK: Heat shock protein genes and their functional significance in fish. Gene. 2002, 295: 173-183. 10.1016/S0378-1119(02)00687-X.View ArticlePubMedGoogle Scholar
- Piferer F, Blázquez M: Aromatase distribution and regulation in fish. Fish Physiol Biochem. 2005, 31: 215-226. 10.1007/s10695-006-0027-0.View ArticleGoogle Scholar
- Nickle B, Robinson PR: The opsins of the vertebrate retina; insights from structural, biochemical and evolutionary studies. Cell Mol Life Sci. 2007, 64: 2917-2932. 10.1007/s00018-007-7253-1.View ArticlePubMedGoogle Scholar
- Pörtner HO, Knust R: Climate change affects marine fishes through the oxygen limitation of thermal tolerance. Science. 2007, 315: 95-97. 10.1126/science.1135471.View ArticlePubMedGoogle Scholar
- Pörtner HO, Bock C, Knust R, Lannig G, Lucassen M, Mark FC, Sartoris FJ: Cod and climate in a latitudinal cline: physiological analyses of climate effects in marine fishes. Climate research. 2008, 37: 253-270. 10.3354/cr00766.View ArticleGoogle Scholar
- Marcil J, Swain DP, Hutchings JA: Genetic and environmental components of phenotypic variation in body shape among populations of Atlantic cod (Gadus morhua). Biol J Linn Soc. 2006, 88: 351-365. 10.1111/j.1095-8312.2006.00656.x.View ArticleGoogle Scholar
- Nissling A, Kryvi H, Vallin L: Variation in egg buoyancy of Baltic cod Gadus morhua and its implications for egg survival in prevailing conditions in the Baltic Sea. Mar Ecol-Prog Ser. 1994, 110: 67-74. 10.3354/meps110067.View ArticleGoogle Scholar
- Schulte PM, Glémet HC, Fiebig AA, Powers DA: Adaptive variation in lactate dehydrogenase-B gene expression: Role of a stress-responsive regulatory element. P Natl Acad Sci USA. 2000, 97: 6597-6602. 10.1073/pnas.97.12.6597.View ArticleGoogle Scholar
- Hemmer-Hansen J, Nielsen EE, Frydenberg J, Loeschcke V: Adaptive divergence in a high gene flow environment: Hsc70 variation in the European flounder (Platichthys flesus L.). Heredity. 2007, 99: 592-600. 10.1038/sj.hdy.6801055.View ArticlePubMedGoogle Scholar
- Colosimo PF, Peichel CL, Nereng K, Blackman BK, Shapiro MD, Schluter D, Kingsley DM: The genetic architecture of parallel armor plate reduction in threespine sticklebacks. PloS Biol. 2004, 2: 635-641. 10.1371/journal.pbio.0020109.View ArticleGoogle Scholar
- Rogers SM, Bernatchez L: The genetic architecture of ecological speciation and the association with signatures of selection in natural lake whitefish (Coregonas sp Salmonidae) species pairs. Mol Biol Evol. 2007, 24: 1423-1438. 10.1093/molbev/msm066.View ArticlePubMedGoogle Scholar
- Árnason E, Hernandez UB, Kristinsson K: Intense habitat-specific fisheries-induced selection at the molecular Pan I locus predicts imminent collapse of a major cod fishery. PLoS ONE. 2009, 4: e5529-10.1371/journal.pone.0005529.PubMed CentralView ArticlePubMedGoogle Scholar
- O'Brien CM, Fox CJ, Planque B, Casey J: Climate variability and North Sea cod. Nature. 2000, 404: 142-10.1038/35004654.View ArticlePubMedGoogle Scholar
- Estoup A, Largiader CR, Perrot E, Chourrout D: Rapid one-tube DNA extraction for reliable PCR detection of fish polymorphic markers and transgenes. Mol Mar Biol Biotech. 1996, 5: 295-298.Google Scholar
- Guo S, Thompson E: Performing the exact test of Hardy-Weinberg proportion for multiple alleles. Biometrics. 1992, 48: 361-372. 10.2307/2532296.View ArticlePubMedGoogle Scholar
- Excoffier L, Laval G, Schneider S: Arlequin ver. 3.0: An integrated software package for population genetics data analysis. Evol Bioinform Online. 2005, 1: 47-50.PubMed CentralGoogle Scholar
- Weir BS, Cockerham CC: Estimating F-statistics for the analysis of population structure. Evolution. 1984, 38: 1358-1370. 10.2307/2408641.View ArticleGoogle Scholar
- Raymond M, Rousset F: GENEPOP (version 1.2): population genetics software for exact tests and ecumenicism. J Hered. 1995, 86: 248-249.Google Scholar
- Young FW: ViSta: the visual statistics system. 1996, Research Memorandum 94-1(b). L.L. Thursone Psychometric Laboratory, University of North Carolina, Chaper Hill, NC, 2Google Scholar
- Beaumont MA, Nichols RA: Evaluating loci for use in the genetic analysis of population structure. Proc R Soc Lond B Bio. 1996, 263: 1619-1626. 10.1098/rspb.1996.0237.View ArticleGoogle Scholar
- Beaumont MA, Balding DJ: Identifying adaptive genetic divergence among populations from genome scans. Mol Ecol. 2004, 13: 969-980. 10.1111/j.1365-294X.2004.02125.x.View ArticlePubMedGoogle Scholar
- Robertson A: Letters to the editor: remarks on the Lewontin-Krakauer test. Genetics. 1975, 80: 396-PubMed CentralPubMedGoogle Scholar
- Excoffier L, Hofer T, Foll M: Detecting loci under selection in a hierarchically structured population. Heredity. 2009, 103: 285-298. 10.1038/hdy.2009.74.View ArticlePubMedGoogle Scholar
- Joost S, Kalbermatten M, Bonin A: Spatial analysis method (SAM): a software tool combining molecular and environmental data to identify candidate loci for selection. Mol Ecol Res. 2008, 8: 957-960. 10.1111/j.1755-0998.2008.02162.x.View ArticleGoogle Scholar
- Case RAJ, Hutchinson WF, Hauser L, Van Oosterhout C, Carvalho GR: Macro- and micro-geographic variation in pantophysin (Pan I) allele frequencies in NE Atlantic cod Gadus morhua. Mar Ecol-Prog Ser. 2005, 301: 267-278. 10.3354/meps301267.View ArticleGoogle Scholar
- Ashburner M, Ball CA, Blake JA, Botstein D, Butler H, Cherry JM, Davis AP, Dolinski K, Dwight SS, Eppig JT, Harris MA, Hill DP, Issel-Tarver L, Kasarskis A, Lewis S, Matese JC, Richardson JE, Ringwald M, Rubin GM, Sherlock G: Gene Ontology: tool for the unification of biology. Nat Genet. 2000, 25: 25-29. 10.1038/75556.PubMed CentralView ArticlePubMedGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.