Hybridization between multiple fence lizard lineages in an ecotone: locally discordant variation in mitochondrial DNA, chromosomes, and morphology

We investigated a hybrid zone between two major lineages of fence lizards (Sceloporus cowlesi and Sceloporus tristichus) in the Sceloporus undulatus species complex in eastern Arizona. This zone occurs in an ecotone between Great Basin Grassland and Conifer Woodland habitats. We analysed spatial variation in mtDNA (N=401; 969 bp), chromosomes (N=217), and morphology (N=312; 11 characters) to characterize the hybrid zone and assess species limits. A fine-scale population level phylogenetic analysis refined the boundaries between these species and indicated that four nonsister mtDNA clades (three belonging to S. tristichus and one to S. cowlesi) are sympatric at the centre of the zone. Estimates of cytonuclear disequilibria in the population closest to the centre of the hybrid zone suggest that the S. tristichus clades are randomly mating, but that the S. cowlesi haplotype has a significant nonrandom association with nuclear alleles. Maximum-likelihood cline-fitting analyses suggest that the karyotype, morphology, and dorsal colour pattern clines are all coincident, but the mtDNA cline is skewed significantly to the south. A temporal comparison of cline centres utilizing karyotype data collected in the early 1970s and in 2002 suggests that the cline may have shifted by approximately 1.5 km to the north over a 30-year period. The recent northward expansion of juniper trees into the Little Colorado River Basin resulting from intense cattle overgrazing provides a plausible mechanism for a shifting hybrid zone and the introgression of the mtDNA haplotypes, which appear to be selectively neutral. It is clear that complex interactions are operating simultaneously in this contact zone, including the formation of hybrids between populations within S. tristichus having diagnostic mtDNA, morphology, karyotypes, and dorsal colour patterns, and secondary contact between these and a distantly related yet morphologically cryptic mtDNA lineage (S. cowlesi).     

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.     

Is Gene Flow Promoting the Reversal of Pleistocene Divergence in the Mountain Chickadee (Poecile gambeli)?

The Pleistocene glacial cycles left a genetic legacy on taxa throughout the world; however, the persistence of genetic lineages that diverged during these cycles is dependent upon levels of gene flow and introgression. The consequences of secondary contact among taxa may reveal new insights into the history of the Pleistocene’s genetic legacy. Here, we use phylogeographic methods, using 20 nuclear loci from regional populations, to infer the consequences of secondary contact following divergence in the Mountain Chickadee (Poecile gambeli). Analysis of nuclear data identified two geographically-structured genetic groups, largely concordant with results from a previous mitochondrial DNA (mtDNA) study. Additionally, the estimated multilocus divergence times indicate a Pleistocene divergence, and are highly concordant with mtDNA. The previous mtDNA study showed a paucity of sympatry between clades, while nuclear patterns of gene flow show highly varied patterns between populations. The observed pattern of gene flow, from coalescent-based analyses, indicates southern populations in both clades exhibit little gene flow within or between clades, while northern populations are experiencing higher gene flow within and between clades. If this pattern were to persist, it is possible the historical legacy of Pleistocene divergence may be preserved in the southern populations only, and the northern populations would become a genetically diverse hybrid species.     

Mito‐nuclear discord in six congeneric lineages of Holarctic ducks (genus Anas)

Many species have Holarctic distributions that extend across Europe, Asia and North America. Most genetics research on these species has examined only mitochondrial (mt) DNA, which has revealed wide variance in divergence between Old World (OW) and New World (NW) populations, ranging from shallow, unstructured genealogies to deeply divergent lineages. In this study, we sequenced 20 nuclear introns to test for concordant patterns of OW–NW differentiation between mtDNA and nuclear (nu) DNA for six lineages of Holarctic ducks (genus Anas). Genetic differentiation for both marker types varied widely among these lineages (idiosyncratic population histories), but mtDNA and nuDNA divergence within lineages was not significantly correlated. Moreover, compared with the association between mtDNA and nuDNA divergence observed among different species, OW–NW nuDNA differentiation was generally lower than mtDNA divergence, at least for lineages with deeply divergent mtDNA. Furthermore, coalescent estimates indicated significantly higher rates of gene flow for nuDNA than mtDNA for four of the six lineages. Thus, Holarctic ducks show prominent mito‐nuclear discord between OW and NW populations, and we reject differences in sorting rates as the sole cause of the within‐species discord. Male‐mediated intercontinental gene flow is likely a leading contributor to this discord, although selection could also cause increased mtDNA divergence relative to weak nuDNA differentiation. The population genetics of these ducks contribute to growing evidence that mtDNA can be an unreliable indicator of stage of speciation and that more holistic approaches are needed for species delimitation.     

Differential patterns of male and female mtDNA exchange across the Atlantic Ocean in the blue mussel, Mytilus edulis

Comparisons among loci with differing modes of inheritance can reveal unexpected aspects of population history. We employ a multilocus approach to ask whether two types of independently assorting mitochondrial DNAs (maternally and paternally inherited: F- and M-mtDNA) and a nuclear locus (ITS) yield concordant estimates of gene flow and population divergence. The blue mussel, Mytilus edulis, is distributed on both North American and European coastlines and these populations are separated by the waters of the Atlantic Ocean. Gene flow across the Atlantic Ocean differs among loci, with F-mtDNA and ITS showing an imprint of some genetic interchange and M-mtDNA showing no evidence for gene flow. Gene flow of F-mtDNA and ITS causes trans-Atlantic population divergence times to be greatly underestimated for these loci, although a single trans-Atlantic population divergence time (1.2 MYA) can be accommodated by considering all three loci in combination in a coalescent framework. The apparent lack of gene flow for M-mtDNA is not readily explained by different dispersal capacities of male and female mussels. A genetic barrier to M-mtDNA exchange between North American and European mussel populations is likely to explain the observed pattern, perhaps associated with the double uniparental system of mitochondrial DNA inheritance.     

Genetic, morphological, and ecological characterization of a hybrid zone that spans a migratory divide

This study characterizes a hybrid zone that spans a migratory divide between subspecies of the Swainson's thrush (Catharus ustulatus), a long distance migratory songbird, in the Coast Mountains of British Columbia. To assess the potential for a barrier to gene flow between the subspecies, I: (1) analyzed the shape and width of genetic and morphological clines relative to estimates of dispersal distance, (2) assessed the ratio of parental to hybrid genotypes across the hybrid zone, (3) estimated population density across the hybrid zone, and (4) compared the spatial relationship between the hybrid zone and an existing environmental gradient. The results indicate that the hybrid zone is characterized by mostly concordant character clines that are narrow relative to dispersal, the absence of a hybrid swarm, and low population density at the center of the zone. This hybrid zone and additional regions of contact between these subspecies are found on the border between coastal and interior climatic regions throughout the Pacific Northwest. An identified shift in the location, but not the width, of the mtDNA cline relative to the nuclear clines is consistent with asymmetrical hybridization. Neutral diffusion of populations following secondary contact and hybrid superiority within an ecotone are insufficient explanations for the observed patterns. The hypothesis that best fits the data is that the Swainson's thrush hybrid zone is a tension zone maintained by dispersal and ecologically mediated barriers to gene flow.     

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.     

Assessment of gene flow across a hybrid zone in red-tailed chipmunks (Tamias ruficaudus)

The role of hybridization in animal speciation is controversial and recent research has challenged the long-standing criterion of complete reproductive isolation to attain species status. The speciation-with-gene-flow model posits that the genome is semi-permeable and hybridization may be a phase in the process of divergence. Here, we apply these concepts to a previously identified zone of mtDNA introgression between the two strongly morphologically differentiated subspecies of red-tailed chipmunk ( Tamias ruficaudus ) in the US Inland Northwest. Using multilocus genotype data from the southern, older contact zone, we demonstrate that neutral gene flow is unusually low between the subspecies across the Lochsa River. This is geographically congruent with the discontinuity in bacular morphology, indicating that the cline of mitochondrial DNA (mtDNA) haplotypes is displaced. Furthermore, we elucidate the evolutionary forces responsible by testing hypotheses of lineage sorting and hybridization. We determined that introgressive hybridization is the cause of mtDNA/morphology incongruence because there are non-zero levels of migration and gene flow. Although our estimate of the age of the hybrid zone has wide credibility intervals, the hybridization events occurred in the Late Pleistocene and the divergence occurred in the Middle Pleistocene. Finally, we assessed substructure within and adjacent to the hybrid zone and found that the hybrid zone constitutes a set of populations that are genetically differentiated from parental sets of populations; therefore, hybridization in this system is not likely an evolutionary sink, but has generated novel combinations of genotypes.     

Inter‐island divergence within Drosophila mauritiana,a species of the D. simulans complex: Past history and/or speciation in progress?

Speciation with gene flow may be more common than generally thought, which makes detailed understanding of the extent and pattern of genetic divergence between geographically isolated populations useful. Species of the Drosophila simulans complex provide a good model for speciation and evolutionary studies, and hence understanding their population genetic structure will increase our understanding of the context in which speciation has occurred. Here, we describe genetic diversity and genetic differentiation of two distant populations of D. mauritiana (Mauritius and Rodrigues Islands) at mitochondrial and nuclear loci. We surveyed the two populations for their mitochondrial haplotypes, eight nuclear genes and 18 microsatellite loci. A new mitochondrial type is fixed in the Rodrigues population of D. mauritiana. The two populations are highly differentiated, their divergence appears relatively ancient (100,000 years) compared to the origin of the species, around 0.25MYA, and they exhibit very limited gene flow. However, they have similar levels of divergence from their sibling, D. simulans. Both nuclear genes and microsatellites revealed contrasting demographic histories between the two populations, expansion for the Mauritius population and stable population size for the Rodrigues Island population. The discovery of pronounced geographic structure within D. mauritiana combined to genetic structuring and low gene flow between the two island populations illuminates the evolutionary history of the species and clearly merits further attention in the broad context of speciation.     

Huge populations and old species of Costa Rican and Panamanian dirt frogs inferred from mitochondrial and nuclear gene sequences

Molecular genetic data were used to investigate population sizes and ages of Eleutherodactylus (Anura: Leptodactylidae), a species-rich group of small leaf-litter frogs endemic to Central America. Population genetic structure and divergence was investigated for four closely related species surveyed across nine localities in Costa Rica and Panama. DNA sequence data were collected from a mitochondrial gene (ND2) and a nuclear gene (c-myc). Phylogenetic analyses yielded concordant results between loci, with reciprocal monophyly of mitochondrial DNA haplotypes for all species and of c-myc haplotypes for three of the four species. Estimates of genetic differentiation among populations (FST) based upon mitochondrial data were always higher than nuclear-based FST estimates, even after correcting for the expected fourfold lower effective population size (Ne) of the mitochondrial genome. Comparing within-population variation and the relative mutation rates of the two genes revealed that the Ne of the mitochondrial genome was 15-fold lower than the estimate of the nuclear genome based on c-myc. Nuclear FST estimates were approximately 0 for the most proximal pairs of populations, but ranged from 0.5 to 1.0 for all other pairs, even within the same nominal species. The nuclear locus yielded estimates of Ne within localities on the order of 105. This value is two to three orders of magnitude larger than any previous Ne estimate from frogs, but is nonetheless consistent with published demographic data. Applying a molecular clock model suggested that morphologically indistinguishable populations within one species may be 107 years old. These results demonstrate that even a geologically young and dynamic region of the tropics can support very old lineages that harbour great levels of genetic diversity within populations. The association of high nucleotide diversity within populations, large divergence between populations, and high species diversity is also discussed in light of neutral community models.     

Genetic differentiation and history of populations of the Italian treefrog Hyla intermedia: lack of concordance between mitochondrial and nuclear markers

The genetic differentiation among 33 populations of the Italian treefrog, Hyla intermedia (Anura: Hylidae), was investigated using both biparentally (23 allozyme loci) and maternally (partial mitochondrial cytochrome b gene) inherited markers. Two main population groups were evidenced by both markers, located north and south of the northern Apennines. However, the pattern of differentiation between these two groups was much less pronounced at allozymes than at mtDNA, leading to gene flow estimates that were 25 times lower at mitochondrial than at nuclear level. Also, the mtDNA divergence between the two groups was particularly marked for two cospecific lineages of anuran amphibians (the P-distance being on average 9.04%), while their average genetic distance at allozymes was comparatively low (D _NEI = 0.07). This contrasting pattern of nuclear versus mitochondrial genetic variation is discussed in the context of: (1) marker specific selection, (2) secondary contact and sex-biased gene flow and (3) ancestral polymorphism and colonization from north to south. Finally we emphasize how, for population genetic studies, the use of multiple markers having distinct evolutionary properties can help unravel the existence of more complex evolutionary histories.     

Inferring human population sizes, divergence times and rates of gene flow from mitochondrial, X and Y chromosome resequencing data

We estimate parameters of a general isolation-with-migration model using resequence data from mitochondrial DNA (mtDNA), the Y chromosome, and two loci on the X chromosome in samples of 25-50 individuals from each of 10 human populations. Application of a coalescent-based Markov chain Monte Carlo technique allows simultaneous inference of divergence times, rates of gene flow, as well as changes in effective population size. Results from comparisons between sub-Saharan African and Eurasian populations estimate that 1500 individuals founded the ancestral Eurasian population approximately 40 thousand years ago (KYA). Furthermore, these small Eurasian founding populations appear to have grown much more dramatically than either African or Oceanian populations. Analyses of sub-Saharan African populations provide little evidence for a history of population bottlenecks and suggest that they began diverging from one another upward of 50 KYA. We surmise that ancestral African populations had already been geographically structured prior to the founding of ancestral Eurasian populations. African populations are shown to experience low levels of mitochondrial DNA gene flow, but high levels of Y chromosome gene flow. In particular, Y chromosome gene flow appears to be asymmetric, i.e., from the Bantu-speaking population into other African populations. Conversely, mitochondrial gene flow is more extensive between non-African populations, but appears to be absent between European and Asian populations.     

World phylogeography and male‐mediated gene flow in the sandbar shark,Carcharhinus plumbeus

The sandbar shark, Carcharhinus plumbeus, is a large, cosmopolitan, coastal species. Females are thought to show philopatry to nursery grounds while males potentially migrate long distances, creating an opportunity for male‐mediated gene flow that may lead to discordance in patterns revealed by mitochondrial DNA (mtDNA) and nuclear markers. While this dynamic has been investigated in elasmobranchs over small spatial scales, it has not been examined at a global level. We examined patterns of historical phylogeography and contemporary gene flow by genotyping 329 individuals from nine locations throughout the species’ range at eight nuclear microsatellite markers and sequencing the complete mtDNA control region. Pairwise comparisons often resulted in fixation indices and divergence estimates of greater magnitude using mtDNA sequence data than microsatellite data. In addition, multiple methods of estimation suggested fewer populations based on microsatellite loci than on mtDNA sequence data. Coalescent analyses suggest divergence and restricted migration among Hawaii, Taiwan, eastern and western Australia using mtDNA sequence data and no divergence and high migration rates, between Taiwan and both Australian sites using microsatellite data. Evidence of secondary contact was detected between several localities and appears to be discreet in time rather than continuous. Collectively, these data suggest complex spatial/temporal relationships between shark populations that may feature pulses of female dispersal and more continuous male‐mediated gene flow.     

Mitochondrial DNA variation and the evolutionary history of cryptic Gammarus fossarum types

The evolutionary history of the cryptic Gammarus fossarum species complex (Crustacea, Amphipoda) in Central Europe was approached by investigating the genetic variation in populations of a natural contact zone. Nucleotide sequence variation of a 395-bp segment of the mitochondrial 16S rRNA gene was compared to that of six nuclear allozyme loci. Three major mtDNA lineages were found, the eastern clade being consistent with the former allozyme type A. The two western clades (types B and C) were not distinguished previously. Strong sequence divergence and correlation with nuclear genetic isolation in syntopic populations, however, justifies the specific status of the three G. fossarum types. The common speciation event is believed to be very old (Miocene). The within-type mtDNA variation is probably molded by the ice ages, with type B populations being most affected. Moreover, the patch-like distribution of mtDNA type B lineages in an area near the contact zone corroborates the hypothesis of a recent colonization.     

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.     

Divergent character clines across a recent secondary contact zone in a Hispaniolan lizard

Studies of genetic contact zones provide valuable information regarding the processes of population divergence, adaptation and speciation. In this paper, I examine transitions in morphology, mitochondrial DNA (mtDNA) and nuclear DNA (nDNA) haplotypes across a recent secondary contact zone in a Hispaniolan lizard Ameiva chrysolaema . Maximum likelihood cline fitting analyses suggest non-coincidence of cline centers and that the mtDNA cline is significantly displaced to the west of the remaining clines. nDNA and morphological clines are coincident and tend to be associated with the prevailing environmental gradient. The lack of cytonuclear disequilibrium near the center of the contact zone and the non-coincidence of character clines suggest that this zone does not conform to a tension zone model of hybridization; thus, gene flow across the zone does not seem to be impeded. The extremely narrow width of the dorsal scale size cline and the close association of this cline with the steepness of the environmental (precipitation) gradient suggest that this character may be under environmental selection. Taken together, this contact zone appears to be structured by a combination of mtDNA introgression, possibly associated with eastward movement of the zone, and environmental selection on some characters.     

Genetic structure, introgression, and a narrow hybrid zone between northern and California spotted owls (Strix occidentalis)

The northern spotted owl (Strix occidentalis caurina) is a threatened subspecies and the California spotted owl (Strix occidentalis occidentalis) is a subspecies of special concern in the western United States. Concern for their continued viability has arisen because of habitat loss caused by timber harvesting. The taxonomic status of the northern subspecies has been the subject of continuing controversy. We investigated the phylogeographical and population genetic structure of northern and California spotted owls with special reference to their region of contact. Mitochondrial DNA (mtDNA) control region sequences confirmed the existence of two well-differentiated lineages connected by a narrow hybrid zone in a region of low population density in north central California. Maximum-likelihood estimates indicated bidirectional gene flow between the lineages but limited introgression outside the region of contact. The lengths of both the mtDNA hybrid zone and the reduced density patch were similar and slightly exceeded estimates of natal dispersal distances. This suggests that the two subspecies were in secondary contact in a hybrid zone trapped by a population density trough. Consequently, the zone of interaction is expected to be geographically stable. We discovered a third, rare clade of haplotypes, which we interpreted to be a result of incomplete lineage sorting; those haplotypes result in a paraphyletic northern spotted owl with respect to the California spotted owl. A congeneric species, the barred owl (Strix varia), occasionally hybridizes with spotted owls; our results indicated an upper bound for the frequency of barred owl mtDNA haplotypes in northern spotted owl populations of 3%.     

Spatial and temporal genetic structure of the planktonic Sagitta setosa (Chaetognatha) in European seas as revealed by mitochondrial and nuclear DNA markers

Little is known about the spatial and temporal scales at which planktonic organisms are genetically structured. A previous study of mitochondrial DNA (mtDNA) in the holoplanktonic chaetognath Sagitta setosa revealed strong phylogeographic structuring suggesting that Northeast (NE) Atlantic, Mediterranean and Black Sea populations are genetically disjunct. The present study used a higher sampling intensity and a combination of mitochondrial and four microsatellite markers to reveal population structuring between and within basins. Between basins, both marker sets indicated significant differentiation confirming earlier results that gene flow is probably absent between the respective S. setosa populations. At the within-basin scale, we found no evidence of spatial or temporal structuring within the NE Atlantic. In the Mediterranean basin, both marker sets indicated significant structuring, but only the mtDNA data indicated a sharp genetic division between Adriatic and all other Mediterranean populations. Data were inconclusive about population structuring in the Black Sea. The levels of differentiation indicated by the two marker sets differed substantially, with far less pronounced structure detected by microsatellite than mtDNA data. This study also uncovered the presence of highly divergent mitochondrial lineages that were discordant with morphology, geography and nuclear DNA. We thus propose the hypothesis that highly divergent mitochondrial lineages may be present within interbreeding S. setosa populations.     

Genetic population structure of the paper wasp Polistes olivaceus (Hymenoptera: Vespidae) in Bangladesh

Dispersal triggers gene flow, which in turn strongly affects the ensuing genetic population structure of a species. Using nuclear microsatellite loci and mitochondrial DNA (mtDNA), we estimated the genetic population structure of the wasp Polistes olivaceus throughout Bangladesh. The level of population differentiation using nuclear markers (F _ST) appeared to be much lower than that estimated using mtDNA haplotype sequences (Ф_ST), even after correcting for effective population size differences between the two markers. These results suggest a philopatric tendency, in which gynes disperse less than males. We observed no isolation by distance among the study populations at either the nuclear or mtDNA level, suggesting nonequilibrium between gene flow and drift as a result of very frequent interpopulation movement. For the nuclear markers, an individual assignment test showed no genetically and geographically distinct groups. Instead, phylogenetic analyses as well as a minimum spanning network using mtDNA haplotypes consistently revealed two distinct lineages. The distribution of haplotypes indicated western populations with a single lineage and offered clear evidence for restricted gene flow across the Jamuna–Padma–Upper Meghna river system. Mismatch distributions exhibited a unimodal distribution, which along with a starlike haplotype network, suggested a population expansion in lineage I but not in lineage II. Overall, these results suggest that gene flow among populations of P. olivaceus was affected by both female philopatry and a major river system across Bangladesh.