Pervasive selection or is it…? why are FST outliers sometimes so frequent?

It is now common for population geneticists to estimate F_ST for a large number of loci across the genome, before testing for selected loci as being outliers to the F_ST distribution. One surprising result of such F_ST scans is the often high proportion (>1% and sometimes >10%) of outliers detected, and this is often interpreted as evidence for pervasive local adaptation. In this issue of Molecular Ecolog, Fourcade et al. ( 2013 ) observe that a particularly high rate of F_ST outliers has often been found in river organisms, such as fishes or damselflies, despite there being no obvious reason why selection should affect a larger proportion of the genomes of these organisms. Using computer simulations, Fourcade et al. ( 2013 ) show that the strong correlation in co‐ancestry produced in long one‐dimensional landscapes (such as rivers, valleys, peninsulas, oceanic ridges or coastlines) greatly increases the neutral variance in F_ST, especially when the landscape is further reticulated into fractal networks. As a consequence, outlier tests have a high rate of false positives, unless this correlation can be taken into account. Fourcade et al.'s study highlights an extreme case of the general problem, first noticed by Robertson ( 1975a , b ) and Nei & Maruyama ( 1975 ), that correlated co‐ancestry inflates the neutral variance in F_ST when compared to its expectation under an island model of population structure. Similar warnings about the validity of outlier tests have appeared regularly since then but have not been widely cited in the recent genomics literature. We further emphasize that F_ST outliers can arise in many different ways and that outlier tests are not designed for situations where the genetic architecture of local adaptation involves many loci.     

Microgeographic structure of Anopheles gambiae in western Kenya based on mtDNA and microsatellite loci

The population genetic structure of the Anopheles gambiae in western Kenya was studied using length variation at five microsatellite loci and sequence variation in a 648-nt mtDNA fragment. Mosquitoes were collected from houses in villages spanning up to 50 km distance, The following questions were answered, (i) Are mosquitoes in a house more related genetically to each other than mosquitoes between houses? (ii) What degree of genetic differentiation occurs on these geographical scales? (iii) How consistent are the results obtained with both types of genetic markers? At the house level, no differentiation was detected by F_ST and R_ST, and the band sharing index test revealed no significant associations of alleles across loci. Likewise, indices of kinship based on mtDNA haplotypes in houses were even lower than in the pooled sample. Therefore, the hypothesis that mosquitoes in a house are more related genetically was rejected. At increasing geographical scales, microsatellite allele distributions were similar among all population samples and no subdivision of the gene pool was detected by F_ST or R_ST. Likewise, estimates of haplotype divergence of mtDNA between populations were not higher than the within population estimates, and mtDNA-based F_ST values were not significantly different from zero. That sequence variation in mtDNA provided matching results with microsatellite loci (while high genetic variation was observed in all loci), suggested that this pattern represents the whole genome. The minimum area associated with a deme of A. gambiae in western Kenya is therefore larger than 50 km in diameter.     

Adaptive divergence despite strong genetic drift: genomic analysis of the evolutionary mechanisms causing genetic differentiation in the island fox (Urocyon littoralis)

The evolutionary mechanisms generating the tremendous biodiversity of islands have long fascinated evolutionary biologists. Genetic drift and divergent selection are predicted to be strong on islands and both could drive population divergence and speciation. Alternatively, strong genetic drift may preclude adaptation. We conducted a genomic analysis to test the roles of genetic drift and divergent selection in causing genetic differentiation among populations of the island fox (Urocyon littoralis). This species consists of six subspecies, each of which occupies a different California Channel Island. Analysis of 5293 SNP loci generated using Restriction‐site Associated DNA (RAD) sequencing found support for genetic drift as the dominant evolutionary mechanism driving population divergence among island fox populations. In particular, populations had exceptionally low genetic variation, small N_e (range = 2.1–89.7; median = 19.4), and significant genetic signatures of bottlenecks. Moreover, islands with the lowest genetic variation (and, by inference, the strongest historical genetic drift) were most genetically differentiated from mainland grey foxes, and vice versa, indicating genetic drift drives genome‐wide divergence. Nonetheless, outlier tests identified 3.6–6.6% of loci as high F_ST outliers, suggesting that despite strong genetic drift, divergent selection contributes to population divergence. Patterns of similarity among populations based on high F_ST outliers mirrored patterns based on morphology, providing additional evidence that outliers reflect adaptive divergence. Extremely low genetic variation and small N_e in some island fox populations, particularly on San Nicolas Island, suggest that they may be vulnerable to fixation of deleterious alleles, decreased fitness and reduced adaptive potential.     

Genetic differentiation between humpback whales (Megaptera novaeangliae) from Atlantic and Pacific breeding grounds of South America

Humpback whales wintering in tropical waters along the Atlantic and Pacific coasts of the South American continent are thought to represent distinct populations or “stocks.” Here we present the first analysis of genetic differentiation and estimates of gene flow between these breeding stocks, based on both mitochondrial DNA (mtDNA) control region sequences (465 bp) and 16 microsatellite loci from samples collected off Brazil (n = 277) and Colombia (n = 148), as well as feeding areas near the western Antarctic Peninsula (n = 86). We found significant differentiation between Brazilian and Colombian breeding grounds at both mtDNA (F_ST = 0.058) and microsatellite (F_ST = 0.011) markers and corroborated previous studies showing genetic similarity between humpbacks from Colombia and those from Antarctic Peninsula feeding areas. Estimates of long‐term gene flow between Brazil and Colombia were low to moderate, asymmetrical, and mostly mediated by males. Assignment procedures detected some cases of interchange and individuals of admixed ancestry between breeding grounds, indicating limited mixing of individuals between these stocks. Overall, results highlight the differentiation of humpback whale breeding populations with adjacent feeding grounds. This appears to be a remarkable example of fidelity to seasonal habitat in the absence of any contemporary barriers.     

Application of Microsatellite Markers to Population Genetics Studies of Japanese Flounder Paralichthys olivaceus

We examined population genetic structure by means of microsatellite analysis among 7 Japanese flounder (Paralichthys olivaceus) populations collected from coastal sea areas around Japan. As was expected, all of the 11 microsatellite loci examined were variable in all populations (number of alleles per locus, 15.2–18.2; average of expected heterozygosity, 0.74–0.78). Eleven population pairs in 21 possible pairwise comparisons showed significant genetic heterogeneity associated with allele frequency distributions or fixation index (F _ST ). Modified Cavalli-Sforza chord distance (D _A ) and Nei's standard genetic distance (D _ST ) ranged from 0.051 to 0.090, and from 0.000 to 0.025, respectively. There was evidence that the populations assessed in this study were not drawn from a single panmictic population; however, it appears that Japanese flounder populations around Japan are not well-structured, as an estimate of the fixation index value among the 7 localities was very low (F _ST = 0.0025). Received April 27, 2001; accepted June 27, 2001.     

Genetic diversity and differentiation of Acanthoscelides obtectus Say (Coleoptera: Bruchidae) populations in China

相似文献   

High degree of genetic differentiation in marine three‐spined sticklebacks (Gasterosteus aculeatus)

Populations of widespread marine organisms are typically characterized by a low degree of genetic differentiation in neutral genetic markers, but much less is known about differentiation in genes whose functional roles are associated with specific selection regimes. To uncover possible adaptive population divergence and heterogeneous genomic differentiation in marine three‐spined sticklebacks (Gasterosteus aculeatus), we used a candidate gene‐based genome‐scan approach to analyse variability in 138 microsatellite loci located within/close to (<6 kb) functionally important genes in samples collected from ten geographic locations. The degree of genetic differentiation in markers classified as neutral or under balancing selection—as determined with several outlier detection methods—was low (F_ST = 0.033 or 0.011, respectively), whereas average F_ST for directionally selected markers was significantly higher (F_ST = 0.097). Clustering analyses provided support for genomic and geographic heterogeneity in selection: six genetic clusters were identified based on allele frequency differences in the directionally selected loci, whereas four were identified with the neutral loci. Allelic variation in several loci exhibited significant associations with environmental variables, supporting the conjecture that temperature and salinity, but not optic conditions, are important drivers of adaptive divergence among populations. In general, these results suggest that in spite of the high degree of physical connectivity and gene flow as inferred from neutral marker genes, marine stickleback populations are strongly genetically structured in loci associated with functionally relevant genes.     

Genetic variation in two populations of the rough‐skinned newt (Taricha granulosa) assessed using novel tetranucleotide microsatellite loci

Seven novel tetranucleotide microsatellite loci were identified from a partial genomic DNA library, enriched for GATA‐motif microsatellites, from the rough‐skinned newt (Taricha granulosa). All loci were polymorphic, and one displayed a high frequency null allele. A related species, T. rivularis, displays strong site fidelity and detectable population structure over small spatial scales, so we assessed genetic variation in two samples of T. granulosa separated by 16 km. Distributions of allele frequencies differ significantly between our two sites, but small values of F_ST and Rho_ST suggest that the populations are linked by a large amount of gene flow.     

Comparison of quantitative and molecular variation in agroforestry populations of the shea tree (Vitellaria paradoxa C.F. Gaertn) in Mali

In this study we investigated the within- and between-population genetic variation using microsatellite markers and quantitative traits of the shea tree, Vitellaria paradoxa, an important agroforestry tree species of the Sudano–Sahelian region in Africa. Eleven populations were sampled across Mali and in northern Côte d’Ivoire. Leaf size and form and growth traits were measured in a progeny test at the nursery stage. Eight microsatellites were used to assess neutral genetic variation. Low levels of heterozygosity were recorded (1.6–3.0 alleles/locus; H_E = 0.25–0.42) and the fixation index (F_IS = −0.227–0.186) was not significantly different from zero suggesting that Hardy–Weinberg equilibrium is encountered in all populations sampled. Quantitative traits exhibited a strong genetic variation between populations and between families within populations. The degree of population differentiation of the quantitative traits (Q_ST = 0.055–0.283, Q_STmean = 0.189) strongly exceeds that in eight microsatellite loci (F_ST = −0.011–0.142, F_STmean = 0.047). Global and pairwise F_ST values were very low and not significantly different from zero suggesting agroforestry practices are amplifying gene flow (Nm = 5.07). The population means for quantitative traits and the rainfall variable were not correlated, showing variation was not linked with this climatic cline. It is suggested that this marked differentiation for quantitative traits, independent of environmental clines and despite a high gene flow, is a result of local adaptation and human selection of shea trees. This process has induced high linkage disequilibrium between underlying loci of polygenic characters.     

Risso’s dolphins (<Emphasis Type="Italic">Grampus griseus</Emphasis>) in UK waters are differentiated from a population in the Mediterranean Sea and genetically less diverse

The Risso’s dolphin (Grampus griseus) has a worldwide distribution, but little is known about their population genetic structure. Local coastal populations are never known to be abundant, and are sometimes under anthropogenic impact. Therefore the question of regional differentiation by genetic drift and overall diversity levels is of conservation interest. Here we present preliminary data that clearly indicate genetic differentiation and lower genetic diversity of a population sampled in UK waters (primarily from the Western Isles, Scotland) compared to a Mediterranean sample. Significant differentiation was found in comparison with a sample from the Mediterranean for both microsatellite DNA markers (F _ST = 0.0296) and mtDNA sequence data (F _ST = 0.260; ϕ _ST = 0.542). Allelic diversity was lower in the UK for nearly all loci.     

Heterogeneity and concordance in locus‐specific differentiation and introgression between species of towhees

The maintenance or breakdown of reproductive isolation is an observable outcome of secondary contact between species. In cases where hybrids beyond the F1 are formed, the representation of each species' ancestry can vary dramatically among genomic regions. This genomic heterogeneity in ancestry and introgression can offer insight into evolutionary processes, particularly if introgression is compared in multiple hybrid zones. Similarly, considerable heterogeneity exists across the genome in the extent to which populations and species have diverged, reflecting the combined effects of different evolutionary processes on genetic variation. We studied hybridization across two hybrid zones of two phenotypically well‐differentiated bird species in Mexico (Pipilo maculatus and P. ocai), to investigate genomic heterogeneity in differentiation and introgression. Using genotyping‐by‐sequencing (GBS) and hierarchical Bayesian models, we genotyped 460 birds at over 41 000 single nucleotide polymorphism (SNP) loci. We identified loci exhibiting extreme introgression relative to the genome‐wide expectation using a Bayesian genomic cline model. We also estimated locus‐specific F_ST and identified loci with exceptionally high genetic divergence between the parental species. We found some concordance of locus‐specific introgression in the two independent hybrid zones (6–20% of extreme loci shared across zones), reflecting areas of the genome that experience similar gene flow when the species interact. Additionally, heterogeneity in introgression and divergence across the genome revealed another subset of loci under the influence of locally specific factors. These results are consistent with a history in which reproductive isolation has been influenced by a common set of loci in both hybrid zones, but where local environmental and stochastic factors also lead to genomic differentiation.     

Genome-Wide Scan for Signatures of Human Population Differentiation and Their Relationship with Natural Selection,Functional Pathways and Diseases

Genetic differences both between individuals and populations are studied for their evolutionary relevance and for their potential medical applications. Most of the genetic differentiation among populations are caused by random drift that should affect all loci across the genome in a similar manner. When a locus shows extraordinary high or low levels of population differentiation, this may be interpreted as evidence for natural selection. The most used measure of population differentiation was devised by Wright and is known as fixation index, or F_ST. We performed a genome-wide estimation of F_ST on about 4 millions of SNPs from HapMap project data. We demonstrated a heterogeneous distribution of F_ST values between autosomes and heterochromosomes. When we compared the F_ST values obtained in this study with another evolutionary measure obtained by comparative interspecific approach, we found that genes under positive selection appeared to show low levels of population differentiation. We applied a gene set approach, widely used for microarray data analysis, to detect functional pathways under selection. We found that one pathway related to antigen processing and presentation showed low levels of F_ST, while several pathways related to cell signalling, growth and morphogenesis showed high F_ST values. Finally, we detected a signature of selection within genes associated with human complex diseases. These results can help to identify which process occurred during human evolution and adaptation to different environments. They also support the hypothesis that common diseases could have a genetic background shaped by human evolution.     

Climate variables explain neutral and adaptive variation within salmonid metapopulations: the importance of replication in landscape genetics

Understanding how environmental variation influences population genetic structure is important for conservation management because it can reveal how human stressors influence population connectivity, genetic diversity and persistence. We used riverscape genetics modelling to assess whether climatic and habitat variables were related to neutral and adaptive patterns of genetic differentiation (population‐specific and pairwise F_ST) within five metapopulations (79 populations, 4583 individuals) of steelhead trout (Oncorhynchus mykiss) in the Columbia River Basin, USA. Using 151 putatively neutral and 29 candidate adaptive SNP loci, we found that climate‐related variables (winter precipitation, summer maximum temperature, winter highest 5% flow events and summer mean flow) best explained neutral and adaptive patterns of genetic differentiation within metapopulations, suggesting that climatic variation likely influences both demography (neutral variation) and local adaptation (adaptive variation). However, we did not observe consistent relationships between climate variables and F_ST across all metapopulations, underscoring the need for replication when extrapolating results from one scale to another (e.g. basin‐wide to the metapopulation scale). Sensitivity analysis (leave‐one‐population‐out) revealed consistent relationships between climate variables and F_ST within three metapopulations; however, these patterns were not consistent in two metapopulations likely due to small sample sizes (N = 10). These results provide correlative evidence that climatic variation has shaped the genetic structure of steelhead populations and highlight the need for replication and sensitivity analyses in land and riverscape genetics.     

Dissecting the genetic structure of Korean population using genome-wide SNP arrays

Genome-wide SNP arrays have generated unprecedented quantities of data allow the detection of human evolutionary history and dense genome-wide data also enable the identification of distance ancestry among individuals or ethnic groups. To explain wider aspects of the genetic structure of Koreans and the East Asian population, we analyzed 79 individuals from the Korean HapMap project at 555,352 common single-nucleotide polymorphism loci, and compared this data with the worldwide population groups with the 53 ethnic groups from Human Genome Diversity Panel (HGDP-CEPH). Population differentiation (F_ST), Principal Component Analyses, STRUCTURE and ADMIXTURE are examined. In general, all the individual samples studies here were classified into subset of ethnic groups according to their geographical origins. Korean HapMap individuals were grouped together with East Asian populations from HGDP panel. Recently, a sub-population structure within Korean population has been reported. Our result, however, revealed the genetic homogeneity of Korean population. The ADMIXTURE analysis showed that, overall the Korean populations derive 79 % of their genomic ancestry from southern Asia and have relatively little northern Asian ancestry (21 %). The present work, therefore, provide the evidence that the male-biased southern-to-northern migration influenced not only for the genetic make up of the Y chromosome in the Korean population but also, its autosomal composition.     

Fine-scale genetic population structure in a mobile marine mammal: inshore bottlenose dolphins in Moreton Bay, Australia

Highly mobile marine species in areas with no obvious geographic barriers are expected to show low levels of genetic differentiation. However, small‐scale variation in habitat may lead to resource polymorphisms and drive local differentiation by adaptive divergence. Using nuclear microsatellite genotyping at 20 loci, and mitochondrial control region sequencing, we investigated fine‐scale population structuring of inshore bottlenose dolphins (Tursiops aduncus) inhabiting a range of habitats in and around Moreton Bay, Australia. Bayesian structure analysis identified two genetic clusters within Moreton Bay, with evidence of admixture between them (F_ST = 0.05, P = 0.001). There was only weak isolation by distance but one cluster of dolphins was more likely to be found in shallow southern areas and the other in the deeper waters of the central northern bay. In further analysis removing admixed individuals, southern dolphins appeared genetically restricted with lower levels of variation (AR = 3.252, π = 0.003) and high mean relatedness (r = 0.239) between individuals. In contrast, northern dolphins were more diverse (AR = 4.850, π = 0.009) and were mixing with a group of dolphins outside the bay (microsatellite‐based STRUCTURE analysis), which appears to have historically been distinct from the bay dolphins (mtDNA Φ_ST = 0.272, P < 0.001). This study demonstrates the ability of genetic techniques to expose fine‐scale patterns of population structure and explore their origins and mechanisms. A complex variety of inter‐related factors including local habitat variation, differential resource use, social behaviour and learning, and anthropogenic disturbances are likely to have played a role in driving fine‐scale population structure among bottlenose dolphins in Moreton Bay.     

Signatures of adaptation and genetic structure among the mainland populations of Pinus radiata (D. Don) inferred from SNP loci

Insights into the relative contributions of locus specific and genome-wide effects on population genetic diversity can be gained through separation of their resulting genetic signals. Here we explore patterns of adaptive and neutral genetic diversity in the disjunct natural populations of Pinus radiata (D. Don) from mainland California. A first-generation common garden of 447 individuals revealed significant differentiation of wood phenotypes among populations (P _ST), possibly reflecting local adaptation in response to environment. We subsequently screened all trees for genetic diversity at 149 candidate gene single nucleotide polymorphism (SNP) loci for signatures of adaptation. Ten loci were identified as being possible targets of diversifying selection following F _ST outlier tests. Multivariate canonical correlation performed on a data set of 444 individuals identified significant covariance between environment, adaptive phenotypes and outlier SNP diversity, lending support to the case for local adaptation suggested from F _ST and P _ST tests. Covariation among discrete sets of outlier SNPs and adaptive phenotypes (inferred from multivariate loadings) with environment are supported by existing studies of candidate gene function and genotype–phenotype association. Canonical analyses failed to detect significant correlations between environment and 139 non-outlier SNP loci, which were applied to estimate neutral patterns of genetic differentiation among populations (F _ST 4.3 %). Using this data set, significant hierarchical structure was detected, indicating three populations on the mainland. The hierarchical relationships based on neutral SNP markers (and SSR) were in contrast with those inferred from putatively adaptive loci, potentially highlighting the independent action of selection and demography in shaping genetic structure in this species.     

Relative performance of Bayesian clustering software for inferring population substructure and individual assignment at low levels of population differentiation

Traditional methods for characterizing genetic differentiation among populations rely on a priori grouping of individuals. Bayesian clustering methods avoid this limitation by using linkage and Hardy–Weinberg disequilibrium to decompose a sample of individuals into genetically distinct groups. There are several software programs available for Bayesian clustering analyses, all of which describe a decrease in the ability to detect distinct clusters as levels of genetic differentiation among populations decrease. However, no study has yet compared the performance of such methods at low levels of population differentiation, which may be common in species where populations have experienced recent separation or high levels of gene flow. We used simulated data to evaluate the performance of three Bayesian clustering software programs, PARTITION, STRUCTURE, and BAPS, at levels of population differentiation below F _ST=0.1. PARTITION was unable to correctly identify the number of subpopulations until levels of F _ST reached around 0.09. Both STRUCTURE and BAPS performed very well at low levels of population differentiation, and were able to correctly identify the number of subpopulations at F _ST around 0.03. The average proportion of an individual’s genome assigned to its true population of origin increased with increasing F _ST for both programs, reaching over 92% at an F _ST of 0.05. The average number of misassignments (assignments to the incorrect subpopulation) continued to decrease as F _ST increased, and when F _ST was 0.05, fewer than 3% of individuals were misassigned using either program. Both STRUCTURE and BAPS worked extremely well for inferring the number of clusters when clusters were not well-differentiated (F _ST=0.02–0.03), but our results suggest that F _ST must be at least 0.05 to reach an assignment accuracy of greater than 97%.     

Genetic diversity and population differentiation in the yellow drum Nibea albiflora along the coast of the China Sea

Nibea albiflora (yellow drum) is an important seafood fish species in East Asia. We explored the population genetic variation of N. albiflora along the coastal waters of the China Sea using microsatellite markers to facilitate a selective breeding programme that is undertaken in China. A total of 256 alleles were detected at 12 loci in four wild populations. A high level of genetic diversity was observed with the mean number of alleles and the observed and expected heterozygosity in each population ranging from 7.917 to 14.083, 0.701 to 0.764 and 0.765 to 0.841, respectively. Pairwise fixation index (F_ST) analysis indicated significant but weak genetic differentiation among populations from four localities (F_ST?=?0.030, P?<?0.01), which was also confirmed by analysis of molecular variance (AMOVA). Significant genetic differentiation was detected between Ningde and the other populations (F_ST?=?0.047–0.056, P?<?0.01). Structure analysis suggested that N. albiflora within the examined range might be composed of two stocks. The data of the present study revealed high genetic diversity and low genetic differentiation among the N. albiflora populations along the coast of the China Sea. This baseline information could be valuable for future selective breeding programmes of N. albiflora.     

Extensive gene flow along the urban–rural gradient in a migratory colonial bird

Urban colonization by wildlife involves a combination of several different mechanisms, including phenotype or genotype sorting, phenotypic plasticity and microevolutionary adaptation. Combination of these processes can produce a rapid phenotypic, but also genetic divergence of urban versus rural populations. Here, we examined the pattern of genetic differentiation between urban and rural populations of a colonial migratory bird, the black‐headed gull Chroicocephalus ridibundus. To this end, we sampled ca 170 individuals from six (two urban and four rural) colonies in northern Poland, and genotyped them at ten microsatellite loci. Our analysis provided evidence for negligible genetic divergence of urban and rural colonies, as assessed with fixation index F_ST and Nei's unbiased genetic distance D (mean pairwise urban‐rural comparisons: F_ST = 0.003 ± 0.001 [SE] and D = 0.012 ± 0.006 [SE]). Bayesian clustering methods provided support for homogeneous genetic structure across all urban and rural populations. Also, we found no support for reduced allelic diversity in urban versus rural colonies. These results stand in a stark contrast to the previous findings on the genetic consequences of urbanization in birds. We hypothesize that this pattern could possibly be attributed to the important life‐history characters of the black‐headed gull, including coloniality, migratoriness, and high dispersal propensity. Our study provides a novel insight into the urban landscape genetics, underlining large variation in the mechanisms of urban colonization and its genetic consequences in wild animal populations.     

Historical vicariance and male-mediated gene flow in the toad-headed lizards Phrynocephalus przewalskii

Using mitochondrial and microsatellite DNA data and a population genetic approach, we tested male‐mediated gene flow in the toad‐headed lizards Phrynocephalus przewalskii. The mitochondrial DNA (ND2 gene), on the one hand, revealed two major lineages and a strong population genetic structure (F_ST = 0.692; F_ST′ = 0.995). The pairwise differences between the two lineages ranged from 2.1% to 6.4% and the geographical division of the two lineages coincided with a mountain chain consisting of the Helan and Yin Mountains, suggesting a historical vicariant pattern. On the other hand, the nuclear microsatellite DNA revealed a significant but small population genetic structure (F_ST = 0.017; F_ST′ = 0.372). The pairwise F_ST among the nine populations examined with seven microsatellite DNA loci ranged from 0.0062 to 0.0266; the assignment test failed to detect any naturally occurring population clusters. Furthermore, the populations demonstrated a weak isolation by distance and a northeast to southwest clinal variation, rather than a vicariant pattern. A historical vicariant event followed by male‐mediated gene flow appears to be the best explanation for the data. Approximately 2–5 Ma, climatic change may have created an uninhabitable zone along the Helan‐Yin mountain chain and initiated the divergence between the two mitochondrial lineages. With further climatic changes, males were able to disperse across the mountain chain, causing sufficient gene flow that eventually erased the vicariant pattern and drastically reduced the population genetic structure, while females remained philopatric and maintained the mitochondrial DNA (mtDNA) divergence. Although polygyny mating system and female philopatry may partially contribute to the reduced movement of females, other hypotheses, such as female intrasexual aggression, should also be explored.