Evolutionary and statistical properties of three genetic distances

Many genetic distances have been developed to summarize allele frequency differences between populations. I review the evolutionary and statistical properties of three popular genetic distances: DS, DA, and theta;, using computer simulation of two simple evolutionary histories: an isolation model of population divergence and an equilibrium migration model. The effect of effective population size, mutation rate, and mutation mechanism upon the parametric value between pairs of populations in these models explored, and the unique properties of each distance are described. The effect of these evolutionary parameters on study design is also investigated and similar results are found for each genetic distance in each model of evolution: large sample sizes are warranted when populations are relatively genetically similar; and loci with more alleles produce better estimates of genetic distance.     

Craniometric variation,genetic theory,and modern human origins

Recent controversies surrounding models of modern human origins have focused on among-group variation, particularly the reconstruction of phylogenetic trees from mitochondrial DNA (mtDNA) and, the dating of population divergence. Problems in tree estimation have been seen as weakening the case for a replacement model and favoring a multiregional evolution model. There has been less discussion of patterns of within-group variation, although the mtDNA evidence has consistently shown the greatest diversity within African populations. Problems of interpretation abound given the numerous factors that can influence within-group variation, including the possibility of earlier divergence, differences in population size, patterns of population expansion, and variation in migration rates. We present a model of within-group phenotypic variation and apply it to a large set of craniometric data representing major Old World geographic regions (57 measurements for 1,159 cases in four regions: Europe, Sub-Saharan Africa, Australasia, and the Far East). The model predicts a linear relationship between variation within populations (the average within-group variance) and variation between populations (the genetic distance of populations to pooled phenotypic means). On a global level this relationship should hold if the long-term effective population sizes of each region are correctly specified. Other potential effects on withingroup variation are accounted for by the model. Comparison of observed and expected variances under the assumption of equal effective sizes for four regions indicates significantly greater within-group variation in Africa and significantly less within-group variation in Europe. These results suggest that the long-term effective population size was greatest in Africa. Closer examination of the model suggests that the long-term African effective size was roughly three times that of any other geographic region. Using these estimates of relative population size, we present a method for analyzing ancient population structure, which provides estimates of ancient migration. This method allows us to reconstruct migration history between geographic regions after adjustment for the effect of genetic drift on interpopulational distances. Our results show a clear isolation of Africa from other regions. We then present a method that allows direct estimation of the ancient migration matrix, thus providing us with information on the actual extent of interregional migration. These methods also provide estimates of time frames necessary to reach genetic equilibrium. The ultimate goal is extracting as much information from present-day patterns of human variation relevannt to issues of human origins. Our results are in agreement with mismatch distribution analysis of mtDNA, and they support a “weak Garden o Eden” model. In this model, modern-day variation can be explained by divergence from an initial source (perhaps Africa) into a number o small isolated populations, followed by later population expansion throughout our species. The major populationn expansions of Homo sapiens during and after the late Pleistocene have had the effect of “freezing” ancient patterns of population structure. While this is not the only possible scenario, we do note the close agreement with ecent analyses of mtDNA mismatch distibutions. © 1994 Wiley-Liss, Inc.     

Estimation of hominoid ancestral population sizes under bayesian coalescent models incorporating mutation rate variation and sequencing errors

Estimation of population parameters for the common ancestors of humans and the great apes is important in understanding our evolutionary history. In particular, inference of population size for the human-chimpanzee common ancestor may shed light on the process by which the 2 species separated and on whether the human population experienced a severe size reduction in its early evolutionary history. In this study, the Bayesian method of ancestral inference of Rannala and Yang (2003. Bayes estimation of species divergence times and ancestral population sizes using DNA sequences from multiple loci. Genetics. 164:1645-1656) was extended to accommodate variable mutation rates among loci and random species-specific sequencing errors. The model was applied to analyze a genome-wide data set of approximately 15,000 neutral loci (7.4 Mb) aligned for human, chimpanzee, gorilla, orangutan, and macaque. We obtained robust and precise estimates for effective population sizes along the hominoid lineage extending back approximately 30 Myr to the cercopithecoid divergence. The results showed that ancestral populations were 5-10 times larger than modern humans along the entire hominoid lineage. The estimates were robust to the priors used and to model assumptions about recombination. The unusually low X chromosome divergence between human and chimpanzee could not be explained by variation in the male mutation bias or by current models of hybridization and introgression. Instead, our parameter estimates were consistent with a simple instantaneous process for human-chimpanzee speciation but showed a major reduction in X chromosome effective population size peculiar to the human-chimpanzee common ancestor, possibly due to selective sweeps on the X prior to separation of the 2 species.     

Genetic divergence of climatically marginal populations of Vicia pisiformis on the Scandinavian Peninsula

Vicia pisiformis L. is a perennial leguminous plant with a main distribution in broadleaved forest-steppes of eastern Europe. The species is classified as endangered (EN) according to the IUCN red-lists in both Norway and Sweden, due to severe fragmentation, small population sizes and continuing population decline. The populations on the Scandinavian Peninsula constitute the northern limit of the species distribution and are mostly restricted to warm stony slopes with predominantly southern aspects. In this study we used the AFLP method, which is a high-resolution genetic fingerprint method. Samples were collected from 22 Scandinavian populations. The overall genetic structure was analysed in an AMOVA, in a Mantel test and through constrained correspondence analysis (CCA). The ordination scores representing non-geographic genetic divergence were extracted from the CCA and analysed in a linear model using habitat variables and population size as explanatory variables. We found (i) a strong geographic structure, (ii) significant genetic divergence between populations, (iii) that this genetic divergence remained significant even after removing the effect of geography in a partial CCA and (iv) that the remaining non-geographic part of genetic divergence (distance from the ordination centre) was associated with aspect, populations with a northern aspect were more genetically divergent. Aspect explains more variation than population size and is the only variable retained in the minimal adequate model. We suggest that local adaptation has caused this divergence from an expected geographical pattern of genetic variation. This explanation is further supported by the association between aspect and specific AFLP fragments.     

Estimating divergence times from DNA sequences

The patterns of genetic variation within and among individuals and populations can be used to make inferences about the evolutionary forces that generated those patterns. Numerous population genetic approaches have been developed in order to infer evolutionary history. Here, we present the “Two-Two (TT)” and the “Two-Two-outgroup (TTo)” methods; two closely related approaches for estimating divergence time based in coalescent theory. They rely on sequence data from two haploid genomes (or a single diploid individual) from each of two populations. Under a simple population-divergence model, we derive the probabilities of the possible sample configurations. These probabilities form a set of equations that can be solved to obtain estimates of the model parameters, including population split times, directly from the sequence data. This transparent and computationally efficient approach to infer population divergence time makes it possible to estimate time scaled in generations (assuming a mutation rate), and not as a compound parameter of genetic drift. Using simulations under a range of demographic scenarios, we show that the method is relatively robust to migration and that the TTo method can alleviate biases that can appear from drastic ancestral population size changes. We illustrate the utility of the approaches with some examples, including estimating split times for pairs of human populations as well as providing further evidence for the complex relationship among Neandertals and Denisovans and their ancestors.     

Genetic divergence disclosing a rapid prehistorical dispersion of native americans in central and South america

An accurate estimate of the divergence time between Native Americans is important for understanding the initial entry and early dispersion of human beings in the New World. Current methods for estimating the genetic divergence time of populations could seriously depart from a linear relationship with the true divergence for multiple populations of a different population size and significant population expansion. Here, to address this problem, we propose a novel measure to estimate the genetic divergence time of populations. Computer simulation revealed that the new measure maintained an excellent linear correlation with the population divergence time in complicated multi-population scenarios with population expansion. Utilizing the new measure and microsatellite data of 21 Native American populations, we investigated the genetic divergences of the Native American populations. The results indicated that genetic divergences between North American populations are greater than that between Central and South American populations. None of the divergences, however, were large enough to constitute convincing evidence supporting the two-wave or multi-wave migration model for the initial entry of human beings into America. The genetic affinity of the Native American populations was further explored using Neighbor-Net and the genetic divergences suggested that these populations could be categorized into four genetic groups living in four different ecologic zones. The divergence of the population groups suggests that the early dispersion of human beings in America was a multi-step procedure. Further, the divergences suggest the rapid dispersion of Native Americans in Central and South Americas after a long standstill period in North America.     

Inferring Demographic History from a Spectrum of Shared Haplotype Lengths

There has been much recent excitement about the use of genetics to elucidate ancestral history and demography. Whole genome data from humans and other species are revealing complex stories of divergence and admixture that were left undiscovered by previous smaller data sets. A central challenge is to estimate the timing of past admixture and divergence events, for example the time at which Neanderthals exchanged genetic material with humans and the time at which modern humans left Africa. Here, we present a method for using sequence data to jointly estimate the timing and magnitude of past admixture events, along with population divergence times and changes in effective population size. We infer demography from a collection of pairwise sequence alignments by summarizing their length distribution of tracts of identity by state (IBS) and maximizing an analytic composite likelihood derived from a Markovian coalescent approximation. Recent gene flow between populations leaves behind long tracts of identity by descent (IBD), and these tracts give our method power by influencing the distribution of shared IBS tracts. In simulated data, we accurately infer the timing and strength of admixture events, population size changes, and divergence times over a variety of ancient and recent time scales. Using the same technique, we analyze deeply sequenced trio parents from the 1000 Genomes project. The data show evidence of extensive gene flow between Africa and Europe after the time of divergence as well as substructure and gene flow among ancestral hominids. In particular, we infer that recent African-European gene flow and ancient ghost admixture into Europe are both necessary to explain the spectrum of IBS sharing in the trios, rejecting simpler models that contain less population structure.     

On the size distribution of private microsatellite alleles

Private microsatellite alleles tend to be found in the tails rather than in the interior of the allele size distribution. To explain this phenomenon, we have investigated the size distribution of private alleles in a coalescent model of two populations, assuming the symmetric stepwise mutation model as the mode of microsatellite mutation. For the case in which four alleles are sampled, two from each population, we condition on the configuration in which three distinct allele sizes are present, one of which is common to both populations, one of which is private to one population, and the third of which is private to the other population. Conditional on this configuration, we calculate the probability that the two private alleles occupy the two tails of the size distribution. This probability, which increases as a function of mutation rate and divergence time between the two populations, is seen to be greater than the value that would be predicted if there was no relationship between privacy and location in the allele size distribution. In accordance with the prediction of the model, we find that in pairs of human populations, the frequency with which private microsatellite alleles occur in the tails of the allele size distribution increases as a function of genetic differentiation between populations.     

Detection of barriers to dispersal is masked by long lifespans and large population sizes

Population genetic analyses of species inhabiting fragmented landscapes are essential tools for conservation. Occasionally, analyses of fragmented populations find no evidence of isolation, even though a barrier to dispersal is apparent. In some cases, not enough time may have passed to observe divergence due to genetic drift, a problem particularly relevant for long‐lived species with overlapping generations. Failing to consider this quality during population structure analyses could result in incorrect conclusions about the impact of fragmentation on the species. We designed a model to explore how lifespan and population size influence perceived population structure of isolated populations over time. This iterative model tracked how simulated populations of variable lifespan and population size were affected by drift alone, using a freshwater mussel, Quadrula quadrula (mapleleaf), as a model system. In addition to exhibiting dramatic lifespan variability among species, mussels are also highly imperiled and exhibit fragmentation by dams throughout the range of many species. Results indicated that, unless population size was small (<50 individuals) or lifespan short (<22 years), observing genetic divergence among populations was unlikely. Even if wild populations are isolated, observing population structure in long‐lived mussels from modern damming practices is unlikely because it takes longer for population structure to develop in these species than most North American dams have existed. Larger population sizes and longer lifespans increase the time needed for significant divergence to occur. This study helps illuminate the factors that influence genetic responses by populations to isolation and provides a useful model for conservation‐oriented research.     

The Allelic Correlation Structure of Gainj- and Kalam-Speaking People. II. the Genetic Distance between Population Subdivisions

The patterning of allele frequency variability among 18 local groups of Gainj and Kalam speakers of highland Papua New Guinea is investigated using new genetic distance methods. The genetic distances proposed here are obtained by decomposing Sewall Wright's coefficient FST into a set of coefficients corresponding to all pairs of population subdivisions. Two statistical methods are given to estimate these quantities. One method provides estimates weighted by sample sizes, while the other method does not use sample size weighting. Both methods correct for the within-individual and between-individual-within-groups sums of squares. Genetic distances among the Gainj and Kalam subdivisions are analyzed with respect to demographic, geographic, and linguistic variables. We find that a demographic feature, group size, has the greatest demonstrable association with the patterning of genetic distances. The pattern of geographic distances among groups displays a weak congruence with the pattern of genetic distances, and the association of genetic and linguistic diversity is very low. An effect of differences in group size on genetic distances is not surprising, from basic theoretical considerations, but genetic distances have not often been analyzed with respect to these variables in the past. The lack of correspondence between genetic distances and linguistic and geographic differences is an unusual feature that distinguishes the Gainj and Kalam from most other tribal populations.     

Stochastic interplay between mutation and recombination during the acquisition of drug resistance mutations in human immunodeficiency virus type 1

The emergence of drug resistance mutations in human immunodeficiency virus (HIV) has been a major setback in the treatment of infected patients. Besides the high mutation rate, recombination has been conjectured to have an important impact on the emergence of drug resistance. Population genetic theory suggests that in populations limited in size recombination may facilitate the acquisition of beneficial mutations. The viral population in an infected patient may indeed represent such a population limited in size, since current estimates of the effective population size range from 500 to 10(5). To address the effects of limited population size, we therefore expand a previously described deterministic population genetic model of HIV replication by incorporating the stochastic processes that occur in finite populations of infected cells. Using parameter estimates from the literature, we simulate the evolution of drug-resistant viral strains. The simulations show that recombination has only a minor effect on the rate of acquisition of drug resistance mutations in populations with effective population sizes as small as 1,000, since in these populations, viral strains typically fix beneficial mutations sequentially. However, for intermediate effective population sizes (10(4) to 10(5)), recombination can accelerate the evolution of drug resistance by up to 25%. Furthermore, a reduction in population size caused by drug therapy can be overcome by a higher viral mutation rate, leading to a faster evolution of drug resistance.     

Likelihood analysis of ongoing gene flow and historical association

Abstract.— We develop a Monte Carlo-based likelihood method for estimating migration rates and population divergence times from data at unlinked loci at which mutation rates are sufficiently low that, in the recent past, the effects of mutation can be ignored. The method is applicable to restriction fragment length polymorphisms (RFLPs) and single nucleotide polymorphisms (SNPs) sampled from a subdivided population. The method produces joint maximum-likelihood estimates of the migration rate and the time of population divergence, both scaled by population size, and provides a framework in which to test either for no ongoing gene flow or for population divergence in the distant past. We show the method performs well and provides reasonably accurate estimates of parameters even when the assumptions under which those estimates are obtained are not completely satisfied. Furthermore, we show that, provided that the number of polymorphic loci is sufficiently large, there is some power to distinguish between ongoing gene flow and historical association as causes of genetic similarity between pairs of populations.     

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.     

Biased mutations and microsatellite variation

Mutation bias is one of the forces that may constrain the variation at microsatellite loci. Here, we study the dynamics of population statistics and the genetic distance between two populations under multiple stepwise mutations with linear bias and random drift. Expressions are derived for these statistics as functions of time, as well as at mutation-drift equilibrium. Applying these expressions to published data on humans and chimpanzees, the regression coefficient of mutation bias on allele size was estimated to be at least between - 0.0064 and -0.013. The assumption of mutational bias produces larger estimates of divergence times than are obtained in its absence; in particular, the time of split between African and non-African human populations is estimated to be between 183,000 and 222,000 years, assuming one-step mutations and no selection. With multistep mutations, the divergence time is estimated to be lower.      

Population genetic structure,differentiation, and diversity in Tetrix subulata pygmy grasshoppers: roles of population size and immigration

Genetic diversity within and among populations and species is influenced by complex demographic and evolutionary processes. Despite extensive research, there is no consensus regarding how landscape structure, spatial distribution, gene flow, and population dynamics impact genetic composition of natural populations. Here, we used amplified fragment length polymorphisms (AFLPs) to investigate effects of population size, geographic isolation, immigration, and gene flow on genetic structure, divergence, and diversity in populations of Tetrix subulata pygmy grasshoppers (Orthoptera: Tetrigidae) from 20 sampling locations in southern Sweden. Analyses of 1564 AFLP markers revealed low to moderate levels of genetic diversity (PPL = 59.5–90.1; Hj = 0.23–0.32) within and significant divergence among sampling localities. This suggests that evolution of functional traits in response to divergent selection is possible and that gene flow is restricted. Genetic diversity increased with population size and with increasing proportion of long‐winged phenotypes (a proxy of recent immigration) across populations on the island of Öland, but not on the mainland. Our data further suggested that the open water separating Öland from the mainland acts as a dispersal barrier that restricts migration and leads to genetic divergence among regions. Isolation by distance was evident for short interpopulation distances on the mainland, but gradually disappeared as populations separated by longer distances were included. Results illustrate that integrating ecological and molecular data is key to identifying drivers of population genetic structure in natural populations. Our findings also underscore the importance of landscape structure and spatial sampling scheme for conclusions regarding the role of gene flow and isolation by distance.     

Genetic Diversity in Introduced Populations with an Allee Effect

A phenomenon that strongly influences the demography of small introduced populations and thereby potentially their genetic diversity is the demographic Allee effect, a reduction in population growth rates at small population sizes. We take a stochastic modeling approach to investigate levels of genetic diversity in populations that successfully overcame either a strong Allee effect, in which populations smaller than a certain critical size are expected to decline, or a weak Allee effect, in which the population growth rate is reduced at small sizes but not negative. Our results indicate that compared to successful populations without an Allee effect, successful populations with a strong Allee effect tend to (1) derive from larger founder population sizes and thus have a higher initial amount of genetic variation, (2) spend fewer generations at small population sizes where genetic drift is particularly strong, and (3) spend more time around the critical population size and thus experience more genetic drift there. In the case of multiple introduction events, there is an additional increase in diversity because Allee-effect populations tend to derive from a larger number of introduction events than other populations. Altogether, a strong Allee effect can either increase or decrease genetic diversity, depending on the average founder population size. By contrast, a weak Allee effect tends to decrease genetic diversity across the entire range of founder population sizes. Finally, we show that it is possible in principle to infer critical population sizes from genetic data, although this would require information from many independently introduced populations.     

Characterising genetic diversity and effective population size in one reservoir and two riverine populations of the threatened Macquarie perch

Nine polymorphic microsatellite loci have been used to infer population genetic diversity and structure of the threatened Australian freshwater fish, Macquarie perch, across three tributaries of the Murrumbidgee River in south-eastern Australia. This investigation has revealed a high level of divergence among all three populations, along with contrasting patterns of genetic diversity. The Cotter Reservoir, which is a stronghold population for the species, has typically higher diversity and effective population size than nearby riverine populations. This suggests that the reservoir population is unlikely to have undergone a genetic bottleneck during and following dam construction. Genetic diversity estimates were comparable with one riverine site but were significantly higher than a population sampled from the Queanbeyan River. This comparison revealed significantly less heterozygotes in the Queanbeyan River and lower estimates of effective population size. Options and considerations for stock replenishment of this population are discussed.     

Relative influences of historical and contemporary forces shaping the distribution of genetic variation in the Atlantic killifish, Fundulus heteroclitus

A major goal of population genetics research is to identify the relative influences of historical and contemporary processes that serve to structure genetic variation. Most population genetic models assume that populations exist in a state of migration-drift equilibrium. However, in the past this assumption has rarely been verified, and is likely rarely achieved in natural populations. We assessed the equilibrium status at both local and regional scales of the Atlantic killifish, Fundulus heteroclitus . This species is a model organism for the study of adaptive clinal variation, but has also experienced a complicated history of range expansion and secondary contact following allopatric divergence, potentially obscuring the influence of contemporary evolutionary processes. Presumptively neutral genetic markers (microsatellites) demonstrated zones of secondary intergradation among coastal populations centred around northern New Jersey and the Chesapeake Bay region. Analysis of genetic variation indicated isolation by distance among some populations and provided supporting evidence that the Delaware Bay, but not the Chesapeake Bay, has acted as a barrier to dispersal among coastal populations. Bayesian estimates indicated large effective population sizes and low migration rates, and were in good agreement with empirically derived estimates of population and neighbourhood size from mark–recapture studies. These data indicate that populations are not in migration-drift equilibrium at a regional scale, and suggest that contributing factors include large population size combined with relatively low migration rates. These conditions should be considered when interpreting the evolutionary significance of the distribution of genetic variation among F. heteroclitus populations.     

Likelihood and Bayes estimation of ancestral population sizes in hominoids using data from multiple loci

Polymorphisms in an ancestral population can cause conflicts between gene trees and the species tree. Such conflicts can be used to estimate ancestral population sizes when data from multiple loci are available. In this article I extend previous work for estimating ancestral population sizes to analyze sequence data from three species under a finite-site nucleotide substitution model. Both maximum-likelihood (ML) and Bayes methods are implemented for joint estimation of the two speciation dates and the two population size parameters. Both methods account for uncertainties in the gene tree due to few informative sites at each locus and make an efficient use of information in the data. The Bayes algorithm using Markov chain Monte Carlo (MCMC) enjoys a computational advantage over ML and also provides a framework for incorporating prior information about the parameters. The methods are applied to a data set of 53 nuclear noncoding contigs from human, chimpanzee, and gorilla published by Chen and Li. Estimates of the effective population size for the common ancestor of humans and chimpanzees by both ML and Bayes methods are approximately 12,000-21,000, comparable to estimates for modern humans, and do not support the notion of a dramatic size reduction in early human populations. Estimates published previously from the same data are several times larger and appear to be biased due to methodological deficiency. The divergence between humans and chimpanzees is dated at approximately 5.2 million years ago and the gorilla divergence 1.1-1.7 million years earlier. The analysis suggests that typical data sets contain useful information about the ancestral population sizes and that it is advantageous to analyze data of several species simultaneously.     

Effective population sizes and temporal stability of genetic structure in Rana pipiens, the northern leopard frog

Although studies of population genetic structure are very common, whether genetic structure is stable over time has been assessed for very few taxa. The question of stability over time is particularly interesting for frogs because it is not clear to what extent frogs exist in dynamic metapopulations with frequent extinction and recolonization, or in stable patches at equilibrium between drift and gene flow. In this study we collected tissue samples from the same five populations of leopard frogs, Rana pipiens, over a 22-30 year time interval (11-15 generations). Genetic structure among the populations was very stable, suggesting that these populations were not undergoing frequent extinction and colonization. We also estimated the effective size of each population from the change in allele frequencies over time. There exist few estimates of effective size for frog populations, but the data available suggest that ranid frogs may have much larger ratios of effective size (Ne) to census size (Nc) than toads (bufonidae). Our results indicate that R. pipiens populations have effective sizes on the order of hundreds to at most a few thousand frogs, and Ne/Nc ratios in the range of 0.1-1.0. These estimates of Ne/Nc are consistent with those estimated for other Rana species. Finally, we compared the results of three temporal methods for estimating Ne. Moment and pseudolikelihood methods that assume a closed population gave the most similar point estimates, although the moment estimates were consistently two to four times larger. Wang and Whitlock's new method that jointly estimates Ne and the rate of immigration into a population (m) gave much smaller estimates of Ne and implausibly large estimates of m. This method requires knowing allele frequencies in the source of immigrants, but was thought to be insensitive to inexact estimates. In our case the method may have failed because we did not know the true source of immigrants for each population. The method may be more sensitive to choice of source frequencies than was previously appreciated, and so should be used with caution if the most likely source of immigrants cannot be identified clearly.