Two adjacent inversions maintain genomic differentiation between migratory and stationary ecotypes of Atlantic cod

Atlantic cod is composed of multiple migratory and stationary populations widely distributed in the North Atlantic Ocean. The Northeast Arctic cod (NEAC) population in the Barents Sea undertakes annual spawning migrations to the northern Norwegian coast. Although spawning occurs sympatrically with the stationary Norwegian coastal cod (NCC), phenotypic and genetic differences between NEAC and NCC are maintained. In this study, we resolve the enigma by revealing the mechanisms underlying these differences. Extended linkage disequilibrium (LD) and population divergence were demonstrated in a 17.4‐Mb region on linkage group 1 (LG1) based on genotypes of 494 SNPs from 192 parents of farmed families of NEAC, NCC or NEACxNCC crosses. Linkage analyses revealed two adjacent inversions within this region that repress meiotic recombination in NEACxNCC crosses. We identified a NEAC‐specific haplotype consisting of 186 SNPs that was fixed in NEAC sampled from the Barents Sea, but segregating under Hardy–Weinberg equilibrium in eight NCC stocks. Comparative genomic analyses determine the NEAC configuration of the inversions to be the derived state and date it to ~1.6–2.0 Mya. The haplotype block harbours 763 genes, including candidates regulating swim bladder pressure, haem synthesis and skeletal muscle organization conferring adaptation to long‐distance migrations and vertical movements down to large depths. Our results suggest that the migratory ecotype experiences strong directional selection for the two adjacent inversions on LG1. Despite interbreeding between NEAC and NCC, the inversions are maintaining genetic differentiation, and we hypothesize the co‐occurrence of multiple adaptive alleles forming a ‘supergene’ in the NEAC population.     

Pantophysin (Pan I) locus divergence between inshore v. offshore and northern v. southern populations of Atlantic cod in the north‐east Atlantic

More than 6000 cod Gadus morhua , sampled in coastal and offshore waters stretching from the Barents Sea down to the North Sea, were analysed for frequencies of alleles at the scnDNA pantophysin locus ( Pan I)[formerly called synaptophysin ( Syp I)]. The significant allele frequency difference between the two major stocks of cod in Norway, north‐east Arctic cod (NEAC) and Norwegian coastal cod (NCC), was upheld in all years of the investigation (1993 to 2001), and applied both to larval cod and post‐juveniles of various ages. On a north‐south axis, the appearance of a latitudinal cline of post‐juvenile (≥1 year) allele frequencies was exposed. The intermediate allele frequencies in coastal areas of northern Norway, seem to a large extent to be caused by intermingling of the two stocks, although the existence of populations of coastal cod with alternative Pan I frequencies could not be ruled out. The role of selection is yet unresolved. Depth of the sampling location seemed to have an effect on the allele frequencies and their temporal stability, while there was no indication of seasonal variation in the frequencies. Breeding structure was the most likely cause for upholding the extreme divergence in Pan I frequencies between NEAC and NCC.     

Recolonization history and large-scale dispersal in the open sea: the case study of the North Atlantic cod, Gadus morhua L.

Most studies of the genetic structure of Atlantic cod have focused on small geographical scales. In the present study, the genetic structure of cod sampled on spawning grounds in the North Atlantic was examined using eight microsatellite loci and the Pan I locus. A total of 954 cod was collected from nine different regions: the Baltic Sea, the North Sea, the Celtic Sea, the Irish Sea and Icelandic waters during spring 2002 and spring 2003, from Norwegian waters and the Faroe Islands (North and West spawning grounds) in spring 2003, and from Canadian waters in 1998. Temporal stability among spawning grounds was observed in Icelandic waters and the Celtic Sea, and no significant difference was observed between the samples from the Baltic Sea and between the samples from Faroese waters. F -statistics showed significant differences between most populations and a pattern of isolation-by-distance was described with microsatellite loci. The Pan I locus revealed the presence of two genetically distinguishable basins, the North-west Atlantic composed of the Icelandic and Canadian samples and the North-east Atlantic composed of all other samples. Permutation of allele sizes at each microsatellite locus among allelic states supported a mutational component to the genetic differentiation, indicating a historical origin of the observed variation. Estimation of the time of divergence was approximately 3000 generations, which places the origin of current genetic pattern of cod in the North Atlantic in the late Weichselian (Wisconsinian period), at last glacial maximum.  © 2008 The Linnean Society of London, Biological Journal of the Linnean Society , 2008, 94 , 315–329.     

Evidence of microsatellite hitch-hiking selection in Atlantic cod (Gadus morhua L.): implications for inferring population structure in nonmodel organisms

Microsatellites have gained wide application for elucidating population structure in nonmodel organisms. Since they are generally noncoding, neutrality is assumed but rarely tested. In Atlantic cod (Gadus morhua L.), microsatellite studies have revealed highly heterogeneous estimates of genetic differentiation among loci. In particular one locus, Gmo 132, has demonstrated elevated genetic differentiation. We investigated possible hitch-hiking selection at this and other microsatellite loci in Atlantic cod. We employed 11 loci for analysing samples from the Baltic Sea, North Sea, Barents Sea and Newfoundland covering a large part of the species' distributional range. The 'classical' Lewontin-Krakauer test for selection based on variance in estimates of F(ST) and (standardized genetic differentiation) revealed only one significant pairwise test (North Sea-Barents Sea), and the source of the elevated variance could not be ascribed exclusively to Gmo 132. In contrast, different variants of the recently developed ln Rtheta test for selective sweeps at microsatellite loci revealed a high number of significant outcomes of pair-wise tests for Gmo 132. Further, the presence of selection was indicated in at least one other locus. The results suggest that many previous estimates of genetic differentiation in cod based on microsatellites are inflated, and in some cases relationships among populations are obscured by one or more loci being the subject to hitch-hiking selection. Likewise, temporal estimates of effective population sizes in Atlantic cod may be flawed. We recommend, generally, to use a higher number of microsatellite loci to elucidate population structure in marine fishes and other nonmodel species to allow for identification of outlier loci that are subject to selection.     

Differences in Salinity Tolerance and Gene Expression Between Two Populations of Atlantic Cod (Gadus morhua) in Response to Salinity Stress

Populations of marine fish, even from contrasting habitats, generally show low genetic differentiation at neutral genetic markers. Nevertheless, there is increasing evidence for differences in gene expression among populations that may be ascribed to adaptive divergence. Studying variation in salinity tolerance and gene expression among Atlantic cod (Gadus morhua) from two populations distributed across a steep salinity gradient, we observed high mortality (45% North Sea cod and 80% Baltic Sea cod) in a reciprocal common garden setup. Quantitative RT-PCR assays for expression of hsp70 and Na/K-ATPase α genes demonstrated significant differences in gene regulation within and between populations and treatment groups despite low sample sizes. Most interesting are the significant differences observed in expression of the Na/K-ATPase α gene in gill tissue between North Sea and Baltic cod. The findings strongly suggest that Atlantic cod are adapted to local saline conditions, despite relatively low levels of neutral genetic divergence between populations.     

Comparative analysis of genetic variability of white sea cod (Gadus morhua marisalbi) at allozyme and microsatellite markers

Variability of cod spawning and feeding schools from Kandalaksha Bay of the White Sea, was examined at six allozyme and eight microsatellite loci. The degree of genetic differentiation at allozyme loci constituted θ = 0.36% [95% bootstrap interval 0.0458; 0.6743]. The differentiation estimates obtained using microsatellite markers were higher, θ = 1.33% [0.057; 3.11]. It was demonstrated that the level of genetic diversity in the White Sea cod was lower than that established for the Atlantic cod from Barents Sea using the same set of allozyme and microsatellite markers. The genetic data obtained support the opinion that the White Sea cod is a reproductively independent group formed as a result of the Holocene dispersal of Atlantic cod.     

Mixed stock analysis and the power of different classes of molecular markers in discriminating coastal and oceanic Atlantic cod (<Emphasis Type="Italic">Gadus morhua</Emphasis> L.) on the Lofoten spawning grounds,Northern Norway

Atlantic cod (Gadus morhua) encompasses many different populations or stocks, which in part are managed separately. In the northeast Atlantic cod is divided into two main management units; northeast Arctic cod and coastal cod. These two groups have traditionally been separated by otolith classification. In this study, the power of different classes of genetic markers in separating the two cod groups was investigated. The variation in thirteen genetic markers, including allozymes, haemoglobin, the scDNA locus Pantophysin (Pan I) and a number of microsatellites was investigated, and mixed stock analysis and individual assignment tests were performed on samples comprising a mixture of individuals of putative coastal and oceanic type cod. The genetic analyses showed a large genetic differentiation between outer stations and stations located closer to the mainland shore. Mixed stock analysis and individual assignment tests used for estimation of stock proportions gave results similar to those obtained by otolith classification. Guest editors: J. Davenport, G. Burnell, T. Cross, M. Emmerson, R. McAllen, R. Ramsay & E. Rogan Challenges to Marine Ecosystems     

Adaptive divergence in a high gene flow environment: Hsc70 variation in the European flounder (Platichthys flesus L.)

Little is known about local adaptations in marine fishes since population genetic surveys in these species have typically not applied genetic markers subject to selection. In this study, we used a candidate gene approach to investigate adaptive population divergence in the European flounder (Platichthys flesus L.) throughout the northeastern Atlantic. We contrasted patterns of genetic variation in a presumably neutral microsatellite baseline to patterns from a heat-shock cognate protein gene, Hsc70. Using two different neutrality tests we found that the microsatellite data set most likely represented a neutral baseline. In contrast, Hsc70 strongly deviated from neutral expectations. Importantly, when estimating standardized levels of population divergence (F(ST)'), we also found a large discrepancy in the patterns of structuring in the two data sets. Thus, samples grouped according to geographical or historical proximity with regards to microsatellites, but according to environmental similarities with regards to Hsc70. The differences between the data sets were particularly pronounced in pairwise comparisons involving populations in the western and central Baltic Sea. For instance, the genetic differentiation between geographically close Baltic Sea and North Sea populations was found to be 0.02 and 0.45 for microsatellites and Hsc70 respectively. Our results strongly suggest adaptive population divergence and indicate local adaptations at the DNA level in a background of high levels of gene flow, typically found in many marine fish species. Furthermore, this study highlights the usefulness of the candidate gene approach for demonstrating local selection in non-model organisms such as most marine fishes.     

Characterization of 155 EST-derived microsatellites from Atlantic cod (Gadus morhua) and validation for linkage mapping

Microsatellite markers for Atlantic cod (Gadus morhua) were identified from a collection of 30 630 expressed sequence tags. Primers were designed for 395 microsatellites and 155 were successfully amplified. Allele number varied from 1 to 26 (average 6.34). Average observed and expected heterozygosities were 0.50 and 0.54, respectively. A subset of 105 microsatellites tested for Mendelian segregation showed no significant distortion of segregation when correcting for multiple tests. Null alleles were detected at four loci. Significant blastx matches were found for 23 loci. These microsatellites will be used to create a linkage map to enhance genetic selection in commercial cod breeding.     

Variability of microsatellite loci of Greenland cod <Emphasis Type="Italic">Gadus ogac</Emphasis> Richardson 1836: Comparison with other species of <Emphasis Type="Italic">Gadus</Emphasis> genus (Gadidae)

Comparative analysis of variability of seven microsatellite loci—Gmo3, Gmo-G12, Gmo-G18, Gmo19, Gmo34, Gmo35 and Pgmo32—was performed for the Greenland cod Gadus ogac, Pacific cod G. macrocephalus, Atlantic cod G. morhua, and White Sea cod G. morhua marisalbi. High genetic identity was observed between the Greenland cod and Pacific cod (I = 0.9520). Pair analysis of genetic differentiation was performed on the studied microsatellite loci according to θ (analogue of F _ST). The Greenland cod differed significantly from the Pacific, Atlantic, and the White Sea cod; however, the differentiation level varied. The lowest value was observed for the pair Greenland cod-Pacific cod (0.123), and the highest levels were registered for the pairs Greenland cod-Atlantic cod (0.605) and Greenland cod-White Sea cod (0.535).     

Adaptive and Neutral Genetic Variation and Colonization History of Atlantic Salmon, Salmo salar

Synopsis I combined neutral microsatellite markers with the major histocompatibility complex (MHC) class IIB to study genetic differentiation and colonization history in Atlantic salmon, Salmo salar, in the Baltic Sea and in the north-eastern Atlantic. Baltic salmon populations have lower levels of microsatellite genetic variation, in terms of heterozygosity and allelic richness than Atlantic populations, confirming earlier findings with other genetic markers, suggesting that the Baltic Sea populations have been exposed to genetic bottlenecks, most likely at a founding event. On the other hand, the level of MHC variation was similar in the Baltic and in the north-eastern Atlantic, indicating that positive balancing selection has increased the level of MHC-variation. Both microsatellite and MHC class IIB genetic variation give strong support to the hypothesis that the Baltic salmon are of a biphyletic origin, the southern population in this study is strongly differentiated from both the northern Baltic salmon populations and from the north-eastern Atlantic populations. Salmon may have colonized the northern Baltic Sea either from the south, via the so called “N?rke strait” or from the north, via a proposed historical connection between the White Sea and the northern Baltic. At microsatellites, no significant isolation-by distance was found at either colonization route. At the MHC, populations were significantly isolated by distance when assuming that colonization occurred via the “N?rke strait”.     

Evidence for positive selection at the pantophysin (Pan I) locus in walleye pollock, Theragra chalcogramma

Nucleotide polymorphism at the pantophysin (Pan I) locus in walleye pollock, Theragra chalcogramma, was examined using DNA sequence data. Two distinct allelic lineages were detected in pollock, resulting from three amino acid replacement mutations in the first intravesicular domain of the protein. The common Pan I allelic group, comprising 94% of the samples, was less polymorphic (pi = 0.005) than the uncommon group (pi = 0.008), and nucleotide diversity in both was higher than for two allelic lineages in the related Atlantic cod, Gadus morhua. Phylogenetic analyses of Pan I sequences from these two species did not clearly resolve orthology among allelic groups, in part because of recombination that has occurred between the two pollock lineages. Conventional tests of neutrality comparing polymorphisms within and between homologous regions of the Pan I locus in walleye pollock and Atlantic cod did not detect the effects of selection. This result is likely attributed to low levels of synonymous divergence among allelic lineages and a lack of mutation-drift equilibrium inferred from nucleotide mismatch frequency distributions. However, the ratio of nonsynonymous to synonymous substitutions per site (dN/dS) exceeded unity in two intravesicular domains of the protein and the influence of positive selection at multiple codon sites was strongly inferred through the use of maximum-likelihood analyses. In addition, the frequency spectrum of linked neutral variation showed indirect effects of adaptive hitchhiking in pollock resulting from a selective sweep of the common allelic lineage. Recombination between the two allelic classes may have prevented complete loss of the older, more polymorphic lineage. The results suggest that recurrent sweeps driven by positive selection is the principle mode of evolution at the Pan I locus in gadid fishes.     

Disentangling structural genomic and behavioural barriers in a sea of connectivity

Genetic divergence among populations arises through natural selection or drift and is counteracted by connectivity and gene flow. In sympatric populations, isolating mechanisms are thus needed to limit the homogenizing effects of gene flow to allow for adaptation and speciation. Chromosomal inversions act as an important mechanism maintaining isolating barriers, yet their role in sympatric populations and divergence with gene flow is not entirely understood. Here, we revisit the question of whether inversions play a role in the divergence of connected populations of the marine fish Atlantic cod (Gadus morhua), by exploring a unique data set combining whole‐genome sequencing data and behavioural data obtained with acoustic telemetry. Within a confined fjord environment, we find three genetically differentiated Atlantic cod types belonging to the oceanic North Sea population, the western Baltic population and a local fjord‐type cod. Continuous behavioural tracking over 4 year revealed temporally stable sympatry of these types within the fjord. Despite overall weak genetic differentiation consistent with high levels of gene flow, we detected significant frequency shifts of three previously identified inversions, indicating an adaptive barrier to gene flow. In addition, behavioural data indicated that North Sea cod and individuals homozygous for the LG12 inversion had lower fitness in the fjord environment. However, North Sea and fjord‐type cod also occupy different depths, possibly contributing to prezygotic reproductive isolation and representing a behavioural barrier to gene flow. Our results provide the first insights into a complex interplay of genomic and behavioural isolating barriers in Atlantic cod and establish a new model system towards an understanding of the role of genomic structural variants in adaptation and diversification.     

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.     

Postglacial recolonization patterns and genetic relationships among whitefish (Coregonus sp.) populations in Denmark, inferred from mitochondrial DNA and microsatellite markers

The genetic relationships among morphologically and geographically divergent populations of whitefish (genus: Coregonus ) from Denmark and the Baltic Sea region were studied by analysis of microsatellites and polymerase chain reaction–restriction fragment length polymorphism (PCR–RFLP) analysis of mitochondrial DNA (mtDNA) segments. The endangered North Sea houting (classified as C. oxyrhynchus ) differs morphologically and physiologically from other Danish whitefish ( C. lavaretus ). However, limited divergence of North Sea houting was observed both at the level of mtDNA and microsatellites. The implications of these results for the conservation status of North Sea houting are discussed in the light of current definitions of evolutionary significant units. Both mtDNA and microsatellite data indicated that postglacial recolonization by C. lavaretus in Denmark was less likely to have taken place from the Baltic Sea. Instead, the data suggested a recent common origin of all Danish whitefish populations, including North Sea houting, probably by recolonization via the postglacial Elbe River system. Estimates of genetic differentiation among populations based on mtDNA and microsatellites were qualitatively different. In addition, for both classes of markers analyses of genetic differentiation yielded different results, depending on whether molecular distances between alleles or haplotypes were included.     

Genetic variation within and among North Atlantic and Baltic populations of the benthic alga Phycodrys rubens (Rhodophyta)

Genetic variation was examined within and among North Atlantic, North Sea and Baltic populations of the benthic red alga Phycodrys rubens using allozymes and random amplified polymorphic DNA (RAPD) markers. On western and eastern North Atlantic coasts distinct allozyme types were found, with the exception of western Newfoundland where East and West Atlantic types co-occur. Along the European coasts, two genetic groups were distinguished by fixed allelic differences: an outer oceanic group and a North Sea/Baltic group. The two genetic types co-occur in the Skagerrak and Kattegat region. Reproductive isolation between the two types is suggested by the lack of hybrids in the overlap zones, and they may therefore represent sibling species. Unexpectedly, an analysis of RAPD variation was unable to recover the two cryptic species identified using allozymes. Within-population RAPD variation was similar to or greater than between-population variation. The lack of structure in the RAPD data cannot be attributed solely to technical artefacts of the method but appears to reflect real biological variability. Within-population genomic polymorphisms caused by frequent mutational events are discussed, as are high amounts of genetic drift and possible disruptive selection brought about by stressed habitats. Finally, Baltic and extra-Baltic salinity ecotypes are known to exist in P. rubens. However, no correlation between ecotypic variation and allozyme groups was detected.     

Allozyme and microsatellite loci provide discordant estimates of population differentiation in the endangered dusky grouper (Epinephelus marginatus) within the Mediterranean Sea

The dusky grouper, Epinephelus marginatus, inhabits coastal reefs in the Mediterranean Sea and Atlantic Ocean. A decline in the abundance of this long-lived protogynous hermaphrodite has led to its listing as an endangered species in the Mediterranean, and heightened management concerns regarding its genetic variability and population substructure. To address these concerns, we analysed genetic variation at seven microsatellite and 28 allozyme loci in dusky groupers sampled from seven areas (for microsatellites) and three areas (for allozymes) in the west-central Mediterranean. Levels of genetic variability were higher for microsatellites than for allozymes (mean H(E) = 0.78 and 0.07, respectively), but similar to those observed in other marine fishes with comparable markers. Both microsatellites and allozymes revealed significant genetic differentiation among all areas analysed with each class of marker, but the magnitude of differentiation revealed by allozymes over three locales (F(ST) = 0.214) was greater than that detected with microsatellites over seven areas, or over the three areas shared with the allozyme analysis (F(ST) = 0.018 and approximately 0, respectively). A large proportion of the allozyme differentiation was due to a single locus (ADA*) possibly influenced by selection, but allozyme differentiation over the three areas was still highly significant (F(ST) = 0.06, P < 0.0001), and the 95% confidence intervals for allozyme and microsatellite F(ST) did not overlap when this locus was excluded. There was no evidence of isolation by distance with either class of markers. Our results lead us to conclude that dusky groupers are not panmictic in the Mediterranean Sea and suggest that they should be managed on a local basis. However, more work is needed to elucidate genetic relationships among populations.     

Female philopatry in coastal basins and male dispersion across the North Atlantic in a highly mobile marine species, the sperm whale (Physeter macrocephalus)

The mechanisms that determine population structure in highly mobile marine species are poorly understood, but useful towards understanding the evolution of diversity, and essential for effective conservation and management. In this study, we compare putative sperm whale populations located in the Gulf of Mexico, western North Atlantic, Mediterranean Sea and North Sea using mtDNA control region sequence data and 16 polymorphic microsatellite loci. The Gulf of Mexico, western North Atlantic and North Sea populations each possessed similar low levels of haplotype and nucleotide diversity at the mtDNA locus, while the Mediterranean Sea population showed no detectable mtDNA diversity. Mitochondrial DNA results showed significant differentiation between all populations, while microsatellites showed significant differentiation only for comparisons with the Mediterranean Sea, and at a much lower level than seen for mtDNA. Samples from either side of the North Atlantic in coastal waters showed no differentiation for mtDNA, while North Atlantic samples from just outside the Gulf of Mexico (the western North Atlantic sample) were highly differentiated from samples within the Gulf at this locus. Our analyses indicate a previously unknown fidelity of females to coastal basins either side of the North Atlantic, and suggest the movement of males among these populations for breeding.     

Structural and functional connectivity of marine fishes within a semi-enclosed Newfoundland fjord

The interplay between structural connectivity ( i.e. habitat continuity) and functional connectivity ( i.e. dispersal probability) in marine fishes was examined in a coastal fjord (Holyrood Pond, Newfoundland, Canada) that is completely isolated from the North Atlantic Ocean for most of the year. Genetic differentiation was described in three species (rainbow smelt Osmerus mordax , white hake Urophycis tenuis and Atlantic cod Gadus morhua ) with contrasting life histories using seven to 10 microsatellite loci and a protein-coding locus, Pan I ( G. morhua ). Analysis of microsatellite differentiation indicated clear genetic differences between the fjord and coastal regions; however, the magnitude of difference was no more elevated than adjacent bays and was not enhanced by the fjord's isolation. Osmerus mordax was characterized by the highest structure overall with moderate differentiation between the fjord and St Mary's Bay ( F _ST c. 0·047). In contrast, U. tenuis and G. morhua displayed weak differentiation ( F _ST < 0·01). Nonetheless, these populations did demonstrate high rates (< 75%) of Bayesian self-assignment. Furthermore, elevated differentiation was observed at the Pan I locus in G. morhua between the fjord and other coastal locations. Interestingly, locus-specific genetic differentiation and expected heterozygosity were negatively associated in O. mordax , in contrast to the positive associations observed in U. tenuis and G. morhua . Gene flow in these species is apparently unencumbered by limited structural connectivity, yet the observed differentiation suggests that population structuring exists over small scales despite high dispersal potential.     

Twenty-three novel microsatellite markers developed from Atlantic cod Gadus morhua L. expressed sequence tags

Twenty-three novel polymorphic microsatellite markers were developed from c. 2000 expressed sequence tags of Atlantic cod Gadus morhua L. Gene identity was determined at 12 loci, confirming the associated microsatellites as type I markers. These microsatellite markers provide useful tools for studies of population genetics and reproductive ecology and for constructing linkage maps of G. morhua .