Polymorphism at five microsatellite loci were screened to determine the genetic variability and the temporal stability of population structure in natural populations of European hake (Merluccius merluccius, L.) within the Bay of Biscay. In addition, the control region (900 bp) and two protein coding genes (ATPase, subunits 6 and 8, 842 bp and a partial sequence of the ND1, 800 bp) of the mitochondrial DNA (mtDNA) were sequenced from geographically distant populations from the extremes of the species range. One hundred individuals from either side of a supposed stock boundary within the bay were collected in autumn 1997. This sampling strategy was repeated during hake spawning seasons in late spring of 1998 and 1999. Low levels of population subdivision were found between putative populations within years. Similarly, low levels of differentiation were found between autumn 1997 northern samples and spring 1998 southern samples which were collected 7 months later on spawning grounds. These results are discussed in relation to ecological, behavioural and oceanographic information. Sampling effects, which may influence these results, are also discussed. Theta (theta) estimates were significantly different from zero in every other pairwise comparison between geographical areas (north and south of the Bay of Biscay) and between years within the same area (P<0.05). Hierarchical analysis of molecular variance (AMOVA) does not confirm the temporal persistence of population structure. These results are discussed in relation to variance in reproductive success, and temporal spawning patterns, which may exist within the bay. mtDNA variability was very low between geographically distant samples from Norway and the Mediterranean Sea with only 10 variable sites found in a total of 2542 bp of mtDNA, these differences being exclusively in the D-loop.