From the Old World to the New World: a molecular chronicle of the phylogeny and biogeography of hystricognath rodents

The spider genus Hypochilus is currently restricted to cool, moist microhabitats in three widely separated montane regions of North America, providing an opportunity to study both deep (i.e., continental level) and shallow (within montane region) biogeographic history. Members of the genus also retain many plesiomorphic morphological characteristics, inviting the study of comparative rates of morphological evolution. In this paper, Hypochilus phylogeny and associated evolutionary problems are addressed using both new molecular (28S nDNA and CO1 mtDNA) and previously published (K. M. Catley, 1994, Am. Mus. Nov. 3088, 1-27) morphological data. Although the molecular data provide limited resolution of root placement within Hypochilus, most analyses are at least consistent with morphology-supported montane relationships of (Rockies (California, Appalachian)). The monophyly of Hypochilus species distributed in the California mountains is ambiguous, with several analyses indicating that this fauna may be paraphyletic with respect to a monophyletic Appalachian lineage. The montane regions differ in consistent ways in depths of both mitochondrial and nuclear phylogenetic divergence. Molecular clock analyses, in combination with arthropod-based mtDNA rate calibrations, suggest that the regional faunas are of different ages and that speciation in all faunas likely occurred prior to the Pleistocene. Limited intraspecific sampling reveals extraordinarily high levels of mtDNA cytochrome oxidase sequence divergence. These extreme divergences are most consistent with morphological stasis at the species level, despite preliminary evidence that Hypochilus taxa are characterized by fragmented population structures.     

Evidence for long-distance dispersal in a sedentary passerine, Gymnorhina tibicen (Artamidae)

Australian magpies (Gymnorhina tibicen) are group-living birds found across much of mainland Australia. Adults commonly remain in a breeding territory until death. Young of the year either remain on the natal (birth) site or are forced by their parents to disperse. Observational studies in south-eastern Australia suggest that most dispersing juveniles settle within 7 km of their natal territory. Therefore, despite potential for considerable gene flow (via flight), social organization predisposes magpies towards local population structuring. In this study, we measured genetic variation at both nuclear (allozyme) and mitochondrial loci and found evidence of substantial gene flow over very large distances (up to 1599 km). Thus, some juvenile magpies may disperse much greater distances than was previously thought. For mtDNA, geographic and genetic distance were strongly correlated, consistent with a pattern of isolation by distance. Therefore, although female gene flow is substantial it is apparently geographically restricted over large distances, in approximately a stepping-stone fashion. We conclude that a strong relationship between gene flow and geographic distance can develop even over large distances if populations have experienced no major historical disturbances to gene flow.     

Inferring contemporary levels of gene flow and demographic history in a local population of the leaf beetle Gonioctena olivacea from mitochondrial DNA sequence variation

We have studied mitochondrial DNA variation in a local population of the leaf beetle species Gonioctena olivacea, to check whether its apparent low dispersal behaviour affects its pattern of genetic variation at a small geographical scale. We have sampled 10 populations of G. olivacea within a rectangle of 5 x 2 km in the Belgian Ardennes, as well as five populations located approximately along a straight line of 30 km and separated by distances of 3-12 km. For each sampled individual (8-19 per population), a fragment of the mtDNA control region was polymerase chain reaction-amplified and sequenced. Sequence data were analysed to test whether significant genetic differentiation could be detected among populations separated by such relatively short distances. The reconstructed genealogy of the mitochondrial haplotypes was also used to investigate the demographic history of these populations. Computer simulations of the evolution of populations were conducted to assess the minimum amount of gene flow that is necessary to explain the observed pattern of variation in the samples. Results show that migration among populations included in the rectangle of 5 x 2 km is substantial, and probably involves the occurrence of dispersal flights. This appears difficult to reconcile with the results of a previous ecological field study that concluded that most of this species dispersal occurs by walking. While sufficient migration to homogenize genetic diversity occurs among populations separated by distances of a few hundred metres to a few kilometres, distances greater than 5 km results in contrast in strong differentiation among populations, suggesting that migration is drastically reduced on such distances. Finally, the results of coalescent simulations suggest that the star-like genealogy inferred from the mtDNA sequence data is fully compatible with a past demographic expansion. However, a metapopulation structure alone (without the need to invoke a population expansion event) cannot be dismissed as the cause of this star shape.     

Patterns of population subdivision, gene flow and genetic variability in the African wild dog (Lycaon pictus)

African wild dogs are large, highly mobile carnivores that are known to disperse over considerable distances and are rare throughout much of their geographical range. Consequently, genetic variation within and differentiation between geographically separated populations is predicted to be minimal. We determined the genetic diversity of mitochondrial DNA (mtDNA) control region sequences and microsatellite loci in seven populations of African wild dogs. Analysis of mtDNA nucleotide diversity suggests that, historically, wild dog populations have been small relative to other large carnivores. However, population declines due to recent habitat loss have not caused a dramatic reduction in genetic diversity. We found one historical and eight recent mtDNA genotypes in 280 individuals that defined two highly divergent clades. In contrast to a previous, more limited, mtDNA analysis, sequences from these clades are not geographically restricted to eastern or southern African populations. Rather, we found a large admixture zone spanning populations from Botswana, Zimbabwe and south-eastern Tanzania. Mitochondrial and microsatellite differentiation between populations was significant and unique mtDNA genotypes and alleles characterized the populations. However, gene flow estimates (Nm) based on microsatellite data were generally greater than one migrant per generation. In contrast, gene flow estimates based on the mtDNA control region were lower than expected given differences in the mode of inheritance of mitochondrial and nuclear markers which suggests a male bias in long-distance dispersal.     

Mitochondrial genealogy of Western Australian barramundi: applications of inbreeding coefficients and coalescent analysis for separating temporal population processes

In Australian populations of barramundi Lates calcarifer , phylogenetic reconstruction of mtDNA sequences provided evidence of significant historical levels of gene flow, despite the substantial structuring of contemporary populations. The geographical pattern of mtDNA sequences among the populations was not congruent with previous evidence of a major disjunction between western and eastern populations of barramundi in Australia.     

Genetic consequences of the ice ages on nurseries of the bat Myotis myotis: a mitochondrial and nuclear survey

Analyses of mitochondrial DNA (mtDNA) control region polymorphism and of variation at 10 nuclear microsatellite loci were used to investigate the mechanisms and genetic consequences of postglacial expansion of Myotis myotis in Europe. Initial sampling consisted of 480 bats genotyped in 24 nursery colonies arranged along a transect of approximately 3000 km. The phylogeographical survey based on mtDNA sequences revealed the existence of major genetic subdivisions across this area, with several suture zones between haplogroups. Such zones of secondary contact were found in the Alps and Rhodopes, whereas other potential barriers to gene flow, like the Pyrenees, did not coincide with genetic discontinuities. Areas of population admixture increased locally the genetic diversity of colonies, which confounded the northward decrease in nucleotide diversity predicted using classical models of postglacial range expansion. However, when analyses were restricted to a subset of 15 nurseries originating from a single presumed glacial refugium, mtDNA polymorphism did indeed support a northwards decrease in diversity. Populations were also highly structured (PhiST = 0.384). Conversely, the same subset of colonies showed no significant latitudinal decrease in microsatellite diversity and much less population structure (FST = 0.010), but pairwise genetic differentiation at these nuclear markers was strongly correlated with increasing geographical distance. Together, this evidence suggests that alleles carried via male bats have maintained enough nuclear gene flow to counteract the effects of recurrent bottlenecks generally associated with recolonization processes. As females are highly philopatric, we argue that the maternally transmitted mtDNA marker better reflects the situation of past, historical gene flow, whereas current levels of gene flow are better reflected by microsatellite markers.     

Influence of habitat discontinuity, geographical distance, and oceanography on fine-scale population genetic structure of copper rockfish (Sebastes caurinus)

The copper rockfish is a benthic, nonmigratory, temperate rocky reef marine species with pelagic larvae and juveniles. A previous range-wide study of the population-genetic structure of copper rockfish revealed a pattern consistent with isolation-by-distance. This could arise from an intrinsically limited dispersal capability in the species or from regularly-spaced extrinsic barriers that restrict gene flow (offshore jets that advect larvae offshore and/or habitat patchiness). Tissue samples were collected along the West Coast of the contiguous USA between Neah Bay, WA and San Diego, CA, with dense sampling along Oregon. At the whole-coast scale (approximately 2200 km), significant population subdivision (F(ST) = 0.0042), and a significant correlation between genetic and geographical distance were observed based on 11 microsatellite DNA loci. Population divergence was also significant among Oregon collections (approximately 450 km, F(ST) = 0.001). Hierarchical amova identified a weak but significant 130-km habitat break as a possible barrier to gene flow within Oregon, across which we estimated that dispersal (N(e)m) is half that of the coast-wide average. However, individual-based Bayesian analyses failed to identify more than a single population along the Oregon coast. In addition, no correlation between pairwise population genetic and geographical distances was detected at this scale. The offshore jet at Cape Blanco was not a significant barrier to gene flow in this species. These findings are consistent with low larval dispersal distances calculated in previous studies on this species, support a mesoscale dispersal model, and highlight the importance of continuity of habitat and adult population size in maintaining gene flow.     

Islands within an island: phylogeography and conservation genetics of the endangered Hawaiian tree snail Achatinella mustelina

Mitochondrial DNA (mtDNA) sequences were used to evaluate phylogeographic structure within and among populations of three endangered Hawaiian tree snail species (n = 86). The primary focus of this investigation was on setting conservation priorities for Achatinella mustelina. Limited data sets for two additional endangered Hawaiian tree snails, A. livida and A. sowerbyana, were also developed for comparative purposes. Pairwise genetic distance matrices and phylogenetic trees were generated, and an analysis of molecular variance was performed on 675-base pair cytochrome oxidase I gene sequences from multiple populations of Hawaiian tree snails. Sequence data were analysed under distance-based maximum-likelihood, and maximum-parsimony optimality criteria. Within the focal species, A. mustelina, numbers of variable and parsimony informative sites were 90 and 69, respectively. Pairwise intraspecific mtDNA sequence divergence ranged from 0 to 5.3% in A. mustelina, from 0 to 1.0% in A. livida and from 0 to 1.9% in A. sowerbyana. For A. mustelina, population genetic structure and mountain topography were strongly correlated. Maximum genetic distances were observed across deep, largely deforested valleys, and steep mountain peaks, independent of geographical distance. However, in certain areas where forest cover is presently fragmented, little mtDNA sequence divergence exists despite large geographical scales (8 km). Genetic data were used to define evolutionarily significant units for conservation purposes including decisions regarding placement of predator exclusion fences, captive propagation, re-introduction and translocation.     

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.     

Population structure in two sympatric species of the Lake Tanganyika cichlid tribe Eretmodini: evidence for introgression

Patterns of genetic differentiation were analysed and compared in two sympatric species of the endemic Lake Tanganyika cichlid tribe Eretmodini by means of mitochondrial DNA (mtDNA) sequences of the control region and six microsatellite DNA loci. The sample area covers a total of 138 km of mostly uninterrupted rocky shoreline in the Democratic Republic of Congo and includes the entire distribution range of Tanganicodus cf. irsacae that stretches over a distance of 35 km. Both markers detected significant genetic differentiation within and between the two species. T. cf. irsacae contained lower overall genetic variation than Eretmoduscyanostictus, possibly due to its more restricted range of distribution and its smaller effective population sizes. Complete fixation of Tanganicodus mtDNA haplotypes was observed in Eretmodus at two localities, while at two other localities some Tanganicodus individuals possessed Eretmodus mtDNA haplotypes. Taking into account the relatively large average sequence divergence of 6.2% between the two species, as well as the geographical distribution of mtDNA haplotypes in the lake, the observed pattern is more likely to be a consequence of asymmetric introgression than of shared ancestral polymorphism. As there is significant population differentiation between sympatric Tanganicodus and Eretmodus populations, the events of introgressions may have happened after secondary contact, but our data provide no evidence for ongoing gene flow and suggest that both species are reproductively isolated at present time.     

Indo-West Pacific patterns of genetic differentiation in the high-dispersal starfish Linckia laevigata

Genetic variation in four natural populations of the starfish Linckia laevigata from the Indo-West Pacific was examined using restriction fragment analysis of a portion of the mtDNA including the control region. Digestion with seven restriction enzymes identified 47 haplotypes in a sample of 326 individuals. Samples collected from reef sites within each location were not significantly differentiated based on Φ_ST or spatial distribution of haplotypes, indicating that dispersal is high over short to moderate distances. Evidence of gene flow is further supported by the low divergence among haplotypes and the lack of any clear geographical structuring among different haplotypes in the gene phylogeny. However, analysis of molecular variance ( AMOVA ), Φ_ST and contingency χ² analyses of the spatial distribution of haplotypes demonstrate the presence of significant broad scale population genetic structure among the four widespread locations examined. RFLP data are consistent with high gene flow between the Philippines and Western Australia and moderate gene flow between the Great Barrier Reef (GBR) and Fiji, but only limited gene flow between either the Philippines or Western Australia and either the GBR or Fiji. The presence of mtDNA structure contrasts with previous allozyme data which suggest that dispersal among widely separated locations is equivalent to dispersal among populations within the highly connected GBR studies. This discordance between patterns of gene flow inferred from these two markers cannot be fully accounted for by differences in effective population size for mtDNA. This might suggest that while mtDNA variation may represent contemporary patterns of gene flow, allozyme variation among populations is yet to reach equilibrium between drift and migration over the range surveyed.     

Hierarchical genetic structure and gene flow in three sympatric species of Amazonian rodents

The population genetic structure of three species of Amazonian rodents ( Oligoryzomys microtis, Oryzomys capito , and Mesomys hispidus ) is examined for mtDNA sequence haplotypes of the cytochrome b gene by hierarchical analysis of variance and gene flow estimates based on fixation indices ( N _ST) and coalescence methods. Species samples are from the same localities along 1000 km of the Rio Juruá in western Amazonian Brazil, but each species differs in important life history traits such as population size and reproductive rate. Average haplotype differentiation, hierarchical haplotype apportionment, and gene flow estimates are contrasted in discussing the current and past population structure. Two species exhibit isolation by distance patterns wherein gene flow is largely limited to geographically adjacent localities. Mesomys exhibits this pattern throughout its range along the river. More than 75% of haplotype variation is apportioned among localities and regions, and estimates of Nm for pair-wise comparisons are nearly always less than 1. Oligoryzomys shows weak isolation by distance, but only over the largest geographical distances. Nm values for this species are nearly always above 1 and most (about 80%) of haplotype variation is contained within local populations. In contrast, Oryzomys exhibits no genetic structure throughout its entire distribution; Nm values average 17 and nearly 90% of the total haplotype variance is contained within local populations. Although gene flow estimates are high, the pattern of Nm as a function of geographical distance suggests that this species experienced a more recent invasion of the region and is still in genetic disequilibrium under its current demographic conditions.     

Population genetic structure and male-biased dispersal in the queenless ant Diacamma cyaneiventre

In this study we investigated the population genetic structure of the queenless ant Diacamma cyaneiventre. This species, lacking winged queens, is likely to have a restricted female dispersal. We used both mitochondrial and microsatellite markers to assess the consequence of such restricted female dispersal at three geographical scales: within a given locality (< 1 km), between localities within a given region (< 10 km) and between regions (> 36 km). Within a locality, a strong population structure was observed for mitochondrial DNA (mtDNA) whereas weak or nonexistent population genetic structure was observed for the microsatellites (around 5% of the value for mtDNA). Male gene flow was estimated to be about 20-30 times higher than female gene flow at this scale. At a larger spatial scale, very strong genetic differentiation for both markers was observed between localities - even within a single region. Female dispersal is nonexistent at these scales and male dispersal is very restricted, especially between regions. The phylogeographical structure of the mtDNA haplotypes as well as the very low genetic diversity of mtDNA within localities indicate that new sites are colonized by a single migration event from adjacent localities, followed by successive colony fissions. These patterns of genetic variability and differentiation agree with what is theoretically expected when colonization events are kin-structured and when, following colonization, dispersion is mainly performed by males.     

Genetic and ecological data provide incongruent interpretations of population structure and dispersal in naturally subdivided populations of white-tailed ptarmigan (Lagopus leucura)

The dispersal of individuals among populations affects the demographic and adaptive trajectories of animal populations and is fundamental to understanding population dynamics. White-tailed ptarmigan (Lagopus leucura) are a high elevation grouse species that live year-round in patchily distributed alpine areas in western North America. We investigated the patterns of dispersal and identified barriers to gene flow for a threatened subspecies (L. l. saxatilis) endemic to Vancouver Island, Canada. Connectivity among seven sites was examined using nine microsatellite loci (n = 133 individuals, H(O) = 0.62, mean number of alleles = 10) and direct movement observations using radio-telemetry (n = 118 individuals). Average movement distances of individuals measured by radio-telemetry were 0.63-3.23 km and considerably less than the shortest distance between sampling sites (18 km). Furthermore, despite extensive radio-telemetry data, movement was never observed between any of the seven sampling sites. In contrast, genetic results (STRUCTURE, TESS) showed connectivity among most of the seven sampling sites and suggested that genetic variation is best explained by two clusters of individuals which separated the South sampling site from all other areas of Vancouver Island. Analysis of molecular data also showed a generally consistent pattern of isolation by distance (Mantel test r = 0.11, P < 0.01) with large areas of unsuitable low elevation habitat possibly acting as barriers to gene flow. Despite the naturally subdivided distribution of populations, white-tailed ptarmigan do not fit well into any common definition of a metapopulation. We conclude the incongruities between the genetic and radio-telemetry data are best explained by episodic dispersal patterns. In this study, we demonstrated the importance of combining genetic and ecological data in understanding patterns of dispersal and population structure.     

Population differentiation and nuclear gene flow in the Dominican anole (Anolis oculatus)

Allele frequency data from nuclear microsatellite loci were used to investigate patterns of nuclear gene flow and population structure in the morphologically variable Dominican anole (Anolis oculatus). All six loci used proved to be highly polymorphic, with an average of 18.8 alleles per locus. Test for Hardy-Weinberg equilibrium revealed small numbers of heterozygote deficiencies at single loci in single populations and consistent patterns of increasingly significant heterozygote deficiency in global tests across populations and loci. No significant relationship between FST and patristic distances estimated from mitochondrial DNA sequences was detected and estimates of FIS were significantly higher in females than in males, indicating that gene flow may be sex-biased and mediated mainly by male migration. A highly significant correlation between linearized FST and loge (geographical distance) indicates that geographical proximity is a significant factor in the genetic structure of A. oculatus populations. However, levels of gene flow between morphologically differentiated parapatric populations are frequently seen to be relatively high. This supports the hypothesis of natural selection being the driving force behind the development and maintenance of morphological variation and shows that adaptive differentiation may be maintained despite the homogenizing influence of gene flow. Generally, the morphologically variable populations of A. oculatus are seen to be poor candidates for in situ speciation, but an exceptional case on the west coast of Dominica indicates that isolation resulting from vicariant events may lead to rapid differentiation at both mitochondrial and nuclear loci. This provides a possible mechanism for anole speciation on other Caribbean islands.     

Spatial structure of lemming populations (Dicrostonyx groenlandicus) fluctuating in density

The pattern and scale of the genetic structure of populations provides valuable information for the understanding of the spatial ecology of populations, including the spatial aspects of density fluctuations. In the present paper, the genetic structure of periodically fluctuating lemmings (Dicrostonyx groenlandicus) in the Canadian Arctic was analysed using mitochondrial DNA (mtDNA) control region sequences and four nuclear microsatellite loci. Low genetic variability was found in mtDNA, while microsatellite loci were highly variable in all localities, including localities on isolated small islands. For both genetic markers the genetic differentiation was clear among geographical regions but weaker among localities within regions. Such a pattern implies gene flow within regions. Based on theoretical calculations and population census data from a snap-trapping survey, we argue that the observed genetic variability on small islands and the low level of differentiation among these islands cannot be explained without invoking long distance dispersal of lemmings over the sea ice. Such dispersal is unlikely to occur only during population density peaks.     

Land, language, and loci: mtDNA in Native Americans and the genetic history of Peru

Despite a long history of complex societies and despite extensive present-day linguistic and ethnic diversity, relatively few populations in Peru have been sampled for population genetic investigations. In order to address questions about the relationships between South American populations and about the extent of correlation between genetic distance, language, and geography in the region, mitochondrial DNA (mtDNA) hypervariable region I sequences and mtDNA haplogroup markers were examined in 33 individuals from the state of Ancash, Peru. These sequences were compared to those from 19 American Indian populations using diversity estimates, AMOVA tests, mismatch distributions, a multidimensional scaling plot, and regressions. The results show correlations between genetics, linguistics, and geographical affinities, with stronger correlations between genetics and language. Additionally, the results suggest a pattern of differential gene flow and drift in western vs. eastern South America, supporting previous mtDNA and Y chromosome investigations.     

Fine scale gene flow and individual movements among subpopulations of Centrolene prosoblepon (Anura: Centrolenidae)

Dispersal capabilities determine and maintain local gene flow, and this has implications for population persistence and/or recolonization following environmental perturbations (natural or anthropogenic), disease outbreaks, or other demographic collapses. To predict recolonization and understand dispersal capacity in a stream-breeding frog, we examined individual movement patterns and gene flow among four subpopulations of the Neotropical glassfrog, Centrolene prosoblepon, at a mid-elevation cloud forest site at El Copé, Panama. We measured male movement directly during a two year mark-recapture study, and indirectly with gene flow estimates from mitochondrial DNA sequences (mtDNA). Individuals of this species showed strong site fidelity: over two years, male frogs in all four headwater streams moved very little (mean = 2.33 m; mode = 0 m). Nine individuals changed streams within one or two years, moving 675-1,108 m. For those males moving more than 10 m, movement was biased upstream (p < 0.001). Using mtDNA ND1 gene sequences, we quantified gene flow within and among headwater streams at two spatial scales: among headwater streams within two adjacent watersheds (2.5 km2) and among streams within a longitudinal gradient covering 5.0 km2. We found high gene flow among headwater streams (phi(ST) = 0.007, p = 0.325) but gene flow was more limited across greater distances (phi(CT) = 0.322, p = 0.065), even within the same drainage network. Lowland populations of C. prosoblepon potentially act as an important source of colonists for upland populations in this watershed.     

Population structure of the spotted babylon, Babylonia areolata in three wild populations along the Chinese coastline revealed using AFLP fingerprinting

Many marine gastropods are sedentary as adults but have planktonic larvae which can potentially be dispersed over large distances. Consequently larval transport is expected to play a prominent role in facilitating gene flow and determining population structure. The spotted babylon (Babylonia areolata) is a dioecious species possessing an approximately two week planktonic larval stage. We analyzed the population structure of the spotted babylon using amplified fragment length polymorphism (AFLP). One hundred and sixteen AFLP loci were analyzed in 63 individuals from three populations and revealed a high level of genetic diversity, with all individuals harboring a unique banding pattern. AMOVA results and an assignment test revealed that population differentiation was present. PCoA, pairwise F_ST and UPGMA tree all revealed that gene flow might be present only on a small geographic scale (around 160 km), but, over a large distance (around 1000 km), only reduced gene flow occurred. A mantel test indicated a highly significant positive correlation between genetic differentiation and geographical distance.     

Great bustard population structure in central Spain: concordant results from genetic analysis and dispersal study.

We found significant sex differences in the mtDNA genetic structure and dispersal patterns of great bustards in a population of 11 breeding groups, "leks", in central Spain. The analysis of genetic distances showed that the female population was divided into three groups of leks separated by ca. 50 km, whereas male haplotypes were randomly distributed among leks. Genetic distances among pairs of leks were positively correlated with geographical distances in females but not in males. While female haplotype distributions were homogeneous among leks at close distances, differences in male genetic structure were highly variable even between two close leks. These results from genetic analyses were concordant with those from a radiotracking study on natal dispersal. Natal dispersal distances were higher in males than in females. Also, the frequency of movement of a female between two leks was positively correlated with their genetic affinity and geographical proximity. In males, the frequency of movement was correlated with geographical proximity but not with genetic affinity. Males dispersed among genetically unrelated leks, contributing to keep nuclear genetic diversity in the population, whereas females tended to be philopatric. These results suggest that isolation-by-distance influences the distribution of maternal lineages at a regional level.