History, chance and adaptation during biological invasion: separating stochastic phenotypic evolution from response to selection

Introduced species often exhibit changes in genetic variation, population structure, selection regime and phenotypic traits as they colonize and expand into new ranges. For these reasons, species invasions are increasingly recognized as promising systems for studying adaptive evolution over contemporary time scales. However, changes in phenotypic traits during invasion occur under non-equilibrium demographic conditions and may reflect the influences of prior evolutionary history and chance events, as well as selection. We briefly review the evidence for phenotypic evolution and the role of selection during invasion. While there is ample evidence for evolutionary change, it is less clear if selection is the primary mechanism. We then discuss the likelihood that stochastic events shift phenotypic distributions during invasion, and argue that hypotheses of adaptation should be tested against appropriate null models. We suggest two experimental frameworks for separating stochastic evolution from adaptation: statistically accounting for phenotypic variation among putative invasion sources identified by using phylogenetic or assignment methods and by comparing estimates of differentiation within and among ranges for both traits and neutral markers ( Q _ST vs. F _ST). Designs that incorporate a null expectation can reveal the role of history and chance in the evolutionary process, and provide greater insights into evolution during species invasions.     

Geographic patterns of genetic differentiation and plumage colour variation are different in the pied flycatcher (Ficedula hypoleuca)

The pied flycatcher is one of the most phenotypically variable bird species in Europe. The geographic variation in phenotypes has often been attributed to spatial variation in selection regimes that is associated with the presence or absence of the congeneric collared flycatcher. Spatial variation in phenotypes could however also be generated by spatially restricted gene flow and genetic drift. We examined the genetic population structure of pied flycatchers across the breeding range and applied the phenotypic Q _ST ( P _ST)– F _ST approach to detect indirect signals of divergent selection on dorsal plumage colouration in pied flycatcher males. Allelic frequencies at neutral markers were found to significantly differ among populations breeding in central and southern Europe whereas northerly breeding pied flycatchers were found to be one apparently panmictic group of individuals. Pairwise differences between phenotypic ( P _ST) and neutral genetic distances ( F _ST) were positively correlated after removing the most differentiated Spanish and Swiss populations from the analysis, suggesting that genetic drift may have contributed to the observed phenotypic differentiation in some parts of the pied flycatcher breeding range. Differentiation in dorsal plumage colouration however greatly exceeded that observed at neutral genetic markers, which indicates that the observed pattern of phenotypic differentiation is unlikely to be solely maintained by restricted gene flow and genetic drift.     

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.     

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.     

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.     

Scanning the genome for gene single nucleotide polymorphisms involved in adaptive population differentiation in white spruce

Conifers are characterized by a large genome size and a rapid decay of linkage disequilibrium, most often within gene limits. Genome scans based on noncoding markers are less likely to detect molecular adaptation linked to genes in these species. In this study, we assessed the effectiveness of a genome-wide single nucleotide polymorphism (SNP) scan focused on expressed genes in detecting local adaptation in a conifer species. Samples were collected from six natural populations of white spruce ( Picea glauca ) moderately differentiated for several quantitative characters. A total of 534 SNPs representing 345 expressed genes were analysed. Genes potentially under natural selection were identified by estimating the differentiation in SNP frequencies among populations ( F _ST) and identifying outliers, and by estimating local differentiation using a Bayesian approach. Both average expected heterozygosity and population differentiation estimates ( H _E = 0.270 and F _ST = 0.006) were comparable to those obtained with other genetic markers. Of all genes, 5.5% were identified as outliers with F _ST at the 95% confidence level, while 14% were identified as candidates for local adaptation with the Bayesian method. There was some overlap between the two gene sets. More than half of the candidate genes for local adaptation were specific to the warmest population, about 20% to the most arid population, and 15% to the coldest and most humid higher altitude population. These adaptive trends were consistent with the genes' putative functions and the divergence in quantitative traits noted among the populations. The results suggest that an approach separating the locus and population effects is useful to identify genes potentially under selection. These candidates are worth exploring in more details at the physiological and ecological levels.     

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.     

Population divergence in the amphicarpic species Amphicarpaea edgeworthii Benth. (Fabaceae): microsatellite markers and leaf morphology

Comparative analyses of the genetic differentiation in microsatellite markers ( F _ST) and leaf morphology characters ( Q _ST) of Amphicarpaea edgeworthii Benth. were conducted to gain insight into the roles of random processes and natural selection in the population divergence. Simple sequence repeat analyses on 498 individuals of 19 natural populations demonstrate that a significant genetic differentiation occurs among populations (mean F _ST = 0.578), and A. edgeworthii is a highly self-fertilized species (mean selfing rate s  = 0.989). The distribution pattern of genetic diversity in this species shows that central populations possess high genetic diversity (e.g. population WL with H _E = 0.673 and population JG with H _E = 0.663), whereas peripheral ones have a low H _E as in population JD (0.011). The morphological divergence of leaf shape was estimated by the elliptical Fourier analysis on the data from 11 natural and four common garden populations. Leaf morphology analyses indicate the morphological divergence does not show strong correlation with the genetic differentiation ( R  = 0.260, P  = 0.069). By comparing the 95% confidence interval of Q _ST with that of F _ST, Q _ST values for five out of 12 quantitative traits are significantly higher than the average F _ST value over eight microsatellite loci. The comparison of F _ST and Q _ST suggests that two kinds of traits can be driven by different evolutionary forces, and the population divergence in leaf morphology is shaped by local selections.  © 2009 The Linnean Society of London, Biological Journal of the Linnean Society , 2009, 96 , 505–516.     

Landscape genetics, historical isolation and cross-Andean gene flow in the wax palm, Ceroxylon echinulatum (Arecaceae)

Knowledge of the role of landscapes in shaping genetic connectivity and divergence is essential for understanding patterns of biogeography and diversity. This is particularly relevant for the Andes region, a major biodiversity hotspot of relatively recent origin. We examined the phylogeography and landscape genetics of the Andean wax palm Ceroxylon echinulatum (Arecaceae) that occurs in two narrow bands of montane forests on each side of the Andes in Ecuador and northeastern Peru. First, we tested the hypothesis of C. echinulatum being a geographic cline species crossing the Andes in the Amotape–Huancabamba zone (AHZ) of southern Ecuador/northern Peru, as indicated by observations on fruit morphology. Second, we assessed the timeframe of cross-Andean divergence, and third, we investigated the impact of contemporary and historical landscape features on observed spatio-genetic patterns. Individual-based Bayesian clustering (BC) identified a northeastern, southeastern, southwestern, and northwestern cluster, with areas of genetic discontinuity coinciding with the Andes and the Giron–Paute deflection. F -statistics derived from BC suggested an east-to-west dispersal history. Population-based analyses revealed strong genetic structuring at both small and large geographic scales. Interpopulation relationships and Mantel tests strongly supported the cline model with cross-Andean dispersal in the AHZ. Along the cline, gene flow measured as F _ST was mainly limited by distance, with less but significant impact of climatic friction. Coalescent analysis revealed that cross-Andean divergence took place during the Quaternary. Significant historical isolation ( R _ST >  F _ST) was found in the southwestern population. The current study illustrates a joint effect of founder dynamics, divergence by distance and historical isolation on patterns of Andean diversity and distribution.     

Evidence for morphological and adaptive genetic divergence between lake and stream habitats in European minnows (Phoxinus phoxinus, Cyprinidae)

Natural selection drives local adaptation, potentially even at small temporal and spatial scales. As a result, adaptive genetic and phenotypic divergence can occur among populations living in different habitats. We investigated patterns of differentiation between contrasting lake and stream habitats in the cyprinid fish European minnow (Phoxinus phoxinus) at both the morphological and genomic levels using geometric morphometrics and AFLP markers, respectively. We also used a spatial correlative approach to identify AFLP loci associated with environmental variables representing potential selective forces responsible for adaptation to divergent habitats. Our results identified different morphologies between lakes and streams, with lake fish presenting a deeper body and caudal peduncle compared to stream fish. Body shape variation conformed to a priori predictions concerning biomechanics and swimming performance in lakes vs. streams. Moreover, morphological differentiation was found to be associated with several environmental variables, which could impose selection on body and caudal peduncle shape. We found adaptive genetic divergence between these contrasting habitats in the form of 'outlier' loci (2.9%) whose genetic divergence exceeded neutral expectations. We also detected additional loci (6.6%) not associated with habitat type (lake vs. stream), but contributing to genetic divergence between populations. Specific environmental variables related to trophic dynamics, landscape topography and geography were associated with several neutral and outlier loci. These results provide new insights into the morphological divergence and genetic basis of adaptation to differentiated habitats.     

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.     

Parallel changes in genetic diversity and species diversity following a natural disturbance

We examined the spatial and temporal variation of species diversity and genetic diversity in a metacommunity comprising 16 species of freshwater gastropods. We monitored species abundance at five localities of the Ain river floodplain in southeastern France, over a period of four years. Using 190 AFLP loci, we monitored the genetic diversity of Radix balthica , one of the most abundant gastropod species of the metacommunity, twice during that period. An exceptionally intense drought occurred during the last two years and differentially affected the study sites. This allowed us to test the effect of natural disturbances on changes in both genetic and species diversity. Overall, local (alpha) diversity declined as reflected by lower values of gene diversity H _S and evenness. In parallel, the among-sites (beta) diversity increased at both the genetic ( F _ST) and species ( F _STC) levels. These results suggest that disturbances can lead to similar changes in genetic and community structure through the combined effects of selective and neutral processes.     

Ancestry and divergence of subtropical montane forest isolates: molecular biogeography of the genus Abies (Pinaceae) in southern México and Guatemala

The genus Abies has a complex history in southern México and Guatemala. In this region, four closely related species, Abies flinckii , A. guatemalensis , A. hickelii , and A. religiosa , are distributed in fragmented and isolated montane populations. Range-wide genetic variation was investigated across species using cytoplasmic DNA markers with contrasted inheritance. Variation at two maternally inherited mitochondrial DNA markers was low. All species shared two of the nine mitotypes detected, while the remaining seven mitochondrial DNA types were restricted to a few isolated stands. Mitochondrial genetic differentiation across taxa was high ( G _ST = 0.933), it was not related to the taxonomic identity ( amova ; P  > 0.05) of the populations, and it was not phylogeographically structured ( G _ST ≈  N _ST). In contrast, variation at three paternally inherited chloroplast DNA microsatellites was high. Chloroplast genetic differentiation was lower ( G _ST = 0.402; R _ST = 0.547) than for mitochondrial DNA, but it was significantly related to taxonomy ( amova ; P  < 0.001), and exhibited a significant phylogeographical structure ( G _ST <  R _ST). Different analyses of population structure indicated that A. flinckii was the most divergent taxon, while the remaining three species formed a relatively homogeneous group. However, a small number of the populations of these three taxa, all located at the limits of their respective ranges or in the Transverse Volcanic Belt, diverged from this main cluster. These trends suggest that the Mesoamerican Abies share a recent common ancestor and that their divergence and speciation is mainly driven by genetic drift and isolation during the warm interglacial periods.     

Allozyme Diversity in Populations of Cymbidium goeringii (Orchidaceae)

Abstract: Using 14 allozyme loci, we investigated levels of genetic diversity within populations, and degree of genetic divergence among 24 populations of Cymbidium goeringii (Orchidaceae) in Korea and Japan. Cymbidium goeringii maintains high levels of genetic diversity both at population (mean expected heterozygosity, H _e = 0.238) and species levels (0.260). Means of H _e found in 24 populations were not significantly different from each other. About 90 % of the total variation in the species is common to all populations (mean G _ST = 0.108). No unique allele was found in any population. The indirect estimate of gene flow based on the mean G _ST was high ( Nm = 2.06). Nei's genetic identities for pairs of populations had high values (mean = 0.974 [SD = 0.013]). The Mantel-Z test showed a significant correlation between genetic distance and geographic distance. However, the mean G _ST value between 17 populations in Korea and seven Japanese populations was relatively low (0.029), even though the land connection between the southern Korean peninsula and southern Japanese archipelagos has not existed since the middle Pleistocene. Large numbers of small seeds of C. goeringii might travel long distances by wind from populations to populations both in Korea and Japan, increasing genetic diversity within populations and maintaining low genetic differentiation among populations.     

Genetically based phenotype-environment association for swimming behavior in lake whitefish ecotypes (Coregonus clupeaformis Mitchill)

Abstract.— Studies of phenotype-environment associations in adaptive radiation have focused largely on morphological traits related to resource-based phenotypic differences. The genetic basis of adaptive behaviors implicated in population divergence remains poorly understood, as few studies have tested the hypothesis of behavioral phenotype-environment associations. We provide evidence of a phenotype-environment association for differential adaptive swimming behaviors through experiments conducted on dwarf, normal, and hybrid lake whitefish ( Coregonus clupeaformis ). Highly significant differences were observed for depth selection, directional changes, and burst swimming, implicating a genetic basis for these behaviors. Hybrid crosses revealed that depth selection is under additive genetic control, while dominance effects were suggested for directional changes and burst swimming. Estimates for the genetic basis of behavioral differentiation from an animal model were consistent with these observations. Comparative estimates of behavioral differentiation ( Q _ST) against neutral expectations ( F _ST) revealed pronounced departures from neutral expectations in all three behavioral phenotypes, consistent with the hypothesis that directional selection has driven the divergence of behavior in dwarf and normal lake whitefish ecotypes.     

Microsatellite variation in cod Gadus morhua throughout its geographic range

Previous genetic studies using neutral markers such as allozymes, mtDNA and minisatellite loci have demonstrated varying amounts of population structure in cod Gadus morhua throughout the Atlantic. Microsatellite loci, which are potentially the most informative of presently available neutral genetic markers, have been applied extensively within western and eastern Atlantic areas but not on a range-wide basis. In the present study, six microsatellite DNA loci were used to screen cod samples from nine locations throughout the geographic range from the Scotian Shelf in the West Atlantic to the Barents and Baltic Seas in the east. Overall F _ST value was 0·03 ( P = < 0·001) across all samples. Statistically significant population differences over all loci combined were evident between more geographically distant samples, using either heterogeneity tests or F _ST analysis, with at least one locus showing significant differences between all samples (prior to Bonferroni correction). A significant correlation was observed between genetic and geographical distance, suggesting a higher level of historical and contemporary gene flow between adjacent populations than more distant populations. Samples from either end of the geographic range (Scotian Shelf and Baltic Sea) were particularly distinct when analysed using the STRUCTURE programme and also showed a high level of self-assignment when individuals of either the Scotian Shelf or Baltic Sea were tested against the entire data set. The present microsatellite study demonstrates a high level of geographic population structure between the western Atlantic, middle and eastern Atlantic and Baltic Sea, and thus, the findings should be useful in devising overall management and conservation strategies for the species.     

Molecular and quantitative trait variation within and among populations of the intertidal copepod Tigriopus californicus

Abstract While molecular and quantitative trait variation may be theoretically correlated, empirical studies using both approaches frequently reveal discordant patterns, and these discrepancies can contribute to our understanding of evolutionary processes. Here, we assessed genetic variation in six populations of the copepod Tigriopus californicus. Molecular variation was estimated using five polymorphic microsatellite loci, and quantitative variation was measured using 22-life history and morphometric characters. Within populations, no correlation was found between the levels of molecular variation (heterozygosity) and quantitative variation (heritability). Between populations, quantitative subdivision ( Q _ST) was correlated with molecular subdivision when measured as F _ST but not when measured as R _ST. Unlike most taxa studied to date, the overall level of molecular subdivision exceeded the level of quantitative subdivision ( F _ST= 0.80, R_ST = 0.89, Q _ST = 0.30). Factors that could contribute to this pattern include stabilizing or fluctuating selection on quantitative traits or accelerated rates of molecular evolution.     

Thinking outside the barrier: neutral and adaptive divergence in Indo-Pacific coral reef faunas

Population divergence in marine species frequently occurs where there is no conspicuous geographic feature to mark a population boundary. In the absence of a discernible barrier to gene flow, strong conclusions about the processes that drive population divergence in the ocean are often elusive. Because our knowledge of patterns far outpaces our understanding of processes, there is a need to expand hypotheses regarding population divergence in marine species. Here, I posit that certain population-level processes (range expansion, spatial diffusion constraints and metapopulation dynamics) are primary causes of spatial structuring in marine populations rather than auxiliary factors that aid divergence across an incomplete physical barrier to gene flow. In particular, I use examples from Indo-Pacific coral reef organisms to illustrate this point. In connection, spatial patterns of adaptive divergence are unlikely to be correlated with patterns of selectively neutral divergence in many cases, as the processes affecting adaptive and neutral genetic markers are different.     

PATTERNS OF GENETIC ARCHITECTURE FOR LIFE-HISTORY TRAITS AND MOLECULAR MARKERS IN A SUBDIVIDED SPECIES

Abstract Understanding the utility and limitations of molecular markers for predicting the evolutionary potential of natural populations is important for both evolutionary and conservation genetics. To address this issue, the distribution of genetic variation for quantitative traits and molecular markers is estimated within and among 14 permanent lake populations of Daphnia pulicaria representing two regional groups from Oregon. Estimates of population subdivision for molecular and quantitative traits are concordant, with Q _ST generally exceeding G _ST. There is no evidence that microsatellites loci are less informative about subdivision for quantitative traits than are allozyme loci. Character-specific comparison of Q _ST and G _ST support divergent selection pressures among populations for the majority of life-history traits in both coast and mountain regions. The level of within-population variation for molecular markers is uninformative as to the genetic variation maintained for quantitative traits. In D. pulicaria , regional differences in the frequency of sex may contribute to variation in the maintenance of expressed within-population quantitative-genetic variation without substantially impacting diversity at the genic level. These data are compared to an identical dataset for 17 populations of the temporary-pond species, D. pulex .     

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.