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.     

Life-stage-dependent supergene haplotype frequencies and metapopulation neutral genetic patterns of Atlantic cod,Gadus morhua,from Canada's Northern cod stock region and adjacent areas

Among highly migratory fish species, nursery areas occupied by juveniles often differ from adult habitats. To better understand the spatial dynamics of Canada's Northern cod stock, juveniles caught off the east coast of Newfoundland and Labrador were compared to adults from the same region as well as individuals from other areas in Atlantic Canada using double-digest restriction site-associated DNA sequencing–derived single nucleotide polymorphisms. A reduced proportion of homozygotes with a chromosomal inversion located in linkage group 1 (LG1) was detected between juvenile and adult samples in the Northern cod stock region, potentially indicating age-dependent habitat use or ontogenetic selection for attributes associated with the many genes located in LG1. No selectively neutral genetic differences were found between samples from the Northern cod stock; nevertheless, significant differences were found between some of these samples and cod collected from St. Pierre Bank, Bay of Fundy, Browns Bank and the southern Scotian Shelf. Clustering analysis of variants at neutral loci provided evidence for three major genetic units: (a) the Newfoundland Atlantic Coast, (b) eastern and southern Gulf of St. Lawrence and Burgeo Bank and (c) the Bay of Fundy, Browns Bank and southern Scotian Shelf. Both adaptive and neutral population structure within the Northern cod stock should be considered by managers to promote demographic rebuilding of the stock.     

Population genetic substructure in blue whiting based on allozyme data

The genetic population structure of the pelagic gadoid blue whiting Micromesistius poutassou throughout its east Atlantic distribution range was explored using polymorphisms at the tissue enzyme loci IDHP-2 * and PGM-1 *. The study included 5025 individuals from 65 locations in the east Atlantic and the Mediterranean. Significant geographic heterogeneity in allele frequencies was demonstrated at both loci ( IDHP-2 *: P=0·030, PGM-1 *: P=0·005). The degree of genetic differentiation ( F _ST= 1·1%) in blue whiting appears to be at the same level as for the demersal gadoids cod and haddock. Several separate reproductive units were indicated at the fringes of the distribution range, i.e. in the Barents Sea, and in one Norwegian fjord (Romsdalsfjord).     

Bay-scale population structure in coastal Atlantic cod in Labrador and Newfoundland, Canada

Polymorphisms at five microsatellite DNA loci provide evidence that Atlantic cod Gadus morhua inhabiting Gilbert Bay, Labrador are genetically distinguishable from offshore cod on the north-east Newfoundland shelf and from inshore cod in Trinity Bay, Newfoundland. Antifreeze activity in the blood suggests that Gilbert Bay cod overwinter within the Bay. Gilbert Bay cod are also smaller (weight and length) for their age and consequently less fecund for their age, than cod elsewhere within the northern cod complex. The productivity and recruitment potential of coastal cod off Labrador may thus be much lower than that of offshore northern cod or of inshore cod farther south, implying that a more conservative management strategy may be required for cod from coastal Labrador than traditionally practised for northern cod inhabiting less harsh environments. Relatively high F _ST and R _ST measures of population structure suggest that important barriers to gene flow exist among five components that include two inshore (Gilbert and Trinity Bay) and three offshore cod aggregations on the north-east Newfoundland Shelf and the Grand Bank. D _A and D _SW estimates of genetic distance that involve Gilbert Bay cod are approximately three- and 10–fold larger, respectively, than estimates not involving Gilbert Bay cod. The differences between inshore cod from Gilbert Bay and Trinity Bay raise the possibility that other genetically distinguishable coastal populations may exist, or may have existed prior to the northern cod fishery collapse. Harvesting strategies for northern cod should recognize the existence of genetic diversity between inshore and offshore components as well as among coastal components.     

Microsatellite analysis indicates an absence of population structure among Hippoglossus hippoglossus in the north‐west Atlantic

To explore the population structure of Atlantic halibut Hippoglossus hippoglossus , 160 fish from four locations in the north‐west Atlantic (Bay of Fundy, Scotian Shelf, Gulf of St Lawrence and Iceland) were examined for evidence of population structure using 18 microsatellite markers. Pair‐wise F _ST and a model‐based cluster analysis revealed no significant differentiation between samples, although uncertainties surrounding Atlantic halibut reproductive behaviour made it difficult to ascertain that only a single breeding population had been sampled at each location     

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.     

Evidence for independent origin of two spring‐spawning ciscoes (Salmoniformes: Coregonidae) in Germany

Combined analyses of mitochondrial DNA (mtDNA) and microsatellite loci were performed to assess the genetic differentiation of two spring‐spawning ciscoes from each other and from sympatric Coregonus albula in two German lakes. Polymorphism was screened at six microsatellite loci and mtDNA for a total of 247 and 94 ciscoes, respectively. Microsatellite data showed a weak differentiation between spring‐spawning Coregonus fontanae and sympatric C. albula in Lake Stechlin ( F _ST = 0–0·008), whereas a significant differentiation was observed between spring‐spawning Coregonus lucinensis and sympatric C. albula in Lake Breiter Luzin ( F _ST = 0·013–0·039). A more pronounced genetic difference was observed between both spring‐spawning species ( F _ST = 0·05–0·128). Shared mtDNA haplotypes among sympatric species within both Lake Stechlin and Lake Breiter Luzin were observed, whereas no haplotype was shared between C. fontanae and C. lucinensis . These results suggest an independent origin for spring‐spawning ciscoes in each lake. Evidence is also provided for mtDNA introgression of Coregonus sardinella into C. lucinensis and C. albula in Lake Breiter Luzin. Postglacially, this species or at least a population which showed mtDNA introgression has colonized the Baltic Sea basin up to the glacial margin that was located between Lakes Stechlin and Breiter Luzin.     

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.     

Long-term stability and effective population size in North Sea and Baltic Sea cod (Gadus morhua)

DNA from archived otoliths was used to explore the temporal stability of the genetic composition of two cod populations, the Moray Firth (North Sea) sampled in 1965 and 2002, and the Bornholm Basin (Baltic Sea) sampled in 1928 and 1997. We found no significant changes in the allele frequencies for the Moray Firth population, while subtle but significant genetic changes over time were detected for the Bornholm Basin population. Estimates of the effective population size ( N _e ) generally exceeded 500 for both populations when employing a number of varieties of the temporal genetic method. However, confidence intervals were very wide and N _e 's most likely range in the thousands. There was no apparent loss of genetic variability and no evidence of a genetic bottleneck for either of the populations. Calculations of the expected levels of genetic variability under different scenarios of N _e showed that the number of alleles commonly reported at microsatellite loci in Atlantic cod is best explained by N _e 's exceeding thousand. Recent fishery-induced bottlenecks can, however, not be ruled out as an explanation for the apparent discrepancy between high levels of variability and recently reported estimates of N _e  << 1000. From life history traits and estimates of survival rates in the wild, we evaluate the compatibility of the species' biology and extremely low N _e / N ratios. Our data suggest that very small N _e 's are not likely to be of general concern for cod populations and, accordingly, most populations do not face any severe threat of losing evolutionary potential due to genetic drift.     

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.     

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.     

Genetic analysis of fish genomes and populations: allozyme variation within and among Atlantic salmon from Downeast rivers of Maine

Analysis of genetic variation within and among samples of naturally produced Atlantic salmon ( n  = 372) from 7 Maine (U.S.A.) and one Canadian river were conducted based on 54 allozyme loci. Eight of the 54 loci proved polymorphic, and estimated heterozygosities ( H _S) based on all loci ranged from 0·012 to 0·026 (mean = 0·021, s . e . = 0·002). Only one of 56 tests revealed genotypic proportions that deviated significantly from Hardy–Weinberg expectations. Genetic distances ( D ) between samples ranged from 0·002 to 0·022. No obvious association existed between genetic and geographic distances. Cluster analysis of genetic distances revealed the Dennys River sample as the most differentiated when all samples were included in the analysis, though bootstrap support of the cluster analysis was generally weak. G ‐tests revealed significant differences in allele frequencies among samples at five of the polymorphic loci, and the G ‐value summed over all loci also indicated significant differences among samples. F _ST values indicated that 3·4% of the total genetic diversity was due to variability among samples, while 96·6% was due to variability within samples. These results indicate that the Atlantic salmon analyzed in this study had levels of genetic variability and differentiation among samples comparable to native populations from other areas collected across a similar geographic range.     

Morphological and genetic differentiation in anadromous smelt Osmerus mordax (Mitchill): disentangling the effects of geography and morphology on gene flow

Morphological analyses were combined with genetic analyses at nine microsatellite loci to examine the determinants of gene flow at 21 spawning locations of rainbow smelt Osmerus mordax along the east coast of Canada. Associations between morphology, geography and gene flow were examined using a computational geometric approach and partial Mantel tests. Significant barriers to gene flow and discontinuities in morphology were observed between Newfoundland and mainland Canada, as well as within Newfoundland samples. On regional scales, contrasting patterns were present with restricted gene flow between Newfoundland populations ( F _ST= c . 0·11) and high gene flow between mainland populations ( F _ST= c . 0·017). Within Newfoundland populations, geographic distance was significantly associated with gene flow ( r = 0·85, P < 0·001) contrasting mainland samples where gene flow was most associated with phenotypic divergence ( r = 0·33, P < 0·001). At large spatial scales, weak ( r = 0·19, P = 0·02) associations between gene flow and geographic distance were observed, and moderate associations were also observed between gene flow and morphology ( r = 0·28, P < 0·001). The presence of significant genetic isolation by distance in Newfoundland samples and the clear discontinuity associated with the Cabot Strait suggest geography may be the primary determinant of gene flow. Interestingly, the association between genetic and morphological divergence within mainland samples and overall, supports the hypothesis that gene flow may be moderated by morphological divergence at larger spatial scales even in high gene flow environments.     

The effect of natural selection on estimates of genetic divergence among populations of the Atlantic salmon

The effect of natural selection on the _mMEP-2 * locus on measures of genetic divergence among Atlantic salmon populations was investigated by examining the pattern of change in the level of genetic differentiation (F_ST) averaged over loci when data on the _mMEP-2 * locus were either included or excluded. The level of F_ST among populations at various geographic scales was estimated from allele frequencies at up to four loci (_s AAT-4 *, IDDH-1 *, IDHP-3 *, and _mMEP-2 *). At smaller geographic scales (within river systems or limited geographic regions) levels of variance in _mMEP-2 * allele frequencies were reduced relative to mean levels. At larger geographic scales (across continents or the species range) variation in _mMEP-2 * allele frequencies was greater than mean levels. These results suggest an a priori hypothesis for the effect of selection on the _mMEP-2 * locus which may be applied in future studies on variation in protein coding or other (e.g. mini- and microsatellite) loci in the Atlantic salmon. It is recommended that estimates of gene flow among populations of the Atlantic salmon based on mean F _ST estimates which include data on the _mMEP-2 * locus should be viewed with caution.     

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.     

Characterization of a major histocompatibility class II A gene (Clha-DAA) with an embedded microsatellite marker in Atlantic herring (Clupea harengus L.)

An Atlantic herring major histocompatibility class II A ( Clha-DAA ) cDNA sequence has been characterized and was shown to encode a leader peptide, alpha-1 domain, alpha-2 domain, connecting peptide, transmembrane and cytoplasmic region. The Clha-DAA protein sequence has all the characteristics of a teleost class II A protein with conserved cysteines in both the alpha-1 and the alpha-2 domains and two potential N-linked glycosylation sites. Exon 2 sequences encoding the polymorphic alpha-1 domain from different individuals were analysed and revealed the presence of at least two loci. The Clha-DAA gene consists of four exons and three short introns. Four unique intron 3 sequences from multiple individuals were obtained and were shown to contain a (TG)_n microsatellite sequence. Primers were optimized such that only a single microsatellite locus designated Clha-DAA-INTR3 was amplified. Four herring populations from the North Sea and the Baltic Sea were genotyped for Clha-DAA-INTR3 . In total, 16 Clha-DAA-INTR3 alleles were detected; the distribution of the alleles showed no deviation from Hardy–Weinberg expectation. Levels of genetic differentiation among samples were of similar magnitude as have been reported earlier for neutral microsatellite loci between northern North Sea and Baltic Sea herring populations.     

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.     

Reproductive characteristics and population decline of four species of skate (Rajidae) off the eastern coast of Canada

Four of the most common species of skate (Rajidae) were studied off eastern Canada to determine if their reproductive characteristics were linked to their population trajectories. The fecundity of the winter skate Leucoraja ocellata , the little skate Leucoraja erinacea , the thorny skate Amblyraja radiata and the smooth skate Malacoraja senta averaged between 41 and 56 egg cases per year for each species. For all species but L. ocellata , males matured at larger sizes and at later ages than females. Theoretical rates of population increase for non-equilibrium populations of L. ocellata ( c . 0·07), M. senta ( c . 0·14) and L. erinacea and A. radiata ( c . 0·20) were low compared to most fishes, indicating that north-west Atlantic skates are intrinsically unproductive, yet are theoretically capable of supporting low-level fisheries. Nevertheless, the results of 36 years of research surveys indicate that the abundance of mature L. ocellata , A. radiata and M. senta all decreased by >90% since 1970, indicating that past fishing mortality (both directed and undirected) has outstripped the net productivity of the skate populations on the eastern Scotian Shelf. The relationship between maximum age ( t _max) and age of maturity ( t _mat) was a better predictor of population growth rate than was body size, with the species exhibiting the highest ratios of t _mat : t _max ( L. ocellata = 0·68, M. senta = 0·66) having the lowest predicted population growth rates. L. ocellata appears to have the lowest productivity and has experienced the greatest population decline, thus raising concerns over its future status.     

Allozyme and microsatellite DNA markers of toothfish population structure in the Southern Ocean

Toothfish represent a major finfish resource in the Southern Ocean, where there has been a rapid increase in exploitation of the Patagonian toothfish Dissostichus eleginoides. Molecular genetic methods have been applied to address questions relating to stock discrimination in D. eleginoides. Eight microsatellite loci were tested in population samples from the Atlantic, Indian, and Pacific Ocean sectors of the Southern Ocean. Some loci were characterised by a large number of alleles, but no alleles were restricted to specific ocean basins. However there was a significant allelic heterogeneity in the total data indicative of population differentiation, and a high level of genetic sub-division was measured with F _ST and R _ST. Seven polymorphic allozyme loci revealed no significant heterogeneity among Pacific and Indian Ocean sector samples. The microsatellite DNA heterogeneity suggests that there is restricted gene flow through the Southern Ocean, and that the different fishing grounds may support independent stocks.     

Genetic diversity and differentiation between the two remaining populations of the critically endangered Mediterranean monk seal

The Mediterranean monk seal Monachus monachus , is a critically-endangered species of which only two populations, separated by c . 4000 km, remain: the eastern Mediterranean (150–300 individuals) and the Atlantic/western Sahara populations (100–130 individuals). We measured current levels of nuclear genetic variation at 24 microsatellite loci in 12 seals from the eastern Mediterranean and 98 seals from the western Sahara population and assessed differences between them. In both populations, genetic variation was found to be low, with mean allelic richness for the loci polymorphic in the species of 2.09 and 1.96, respectively. For most loci, the observed allele frequency distributions in both populations were discontinuous and the size ranges similar. The eastern Mediterranean population had 14 private alleles and the western Sahara had 18, but with a much larger sample size. Highly significant differences in allele frequencies between the two populations were found for 14 out of 17 loci. F _ST between the two populations was 0.578 and the estimated number of migrants per generation was 0.046, both clearly indicating substantial genetic differentiation. From a conservation perspective, these results suggest that each population may act as a source for introducing additional genetic variation into the other population.