Mitochondrial DNA sequence variation and genetic stock structure of Atlantic cod (Gadus morhua) from bay and offshore locations on the Newfoundland continental shelf

Bay cod, Atlantic cod (Gadus morhua) that over-winter in the deep-water bays of northeastern Newfoundland, have historically been regarded as distinct in migration and spawning behaviour from offshore (Grand Bank) cod stocks. To investigate their genetic relationships, we determined the DNA sequence of a 307-base-pair portion of the mitochondrial cytochrome b gene for 236 adult cod taken from the waters off northeastern Newfoundland, including fish found over-wintering and spawning in Trinity Bay. Although 17 genotypes were found, a single common genotype occurs at a frequency of greater than 80% in all samples, and no alternative genotype occurs at a frequency of greater than 3%. Genotype proportions did not differ significantly among samples. Measures of genetic subdivision among sampling locations are nil. Cod over-wintering in Trinity Bay are not genetically distinct from offshore cod. In combination with tagging and physiological studies, these data suggest that there is sufficient movement of cod between bay and offshore locations to prevent the development or maintenance of independent inshore stocks. Adult cod that over-winter in Trinity Bay appear to represent an assemblage of temporarily nonmigratory fish that have become physiologically acclimated to cold-water inshore environments. The pattern of genetic variation in northern cod suggests a recent population structure characterized by extensive movement of contemporary individuals superimposed on an older structure characterized by a bottleneck in the population size of cod in the north-western Atlantic.     

A nuclear DNA basis for shelf- and bank-scale population structure in northwest Atlantic cod (Gadus morhua): Labrador to Georges Bank

Variation at five microsatellite DNA loci scored in ≈ 1300 individuals provided evidence of genetic structure among 14 cod ( Gadus morhua ) populations spanning the range of the species in the northwest (NW) Atlantic. Using D _A and D _SW measures of genetic distance, as well as F _ST and R _ST measures of population structure, differences are revealed among populations at continental shelf scales (NE Newfoundland Shelf, Grand Banks, Flemish Cap, Scotian Shelf, Georges Bank) where regions are separated by submarine saddles, channels and trenches. However, we also provide evidence of genetic structure at spawning-bank scales consistent with variation in oceanographic features and in the spatiotemporal distribution of spawning, each of which may represent barriers to gene flow among geographically contiguous populations inhabiting a highly advective environment. The differences described are consistent with postdispersal spawning fidelity to natal areas, a behaviour that may be facilitated by topographically induced gyre-like circulations that can act as retention mechanisms. Significant degrees of substructure among neighbouring and contiguous cod populations may be most easily explained by the associated oceanographic features and processes that conceivably form the template for the evolution of the structure. We suggest that bathymetric and hydrodynamic structure represents a rational starting point for developing hypotheses to examine the processes that lead to the genetic structuring of marine fish species.     

Ancient chromosomal rearrangement associated with local adaptation of a postglacially colonized population of Atlantic Cod in the northwest Atlantic

Intraspecific diversity is central to the management and conservation of exploited species, yet knowledge of how this diversity is distributed and maintained in the genome of many marine species is lacking. Recent advances in genomic analyses allow for genome‐wide surveys of intraspecific diversity and offer new opportunities for exploring genomic patterns of divergence. Here, we analysed genome‐wide polymorphisms to measure genetic differentiation between an offshore migratory and a nonmigratory population and to define conservation units of Atlantic Cod (Gadus morhua) in coastal Labrador. A total of 141 individuals, collected from offshore sites and from a coastal site within Gilbert Bay, Labrador, were genotyped using an ~11k single nucleotide polymorphism array. Analyses of population structure revealed strong genetic differentiation between migratory offshore cod and nonmigratory Gilbert Bay cod. Genetic differentiation was elevated for loci within a chromosomal rearrangement found on linkage group 1 (LG1) that coincides with a previously found double inversion associated with migratory and nonmigratory ecotype divergence of cod in the northeast Atlantic. This inverted region includes several genes potentially associated with adaptation to differences in salinity and temperature, as well as influencing migratory behaviour. Our work provides evidence that a chromosomal rearrangement on LG1 is associated with parallel patterns of divergence between migratory and nonmigratory ecotypes on both sides of the Atlantic Ocean.     

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.     

Genetic Structure of Capelin (Mallotus villosus) in the Northwest Atlantic Ocean

Capelin (Mallotus villosus) is a commercially exploited, key forage-fish species found in the boreal waters of the North Pacific and North Atlantic Oceans. We examined the population structure of capelin throughout their range in the Canadian northwest Atlantic Ocean using genetic-based methods. Capelin collected at ten beach and five demersal spawning locations over the period 2002 through 2008 (N = 3,433 fish) were genotyped using six polymorphic microsatellite loci. Temporally distinct samples were identified at three beach spawning locations: Chance Cove, Little Lawn and Straitsview, Newfoundland. Four capelin stocks are assumed for fisheries management in the northwest Atlantic Ocean based on meristics, morphometrics, tag returns, and seasonal distribution patterns. Our results suggested groupings that were somewhat different than the assumed structure, and indicate at least seven genetically defined populations arising from two ancestral populations. The spatial mosaic of capelin from each of the two basal cluster groups explains much of the observed geographic variability amongst neighbouring samples. The genetic-defined populations were resolved at Jost’s D _est ≥ 0.01 and were composed of fish collected 1) in the Gulf of St. Lawrence, 2) along the south and east coasts of Newfoundland, 3) along coastal northern Newfoundland and southern Labrador, 4) along coastal northern Labrador, 5) near the Saguenay River, and at two nearshore demersal spawning sites, 6) one at Grebes Nest off Bellevue Beach on the east coast of Newfoundland, and 7) one off the coast of Labrador at Domino Run. Moreover, the offshore demersal spawners on the Scotian Shelf and Southeast Shoal appeared to be related to the inshore demersal spawners at Grebes Nest and in Domino Run and to beach spawners from the Gulf of St. Lawrence.     

Microsatellite data reveals weak population substructuring in Copadichromis sp.'virginalis kajose', a demersal cichlid from Lake Malawi, Africa

Small but significant differences were found in allele frequencies among five populations (overall F _ST estimate (θ)=0·004, P=0·006; overall R _ST estimate (RHO)=0·019, P <0·00001) of the demersal cichlid Copadichromis sp.'virginalis kajose', collected from five locations in Lake Malawi. Pairwise F _ST estimates revealed significant differences between the most southerly population (Cape Maclear), and the three most northerly populations (Mbamba Bay, Metangula and Chilola). Pairwise R _ST estimates also revealed significant differences between some populations, but no geographical pattern was discernible. There was no evidence of isolation by distance using either the shortest straight-line distance between samples, or the distance around the shoreline following a 50 m depth contour. F _ST estimates were considerably lower than found in previous studies on the mbuna (rock-dwelling species), but higher than those found in a study of three pelagic cichlid species from Lake Malawi. Substructuring in C. sp.'virginalis kajose' appears to be on a similar scale to the Atlantic cod.     

Effect of environmental change on parasites of Atlantic cod (Gadus morhua) as bioindicators of populations in the north-western Atlantic Ocean

Studies conducted in the 1980s suggested that parasites were useful in distinguishing adjacent populations of Atlantic cod Gadus morhua in the north-western Atlantic. During the mid-1980s, cod populations began to decline, culminating in closure of the fishery in 1992. A marked decrease of cod off the coast of Labrador was coincident with climatic shifts, including a decline of bottom sea temperature, and virtual absence of capelin Mallotus villosus, its main food prey. These changes in some components of the ecosystem were independent of anthropogenic influences. An initial survey of cod parasites in 2000 from coastal Labrador revealed a dramatic decline in abundance of metazoans and also a change in dominance of species of digenetic trematodes compared to 1980-1981. Comparison of metazoan parasites in three other populations of cod inhabiting the Grand (northern and southern) and St. Pierre Banks, where there was no apparent evidence of environmental changes but only commercial over-exploitation, revealed no major alterations in abundance or dominance between 1980-1981 and 2003. The results of this current study suggest that environmental changes in the food web affected the abundance of metazoan parasites only off coastal Labrador but did not influence the use of the parasites to distinguish between cod populations inhabiting the Grand or St. Pierre Banks.     

Regional differentiation of North American Atlantic salmon at allozyme loci

Allozyme variation was characterised by starch gel electrophoresis at 23 enzyme coding loci and one regulatory locus in Atlantic salmon from 53 rivers in Eastern Canada, encompassing the majority of the species' North American range. Variation among rivers was highly heterogeneous and eight of the 15 polymorphisms showed regionally restricted distributions. Nearest neighbour joining (NJ) analysis and multi‐dimensional scaling suggest six distinct regional groups; Labrador/Ungava, Gulf of Saint Lawrence, Newfoundland (excluding Gulf rivers), the Atlantic shore/Southern Uplands of Nova Scotia, the inner Bay of Fundy, and the outer Bay of Fundy. Approximately 25% of observed genetic variation was distributed among these regions with a weak though significant overall correlation of genetic and geographic distance (Mantel Test, r  = 0·255, P  = 0·005). Collectively, the rivers showed consistent divergence from European populations with strong bootstrap support for the two clusters across loci in the NJ analysis. Mean heterozygosity was 0·061 for both continental groups, but the European population showed more than twice the variation among populations. F _ST values were 0·076 and 0·176 for North America and Europe, respectively, with an overall F _ST of 0·330.     

Microsatellite loci reveal highly significant genetic differentiation among Atlantic salmon (Salmo salar L.) stocks from the east coast of Canada

Allele frequency data from eight microsatellite loci provide evidence of highly significant genetic differentiation among stocks of Atlantic salmon Salmo salar L. from the Bay of Fundy, eastern and north-western Nova Scotia and Newfoundland. Estimates of genetic structure ( R _ST and θ) were significant both among all samples taken from the different geographical locations and among samples from geographical regions for which more than one stock was sampled. Samples from the Bay of Fundy taken from stocks which are phenotypically and behaviourally diverse showed particularly high levels of genetic structure. Rogers', allele sharing and (δμ)² distances also revealed significant differences among stock samples and were significantly correlated [Rogers' and (δμ)²] with sea distance between rivers. Results suggest that stocks of Atlantic salmon in eastern Canada are highly diverse genetically and that this should be an important consideration in any management programme for stocks in the area.     

POPULATION STRUCTURE IN AN INSHORE CETACEAN REVEALED BY MICROSATELLITE AND mtDNA ANALYSIS: BOTTLENOSE DOLPHINS (TURSIOPS SP.) IN SHARK BAY, WESTERN AUSTRALIA

We examined population substructure of bottlenose dolphins ( Tursiops sp). in Shark Bay, Western Australia, using 10 highly polymorphic microsatellite loci, and mitochondrial DNA (mtDNA). For microsatellite analysis, 302 different animals were sampled from seven localities throughout the bay. Analysis of genetic differentiation between sampling localities showed a significant correlation between the number of migrants ( Nm ) calculated from F _ST, R _ST and private alleles, and distance between localities–a pattern of isolation-by-distance. For mtDNA, 220 individuals from all seven localities were sequenced for a 351 base pair fragment of the control region, resulting in eight haplotypes, with two distinct clusters of haplotypes. Values of F _ST and (φ)_ST for mtDNA yielded statistically significant differences, mostly between localities that were not adjacent to each other, suggesting female gene flow over a scale larger than the sampled localities. We also observed a significant correlation between the number of female migrants calculated from F _ST and φ_ST and the distance of sampling localities. Our results indicate that dispersal in female dolphins in Shark Bay is more restricted than that of males.     

Mitochondrial and nuclear DNA analyses reveal fine scale geographic structure in bottlenose dolphins (<Emphasis Type="Italic">Tursiops truncatus</Emphasis>) in the Gulf of Mexico

There is a need for biological information to support current stock designations of bottlenose dolphins (Tursiops truncatus) in the Gulf of Mexico. The existence of many inshore, resident “communities” raises questions as to the relationship these dolphins may hold with dolphins inhabiting neighboring inshore and coastal areas. In this study, population subdivision was examined among four resident, inshore bottlenose dolphin stocks (Sarasota Bay, FL, Tampa Bay, FL, Charlotte Harbor, FL and Matagorda Bay, TX) and one coastal stock (1–12 km offshore) in the Gulf of Mexico. Evidence of significant population structure among all areas was found on the basis of both mitochondrial DNA (mtDNA) control region sequence data and nine nuclear microsatellite loci. Estimates of relatedness showed no population contained a significantly high number of related individuals, while separate AMOVAs for males and females indicated that both sexes exhibit a significant level of site philopatry. Results presented here provide the first genetic evidence of population subdivision between the coastal Gulf of Mexico and adjacent inshore areas along the central west coast of Florida. Such strong genetic subdivision is surprising given the short geographical distance between many of these areas and the lack of obvious geographic barriers to prevent gene flow. These findings support the current, separate identification of stocks for bottlenose dolphins inhabiting the eastern coastal and inshore areas of the Gulf of Mexico.     

Q ST MEETS THE G MATRIX: THE DIMENSIONALITY OF ADAPTIVE DIVERGENCE IN MULTIPLE CORRELATED QUANTITATIVE TRAITS

The Q _ST– F _ST comparison has become an increasingly common method for inferring adaptive quantitative trait divergence among populations. For cases in which there is divergence in multiple traits, most studies have applied the method by performing multiple univariate Q _ST– F _ST comparisons. However, because traits are often genetically correlated, such univariate analyses are likely to paint a simplified picture of adaptive divergence. Here we show how the multivariate analogue of Q _ST, F_STq, which accounts for genetic correlations among traits, can be used to supply a more detailed picture of multitrait divergence. We apply the method to naturally occurring genetic variation for a suite of sexually selected display traits in Drosophila serrata . The analyses suggest the operation of divergent multivariate selection that has influenced multiple independent axes of genetic variance in a sex-specific manner. Finally, we show how a comparison of the components of F_STq, the average within and among population genetic variance–covariance matrices, G_W and G_B, can be used as an additional test of the null expectation of neutral divergence, and allows for an investigation of whether natural populations have diverged along major or minor axes of genetic variance.     

Comparison of quantitative and molecular genetic variation of native vs. invasive populations of purple loosestrife (Lythrum salicaria L., Lythraceae)

Study of adaptive evolutionary changes in populations of invasive species can be advanced through the joint application of quantitative and population genetic methods. Using purple loosestrife as a model system, we investigated the relative roles of natural selection, genetic drift and gene flow in the invasive process by contrasting phenotypical and neutral genetic differentiation among native European and invasive North American populations ( Q _ST −  F _ST analysis). Our results indicate that invasive and native populations harbour comparable levels of amplified fragment length polymorphism variation, a pattern consistent with multiple independent introductions from a diverse European gene pool. However, it was observed that the genetic variation reduced during subsequent invasion, perhaps by founder effects and genetic drift. Comparison of genetically based quantitative trait differentiation ( Q _ST) with its expectation under neutrality ( F _ST) revealed no evidence of disruptive selection ( Q _ST >  F _ST) or stabilizing selection ( Q _ST <  F _ST). One exception was found for only one trait (the number of stems) showing significant sign of stabilizing selection across all populations. This suggests that there are difficulties in distinguishing the effects of nonadaptive population processes and natural selection. Multiple introductions of purple loosestrife may have created a genetic mixture from diverse source populations and increased population genetic diversity, but its link to the adaptive differentiation of invasive North American populations needs further research.     

Assessing changes in numbers and distribution of large whale entanglements in Newfoundland and Labrador,Canada1

Entanglements of large whales in commercial fisheries in Newfoundland and Labrador, Canada, have been consistently recorded since 1979, as part of a program aimed at releasing captured animals and reducing costs to fishermen. This data set represented an opportunity to identify fisheries posing particular entanglement risks to local whale populations. Data were assessed over the periods 1979–1992 and 1993–2008, corresponding to distinct phases in fisheries distribution and intensity. Between 1979 and 2008, 1,209 large whale entanglements were recorded in Newfoundland and Labrador. These were mostly humpback whales (Megaptera novaeangliae; 80%) and minke whales (Balaenoptera acutorostrata; 15%). Dramatic declines in reported inshore whale entanglement rates were observed following the 1992 moratorium on Atlantic cod (Gadus morhua) fisheries. Recently, more entanglements have been reported further offshore, largely due to expansion of fisheries targeting snow crab (Chionoecetes opilio). For all whale species, entanglement rates and associated mortality rates varied considerably in different fishing gear. Fractions of humpback and minke whales found dead in different fishing gear differed substantially, with minke whales far more likely to be found dead than humpback whales.     

Seasonal growth of Atlantic cod: effects of temperature, feeding and reproduction

Growth of 2659 Atlantic cod Gadus morhua aged 4 to 9 years examined in Placentia Bay, Newfoundland, peaked in most cases in June and was at a minimum in October or November. Water temperature, partial fullness index ( I _P) and gonado‐somatic index ( I _G) explained between 31 and 52% of the monthly variability in growth. Temperature and I _P of capelin Mallotus villosus had significant effects on growth of all age groups and explained most of the variance for ages 6–8 and 4–5 years, respectively. The I _P of large invertebrates (ages 4 to 7 years), sandlance ( Ammodytes sp. age 6 years) and demersal fishes (age 9 years) had age‐specific effects in the model. Overall, amphipods, decapods and echinoderms dominated the Atlantic cod diet in most seasons, but fish consumption by Atlantic cod was high in June and July, particularly on capelin. The rapid increase in somatic mass during June and July occurred despite cold water temperatures ( < 3° C at 50 m) and moderate to high gonado‐somatic index. The findings of this study suggest that when food was not a limiting factor, growth tended to increase even when Atlantic cod occupied colder waters, but when food was limiting, the opposite may have occured.     

Population genetics ofAnisakis simplex (Nematoda: Ascaridoidea) in Atlantic salmon (Salmo salar) and their use as biological indicators of host stocks

Synopsis The extent to which data concerning the population genetics of a parasite — namely larvalAnisakis simplex, may be used as an auxilliary source of information on various aspects of the marine migrations of Atlantic salmon (Salmo salar) is examined.Frequencies of the various acid phosphatase alleles could not be used in distinguishing between or estimating the proportions of North American and European salmon occurring off west Greenland. The genetic structure of larvae from salmon taken in Scotland was different from that taken in Ireland and England confirming an earlier conclusion that the feeding habits of the fish from Scotland were significantly different. In the western Atlantic, salmon from northern Newfoundland and also salmon from the Bay of Fundy had apparently ingested different populations (or different proportions of the same nematode population) than salmon from Miramachi and Chaleur Bay areas. Samples from northern Newfoundland and Labrador were also genetically different. These findings are discussed in relation to the available literature concerning salmon migrations, explanations are put forward and areas of future research suggested.     

Microsatellite analysis of genetic population structure in an endangered salmonid: the coastal cutthroat trout (Oncorhynchus clarki clarki)

The genetic population structure of coastal cutthroat trout ( Oncorhynchus clarki clarki ) in Washington state was investigated by analysis of variation in allele frequencies at six highly polymorphic microsatellite loci for 13 anadromous populations, along with one outgroup population from the Yellowstone subspecies ( O. clarki bouvieri) (mean heterozygosity = 67%; average number of alleles per locus = 24). Tests for genetic differentiation revealed highly significant differences in genotypic frequencies for pairwise comparisons between all populations within geographical regions and overall population subdivision was substantial ( F _ST = 0.121, R _ST = 0.093), with 44.6% and 55.4% of the among-population diversity being attributable to differences between streams ( F _SR = 0.054) and between regions ( F _RT = 0.067), respectively. Analysis of genetic distances and geographical distances did not support a simple model of isolation by distance for these populations. With one exception, neighbour-joining dendrograms from the Cavalli-Sforza and Edwards' chord distances and maximum likelihood algorithms clustered populations by physiogeographic region, although overall bootstrap support was relatively low (53%). Our results suggest that coastal cutthroat trout populations are ultimately structured genetically at the level of individual streams. It appears that the dynamic balance between gene flow and genetic drift in the subspecies favours a high degree of genetic differentiation and population subdivision with the simultaneous maintenance of high heterozygosity levels within local populations. Results are discussed in terms of coastal cutthroat trout ecology along with implications for the designation of evolutionarily significant units pursuant to the US Endangered Species Act of 1973 and analogous conservation units.     

Population genetic structure of the butterfly Melitaea didyma (Nymphalidae) along a northern distribution range border

The population genetic structure of the butterfly Melitaea didyma was studied along the northern distribution range border in Central Germany by means of allozyme electrophoresis. Individuals were sampled from a total of 21 habitat patches from four regions, and two provinces. Sampling was designed to estimate local vs. regional differentiation. High levels of variability were found, H _e= 0.14–0.21. The mean expected sample heterozygosity from one region, Mosel, was significantly lower than from the Hammelburg region, H _e= 0.17 and 0.19, respectively. Two hierarchical levels of genetic differentiation were found. Within regions individuals sampled from different patches behaved as belonging to one population with high levels of gene flow (Hammelburg F _ST= 0.015, Mosel F _ST= 0.044), though local isolation barriers did create a substructuring of these populations. The inbreeding coefficients, F _IS, were constant over all sample levels, suggesting a similar distribution of habitat patches within regions. Between regions gene flow was limited. An isolation by distance analysis indicated that the hierarchical structure, at the provincial level, may be breaking down due to isolation of regional populations. A more general observation was that the sampling design may greatly have influenced the estimation of genetic differentiation. Depending on which samples were included, overall F _ST estimates ranged from 0.059–0.090.     

Allozyme analyses of the genus Trisopterus: taxonomic status and population structure of the poor cod

Genetic analysis of the four Trisopterus (Gadidae) taxa suggests that the interrelationships of the two morphs of poor cod ( T. minutus minutus in the Atlantic and T. minutus capelanus in the Mediterranean) should be reconsidered. The Mediterranean poor cod T. m. capelanus is more closely related to bib T. luscus than to the Atlantic poor cod, so the population structure in the Atlantic and Mediterranean poor cod must be considered separately. Among 635 Atlantic individuals there was some evidence of poor cod population differentiation (allele frequency heterogeneity test P <0·0005; F _ST=0·0135, P ≤0·0005). Levels of genetic variation were similar to those reported for related gadoid species. Some differentiation was present on the Norwegian coast (samples from Trondheimsfjord) and between the Faeroe Islands (Faeroe Bank) and the adjacent European coastal location. In contrast no statistically significant population differentiation was evident in Mediterranean poor cod, but fewer samples and individuals were screened.     

Subtle population structuring within a highly vagile marine invertebrate, the veined squid Loligo forbesi, demonstrated with microsatellite DNA markers

Microsatellite DNA markers were applied for the first time in a population genetic study of a cephalopod and compared with previous estimates of genetic differentiation obtained using allozyme and mitochondrial DNA (mtDNA) markers. Levels of genetic variation detected with microsatellites were much higher than found with previous markers (mean number of alleles per locus=10.6, mean expected heterozygosity ( H _E)=0.79; allozyme H _E=0.08; mtDNA restriction fragment length polymorphism (RFLP) H _E=0.16). In agreement with previous studies, microsatellites demonstrated genetic uniformity across the population occupying the European shelf seas of the North East Atlantic, and extreme genetic differentiation of the Azores population ( R _ST/ F _ST=0.252/0.245; allozyme F _ST=0.536; mtDNA F _ST=0.789). In contrast to other markers, microsatellites detected more subtle, and significant, levels of differentiation between the populations of the North East Atlantic offshore banks (Rockall and Faroes) and the shelf population ( R _ST=0.048 and 0.057). Breakdown of extensive gene flow among these populations is indicated, with hydrographic (water depth) and hydrodynamic (isolating current regimes) factors suggested as possible barriers to migration. The demonstration of genetic subdivision in an abundant, highly mobile marine invertebrate has implications for the interpretation of dispersal and population dynamics, and consequent management, of such a commercially exploited species. Relative levels of differentiation indicated by the three different marker systems, and the use of measures of differentiation (assuming different mutation models), are discussed.