Genetic similarity of disjunct populations of the giant sea bass Stereolepis gigas

Sequence variation in nuclear and mitochondrial genes of the giant sea bass Stereolepis gigas collected from the Pacific coast and the northern Sea of Cortez was examined. Restriction fragment length polymorphism analysis and direct sequencing showed extremely low mtDNA sequence diversity (13 closely related haplotypes with no evidence of geographical population subdivision). The mitochondrial haplotype mismatch distribution is consistent with a population expansion following the Last Pleistocene glaciation. Differences in single nucleotide polymorphism frequencies between Pacific and Sea of Cortez populations were detected at two of four nuclear loci, which may reflect natural selection or genetic drift in populations with low effective numbers of males. Although Pacific coast and Sea of Cortez populations of giant sea bass do not exhibit the mitochondrial phylogenetic break characteristic of many species with disjunct Pacific and Gulf populations, the possibility of genetic differentiation at nuclear loci suggests that a cautious approach to broodstock selection for captive breeding and restoration programmes be exercised.     

Mitochondrial DNA under siege in avian phylogeography

Mitochondrial DNA (mtDNA) has been the workhorse of research in phylogeography for almost two decades. However, concerns with basing evolutionary interpretations on mtDNA results alone have been voiced since the inception of such studies. Recently, some authors have suggested that the potential problems with mtDNA are so great that inferences about population structure and species limits are unwarranted unless corroborated by other evidence, usually in the form of nuclear gene data. Here we review the relative merits of mitochondrial and nuclear phylogeographical studies, using birds as an exemplar class of organisms. A review of population demographic and genetic theory indicates that mitochondrial and nuclear phylogeographical results ought to concur for both geographically unstructured populations and for populations that have long histories of isolation. However, a relatively common occurrence will be shallow, but geographically structured mtDNA trees--without nuclear gene corroboration--for populations with relatively shorter periods of isolation. This is expected because of the longer coalescence times of nuclear genes (approximately four times that of mtDNA); such cases do not contradict the mtDNA inference of recent isolation and evolutionary divergence. Rather, the nuclear markers are more lagging indicators of changes in population structure. A review of the recent literature on birds reveals the existence of relatively few cases in which nuclear markers contradict mitochondrial markers in a fashion not consistent with coalescent theory. Preliminary information from nuclear genes suggests that mtDNA patterns will prove to be robust indicators of patterns of population history and species limits. At equilibrium, mitochondrial loci are generally a more sensitive indicator of population structure than are nuclear loci, and mitochondrial estimates of F(ST)-like statistics are generally expected to exceed nuclear ones. Hence, invoking behavioural or ecological explanations of such differences is not parsimonious. Nuclear genes will prove important for quantitative estimates of the depths of haplotype trees, rates of population growth and values of gene flow.     

Inferring population history and demography using microsatellites, mitochondrial DNA, and major histocompatibility complex (MHC) genes

Microsatellites and mitochondrial DNA (mtDNA) have traditionally been used in population genetics because of their variability and presumed neutrality, whereas genes of the major histocompatibility complex (MHC) are increasingly of interest because strong selective pressures shape their standing variation. Despite the potential for MHC genes, microsatellites, and mtDNA sequences to complement one another in deciphering population history and demography, the three are rarely used in tandem. Here we report on MHC, microsatellite, and mtDNA variability in a single large population of the eastern tiger salamander (Ambystoma tigrinum tigrinum). We use the mtDNA mismatch distribution and, on microsatellite data, the imbalance index and bottleneck tests to infer aspects of population history and demography. Haplotype and allelic variation was high at all loci surveyed, and heterozygosity was high at the nuclear loci. We find concordance among neutral molecular markers that suggests our study population originated from post-Pleistocene expansions of multiple, fragmented sources that shared few migrants. Differences in N(e) estimates derived from haploid and diploid genetic markers are potentially attributable to secondary contact among source populations that experienced rapid mtDNA divergence and comparatively low levels of nuclear DNA divergence. We find strong evidence of natural selection acting on MHC genes and estimate long-term effective population sizes (N(e)) that are very large, making small selection intensities significant evolutionary forces in this population.     

Steep,coincident, and concordant clines in mitochondrial and nuclear‐encoded genes in a hybrid zone between subspecies of Atlantic killifish,Fundulus heteroclitus

Steep genetic clines resulting from recent secondary contact between previously isolated taxa can either gradually erode over time or be stabilized by factors such as ecological selection or selection against hybrids. We used patterns of variation in 30 nuclear and two mitochondrial SNPs to examine the factors that could be involved in stabilizing clines across a hybrid zone between two subspecies of the Atlantic killifish, Fundulus heteroclitus. Increased heterozygote deficit and cytonuclear disequilibrium in populations near the center of the mtDNA cline suggest that some form of reproductive isolation such as assortative mating or selection against hybrids may be acting in this hybrid zone. However, only a small number of loci exhibited these signatures, suggesting locus‐specific, rather than genomewide, factors. Fourteen of the 32 loci surveyed had cline widths inconsistent with neutral expectations, with two SNPs in the mitochondrial genome exhibiting the steepest clines. Seven of the 12 putatively non‐neutral nuclear clines were for SNPs in genes related to oxidative metabolism. Among these putatively non‐neutral nuclear clines, SNPs in two nuclear‐encoded mitochondrial genes (SLC25A3 and HDDC2), as well as SNPs in the myoglobin, 40S ribosomal protein S17, and actin‐binding LIM protein genes, had clines that were coincident and concordant with the mitochondrial clines. When hybrid index was calculated using this subset of loci, the frequency distribution of hybrid indices for a population located at the mtDNA cline center was non‐unimodal, suggesting selection against advanced‐generation hybrids, possibly due to effects on processes involved in oxidative metabolism.     

Population genetic structure of savannah elephants in Kenya: conservation and management implications

We investigated population genetic structure and regional differentiation among African savannah elephants in Kenya using mitochondrial and microsatellite markers. We observed mitochondrial DNA (mtDNA) nucleotide diversity of 1.68% and microsatellite variation in terms of average number of alleles, expected and observed heterozygosities in the total study population of 10.20, 0.75, and 0.69, respectively. Hierarchical analysis of molecular variance of mtDNA variation revealed significant differentiation among the 3 geographical regions studied (F(CT) = 0.264; P < 0.05) and a relatively lower differentiation among populations within regions (F(SC) = 0.218; P < 0.0001). Microsatellite variation significantly differentiated among populations within regions (F(SC) = 0.019; P < 0.0001) but not at the regional levels (F(CT) = 0.000; P > 0.500). We attribute the high differentiation at the mitochondrial genome to the matrilineal social structure of elephant populations, female natal philopatry, and probably ancient vicariance. Lack of significant regional differentiation at the nuclear loci vis-a-vis strong differences at mtDNA loci between regions is likely the effect of subsequent homogenization through male-mediated gene flow. Our results depicting 3 broad regional mtDNA groups and the observed population genetic differentiation as well as connectivity patterns should be incorporated in the planning of future management activities such as translocations.     

Large discrepancies in differentiation of allozymes, nuclear and mitochondrial DNA loci in recently founded Pacific populations of the pearl oyster Pinctada margaritifera

This study presents a comparative analysis of population structure applied to the pearl oyster (Pinctada margaritifera) from the Central Pacific islands using three classes of molecular markers: two mitochondrial genes (mtDNA), five anonymous nuclear loci (anDNA), and eight polymorphic allozymes. Very low levels of haplotype diversity and nucleotidic divergence detected for mtDNA validate the hypothesis of a recent (re)colonization of Polynesian lagoons after their exondation during the last glaciations. Some nuclear loci, however, showed highly significant FST values, indicating a reduced amount of larval exchange between archipelagos at present. A large interlocus variance of FST was nevertheless observed. We discuss whether this pattern is inherent to the stochasticity of the drift process since recolonization, or if it could result from balancing selection acting on certain loci. This study illustrates once more the need to combine the analysis of several kinds of loci when unrelated phenomena are likely to leave their footprints on genetic structure.     

Contrasting population structure from nuclear intron sequences and mtDNA of humpback whales

Powerful analyses of population structure require information from multiple genetic loci. To help develop a molecular toolbox for obtaining this information, we have designed universal oligonucleotide primers that span conserved intron-exon junctions in a wide variety of animal phyla. We test the utility of exon-primed, intron-crossing amplifications by analyzing the variability of actin intron sequences from humpback, blue, and bowhead whales and comparing the results with mitochondrial DNA (mtDNA) haplotype data. Humpback actin introns fall into two major clades that exist in different frequencies in different oceanic populations. It is surprising that Hawaii and California populations, which are very distinct in mtDNAs, are similar in actin intron alleles. This discrepancy between mtDNA and nuclear DNA results may be due either to differences in genetic drift in mitochondrial and nuclear genes or to preferential movement of males, which do not transmit mtDNA to offspring, between separate breeding grounds. Opposing mtDNA and nuclear DNA results can help clarify otherwise hidden patterns of structure in natural populations.      

A tight balance between natural selection and gene flow in a southern African arid-zone endemic bird

Gene flow is traditionally thought to be antagonistic to population differentiation and local adaptation. However, recent studies have demonstrated that local adaptation can proceed provided that selection is greater than the homogenizing effects of gene flow. We extend these initial studies by combining ecology (climate), phenotype (body size), physiological genetics (oxidative phosphorylation genes), and neutral loci (nuclear microsatellites and introns) to test whether selection can counter-balance gene flow and hence promote local adaptation in a bird whose distribution spans an aridity gradient. Our results show that the Karoo scrub-robin's climatic niche is spatially structured, providing the potential for local adaptation to develop. We found remarkably discordant patterns of divergence among mtDNA, morphology, and neutral loci. For the mitochondrial genes, two amino acid replacements, strong population structure and reduced gene flow were associated with the environmental gradient separating western coastal sites from the interior of southern Africa. In contrast, morphology and the neutral loci exhibited variation independent of environmental variables, and revealed extensive levels of gene flow across the aridity gradient, 50 times larger than the estimates for mitochondrial genes. Together, our results suggest that selective pressures on physiology, mediated by the mitochondrial genome, may well be a common mechanism for facilitating local adaptation to new climatic conditions.     

Microsatellite analysis supports mitochondrial phylogeography of the hellbender (Cryptobranchus alleganiensis)

We investigated genetic diversity of the hellbender (Cryptobranchus alleganiensis) throughout its range in the eastern US using nuclear markers and compared our results to a previously published mitochondrial analysis. A variety of nuclear markers, including protein-coding gene introns and microsatellites were tested but only microsatellites were variable enough for population level analysis. Microsatellite loci showed moderate among population sharing of alleles, in contrast to the reciprocal monophyly exhibited by mitochondrial DNA. However, analyses using F-statistics and Bayesian clustering algorithms showed considerable population subdivision and clustered hellbender populations into the same major groups as the mtDNA. The microsatellites combined with the mtDNA data suggest that gene flow is severely restricted or non-existent among eight major groups, and potentially among populations (rivers) within groups. The combined mtDNA and microsatellite data suggest that the currently recognized hellbender subspecies are paraphyletic. We suggest that the eight independent groups identified in our study should be managed as such, rather than basing conservation decisions on the two named subspecies of hellbender.     

A quantitative genetic analysis of nuclear-cytoplasmic male sterility in structured populations of Silene vulgaris

Gynodioecy, the coexistence of functionally female and hermaphroditic morphs within plant populations, often has a complicated genetic basis involving several cytoplasmic male-sterility factors and nuclear restorers. This complexity has made it difficult to study the genetics and evolution of gynodioecy in natural populations. We use a quantitative genetic analysis of crosses within and among populations of Silene vulgaris to partition genetic variance for sex expression into nuclear and cytoplasmic components. We also use mitochondrial markers to determine whether cytoplasmic effects on sex expression can be traced to mitochondrial variance. Cytoplasmic variation and epistatic interactions between nuclear and cytoplasmic loci accounted for a significant portion of the variation in sex expression among the crosses. Source population also accounted for a significant portion of the sex ratio variation. Crosses among populations greatly enhanced the dam (cytoplasmic) effect, indicating that most among-population variance was at cytoplasmic loci. This is supported by the large among-population variance in the frequency of mitochondrial haplotypes, which also accounted for a significant portion of the sex ratio variance in our data. We discuss the similarities between the population structure we observed at loci that influence sex expression and previous work on putatively neutral loci, as well as the implications this has for what mechanisms may create and maintain population structure at loci that are influenced by natural selection.     

Discrepancies in population differentiation at microsatellites, mitochondrial DNA and plumage colour in the pied flycatcher--inferring evolutionary processes

Genetic differentiation between three populations of the pied flycatcher Ficedula hypoleuca (Norway, Czech Republic and Spain, respectively) was investigated at microsatellite loci and mitochondrial DNA (mtDNA) sequences and compared with the pattern of differentiation of male plumage colour. The Czech population lives sympatrically with the closely related collared flycatcher (F. albicollis) whereas the other two are allopatric. Allopatric populations are on average more conspicuously coloured than sympatric ones, a pattern that has been explained by sexual selection for conspicuous colour in allopatry and a character displacement on breeding plumage colour in sympatry that reduces the rate of hybridization with the collared flycatcher. The Czech population was genetically indistinguishable from the Norwegian population at microsatellite loci and mtDNA sequences. Recent isolation and/or gene flow may explain the lack of genetic differentiation. Accordingly, different selection on plumage colour in the two populations is either sufficiently strong so that gene flow has little impact on the pattern of colour variation, or differentiation of plumage colour occurred so recently that the (presumably) neutral, fast evolving markers employed here are unable to reflect the differentiation. Genetically, the Spanish population was significantly differentiated from the other populations, but the divergence was much more pronounced at mtDNA compared to microsatellites. This may reflect increased rate of differentiation by genetic drift at the mitochondrial, compared with the nuclear genome, caused by the smaller effective population size of the former genome. In accordance with this interpretation, a genetic pattern consistent with effects of small population size in the Spanish population (genetic drift and inbreeding) were also apparent at the microsatellites, namely reduced allelic diversity and heterozygous deficiency.     

Differences in male and female macaque dispersal lead to contrasting distributions of nuclear and mitochondrial DNA variation

Male macaques typically leave their natal group before sexual maturity, while females remain for life. Thus genes flow between groups and populations almost solely through male transfer. This asymmetrical dispersal pattern, affects the distribution of variation in the nuclear and mitochondrial genomes differently Nuclear genetic variation, measured by allozyme polymorphisms, is relatively evenly distributed throughout the populations of a macaque species, provided there are no major geographical barriers. Conversely the distribution of maternally inherited mitochondrial DNA (mtDNA) diversity is characterized by local homogeneity and large interpopulational differences. Because of differences in inheritance, dispersal, and population structure, the information contained in nuclear and mitochondrial genomes is best used to address different types of behavioral, genetic, and conservation questions.     

PREDICTING NUCLEAR GENE COALESCENCE FROM MITOCHONDRIAL DATA: THE THREE-TIMES RULE

Abstract.— Coalescence theory predicts when genetic drift at nuclear loci will result in fixation of sequence differences to produce monophyletic gene trees. However, the theory is difficult to apply to particular taxa because it hinges on genetically effective population size, which is generally unknown. Neutral theory also predicts that evolution of monophyly will be four times slower in nuclear than in mitochondrial genes primarily because genetic drift is slower at nuclear loci. Variation in mitochondrial DNA (mtDNA) within and between species has been studied extensively, but can these mtDNA data be used to predict coalescence in nuclear loci? Comparison of neutral theories of coalescence of mitochondrial and nuclear loci suggests a simple rule of thumb. The “three‐times rule” states that, on average, most nuclear loci will be monophyletic when the branch length leading to the mtDNA sequences of a species is three times longer than the average mtDNA sequence diversity observed within that species. A test using mitochondrial and nuclear intron data from seven species of whales and dolphins suggests general agreement with predictions of the three‐times rule. We define the coalescence ratio as the mitochondrial branch length for a species divided by intraspecific mtDNA diversity. We show that species with high coalescence ratios show nuclear monophyly, whereas species with low ratios have polyphyletic nuclear gene trees. As expected, species with intermediate coalescence ratios show a variety of patterns. Especially at very high or low coalescence ratios, the three‐times rule predicts nuclear gene patterns that can help detect the action of selection. The three‐times rule may be useful as an empirical benchmark for evaluating evolutionary processes occurring at multiple loci.     

Within‐population genetic effects of mtDNA on metabolic rate in Drosophila subobscura

A growing body of research supports the view that within‐species sequence variation in the mitochondrial genome (mtDNA) is functional, in the sense that it has important phenotypic effects. However, most of this empirical foundation is based on comparisons across populations, and few studies have addressed the functional significance of mtDNA polymorphism within populations. Here, using mitonuclear introgression lines, we assess differences in whole‐organism metabolic rate of adult Drosophila subobscura fruit flies carrying either of three different sympatric mtDNA haplotypes. We document sizeable, up to 20%, differences in metabolic rate across these mtDNA haplotypes. Further, these mtDNA effects are to some extent sex specific. We found no significant nuclear or mitonuclear genetic effects on metabolic rate, consistent with a low degree of linkage disequilibrium between mitochondrial and nuclear genes within populations. The fact that mtDNA haplotype variation within a natural population affects metabolic rate, which is a key physiological trait with important effects on life‐history traits, adds weight to the emergent view that mtDNA haplotype variation is under natural selection and it revitalizes the question as to what processes act to maintain functional mtDNA polymorphism within populations.     

Contrasting levels of genetic differentiation among populations of wolverines (<Emphasis Type="Italic">Gulo gulo</Emphasis>) from northern Canada revealed by nuclear and mitochondrial loci

Habitat loss, fragmentation, overharvest, and other anthropogenic factors have resulted in population and distribution declines for North American wolverines (Gulo gulo). Currently, wolverines east of the Hudson Bay are endangered and possibly extinct, whereas the status of wolverines throughout the remaining Holarctic is vulnerable. Three previous studies using nuclear loci have detected little to no significant structuring among wolverines sampled across northern Canada. Based on these results it has been suggested that wolverines in northern Canada represent a single, panmictic population. However, as has been shown in numerous studies, in cases of female site fidelity, it is possible to have demographically autonomous populations even with male-biased gene flow. To better assess the genetic structure of wolverines in northern Canada, we examined nine microsatellite loci and DNA sequence variation from a 200 bp fragment of the mitochondrial (mtDNA) control region for 270 wolverines from nine collecting areas representing three regions of northern Canada. In agreement with previous studies, microsatellite analyses revealed a lack of significant population substructure (F _ST=0.0004). However, analysis of molecular variance, comparisons of pairwise F _ST values and nested-clade analysis of the mtDNA data revealed considerable genetic structuring among samples of wolverines from these three regions of northern Canada. These mitochondrial data provide evidence that wolverines in Canada are genetically structured due to female philopatry. The contrasting patterns of genetic differentiation based on nuclear and mitochondrial data highlight the importance of examining both nuclear and mitochondrial loci when attempting to elucidate patterns of genetic structure.     

The ghost of hybrids past: fixation of arctic charr (Salvelinus alpinus) mitochondrial DNA in an introgressed population of lake trout (S. namaycush)

Complete fixation of arctic charr ( Salvelinus alpinus ) mitochondrial DNA (mtDNA) was observed in a southern Québec population of lake trout ( S. namaycush ). This introgressed population otherwise appeared to be normal with regard to lake trout morphology and three species-diagnostic microsatellite loci. Arctic charr do not occur in the area, suggesting that the hybridization event was prehistoric. Of several possible hypotheses, the most plausible explanation for this aberrant population is that hybridization occurred in situ soon after deglaciation, with repeated backcrossing of hybrids with lake trout. Fixation of S. alpinus mtDNA in the population may have occurred either by chance (drift) or selection, although indirect evidence and data from similarly introgressed brook trout ( S. fontinalis ) populations in the region suggest that selection favouring the S. alpinus mitochondrial type and/or associated nuclear genes may have been involved.     

Near panmixia at the distribution‐wide scale but evidence of genetic differentiation in a geographically isolated population of the Terek Sandpiper Xenus cinereus

Populations from different parts of a species range may vary in their genetic structure, variation and dynamics. Geographically isolated populations or those located at the periphery of the range may differ from those located in the core of the range. Such peripheral populations may harbour genetic variation important for the adaptive potential of the species. We studied the distribution‐wide population genetic structure of the Terek Sandpiper Xenus cinereus using 13 microsatellite loci and the mitochondrial DNA (mtDNA) control region. In addition, we estimated whether genetic variation changes from the core towards the edge of the breeding range. We used the results to evaluate the management needs of the sampled populations. Distribution‐wide genetic structure was negligible; the only population that showed significant genetic differentiation was the geographically isolated Dnieper River basin population in Eastern Europe. The genetic variation of microsatellites decreased towards the edge of the distribution, supporting the abundant‐centre hypotheses in which the core area of the distribution preserves the most genetic variation; however, no such trend could be seen with mtDNA. Overall genetic variation was low and there were signs of past population contractions followed by expansion; this pattern is found in most northern waders. The current effective population size (N_e) is large, and therefore global conservation measures are not necessary. However, the marginal Dnieper River population needs to be considered its own management unit. In addition, the Finnish population warrants conservation actions due to its extremely small size and degree of isolation from the main range, which makes it vulnerable to genetic depletion.     

Microsatellite analysis of the phylogeography, Pleistocene history and secondary contact hypotheses for the killifish, Fundulus heteroclitus

The mummichog, Fundulus heteroclitus, exhibits extensive latitudinal clinal variation in a number of physiological and biochemical traits, coupled with phylogeographical patterns at mitochondrial and nuclear DNA loci that suggest a complicated history of spatially variable selection and secondary intergradation. This species continues to serve as a model for understanding local and regional adaptation to variable environments. Resolving the influences of historical processes on the distribution of genetic variation within and among extant populations of F. heteroclitus is crucial to a better understanding of how populations evolve in the context of contemporary environments. In this study, we analysed geographical patterns of genetic variation at eight microsatellite loci among 15 populations of F. heteroclitus distributed throughout the North American range of the species from Nova Scotia to Georgia. Genetic variation in Northern populations was lower than in Southern populations and was strongly correlated with latitude throughout the species range. The most common Northern alleles at all eight loci exhibited concordant latitudinal clinal patterns, and the existence of an abrupt transition zone in allele frequencies between Northern and Southern populations was similar to that observed for mitochondrial DNA and allozyme loci. A significant pattern of isolation by distance was observed both within and between northern and southern regions. This pattern was unexpected, particularly for northern populations, given the recent colonization history of post-Pleistocene habitats, and was inconsistent with either a recent northward population expansion or a geographically restricted northern Pleistocene refugium. The data provided no evidence for recent population bottlenecks, and estimates of historical effective population sizes suggest that post-Pleistocene populations have been large throughout the species distribution. These results suggest that F. heteroclitus was broadly distributed throughout most of its current range during the last glacial event and that the abrupt transition in allele frequencies that separate Northern and Southern populations may reflect regional disequilibrium conditions associated with the post-Pleistocene colonization history of habitats in that region.     

Molecular signatures of divergence and selection in closely related pine taxa

Efforts to detect loci under selection in plants have mostly focussed on single species. However, assuming that intraspecific divergence may lead to speciation, comparisons of genetic variation within and among recently diverged taxa can help to locate such genes. In this study, coalescent and outlier detection methods were used to assess nucleotide polymorphism and divergence at 79 nuclear gene fragments (1212 SNPs) in 16 populations (153 individuals) of the closely related, but phenotypically and ecologically distinct, pine taxa Pinus mugo, P. uliginosa and P. uncinata across their European distributions. Simultaneously, mitochondrial DNA markers, which are maternally inherited in pines and distributed by seeds at short geographic distance, were used to assess genetic relationships of the focal populations and taxa. The majority of nuclear loci showed homogenous patterns of variation between the taxa due to a high number of shared SNPs and haplotypes, similar levels of polymorphism, and low net divergence. However, against this common genetic background and an overall low population structure within taxa at mitochondrial markers, we identified several genes showing signatures of selection, accompanied by significant intra- and interspecific divergence. Our results indicate that loci involved in species divergence may be involved in intraspecific local adaptation.     

Population genetic structure of the African elephant in Uganda based on variation at mitochondrial and nuclear loci: evidence for male-biased gene flow

A drastic decline has occurred in the size of the Uganda elephant population in the last 40 years, exacerbated by two main factors; an increase in the size of the human population and poaching for ivory. One of the attendant consequences of such a decline is a reduction in the amount of genetic diversity in the surviving populations due to increased effects of random genetic drift. Information about the amount of genetic variation within and between the remaining populations is vital for their future conservation and management. The genetic structure of the African elephant in Uganda was examined using nucleotide variation of mitochondrial control region sequences and four nuclear microsatellite loci in 72 individuals from three localities. Eleven mitochondrial DNA (mtDNA) haplotypes were observed, nine of which were geographically localized. We found significant genetic differentiation between the three populations at the mitochondrial locus while three out of the four microsatellite loci differentiated KV and QE, one locus differentiated KV and MF and no loci differentiated MF and QE. Expected heterozygosity at the four loci varied between 0.51 and 0.84 while nucleotide diversity at the mitochondrial locus was 1.4%. Incongruent patterns of genetic variation within and between populations were revealed by the two genetic systems, and we have explained these in terms of the differences in the effective population sizes of the two genomes and male-biased gene flow between populations.