Natural selection on protein polymorphism in the rodent genus Peromyscus: evidence from interlocus contrasts

The effects of natural selection are generally locus-specific, whereas migration, drift, and inbreeding are expected to have relatively uniform effects across the entire genome. This suggests that multilocus surveys of multiple populations can be used to distinguish selection from demographic effects. The purpose of this study was to test for evidence of selection on protein polymorphism in natural populations of mice in the genus Peromyscus. We analyzed published data from geographic surveys of allozyme variation and used a coalescent-based simulation model to identify specific loci that deviated from neutral expectations. Observed F(ST) values generally exhibited a remarkably close fit to the expected neutral distributions, indicating that the majority of loci are simply tracking stochastic demographic processes. A smaller number of loci exhibited highly significant departures from the expectations of the neutral model and thus appear to be tracking the direct or indirect effects of selection. Most departures from neutrality were characterized by F(ST) values that far exceeded neutral expectations and were therefore attributable to spatially varying selection. Interestingly, the albumin locus was implicated as a candidate gene for local adaptation in four different species of Peromyscus. The results also demonstrate that selection can severely bias marker-based estimates of neutral parameters.     

Multilocus self-recognition systems in fungi as a cause of trans-species polymorphism

Trans-species polymorphism, meaning the presence of alleles in different species that are more similar to each other than they are to alleles in the same species, has been found at loci associated with vegetative incompatibility in filamentous fungi. If individuals differ at one or more of these loci (termed het for heterokaryon), they cannot form stable heterokaryons after vegetative fusion. At the het-c locus in Neurospora crassa and related species there is clear evidence of trans-species polymorphism: three alleles have persisted for approximately 30 million years. We analyze a population genetic model of multilocus vegetative incompatibility and find the conditions under which trans-species polymorphism will occur. In the model, several unlinked loci determine the vegetative compatibility group (VCG) of an individual. Individuals of different VCGs fail to form productive heterokaryons, while those of the same VCG form viable heterokaryons. However, viable heterokaryon formation between individuals of the same VCG results in a loss in fitness, presumably via transfer of infectious agents by hyphal fusion or exploitation by aggressive genotypes. The result is a form of balancing selection on all loci affecting an individual's VCG. We analyze this model by making use of a Markov chain/strong selection, weak mutation (SSWM) approximation. We find that trans-species polymorphism of the type that has been found at the het-c locus is expected to occur only when the appearance of new incompatibility alleles is strongly constrained, because the rate of mutation to such alleles is very low, because the number of possible incompatibility alleles at each locus is restricted, or because the number of incompatibility loci is limited.     

Population genetics of speciation in two closely related wild tomatoes (Solanum section Lycopersicon)

We present a multilocus sequencing study to assess patterns of polymorphism and divergence in the closely related wild tomato species, Solanum peruvianum and S. chilense (Solanum section Lycopersicon, Solanaceae). The data set comprises seven mapped nuclear loci (approximately 9.3 kb of analyzed sequence across loci) and four local population samples per species that cover much of the species' range (between 80 and 88 sequenced alleles across both species). We employ the analytical framework of divergence population genetics (DPG) in evaluating the utility of the "isolation" model of speciation to explain observed patterns of polymorphism and divergence. Whereas the isolation model is not rejected by goodness-of-fit criteria established via coalescent simulations, patterns of intragenic linkage disequilibrium provide evidence for postdivergence gene flow at two of the seven loci. These results suggest that speciation occurred under residual gene flow, implying that natural selection is one of the evolutionary forces driving the divergence of these tomato species. This inference is fully consistent with their recent divergence, conservatively estimated to be 相似文献   

Combined linkage disequilibrium and linkage mapping: Bayesian multilocus approach

Quantitative trait loci (QTL) affecting the phenotype of interest can be detected using linkage analysis (LA), linkage disequilibrium (LD) mapping or a combination of both (LDLA). The LA approach uses information from recombination events within the observed pedigree and LD mapping from the historical recombinations within the unobserved pedigree. We propose the Bayesian variable selection approach for combined LDLA analysis for single-nucleotide polymorphism (SNP) data. The novel approach uses both sources of information simultaneously as is commonly done in plant and animal genetics, but it makes fewer assumptions about population demography than previous LDLA methods. This differs from approaches in human genetics, where LDLA methods use LA information conditional on LD information or the other way round. We argue that the multilocus LDLA model is more powerful for the detection of phenotype–genotype associations than single-locus LDLA analysis. To illustrate the performance of the Bayesian multilocus LDLA method, we analyzed simulation replicates based on real SNP genotype data from small three-generational CEPH families and compared the results with commonly used quantitative transmission disequilibrium test (QTDT). This paper is intended to be conceptual in the sense that it is not meant to be a practical method for analyzing high-density SNP data, which is more common. Our aim was to test whether this approach can function in principle.     

Multilocus analysis of variation using a large empirical data set: phenylpropanoid pathway genes in Arabidopsis thaliana

Detecting the signature of adaptation on nucleotide variation is often difficult in species that like Arabidopsis thaliana might have a complex demographic history. Recent re-sequencing surveys in this species provided genome-wide information that would mainly reflect its demographic history. We have used a large empirical data set (LED) as well as multilocus coalescent simulations to analyse sequence variation at loci involved in the phenylpropanoid pathway of this species. We surveyed and examined DNA sequence variation at nine of these loci (about 19.7 kb) in 23 accessions of A. thaliana and one accession of its closely related species Arabidopsis lyrata . Nucleotide variation was lower at nonsynonymous sites than at silent sites in all loci, indicating generalized functional constraint at the protein level. No association between variation and position in the metabolic pathway was detected. When the data were contrasted against the standard neutral model, significant deviations for silent variation were detected with Tajima's D , Fu's F_S and Fay and Wu's H multilocus test statistics. These deviations were in the same direction than in previous large-scale multilocus analyses, suggesting a genome-wide effect. When the nine-locus data set was contrasted against the large empirical data set, the level (Watterson's θ) and pattern of variation (Tajima's D ) detected in these loci did not deviate either at the single-locus or multilocus level from the corresponding empirical distributions. These results would support an important role of the demographic history of A. thaliana in shaping nucleotide variation at the nine studied phenylpropanoid loci. The potential and limitations of the empirical distribution approach are discussed.     

Genome scans of DNA variability in humans reveal evidence for selective sweeps outside of Africa

The last 50,000-150,000 years of human history have been characterized by rapid demographic expansions and the colonization of novel environments outside of sub-Saharan Africa. Mass migrations outside the ancestral species range likely entailed many new selection pressures, suggesting that genetic adaptation to local environmental conditions may have been more prevalent in colonizing populations outside of sub-Saharan Africa. Here we report a test of this hypothesis using genome-wide patterns of DNA polymorphism. We conducted a multilocus scan of microsatellite variability to identify regions of the human genome that may have been subject to continent-specific hitchhiking events. Using published polymorphism data for a total of 624 autosomal loci in multiple populations of humans, we used coalescent simulations to identify candidate loci for geographically restricted selective sweeps. We identified a total of 13 loci that appeared as outliers in replicated population comparisons involving different reference samples for Africa. A disproportionate number of these loci exhibited reduced levels of relative variability in non-African populations alone, suggesting that recent episodes of positive selection have been more prevalent outside of sub-Saharan Africa.     

Evidence for a slightly deleterious effect of intron polymorphisms at the EF1alpha gene in the deep-sea hydrothermal vent bivalve Bathymodiolus

A multilocus analysis was initiated in order to infer the general effect of demography and the indirect effect of positive selection on some chromosome segments in Bathymodiolus. Mussels of the genus Bathymodiolus inhabit the very hostile, fragmented and variable environment of deep-sea hydrothermal vents which is thought to cause recurrent population bottlenecks via extinction/colonisation processes and adaptation to new environmental conditions. In the course of this work we discovered that the assumption of neutrality of non-coding polymorphisms usually made in genome scan experiments was likely to be violated at one of the loci we analysed. The direct effect of slight purifying selection on non-coding polymorphisms shares many resemblances with the indirect effect of positive selection through genetic hitchhiking. Combining polymorphism with divergence data for several closely related species allowed us to obtain different expectations for the direct effect of negative selection and the indirect effect of positive selection. We observed a strong excess of rare non-coding polymorphisms at the second intron of the EF1alpha gene in the two species Bathymodiolus azoricus and Bathymodiolus thermophilus, while two other loci, the mitochondrial COI gene and an intron of the Lysozyme gene, did not exhibit such a deviation. In addition, the divergence rate of the EF1alpha intron was estimated to be unexpectedly low when calibrated using the closure of the Panama Isthmus that interrupted gene flow between the two species. The polymorphism to divergence ratio was similar to the one observed for the other two loci, in accordance to the hypothesis of purifying selection. We conclude that slight purifying selection is likely to act on polymorphic intronic mutations of the EF1alpha second intron and discuss the possible relationship with the specific biology of Bathymodiolus mussels.     

Perspective: gene divergence, population divergence, and the variance in coalescence time in phylogeographic studies

Molecular methods as applied to the biogeography of single species (phylogeography) or multiple codistributed species (comparative phylogeography) have been productively and extensively used to elucidate common historical features in the diversification of the Earth's biota. However, only recently have methods for estimating population divergence times or their confidence limits while taking into account the critical effects of genetic polymorphism in ancestral species become available, and earlier methods for doing so are underutilized. We review models that address the crucial distinction between the gene divergence, the parameter that is typically recovered in molecular phylogeographic studies, and the population divergence, which is in most cases the parameter of interest and will almost always postdate the gene divergence. Assuming that population sizes of ancestral species are distributed similarly to those of extant species, we show that phylogeographic studies in vertebrates suggest that divergence of alleles in ancestral species can comprise from less than 10% to over 50% of the total divergence between sister species, suggesting that the problem of ancestral polymorphism in dating population divergence can be substantial. The variance in the number of substitutions (among loci for a given species or among species for a given gene) resulting from the stochastic nature of DNA change is generally smaller than the variance due to substitutions along allelic lines whose coalescence times vary due to genetic drift in the ancestral population. Whereas the former variance can be reduced by further DNA sequencing at a single locus, the latter cannot. Contrary to phylogeographic intuition, dating population divergence times when allelic lines have achieved reciprocal monophyly is in some ways more challenging than when allelic lines have not achieved monophyly, because in the former case critical data on ancestral population size provided by residual ancestral polymorphism is lost. In the former case differences in coalescence time between species pairs can in principle be explained entirely by differences in ancestral population size without resorting to explanations involving differences in divergence time. Furthermore, the confidence limits on population divergence times are severely underestimated when those for number of substitutions per site in the DNA sequences examined are used as a proxy. This uncertainty highlights the importance of multilocus data in estimating population divergence times; multilocus data can in principle distinguish differences in coalescence time (T) resulting from differences in population divergence time and differences in T due to differences in ancestral population sizes and will reduce the confidence limits on the estimates. We analyze the contribution of ancestral population size (theta) to T and the effect of uncertainty in theta on estimates of population divergence (tau) for single loci under reciprocal monophyly using a simple Bayesian extension of Takahata and Satta's and Yang's recent coalescent methods. The confidence limits on tau decrease when the range over which ancestral population size theta is assumed to be distributed decreases and when tau increases; they generally exclude zero when tau/(4Ne) > 1. We also apply a maximum-likelihood method to several single and multilocus data sets. With multilocus data, the criterion for excluding tau = 0 is roughly that l tau/(4Ne) > 1, where l is the number of loci. Our analyses corroborate recent suggestions that increasing the number of loci is critical to decreasing the uncertainty in estimates of population divergence time.     

Multilocus methods for estimating population sizes, migration rates and divergence time, with applications to the divergence of Drosophila pseudoobscura and D. persimilis

The genetic study of diverging, closely related populations is required for basic questions on demography and speciation, as well as for biodiversity and conservation research. However, it is often unclear whether divergence is due simply to separation or whether populations have also experienced gene flow. These questions can be addressed with a full model of population separation with gene flow, by applying a Markov chain Monte Carlo method for estimating the posterior probability distribution of model parameters. We have generalized this method and made it applicable to data from multiple unlinked loci. These loci can vary in their modes of inheritance, and inheritance scalars can be implemented either as constants or as parameters to be estimated. By treating inheritance scalars as parameters it is also possible to address variation among loci in the impact via linkage of recurrent selective sweeps or background selection. These methods are applied to a large multilocus data set from Drosophila pseudoobscura and D. persimilis. The species are estimated to have diverged approximately 500,000 years ago. Several loci have nonzero estimates of gene flow since the initial separation of the species, with considerable variation in gene flow estimates among loci, in both directions between the species.     

MULTIPLE EQUILIBRIA AND MAINTENANCE OF ADDITIVE GENETIC VARIANCE IN A MODEL OF PLEIOTROPY

We describe a multilocus model that incorporates pleiotropic stabilizing selection on a large number of characters. We find many different stable equilibria with different levels of polymorphism and additive genetic variability. The results lend support to Wright's concept of a complex adaptive surface with many peaks of different heights. The model assumes that alleles contribute additively to the characters. We analyze the multilocus model by first considering a two-locus model. The two-locus model depends critically on having loci of different effect and on having the optimum phenotype not be that of a completely heterozygous individual. The effects of different loci need to differ only by less than a factor of two. For the multilocus, multicharacter model, we assume that completely heterozygous individuals do not have the optimum phenotype. By restricting attention to a two-allele model, we also assume that there are no alleles that can affect all characters in all possible combinations of directions.     

Multilocus phylogenetics of a rapid radiation in the genus Thomomys (Rodentia: Geomyidae)

Species complexes undergoing rapid radiation present a challenge in molecular systematics because of the possibility that ancestral polymorphism is retained in component gene trees. Coalescent theory has demonstrated that gene trees often fail to match lineage trees when taxon divergence times are less than the ancestral effective population sizes. Suggestions to increase the number of loci and the number of individuals per taxon have been proposed; however, phylogenetic methods to adequately analyze these data in a coalescent framework are scarce. We compare two approaches to estimating lineage (species) trees using multiple individuals and multiple loci: the commonly used partitioned Bayesian analysis of concatenated sequences and a modification of a newly developed hierarchical Bayesian method (BEST) that simultaneously estimates gene trees and species trees from multilocus data. We test these approaches on a phylogeny of rapidly radiating species wherein divergence times are likely to be smaller than effective population sizes, and incomplete lineage sorting is known, in the rodent genus, Thomomys. We use seven independent noncoding nuclear sequence loci (total approximately 4300 bp) and between 1 and 12 individuals per taxon to construct a phylogenetic hypothesis for eight Thomomys species. The majority-rule consensus tree from the partitioned concatenated analysis included 14 strongly supported bipartitions, corroborating monophyletic species status of five of the eight named species. The BEST tree strongly supported only the split between the two subgenera and showed very low support for any other clade. Comparison of both lineage trees to individual gene trees revealed that the concatenation method appears to ignore conflicting signals among gene trees, whereas the BEST tree considers conflicting signals and downweights support for those nodes. Bayes factor analysis of posterior tree distributions from both analyses strongly favor the model underlying the BEST analysis. This comparison underscores the risks of overreliance on results from concatenation, and ignoring the properties of coalescence, especially in cases of recent, rapid radiations.     

Global convergence properties in multilocus viability selection models: the additive model and the Hardy-Weinberg law

A natural coordinate system is introduced for the analysis of the global stability of the Hardy-Weinberg (HW) polymorphism under the general multilocus additive viability model. A global convergence criterion is developed and used to prove that the HW polymorphism is globally stable when each of the loci is diallelic, provided the loci are overdominant and the multilocus recombination is positive. As a corollary the multilocus Hardy-Weinberg law for neutral selection is derived.Research supported in part by NIH grants GM 39907-01, GM 10452-26 and NSF Grant DMS 86-06244Research supported in part by a US-Israel Binational Science Foundation grant 85-00021 and NIH grant GM 28016     

Evolution and polymorphism in the multilocus Levene model with no or weak epistasis

Evolution and the maintenance of polymorphism under the multilocus Levene model with soft selection are studied. The number of loci and alleles, the number of demes, the linkage map, and the degree of dominance are arbitrary, but epistasis is absent or weak. We prove that, without epistasis and under mild, generic conditions, every trajectory converges to a stationary point in linkage equilibrium. Consequently, the equilibrium and stability structure can be determined by investigating the much simpler gene-frequency dynamics on the linkage-equilibrium manifold. For a haploid species an analogous result is shown. For weak epistasis, global convergence to quasi-linkage equilibrium is established. As an application, the maintenance of multilocus polymorphism is explored if the degree of dominance is intermediate at every locus and epistasis is absent or weak. If there are at least two demes, then arbitrarily many multiallelic loci can be maintained polymorphic at a globally asymptotically stable equilibrium. Because this holds for an open set of parameters, such equilibria are structurally stable. If the degree of dominance is not only intermediate but also deme independent, and loci are diallelic, an open set of parameters yielding an internal equilibrium exists only if the number of loci is strictly less than the number of demes. Otherwise, a fully polymorphic equilibrium exists only nongenerically, and if it exists, it consists of a manifold of equilibria. Its dimension is determined. In the absence of genotype-by-environment interaction, however, a manifold of equilibria occurs for an open set of parameters. In this case, the equilibrium structure is not robust to small deviations from no genotype-by-environment interaction. In a quantitative-genetic setting, the assumptions of no epistasis and intermediate dominance are equivalent to assuming that in every deme directional selection acts on a trait that is determined additively, i.e., by nonepistatic loci with dominance. Some of our results are exemplified in this quantitative-genetic context.     

Central equilibria in multilocus systems. I. Generalized nonepistatic selection regimes

The generalized nonepistatic selection regime encompasses combinations of multiplicative and neutral viability effects distributed across a set of loci. These subsume, in particular, mixtures of the classical modes of multiplicative and additive fitness evaluations for multilocus traits. Exact analytic conditions for existence and stability of a multilocus Hardy-Weinberg (H-W) polymorphic equilibrium configuration are ascertained. It is established that the central H-W polymorphism is stable only if the component loci are "over-dominant" and sufficient recombination is in force. The H-W central equilibrium is never stable for tight linkage whenever some multiplicative selection effects are contributed by at least two of the loci involved. In the case of additive selection expression and individual overdominant loci, the H-W polymorphism is stable independently of the level of recombination. In the context of "natural" recombination schemes, "more recombination" enhances the stability of the H-W polymorphic equilibrium.     

Extraordinary multilocus genetic organization in mole crickets, Gryllotalpidae

Abstract.— Allozymic diversity at 21 loci was analyzed in 470 individuals of three species of mole cricket superspecies, Gryllotalpa gryllotalpa (two new chromosomal species, G. tali and G. marismortui ) and G. africana in Israel, which are distributed along a southward transect of increasing aridity. Two outstanding findings emerged in G. tali and G. marismortui : (1) genetic polymorphism was high but heterozygosity very low, indicating significant deviations from Hardy-Weinberg expectations; and (2) significant linkage disequilibria at an unprecedented level for outbreeders and remarkable intersite differences. The results may characterize subterranean gryllotalpids worldwide because a single sample of Neocurtilla hexadactyla from Tefé, Amazonia, shows the same features. Significant variation of heterozygote paucity among loci , combined with the biology of the species, rejects the simple explanation of inbreeding or any other single explanatory model. Likewise, direct selection against heterozygotes or specific multilocus associations can explain, but is not necessary nor likely to explain, the observed results in mole crickets. To explain these results, we developed a multiple-factor mathematical model combining niche viability selection, niche choice, and positive assortative mating. This model involves a special case of Wahlund effect and inbreeding. Simulations based on this model showed that a combination of these three mechanisms may produce the observed distribution of alleles, via selection on a few loci, to affect the entire genome organization.     

Estimating species trees using multiple-allele DNA sequence data

Several techniques, such as concatenation and consensus methods, are available for combining data from multiple loci to produce a single statement of phylogenetic relationships. However, when multiple alleles are sampled from individual species, it becomes more challenging to estimate relationships at the level of species, either because concatenation becomes inappropriate due to conflicts among individual gene trees, or because the species from which multiple alleles have been sampled may not form monophyletic groups in the estimated tree. We propose a Bayesian hierarchical model to reconstruct species trees from multiple-allele, multilocus sequence data, building on a recently proposed method for estimating species trees from single allele multilocus data. A two-step Markov Chain Monte Carlo (MCMC) algorithm is adopted to estimate the posterior distribution of the species tree. The model is applied to estimate the posterior distribution of species trees for two multiple-allele datasets--yeast (Saccharomyces) and birds (Manacus-manakins). The estimates of the species trees using our method are consistent with those inferred from other methods and genetic markers, but in contrast to other species tree methods, it provides credible regions for the species tree. The Bayesian approach described here provides a powerful framework for statistical testing and integration of population genetics and phylogenetics.     

The long-term evolution of multilocus traits under frequency-dependent disruptive selection

Frequency-dependent disruptive selection is widely recognized as an important source of genetic variation. Its evolutionary consequences have been extensively studied using phenotypic evolutionary models, based on quantitative genetics, game theory, or adaptive dynamics. However, the genetic assumptions underlying these approaches are highly idealized and, even worse, predict different consequences of frequency-dependent disruptive selection. Population genetic models, by contrast, enable genotypic evolutionary models, but traditionally assume constant fitness values. Only a minority of these models thus addresses frequency-dependent selection, and only a few of these do so in a multilocus context. An inherent limitation of these remaining studies is that they only investigate the short-term maintenance of genetic variation. Consequently, the long-term evolution of multilocus characters under frequency-dependent disruptive selection remains poorly understood. We aim to bridge this gap between phenotypic and genotypic models by studying a multilocus version of Levene's soft-selection model. Individual-based simulations and deterministic approximations based on adaptive dynamics theory provide insights into the underlying evolutionary dynamics. Our analysis uncovers a general pattern of polymorphism formation and collapse, likely to apply to a wide variety of genetic systems: after convergence to a fitness minimum and the subsequent establishment of genetic polymorphism at multiple loci, genetic variation becomes increasingly concentrated on a few loci, until eventually only a single polymorphic locus remains. This evolutionary process combines features observed in quantitative genetics and adaptive dynamics models, and it can be explained as a consequence of changes in the selection regime that are inherent to frequency-dependent disruptive selection. Our findings demonstrate that the potential of frequency-dependent disruptive selection to maintain polygenic variation is considerably smaller than previously expected.     

A multilocus sequence survey in Arabidopsis thaliana reveals a genome-wide departure from a neutral model of DNA sequence polymorphism

The simultaneous analysis of multiple genomic loci is a powerful approach to studying the effects of population history and natural selection on patterns of genetic variation of a species. By surveying nucleotide sequence polymorphism at 334 randomly distributed genomic regions in 12 accessions of Arabidopsis thaliana, we examined whether a standard neutral model of nucleotide sequence polymorphism is consistent with observed data. The average nucleotide diversity was 0.0071 for total sites and 0.0083 for silent sites. Although levels of diversity are variable among loci, no correlation with local recombination rate was observed, but polymorphism levels were correlated for physically linked loci (<250 kb). We found that observed distributions of Tajima's D- and D/D(min)- and of Fu and Li's D-, D*- and F-, F*-statistics differed significantly from the expected distributions under a standard neutral model due to an excess of rare polymorphisms and high variances. Observed and expected distributions of Fay and Wu's H were not different, suggesting that demographic processes and not selection at multiple loci are responsible for the deviation from a neutral model. Maximum-likelihood comparisons of alternative demographic models like logistic population growth, glacial refugia, or past bottlenecks did not produce parameter estimates that were more consistent with observed patterns. However, exclusion of highly polymorphic "outlier loci" resulted in a fit to the logistic growth model. Various tests of neutrality revealed a set of candidate loci that may evolve under selection.     

A microsatellite-based multilocus screen for the identification of local selective sweeps

With the availability of completely sequenced genomes, multilocus scans of natural variability have become a feasible approach for the identification of genomic regions subjected to natural and artificial selection. Here, I introduce a new multilocus test statistic, ln RV, which is based on the ratio of observed variances in repeat number at a set of microsatellite loci in two groups of populations. The distribution of ln RV values captures demographic history of the populations as well as variation in microsatellite mutation among loci. Given that microsatellite loci associated with a recent selective sweep differ from the remainder of the genome, they are expected to fall outside of the distribution of neutral ln RV values. The ln RV test statistic is applied to a data set of 94 loci typed in eight non-African and two African human populations.     

Population subdivision and the Hudson-Kreitman-Aguade test: testing for deviations from the neutral model in organelle genomes

The Hudson-Kreitman-Aguade (HKA) test is based on the prediction from the neutral theory that levels of polymorphism within a species and the divergence between two closely related species should be correlated. Population subdivision has been shown to alter both the amounts of polymorphism segregating within species and the rate of divergence between species, meaning that genomic regions with different population structures also differ in their divergence to polymorphism ratios. Population subdivision may hence hamper the utility of the HKA test for detecting deviations from the standard neutral model, especially for organelle genomes that often have different patterns of population structure compared with nuclear genes. In this paper, I show that population subdivision inflates the number of instances where the HKA test detects deviations from the neutral model. Using coalescent simulations I show that this bias is most apparent when population subdivision is strong and differs substantially between the loci included. However, if divergence time is large and population structure substantial even changes in the levels of polymorphism and divergence associated with differences in the effective population size between two loci is enough to substantially alter the number of significant outcomes of the HKA test. A dataset on cytoplasmic diversity in Sileine vulgaris and S. latifolia (Ingvarsson & Taylor, 2002) is also reanalysed. The previous study had shown a marked excess of intraspecific polymorphism in both species. However, when effects of population subdivision were removed, ad hoc, levels of intraspecific polymorphism were no longer significantly different from neutral expectations, suggesting that population subdivision contributed to the observed excess of intraspecific polymorphism seen in both species of Silene.