Fine-scale genetic structure, estuarine colonization and incipient speciation in the marine silverside fish Odontesthes argentinensis

The identification of incipient ecological species represents an opportunity to investigate current evolutionary process where adaptive divergence and reproductive isolation are associated. In this study we analysed the genetic structure of marine and estuarine populations of the silverside fish Odontesthes argentinensis using nine microsatellite loci and 396 bp of the mitochondrial DNA (mtDNA) control region. Our main objective was to investigate the relationship among estuarine colonization, divergent selection and speciation in silversides. Significant genetic structure was detected among all marine and estuarine populations. Despite the low phylogeographic structure in mtDNA haplotypes, there was clear signal of local radiations of haplotypes in more ancient populations. Divergence among marine populations was interpreted as a combined result of homing behaviour, isolation by distance and drift. On the other hand, ecological shifts due to the colonization of estuarine habitats seem to have promoted rapid adaptive divergence and reproductive isolation in estuarine populations, which were considered as incipient ecological species. This conclusion is supported by the existence of a set of environmental factors required for successful reproduction of estuarine ecotypes. The pattern of genetic structure indicates that phenotypic and reproductive divergence evolved in the face of potential gene flow between populations. We suggest that the 'divergence-with-gene-flow' model of speciation may account for the diversification of estuarine populations. The approach used can potentially identify 'incipient estuarine species', being relevant to the investigation of the evolutionary relationships of silversides in several coastal regions of the world.     

Imprints from genetic drift and mutation imply relative divergence times across marine transition zones in a pan-European small pelagic fish (Sprattus sprattus)

Geographic distributions of most temperate marine fishes are affected by postglacial recolonisation events, which have left complex genetic imprints on populations of marine species. This study investigated population structure and demographic history of European sprat (Sprattus sprattus L.) by combining inference from both mtDNA and microsatellite genetic markers throughout the species' distribution. We compared effects from genetic drift and mutation for both genetic markers in shaping genetic differentiation across four transition zones. Microsatellite markers revealed significant isolation by distance and a complex population structure across the species' distribution (overall θ(ST)=0.038, P<0.01). Across transition zones markers indicated larger effects of genetic drift over mutations in the northern distribution of sprat contrasting a stronger relative impact of mutation in the species' southern distribution in the Mediterranean region. These results were interpreted to reflect more recent divergence times between northern populations in accordance with previous findings. This study demonstrates the usefulness of comparing inference from different markers and estimators of divergence for phylogeographic and population genetic studies in species with weak genetic structure, as is the case in many marine species.     

Genetic structure of sigmodontine rodents (Cricetidae) along an altitudinal gradient of the Atlantic Rain Forest in southern Brazil

The population genetic structure of two sympatric species of sigmodontine rodents (Oligoryzomys nigripes and Euryoryzomys russatus) was examined for mitochondrial DNA (mtDNA) sequence haplotypes of the control region. Samples were taken from three localities in the Atlantic Rain Forest in southern Brazil, along an altitudinal gradient with different types of habitat. In both species there was no genetic structure throughout their distribution, although levels of genetic variability and gene flow were high.     

陕西省林麝mtDNA D-loop区序列结构和种群遗传多样性

林麝(Moschus berezovskii)曾广泛分布于中国,由于盗猎和栖息地缩小,秦岭地区野生种群数量迅速下降,圈养繁殖种群已成立了几十年,但大多数圈养种群的遗传背景不清,种群规模增长非常缓慢。为了给这一物种的保护和管理提供有用的信息,调查了陕西省林麝1个圈养种群3个野生种群线粒体DNA(mt DNA)D-Loop 632 bp片段的遗传多样性和种群结构。在69个个体中其碱基组成为A+T的平均含量63.2%高于G+C含量36.8%,共检测到变异位点171个(约占总位点数的27.05%)。核苷酸多样性(Pi)为0.04424,平均核苷酸差异数(K)为19.908。69个个体分属32个单倍型,单倍型间的平均遗传距离(P)为0.070。32个单倍型构建的NJ系统树聚为3个分支,4个林麝群体中的单倍型是随机分布的。4个群体的平均遗传距离为0.043(标准误SE为0.005),凤县养殖场群体与留坝和陇县群体的亲缘关系较远。单倍型间的平均遗传距离为0.043,可见其遗传分化尚未达到种群分化的水平。结果表明,陕西省林麝群体mt DNA D-loop区序列存在着较丰富的变异和遗传多样性,凤县野生群体和凤县养殖场群体的核苷酸多样性和单倍型多样较高,养殖场种群没有出现近亲繁殖及遗传多样性下降的情况。凤县野生群体和凤县养殖场群体两者遗传分化较小,存在着较高的基因流水平。     

Mitogenomic analysis of the Australian lungfish (<Emphasis Type="Italic">Neoceratodus forsteri</Emphasis>) reveals structuring of indigenous riverine populations and late Pleistocene movement between drainage basins

Neoceratodus forsteri: is a freshwater species of Dipnoan currently listed as ‘vulnerable to extinction’ under Australian legislation. The species is restricted to at least two indigenous riverine populations in southeastern Queensland, and several other putatively translocated populations. Current understanding of genetic relationships among populations is based on studies of allozymes, microsatellites and mitochondrial DNA (mtDNA) fragments. A notable feature of all these datasets was low genetic variability. Here we sequence the complete mitogenome of 71 N. forsteri individuals from five populations to improve resolution of mtDNA diversity, examine relationships among populations, and evaluate recent demographic history. We recorded 137 variable positions forming 41 haplotypes in the 16,573 bp mitogenome alignment. Strong genetic structure was observed among riverine samples (global Φ_ST?=?0.342) in a pattern consistent with translocation history. Tinana Creek was confirmed as an isolated and genetically unique subpopulation that should be recognized as a distinct management unit. Two previously unreported mtDNA clades (0.46% mean divergence) were found and suggest that genetic exchange among coastal catchments may have been facilitated by riverine connections on the exposed continental shelf during the late Pleistocene. Extended Bayesian skyline analysis showed no evidence for recent historical change in female effective population size, and codon-based selection tests found no evidence for positive selection in coding genes. Overall, our results emphasise the utility of the full mtDNA molecule for capturing population structure in taxa with low genetic diversity. In such cases, informative variation may be scattered across disparate parts of the mitogenome. Surveying relatively short fragments of mtDNA may lead to significant underestimates of population structure when applied to threatened species with low genetic diversity.     

Mitochondrial DNA variation in natural populations of endangered Indian Feather-Back Fish, Chitala chitala

Genetic variation at mitochondrial cytochrome b (cyt b) and D-loop region reveals the evidence of population sub-structuring in Indian populations of highly endangered primitive feather-back fish Chitala chitala. Samples collected through commercial catches from eight riverine populations from different geographical locations of India were analyzed for cyt b region (307 bp) and D-loop region (636–716 bp). The sequences of the both the mitochondrial regions revealed high haplotype diversity and low nucleotide diversity. The patterns of genetic diversity, haplotypes networks clearly indicated two distinct mitochondrial lineages and mismatch distribution strongly suggest a historical influence on the genetic structure of C. chitala populations. The baseline information on genetic variation and the evidence of population sub-structuring generated from this study would be useful for planning effective strategies for conservation and rehabilitation of this highly endangered species.     

Effects of hydrographic barriers on population genetic structure of the sea star Coscinasterias muricata (Echinodermata, Asteroidea) in the New Zealand fiords

New Zealand's 14 deep-water fiords possess persistent salinity stratification and mean estuarine circulation that may serve to isolate populations of marine organisms that have a dispersal larval phase. In order to investigate this idea, we analysed the population structure of the sea star Coscinasterias muricata using a mitochondrial DNA marker. Genetic differentiation among populations of C. muricata was analysed using 366 base pairs of mtDNA D-loop. We compared populations from the fiords with several others sampled from around New Zealand. At a macro-geographical scale (> 1000 km), restricted gene flow between the North and South Islands was observed. At a meso-geographical scale (10-200 km), significant population structure was found among fiords and between fiords and open coast. The pattern of population genetic structure among the fiords suggests a secondary contact between a northern population and a southern one, separated by a contact or mixing zone. These populations may have diverged by the effects of random genetic drift and population isolation as a consequence of the influence of estuarine circulation on dispersal. In northern Fiordland, genetic structure approximated an isolation by distance model. However, the pattern in genetic differences suggests that distance alone cannot explain the most divergent populations and that fiord hydrography may increase the effect of genetic drift within populations in the fiords. Finally, our study indicates that populations within the fiords underwent recent rapid expansion, followed most probably by genetic drift due to a lack of gene flow among the fiords.     

A Genetic Discontinuity in a Continuously Distributed Species: Mitochondrial DNA in the American Oyster, Crassostrea Virginica

Restriction site variation in mitochondrial DNA (mtDNA) of the American oyster (Crassostrea virginica) was surveyed in continuously distributed populations sampled from the Gulf of St. Lawrence, Canada, to Brownsville, Texas. mtDNA clonal diversity was high, with 82 different haplotypes revealed among 212 oysters with 13 endonucleases. The mtDNA clones grouped into two distinct genetic arrays (estimated to differ by about 2.6% in nucleotide sequence) that characterized oysters collected north vs. south of a region on the Atlantic mid-coast of Florida. The population genetic "break" in mtDNA contrasts with previous reports of near uniformity of nuclear (allozyme) allele frequencies throughout the range of the species, but agrees closely with the magnitude and pattern of mtDNA differentiation reported in other estuarine species in the southeastern United States. This concordance of mtDNA phylogenetic pattern across independently evolving species provides strong evidence for vicariant biogeographic processes in initiating intraspecific population structure. The post-Miocene ecological history of the region suggests that reduced precipitation levels in an enlarged Floridian peninsula may have created discontinuities in suitable estuarine habitat for oysters during glacial periods, and that today such population separations are maintained by the combined influence of ecological gradients and oceanic currents on larval dispersal. The results are consistent with the hypothesis that historical vicariant events, in conjunction with contemporary environmental influences on gene flow, can result in genetic discontinuities in continuously distributed species with high dispersal capability.     

Phylogeographic and population genetic analyses reveal Pleistocene isolation followed by high gene flow in a wide ranging,but endangered,freshwater mussel

Freshwater organisms of North America have had their contemporary genetic structure shaped by vicariant events, especially Pleistocene glaciations. Life history traits promoting dispersal and gene flow continue to shape population genetic structure. Cumberlandia monodonta, a widespread but imperiled (IUCN listed as endangered) freshwater mussel, was examined to determine genetic diversity and population genetic structure throughout its range. Mitochondrial DNA sequences and microsatellite loci were used to measure genetic diversity and simulate demographic events during the Pleistocene using approximate Bayesian computation (ABC) to test explicit hypotheses explaining the evolutionary history of current populations. A phylogeny and molecular clock suggested past isolation created two mtDNA lineages during the Pleistocene that are now widespread. Two distinct groups were also detected with microsatellites. ABC simulations indicated the presence of two glacial refugia and post-glacial admixture of them followed by simultaneous dispersal throughout the current range of the species. The Ouachita population is distinct from others and has the lowest genetic diversity, indicating that this is a peripheral population of the species. Gene flow within this species has maintained high levels of genetic diversity in most populations; however, all populations have experienced fragmentation. Extirpation from the center of its range likely has isolated remaining populations due to the geographic distances among them.     

Phylogeography of Varicorhinus barbatulus (Cyprinidae) in Taiwan based on nucleotide variation of mtDNA and allozymes

The phylogeographical patterns and population genetic structures of Varicorhinus barbatulus in Taiwan were investigated based on genetic diversity of 34 allozyme loci and nucleotide sequences of 3' end of the cytochrome b gene, tRNA genes, D-loop control region, and the 5' end of the 12S rRNA of mtDNA. Allozyme and mtDNA analyses revealing evident geographical structuring suggest limited gene flow between populations (F(ST)=0.511 and 0.791, respectively). Low genetic variability within populations (P=5.56%; He=0.018) based on allozymes and significantly negative Tajima's D statistics based on mtDNA suggest that most populations in Taiwan may have originated from a small number of founders followed by demographic expansion. The gene genealogy of mtDNA identified six lineages corresponding to major drainages that were separated by the geological barriers due to vicariant events. A minimum spanning network based on nucleotide substitutions reflects divergence from populations of the Miao-li Plateau to northern and southern regions of the island. In contrast to a previous hypothesis that suggests an early invasion to eastern part of Taiwan prior to the lifting of central mountain range some one million years ago, the mtDNA genealogy and molecular dating reveal very recent colonization of the eastern population. Nested clade analyses revealing significant associations between genetic structure and geographical division identify past fragmentation and range expansion as major phylogeographical events that shaped the geographical distribution of this species in Taiwan.     

Genomic population structure of freshwater‐resident and anadromous ide (Leuciscus idus) in north‐western Europe

Climate change experts largely agree that future climate change and associated rises in oceanic water levels over the upcoming decades, will affect marine salinity levels. The subsequent effects on fish communities in estuarine ecosystems however, are less clear. One species that is likely to become increasingly affected by changes in salinity is the ide (Leuciscus idus). The ide is a stenohaline freshwater fish that primarily inhabits rivers, with frequent anadromous behavior when sea salinity does not exceed 15%. Unlike most other anadromous Baltic Sea fish species, the ide has yet to be subjected to large‐scale stocking programs, and thus provides an excellent opportunity for studying the natural population structure across the current salinity gradient in the Danish Belts. To explore this, we used Genotyping‐by‐Sequencing to determine genomic population structure of both freshwater resident and anadromous ide populations in the western Baltic Sea region, and relate the results to the current salinity gradient and the demographic history of ide in the region. The sample sites separate into four clusters, with all anadromous populations in one cluster and the freshwater resident populations in the remaining three. Results demonstrate high level of differentiation between sites hosting freshwater resident populations, but little differentiation among anadromous populations. Thus ide exhibit the genomic population structure of both a typical freshwater species, and a typical anadromous species. In addition to providing a first insight into the population structure of north‐western European ide, our data also (1) provide indications of a single illegal introduction by man; (2) suggest limited genetic effects of heavy pollution in the past; and (3) indicate possible historical anadromous behavior in a now isolated freshwater population.     

Phylogeography of Hemibarbus labeo (Cyprinidae): secondary contact of ancient lineages of mtDNA

Lin, C.-J., Lin, H.-D., Wang, J.-P., Chao, S.-C. & Chiang, T.-Y. (2007). Phylogeography of Hemibarbus labeo (Cyprinidae): secondary contact of ancient lineages of mtDNA. —Zoologica Scripta, 39, 23–35. Nucleotide sequences of the D-loop control region of mtDNA were used to assess the genetic structure and phylogeography of Hemibarbus labeo in Mainland China and Taiwan. A hierarchical analysis of molecular variance of populations in 11 major streams from three geographical regions revealed significant structuring among populations and geographical regions. High levels of nucleotide diversity (π = 1.88%) and haplotype diversity (h = 96 ± 0.009%) suggest a large effective population size. A maximum likelihood tree based on mtDNA variation identified two ancient mtDNA lineages, which split approximately 3.39 million years ago. Overlapping distribution of the major lineages displayed low correspondence with geographical regions and reflects a scenario of secondary mergence after long isolation. Gene genealogy further revealed a unidirectional migration. Nevertheless, there existed a phylogeographical structure that mostly agreed with a biogeographical hypothesis. That is, within each lineage, a close phylogeny between populations of the River-Campagna and East-Pacific regions was supported by the mtDNA gene genealogy, although monophyly of each geographical region was not supported. The degree of genetic differentiation was correlated with geographical distances between populations, displaying a pattern of ‘isolation by distance’. Gene genealogy of mtDNA revealed that Yangtzejiang population may act as a divergence centre of H. labeo. In addition, Taiwan population was colonized via a recent a founder event, likely from population Yangtzejiang River about 131 000 years before present. Low haplotype number and genetic variability also suggested possible bottleneck events in the Rongjiang and Dazhangjiang populations.     

Intraspecific DNA sequence variation of the mitochondrial control region of white sturgeon (Acipenser transmontanus)

Intraspecific sequence variation in the D-loop region of mtDNA in white sturgeon (Acipenser transmontanus), a relict North American fish species, was examined in 27 individuals from populations of the Columbia and Fraser rivers. Thirty-three varied nucleotide positions were present in a 462-nucleotide D-loop sequence, amplified using the polymerase chain reaction. Bootstrapped neighbor-joining and maximum- parsimony trees of sequences from 19 haplotypes suggest that the two populations have recently diverged. This is consistent with the hypothesis that the Columbia River, a Pleistocene refugium habitat, was the source of founders for the Fraser River after the last glacial recession. On the basis of a divergence time of 10-12 thousand years ago, the estimated substitution rate of the white sturgeon D-loop region is 1.1-1.3 x 10(-7) nucleotides/site/year, which is comparable to rates for hypervariable sequences in the human D-loop region. Furthermore, the ratio of mean percent nucleotide differences in the D- loop (2.27%) to that in whole mtDNA (0.54%, as estimated from restriction-enzyme data) is 4.3, which is similar to the fourfold-to- fivefold-higher substitution rate estimated for the human D-loop. The high nucleotide substitution rate of the hypervariable region indicates that the vertebrate D-loop has potential as a genetic marker in molecular population studies.      

舟山小黄鱼线粒体DNA D-loop区序列变异的遗传多样性分析

小黄鱼(Pseudosciaena polyactis)为我国重要海产经济鱼类之一,过度捕捞和环境污染等因素造成其资源日益衰退。研究小黄鱼种群遗传结构对其资源的保护及其可持续利用有十分重要的意义。该研究采用聚合酶链式反应(PCR)技术对浙江舟山附近海域小黄鱼种群53个个体的mtDNA D-loop区全序列进行扩增,序列长度在795~801bp之间,长度差异不大。采用ClustalX1.83、MEGA3.1、DnaSP4.0等生物信息学软件进行遗传多样性分析,结果显示:53条小黄鱼线粒体DNA D-loop区的T、C、A和G碱基平均含量分别为30.3%、23.1%、32.3%和14.3%,排除13处核苷酸的插入或缺失后,共检测到93处转换和颠换位点,约占分析序列总长度的11.6%,其中包括53个单一多态位点和40个简约信息位点,共确定了52种单倍型,单倍型多样性(hd)为0.9993,单倍型间的平均遗传距离为0.012,转换/颠换平均值为4.305,平均核苷酸差异数(k)为9.73875,核苷酸多样性(π)为0.01233,表明舟山小黄鱼遗传多样性处于中等水平。     

The effect of landscape processes upon gene flow and genetic diversity in an Australian freshwater fish, Neosilurus hyrtlii

1. Flow regime and riverine architecture are two important landscape characteristics that influence genetic diversity and gene flow in riverine species.
2. Using population genetic markers (mtDNA, microsatellites and allozymes), this study aimed to investigate genetic diversity and gene flow in the freshwater fish, Neosilurus hyrtlii , across two major drainage divisions in northern and central Australia (the Gulf of Carpentaria and Lake Eyre basins). These basins lie adjacent to each other and differ in their hydrological inputs and riverine structure, providing an ideal opportunity to identify the impact of landscape processes upon population dynamics of freshwater fish.
3. Populations were strongly structured among basins, among catchments within basins and were weakly structured within catchments in the Lake Eyre Basin, providing support for the Stream Hierarchy Model.
4. Interestingly, mtDNA and microsatellite diversity was much higher in the Gulf of Carpentaria Basin compared to the Lake Eyre Basin. It was concluded that this difference was due to the extreme hydrological variability in this basin and boom-bust population cycles resulting in smaller effective population sizes in the Lake Eyre Basin.     

Single genetic stock of kawakawa Euthynnus affinis (Cantor, 1849) along the Indian coast inferred from sequence analyses of mitochondrial DNA D-loop region

Kawakawa Euthynnus affinis is an epipelagic migratory tuna species, widely distributed in the tropical and subtropical waters of the Indo-Pacific region. Kawakawa constitutes the largest tuna fishery in the Indian waters. In the present study, genetic variation was assessed using sequence analyses of Mitochondrial DNA (mtDNA) D-loop region. A 500 bp segment of D-loop region was sequenced in 400 samples collected from eight localities (Veraval (VE), Ratnagiri (RA), Kochi (KO), Kavaratti (KA), Port-Blair (PB), Tuticorin (TU), Pondicherry (PO), and Vizag (VI)) along the Indian coast. Analysis of molecular variance of mtDNA data revealed no significant genetic differentiation among sites the (Φ _ST ?=?0.0028, P?=?0.20723) indicating a single population along the Indian coast. Phylogenetic analysis revealed no obvious phylogeographic pattern separating the eight samples of kawakawa. However, the genealogical relationships demonstrated that mtDNA D-loop sequences belong to two different clades (clade I and clade II). Clade I is the major clade which consists of more than 98?% specimens from each regional population while clade II has individuals from only three populations (VE, PO, and VI). Results of genetic analyses of the present study support a single stock management of kawakawa along the Indian coast.     

The history of the introduction of the giant river prawn, Macrobrachium cf. rosenbergii (Decapoda, Palaemonidae), in Brazil: New insights from molecular data

The giant river prawn, Macrobrachium cf. rosenbergii, is one of the most cultivated freshwater prawns in the world and has been introduced into more than 40 countries. In some countries, this prawn is considered an invasive species that requires close monitoring. Recent changes in the taxonomy of this species (separation of M. rosenbergii and M. dacqueti) require a re-evaluation of introduced taxa. In this work, molecular analyses were used to determine which of these two species was introduced into Brazil and to establish the geographic origin of the introduced populations that have invaded Amazonian coastal waters. The species introduced into Brazil was M. dacqueti through two introduction events involving prawns originating from Vietnam and either Bangladesh or Thailand. These origins differ from historical reports of the introductions and underline the need to confirm the origin of other exotic populations around the world. The invading populations in Amazonia require monitoring not only because the biodiversity of this region may be affected by the introduction, but also because admixture of different native haplotypes can increase the genetic variability and the likelihood of persistence of the invading species in new habitats.     

Application of mtDNA d-loop region for the study of Russian sturgeon population structure from Iranian coastline of the Caspian Sea

The population structure of the Russian sturgeon ( Acipenser gueldenstaedtii ) was investigated using PCR amplification of the mtDNA D-loop region. Seven composite haplotypes were detected and average nucleotide and haplotype diversity over all samples were found to be 0.05 ± 0.00 and 0.75 ± 0.00 (mean ± SE) respectively. Restriction digest of the mtDNA D-loop region detected two genotypes A and B with a relative high frequencies of 0.5 and 0.36 respectively. These two genotypes (A and B) can be considered as potential genetic markers for different biological groups, stocks or seasonal races of Russian sturgeon. A maximum of 10.8% nucleotide diversity was found between haplotypes AAAAB and BBBBF. 100% heteroplasmy was found in Russian sturgeon and restriction digest of the mtDNA D-loop region also exhibited site heteroplasmy.     

Analysis of mitochondrial DNA in Bolivian llama,alpaca and vicuna populations: a contribution to the phylogeny of the South American camelids

The objectives of this work were to assess the mtDNA diversity of Bolivian South American camelid (SAC) populations and to shed light on the evolutionary relationships between the Bolivian camelids and other populations of SACs. We have analysed two different mtDNA regions: the complete coding region of the MT‐CYB gene and 513 bp of the D‐loop region. The populations sampled included Bolivian llamas, alpacas and vicunas, and Chilean guanacos. High levels of genetic diversity were observed in the studied populations. In general, MT‐CYB was more variable than D‐loop. On a species level, the vicunas showed the lowest genetic variability, followed by the guanacos, alpacas and llamas. Phylogenetic analyses performed by including additional available mtDNA sequences from the studied species confirmed the existence of the two monophyletic clades previously described by other authors for guanacos (G) and vicunas (V). Significant levels of mtDNA hybridization were found in the domestic species. Our sequence analyses revealed significant sequence divergence within clade G, and some of the Bolivian llamas grouped with the majority of the southern guanacos. This finding supports the existence of more than the one llama domestication centre in South America previously suggested on the basis of archaeozoological evidence. Additionally, analysis of D‐loop sequences revealed two new matrilineal lineages that are distinct from the previously reported G and V clades. The results presented here represent the first report on the population structure and genetic variability of Bolivian camelids and may help to elucidate the complex and dynamic domestication process of SAC populations.     

Temporal and spatial genetic variation in spawning grounds of European hake (Merluccius merluccius) in the Bay of Biscay

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.